A secreted LysM effector protects fungal hyphae through chitin-dependent homodimer polymerization
Authors:
Andrea Sánchez-Vallet aff001; Hui Tian aff001; Luis Rodriguez-Moreno aff001; Dirk-Jan Valkenburg aff001; Raspudin Saleem-Batcha aff003; Stephan Wawra aff004; Anja Kombrink aff001; Leonie Verhage aff001; Ronnie de Jonge aff001; H. Peter van Esse aff001; Alga Zuccaro aff004; Daniel Croll aff005; Jeroen R. Mesters aff003; Bart P. H. J. Thomma aff001
Authors place of work:
Laboratory of Phytopathology, Wageningen University& Research, Wageningen, The Netherlands
aff001; Plant Pathology Group, Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland
aff002; Institute of Biochemistry, Center for Structural and Cell Biology in Medicine, University of Lübeck, Lübeck, Germany
aff003; University of Cologne, Institute for Plant Sciences, Cluster of Excellence on Plant Sciences (CEPLAS), Cologne, Germany
aff004; Institute of Biology, University of Neuchâtel, Neuchâtel, Switzerland
aff005
Published in the journal:
A secreted LysM effector protects fungal hyphae through chitin-dependent homodimer polymerization. PLoS Pathog 16(6): e32767. doi:10.1371/journal.ppat.1008652
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1008652
Summary
Plants trigger immune responses upon recognition of fungal cell wall chitin, followed by the release of various antimicrobials, including chitinase enzymes that hydrolyze chitin. In turn, many fungal pathogens secrete LysM effectors that prevent chitin recognition by the host through scavenging of chitin oligomers. We previously showed that intrachain LysM dimerization of the Cladosporium fulvum effector Ecp6 confers an ultrahigh-affinity binding groove that competitively sequesters chitin oligomers from host immune receptors. Additionally, particular LysM effectors are found to protect fungal hyphae against chitinase hydrolysis during host colonization. However, the molecular basis for the protection of fungal cell walls against hydrolysis remained unclear. Here, we determined a crystal structure of the single LysM domain-containing effector Mg1LysM of the wheat pathogen Zymoseptoria tritici and reveal that Mg1LysM is involved in the formation of two kinds of dimers; a chitin-dependent dimer as well as a chitin-independent homodimer. In this manner, Mg1LysM gains the capacity to form a supramolecular structure by chitin-induced oligomerization of chitin-independent Mg1LysM homodimers, a property that confers protection to fungal cell walls against host chitinases.
Keywords:
Plant fungal pathogens – Fungal structure – Crystal structure – Plant cell walls – Dimerization – Hydrolysis – Polymerization – chitin
Introduction
Fungi constitute an evolutionarily and ecologically diverse group of microorganisms that are characterized by the presence of chitin, an N-acetyl-D-glucosamine (GlcNAc) homopolymer, in their cell walls. In addition to providing strength, shape, rigidity and protection against environmental hazards, chitin is also a well-known inducer of plant immune responses [1,2,3]. A major mechanism of plant defense against fungal invasion includes the secretion of microbial cell wall-degrading enzymes that include chitin-degrading enzymes, known as chitinases, to hinder fungal pathogen ingress [4,5]. Plant chitinases are diverse in nature, and grouped into six different classes that belong to glycosyl hydrolase families 18 and 19 [6,7]. Although many chitinases are specifically produced upon pathogen invasion, others are constitutively expressed [4]
Chitin hydrolysis on the one hand inhibits fungal growth, and on the other hand releases chitooligosacharides [8,9](Kasprzewska, 2003; Liu et al, 2014) that are recognized by cell surface receptors of host cells to mount an appropriate immune response [3,10]. In plants, chitin is recognized in the extracellular space through membrane-exposed Lysin motif (LysM)-containing receptor molecules [1,11,12,13]. In turn, many successful fungal pathogens have evolved effector molecules that either protect their cell walls against plant chitinases or prevent or perturb the elicitation of chitin-triggered host immunity [3,10,14,15,16].
Since decades, the interaction between the foliar fungal pathogen Cladosporium fulvum and its only host tomato has been studied to unravel the role of pathogen virulence and host defense mechanisms, including mechanisms that revolve around chitin [17]. After leaf penetration, C. fulvum secretes an arsenal of apoplastic effector proteins, including the chitin-binding effector proteins Avr4, which protects fungal hyphae against hydrolysis by plant chitinases [18,19], and Ecp6, which perturbs the activation of chitin-triggered host immunity [14,20]. Whereas Avr4 binds chitin through an invertebrate chitin-binding domain [21,22], Ecp6 utilizes LysM domains for chitin binding [14,20]. Previous biochemical analysis revealed that Avr4 monomers require a stretch of at least three exposed GlcNAc residues for binding, and positive allosteric interactions among Avr4 molecules occur during chitin binding to facilitate the shielding of cell wall chitin against host chitinases [21]. Based on X-ray crystallography it was recently shown that two Avr4 molecules interact through their chitohexaose ligand to form a three-dimensional molecular sandwich that encapsulates two chitohexaose molecules within the dimeric assembly [23]. A crystal structure of Ecp6 revealed chitin-induced dimerization of two of the three LysM domains, resulting in the formation of an ultrahigh affinity (pM) chitin-binding groove, conferring the capacity to outcompete plant receptors for chitin binding [20]. Interestingly, whereas Avr4 homologs are found in other, C. fulvum-related, Dothideomycete plant pathogens [24], LysM effectors are widespread in the fungal kingdom [25,26,27]. In several plant pathogenic fungi, including the Dothideomycete Zymoseptoria tritici and the Sodariomycetes Magnaporthe oryzae, Colletotrichum higginsianum and Verticillium dahliae, LysM effectors have been shown to contribute to virulence through chitin binding [15,28,29,30].
The LysM effectors Mg1LysM and Mg3LysM, with one and three LysM domains, respectively, have been characterized from the Septoria tritici blotch pathogen of wheat, Z. tritici [15]. Functional characterization has revealed that Mg3LysM can suppress chitin-induced immunity in a similar fashion as C. fulvum Ecp6. Surprisingly, in contrast to Ecp6 and similar to Avr4, Mg3LysM was additionally shown to have the ability to protect fungal hyphae against chitinase hydrolysis. As expected, based on the presence of a single LysM domain only, a role in suppression of chitin-triggered immunity could not be demonstrated for Mg1LysM [15]. Intriguingly, however, Mg1LysM was characterized as a functional homolog of Avr4 that protects hyphae against hydrolysis by host chitinases [15]. In order to understand how a LysM effector that is composed from little more than only a single LysM domain is able to confer protection of cell wall chitin from hydrolysis by plant enzymes, we aimed to obtain a crystal structure of the Z. tritici effector Mg1LysM in this study. Surprisingly, we discovered that Mg1LysM has the ability to simultaneously undergo ligand-mediated dimerization as well as ligand-independent homodimerization, allowing the formation of a contiguous oligomeric structure that anchors to the fungal cell wall through chitin to confer its protection ability.
Results
Crystal structure of Mg1LysM reveals ligand-dependent and -independent intermolecular dimerization
In order to understand LysM effector functionality, and particularly how Mg1LysM is able to protect chitin against chitinase hydrolysis, a crystal structure of Mg1LysM was determined. To this end, Mg1LysM was heterologously produced in the yeast Pichia pastoris and purified based on the presence of a His-FLAG affinity tag. The large Mg1LysM crystals that were finally obtained by micro-seeding techniques [31]belonged to the space group P 61 2 2. Some crystals were soaked with the Ta6Br14 cluster and initial phases were determined by the single-wavelength anomalous dispersion technique (SAD; Table 1). The initial phases were further improved with the help of an I3C soaked crystal. A native dataset was finally refined to a resolution of 2.41 Å with an Rwork and Rfree of 17.96% and 22.03%, respectively (Table 1). The structure model comprises in total four copies of the complete mature protein sequence except for the first amino acid after the signal peptide for chains A to C and the first two amino acids after the signal peptide for chain D, and one carbohydrate molecule per asymmetric unit (a.u.).
As expected, the tertiary structure of the LysM domain of an Mg1LysM monomer is similar to that of previously described LysM domains [13,20,32,33,34,35]with a conserved βααβ-fold in which the antiparallel β-sheet lies adjacent to two α-helices (Fig 1 and S1 Fig). The compact LysM structure is stabilized by two disulfide bridges between Cys44 and Cys78, and between Cys13 and Cys70 (Fig 1 and S1 Fig). In addition to the single LysM domain, Mg1LysM comprises a relatively long N-terminal tail that contains a short β-strand (Fig 1).
The four Mg1LysM monomers within the a.u. form two pairs of homodimers that are each very tightly packed. The large homodimerization interface that occurs between two monomers was calculated to be 1113 Å2 using PISA (Protein Interfaces, Surfaces and Assemblies; http://www.ebi.ac.uk/pdbe/prot_int/pistart.html) [36]and is stabilized by a total of 25 hydrogen bonds between residues of each of the monomers. In addition, the crystal structure revealed that the N-terminal 12 amino acid tails of the homodimer run anti-parallel and form a small but stable β-sheet (5ITI7 of each chain) that is stabilized by clustering of all four isoleucine side-chain residues and by threonine-threonine sidechain hydrogen bonding, potentially strengthening the homodimer (Fig 1). The latter hypothesis is further supported by the formation of two additional salt bridges formed between Arg2 of one subunit and Asp12 of the other one (Fig 1). The root mean square deviations (rmsd) between the Cα atoms of the two homodimers of the a.u. is 0.267 Å as calculated with Lsqkab of the CCP4 suite [37].
Surprisingly, when determining the crystal structure for the C. fulvum LysM effector Ecp6 in the absence of exogenously added chitin we found a chitin tetramer in a large interdomain groove between two of the three intrachain LysM domains that appeared to constitute an ultra-high affinity binding site, while no chitin binding was observed to the remaining, third LysM domain of Ecp6 [20]. Unexpectedly, the calculated 2|F0|–|Fc| electron density map of the Mg1LysM crystal structure assembly similarly revealed well-defined electron density for a single chitin trimer bound to one monomer of the a.u. only (S2 Fig). Inspection of the crystal packing interactions unveiled the presence of a chitin binding pocket formed between two Mg1LysM protomers of neighbouring homodimers (Fig 1). Since protein purification and crystallization was performed without exogenous addition of chitin in this case as well, we concluded that the chitin once more was derived from the cell wall of the heterologous protein production host P. pastoris. Potentially, this finding indicates that Mg1LysM displays an increased affinity for chitin (low micromolar range) when compared with other, single-acting, LysM domains. The chitin binding site is formed by the loops between the first ß-strand and the first α-helix and between the second α-helix and the second ß-strand of Mg1LysM, encompassing the residues 26GDTLT30 and 56NRI58 that are conserved in many other LysM domains including those of Ecp6 [20](S1 Fig). Remarkably, besides the ligand-independent Mg1LysM homodimerization described above, the crystal structure revealed that chitin induces a dimerization of homodimers and, consequently, that a chitin-binding groove is formed by two LysM domains from two independent protomers (Fig 1, Fig 2). In addition to the amino acids that are in direct contact with the chitin trimer, the ligand-induced dimerization is strengthened by several hydrogen bonds that occur between residues from the two protomers involved. One salt bridge between residues K31 and D54 of the different protomers stabilizes the binding of the single chitin molecule and adds further strength to the dimerization, resulting in a tight binding pocket in which the chitin trimer is strongly bound (Fig 2). Arguably, we would expect an increased chitin-binding affinity of Mg1LysM when compared with a single-acting LysM domain, which can explain in turn why the chitin remained adhered to the Mg1LysM protein during the protein purification procedure.
In order to confirm the chitin-binding activity and determine the Mg1LysM chitin-binding affinity, a polysaccharide affinity precipitation assay and isothermal titration calorimetry (ITC) analysis were used, respectively. Since the crystal structure revealed that a portion of the Mg1LysM binding sites were occupied by chitin in the P. pastoris-produced Mg1LysM preparation, we pursued production of Mg1LysM in the bacterium Escherichia coli as a heterologous system that is devoid of chitin, in order to obtain chitin-free protein. The precipitation assay confirmed the binding of Mg1LysM to chitin and not to other insoluble polysaccharides [15] (S3 Fig). Subsequent ITC analysis based on chitohexaose (GlcNAc)6 titration revealed that this protein preparation bound chitin with a binding affinity of 4.36 μM (Fig 3). As previously demonstrated [20], P. pastoris-produced Ecp6 monomers bind chitin with a stoichiometry of 1:1 (S3 Fig). In contrast, a stoichiometry of 1:2 was observed (n = 0.504) for Mg1LysM based on a single-binding-site model, analogous to the observation that two Mg1LysM protomers originating from two Mg1LysM homodimers bind a single chitin trimer as disclosed in the crystal structure model. Obviously, this ratio also implies a polymerisation reaction in solution upon addition of the ligand chitohexaose. Thus, this finding can be interpreted as an independent confirmation of the ligand-induced Mg1LysM polymerisation as observed in the crystal structure.
Mg1LysM sequence conservation in a world-wide collection of Z. tritici isolates
In order to evaluate Mg1LysM conservation in Z. tritici, the occurrence of sequence polymorphisms was evaluated in a collection of 149 isolates from four different populations collected in Switzerland, Australia, Israel and the USA [38]. This analysis revealed that the Mg1LysM protein sequence is highly conserved (S1 Table). Only five non-synonymous mutations were identified in the full length Mg1LysM protein, three of which were previously identified [15]. Interestingly, none of these non-synonymous SNPs localized within the signal peptide, the homodimerization surface, the chitin-binding groove or concerned the residues involved in disulfide or salt bridge formation (Fig 4A), pointing towards the relevance of these sites for the functionality of Mg1LysM. To test the impact of these polymorphisms on chitin binding, we heterologously produced two allelic variants of Mg1LysM in E. coli, namely the variants from the Z. tritici isolates ST99-CH1E4 and ST99-CH3F4 that, collectively, carry the five non-synonymous mutations (Fig 4C; S1 Table). Whereas both strains share the polymorphisms N3Q and Q34K, strain ST99-CH1E4 additionally carries R24Q, while strain ST99-CH3F4 additionally carries the polymorphisms R48K and Q20T (Fig 4C; S1 Table). Interestingly, a polysaccharide affinity precipitation assay revealed that, like the wild-type protein Mg1LysM, the two variants Mg1LysM_1E4 and Mg1LysM_3F4 still bind chitin (Fig 4D), suggesting that the allelic variants have retained their biological activity.
Chitin-induced Mg1LysM polymerization is crucial for protection of hyphae against the hydrolytic activity of plant chitinases
We previously showed that Mg1LysM is able to protect chitin against chitinase hydrolysis [15]. Localization experiments making use of BODIPY-labelled Mg1LysM protein exogenously applied to Trichoderma viride, a fungal species that exposes cell wall chitin during growth in vitro, revealed that Mg1LysM binds to fungal cell walls (Fig 5A). Next, we attempted to evaluate the contribution of the ligand-independent Mg1LysM homodimerization to hyphal protection against chitinases. To this end, we pursued to produce an Mg1LysM mutant that lacked the 12-amino acid tail that is, besides the large protomer-protomer interface, responsible for ligand-independent homodimerization. Unfortunately, production of this mutant in the heterologous host P. pastoris was not successful as hardly any protein could be detected. The protein is apparently degraded either due to exposure of the hydrophobic residues (V40 and I68) located at the centre of the large (1113 Å2) homodimerization interface, or homodimerization is stringently required for proper folding of the protein.
Subsequently, we evaluated the role of chitin-induced Mg1LysM homodimerization in the protection of fungal hyphae against chitinases by generating three mutant proteins. The T28 residue that makes direct contact with the chitin substrate in the binding site that was previously shown to be essential for chitin binding in the C. fulvum LysM effector Ecp6 was substituted by arginine. In addition, the two residues involved in the formation of the intermolecular salt bridge near the chitin binding site (K31 and D54) were substituted by alanines, respectively. In order to obtain chitin-free proteins, production in E. coli was pursued.
Based on previous findings for Ecp6 [20], we predicted that the T28R mutant was incapable of binding chitin, but the chitin binding capacity of the mutants impaired in the intermolecular salt bridge formation remained to be elucidated [20]. ITC analysis with the mutant T28R revealed a significantly reduced binding affinity of 88.5 μM, which is twenty times weaker when compared with wild-type Mg1LysM protein produced in E. coli (4.36 μM; Fig 3). However, also the binding affinity of the mutants K31A and D54A decreased, to 338.9 μM and 46.3 μM, respectively (Fig 3). Furthermore, besides a lower chitin-binding capacity, the stoichiometry calculated for the mutants K31A and D54A changed from 1:2 as observed for the wild-type Mg1LysM protein to 1:1. This finding implies that a single monomer of K31A or D54A binds a single chitohexaose in solution, whereas a single chitohexaose is bound by two wild-type Mg1LysM protomers, supporting the hypothesis that the chitin-induced dimerization is impaired in K31A and D54A by disruption of the intermolecular salt bridge (Fig 3).
Subsequently, we tested the ability of the Mg1LysM mutants to protect fungal hyphae against the hydrolytic activity of plant chitinases. To this end, spores of Fusarium oxysporum and Trichoderma viride were germinated, incubated with a plant extract containing hydrolytic enzymes including chitinases, and supplemented with wild-type or mutant Mg1LysM. As expected, wild-type Mg1LysM protein produced in E. coli prevented the hydrolysis of hyphae of Fusarium oxysporum f. sp. lycopersici (Fig 5B) and Trichoderma viride (S4 Fig). Furthermore, mutant T28R that is mutated in the substrate-binding loop did not protect F. oxysporum and T. viride hyphae against these hydrolytic enzymes (Fig 5B and S4 Fig), confirming that chitin-binding of Mg1LysM is required to confer protection of cell walls against hydrolysis by plant enzymes. Considering the even lower chitin-binding activity, it is not surprising that also mutant K31A did not protect cell walls against enzymatic hydrolysis. However, also mutant D54A no longer protected cell walls, suggesting that a ten-fold reduction of chitin-binding affinity is sufficient to disrupt the protective activity of Mg1LysM. Unfortunately, based on these findings it is impossible to determine the contribution of the dimerization to the protection activity of Mg1LysM.
We previously determined that LysM effector Ecp6 has two sites that bind chitin with 1:1 stoichiometry, one with ultra-high affinity (kd = 280 pM) and one in the range with which Mg1LysM binds chitin (kd = 1.70 μM) [13,20,32,33,34,35], but both of which bind chitin with higher affinity than Mg1LysM. Intriguingly, Ecp6 fails to protect hyphae against hydrolysis by chitinases [14,20]. Nevertheless, localization experiments making use of constitutive expression of C-terminally GFP-tagged Ecp6 in Verticillium dahliae, and of BODIPY-labelled Ecp6 protein exogenously applied to Botrytis cinerea, two fungal species that expose chitin on their cell walls during growth in vitro, revealed that Ecp6 can bind to fungal cell walls (S5 Fig) in a similar fashion as Cladosporium fulvum effector protein Avr4 that protects fungal cell walls against hydrolysis by chitin binding through an invertebrate chitin-binding domain [18,19]. These findings suggest that binding of a LysM effector to cell wall chitin with high affinity is not sufficient to mediate protection against hydrolytic enzymes. Moreover, from these observations we infer that chitin-induced dimerization of Mg1LysM may be crucial for hyphal protection against plant enzymatic hydrolysis.
Considering that Mg1LysM homodimers possess two chitin-binding sites on opposite sides of the protein complex (Fig 1), combined with the observed chitin-induced dimerization that may be responsible for the protective activity, we hypothesized that Mg1LysM will form highly oligomeric super-complexes in which ligand-independent Mg1LysM homodimers dimerize on both ends in a chitin-dependent manner (Fig 6A). Moreover, we hypothesized that LysM effectors that do not protect hyphae against chitinase hydrolysis would not display such oligomerisation. To test these hypotheses, we first assessed whether we could alter the particle size distribution of Mg1LysM in solutions by adding chitohexaose. Using dynamic light scattering (DLS) we observed that, upon chitin addition at a molar ratio of 1:2 the radius distribution of Mg1LysM particles shifted from around 10 nm in the absence of chitin to 100 nm in the presence of chitin. Moreover, further increase of the chitin concentration to a 1:5 ratio induced a strong signal appearing at 100 μM, demonstrating clear ligand-induced polymerisation of Mg1LysM protein to large protein complexes. Next, we assessed the effect of chitohexaose on the distribution of Ecp6 particles in solution. Interestingly, although the addition of chitohexaose smoothened the Ecp6 particle size distribution, suggesting the stabilization of Ecp6 molecules, chitin addition did not lead to an increased particle size. Thus, in contrast to Mg1LysM, Ecp6 does not display chitin-induced polymerization.
To confirm our observations with respect to the chitin-induced polymerization of Mg1LysM, we pursued an independent validation. To this end, we reasoned that if Mg1LysM indeed polymerizes in the presence of chitin, we should ultimately be able to precipitate such polymeric complexes during centrifugation into a pellet, whereas such a pellet cannot occur in case polymerization does not take place. Thus, we incubated Mg1LysM overnight with chitohexaose and subsequently centrifuged the sample at 20,000 g in the presence of 0.002% methylene blue to visualize the protein. A clear pellet appeared when Mg1LysM was incubated with chitin, but not in the control treatment where no chitin was added to the Mg1LysM protein, nor in the control treatment with chitohexaose only in the absence of Mg1LysM (Fig 7). Collectively, these data confirm the occurrence of chitin-induced Mg1LysM effector polymers. In contrast, similar treatments of Ecp6 with chitin did not result in a pellet after centrifugation, confirming that Ecp6 does not undergo chitin-induced polymerization (Fig 7).
Recently, it was demonstrated that the arbuscular mycorrhizal (AM) fungus Rhizophagus irregularis secretes the LysM effector RiSLM to facilitate its mutualistic symbiosis [39]. Interestingly, like Mg1LysM, RiSLM was shown to be composed of a single LysM only, and furthermore to protect hyphae against hydrolysis by chitinases. Thus, we used RiSLM to test whether chitin-induced polymerisation is restricted to Mg1LysM only, or similarly occurs for RiSLM as well. Interestingly, like with Mg1LysM, a chitin-induced particle size shift was observed in the DLS assay with RiSLM (Fig 6B). Moreover, upon overnight incubation of RiSLM with chitohexaose a clear pellet could be obtained after centrifugation, demonstrating that polymerization occurred (Fig 7). Collectively, these data suggests that chitin-induced polymerization is a common phenomenon that occurs not only with Mg1LysM, but also with other LysM effectors that protect fungal hyphae.
Discussion
Studies on many plant pathogenic fungi have shown that the perception of microbial cell wall-derived glycans by plant hosts plays a central role in microbe–host interactions [10]. Among these glycans, fungal cell wall chitin has emerged as one of the most potent fungal elicitors of host immune responses [3,10]. The widespread glycan perception capacity in plants has spurred the evolution of various fungal strategies to evade glycan perception [3,10]. Many fungal pathogens secrete LysM effectors to perturb the induction of chitin-triggered immunity. Structural analysis of the C. fulvum LysM effector Ecp6 has revealed that this activity could be attributed to the presence of an ultra-high chitin binding affinity site in the LysM effector protein that is established by intramolecular LysM domain dimerization [20]. However, some LysM effectors rather, or additionally, are able to prevent the hydrolysis of fungal cell wall chitin by plant chitinases. Moreover, functional characterization of Mg1LysM, a LysM effector that is merely composed of a single LysM domain and a few additional amino acids, suggested that the ability to protect cell walls is conferred simply by the chitin-binding ability of the LysM domain [15]. Yet, the observation that various other LysM effectors, including C. fulvum Ecp6, M. oryzae Slp1 and C. higginsianum ELP1 and ELP2, are not able to protect hyphae challenged this hypothesis [14,29,30]. Thus, the mechanism by which some LysM effectors are able to protect fungal cell walls remained to be characterized. The crystal structure model that was generated in this study has revealed that Mg1LysM is able to undergo chitin-mediated dimerization such that a chitin molecule is deeply buried in the protein dimer. Nevertheless, the structure of the dimer allows a chitin chain to protrude into the solvent on either side of the binding groove, such that it is conceivable that the dimer can be formed on long-chain chitin polymers of any length, including polymeric cell wall chitin. In addition to several noncovalent bounds between the two Mg1LysM protomers and between the protein and the chitin, a salt bridge between the two Mg1LysM protomers strengthens the chitin-binding affinity and thus supports the chitin-induced dimerization by stabilizing the chitin binding groove. Arguably, it is this particular trait that confers the ability to protect hyphae against plant chitinases, as disruption of the ion bond in the K31A mutant abolished hyphal protection and chitin binding by itself is not sufficient to confer cell wall protection. The crystal structure further revealed that Mg1LysM undergoes ligand-independent homodimerization whereby a large dimerization interface of two Mg1LysM monomers is stabilized by several noncovalent bounds and further strengthened by two salt bridges between the interlaced N-terminal regions of the protein. Despite various efforts, we have not been able to produce monomeric Mg1LysM protein, suggesting that ligand-independent homodimerization is required for proper folding of the protein. Consequently, Mg1LysM homodimers are released and possess two chitin-binding sites on opposite sides of the protein (Fig 1). Combined with the observed chitin-induced dimerization, we postulate that this provides Mg1LysM with the ability to form highly oligomeric super-complexes in the fungal cell wall, in which ligand-independent Mg1LysM homodimers dimerize on both ends in a chitin-dependent manner, leading to the formation of a contiguous structure throughout the cell wall (Fig 6A). Indeed, a chitin-inducible shift in particle size could be demonstrated for Mg1LysM in DLS experiments, confirming the occurrence of polymerisation, whereas such shift was not observed for Ecp6. Moreover, overnight incubation of Mg1LysM with chitohexaose led to the formation of polymers that could be pelleted during centrifugation, demonstrating that the polymers grew to relatively large particles. Also this phenomenon was not observed for Ecp6. Intriguingly, the recently characterized LysM effector RiSLM from the mutualistic fungus Rhizophagus irregularis [39] that, like Mg1LysM, is composed of a single LysM only and protects hyphae against chitinases also showed the ability to polymerize in a chitin-dependent manner, suggesting that polymerization of into oligomeric super-complexes in the fungal cell wall is a common phenomenon among LysM effectors that protect fungal hyphae. Possibly, it is such contiguous structure that provides steric hindrance that renders fungal cell wall chitin inaccessible to chitinase enzymes. Thus, both ligand-independent homodimerization as well as ligand-induced dimerization of Mg1LysM appear to be required for its cell wall protective function. Accordingly, residues shaping these regions are fully conserved in all Mg1LysM isoforms that have been identified to date.
Methods
Protein production and purification
Mg1LysM was previously produced in Pichia pastoris strain GS115 as described [15]. Purification was performed by gel filtration chromatography (Superdex 75: GE Healthcare, Chicago, IL, US) in 20 mM HEPES, pH 7.0, and 50 mM NaCl. P. pastoris produced protein (6–10 mg/mL) was used for protein crystallization. Escherichia coli protein production was performed using pET-SUMO (Thermo Fisher, Waltham, MA, USA). Mg1LysM was cloned into pGEM-T (Promega, Madison, WI, US) using specific primers (S2 Table). Mutants were obtained by PCR using overlapping primers with the corresponding mismatch (S2 Table) followed by digestion of the template by DpnI. For the cloning of ST99_CH3F4 and ST99_CH1E4 versions of Mg1LysM, the sequence was commercially synthesized (Eurofins Genomics, Ebersberg, Germany). After digestion of the resulting vector with SacI and HindIII (Promega, Madison, WI, US), Mg1LysM was cloned in the final vector pET-SUMO. The expression system E. coli ORIGAMI (DE3, Merck, Darmstadt, Germany) cells was used to express the protein following the manufacturer´s instructions. Transformants were selected and grown in Luria broth (LB) medium until an optical density of 0.8 at 600 nm was reached. Protein production was induced with the addition of 0.05 mM Isopropyl ß-D-1-thiogalactopyranoside (IPTG) at 28°C. Cells were harvested by centrifugation ~20 h after induction, the cell pellets were dissolved and lysed using lysozyme from chicken egg (Sigma-Aldrich, St. Louis, MO, US) and the 6xHis-SUMO tagged proteins were purified from the soluble protein fraction after centrifugation using an Ni2+-NTA Superflow column (Qiagen, Venlo, Netherlands). Next, purified proteins were incubated with the SUMO protease ULP1 from Saccharomyces cerevisiae (Thermo Fisher, Waltham, MA, USA), dialysed over night against 200 mM NaCl at 4°C, and again passed through a Ni2+-NTA Superflow column. Native proteins were finally dialysed against 50 mM NaH2PO4, 300 mM NaCl at pH 8.0, concentrated to 0.6 mg/mL over Amicon ultracentrifugal filter units (Sigma-Aldrich, St. Louis, MO, USA) and used for subsequent assays. RiSLM was produced as described previously [39].
Crystallization conditions and structure determination
First crystal hits with 1,4-dioxane as the reservoir solution were obtained overnight in a small initial vapor-diffusion crystallization screening campaign using the Phoenix robot (Art Robbins Instrument LLC, Sunnyvale, CA, USA) with 96-well Intelli Plates (Dunn Labortechnik GmbH, Asbach, Germany) and several different commercial screens (Hampton Research, Aliso Viejo, CA, USA; Molecular Dimensions, Newmarket, Suffolk, UK) [40]. Conditions were further optimized and useful crystals were finally obtained by micro-seeding techniques using 0.1 M sodium citrate pH 5.6, 5%-20% PEG4000 and 5% isopropanol as the reservoir solution [31]. 0.2 M sodium acetate pH 4.6 with 20% ethylene glycol was used as the crystal cryo-buffer. Several crystals were soaked with either I3C (Jena Bioscience GmbH, Jena, Germany), 2 mM in cryo-buffer, quick soak, or Ta6Br14 (Jena Bioscience GmbH, Jena, Germany), 1 mM in cryo-buffer, 1 hr soak, brief wash and prolonged back soak. X-ray diffraction data were collected on BL14.1 at the BESSY II electron storage ring operated by the Helmholz-Zentrum Berlin [41]. Using the Phenix AutoSol wizard [42], initial phases were obtained from the Ta6Br122- derivatized crystals by single-wavelength anomalous dispersion techniques (SAD) that were improved by phase information from the I3C derivatized crystals by single isomorphous replacement with anomalous scattering (SIRAS).
The structure was refined using REFMAC5 [43] and phenix [42] and manually built using Coot [44]. All figures showing structural representations were prepared with the program PyMOL (The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC, DeLano Scientific, Palo Alto, CA, USA]. The quality of the final model was validated with MolProbity [45]. Refinement and phasing statistics are summarized in Table 1.
Chitinase-protection assay
In-vitro chitinase protection assays were performed as described previously [21]. Essentially, ~103 conidiospores of Fusarium oxysporum f. sp. lycopersici or Trichoderma viride were incubated overnight at room temperature in 40 μL of half-strength potato dextrose broth (PDB; Becton Dickinson, Franklin Lakes, NJ, USA) in a 96-well microtiter plate. Subsequently, wild-type or mutant Mg1LysM protein was added at a final concentration of 20 μM. After a 2 h incubation period, 10 μL of tomato extract containing hydrolytic enzymes was added [21]. Fungal growth was assessed microscopically after 4 h of incubation at room temperature.
Polysaccharide precipitation assay
The polysaccharide precipitation assay was performed as described [15]. 800 μL of Mg1LysM (30 μg/mL) was incubated with 50 μL of chitin beads (NEB, Massachusetts, USA), 10 mg shrimp chitin, chitosan, cellulose or xylan (all from Sigma-ALRICH, Missouri, USA) for 6 h at 4°C. The insoluble fraction was pelleted by centrifugation (13,000 g, 5 min), resuspended in 100 μL of water and boiled at 95°C for 10 minutes. The supernatant was concentrated using Microcon Ultracel YM-10 tubes (Merck, Darmstadt, Germany) to 80 μL, boiled at 95°C for 10 minutes with 30 μL of protein loading buffer (4×). The presence of proteins in pellet and supernatant was examined on a Mini-PROTEAN TGX Stain-Free Gel (Bio-Rad, California, USA) followed by Coomassie Brilliant Blue staining.
Isothermal titration calorimetry
Isothermal titration calorimetry (ITC) experiments were performed at 20°C following standard procedures using a Microcal VP-ITC calorimeter (GE Healthcare, Chicago, IL, US). The E. coli-produced wild-type Mg1LysM (20 μM) and the mutants T28R (15 μM), K31A (30 μM) and D54A (30 μM) were titrated with a single injection of 2 μL, followed by 26 injections of 10 μL of (GlcNAc)6 (Isosep AB, Tullinge, Sweden) at 200 μM. Ecp6 (15 μM) was titrated with (GlcNAc)6 at 400 μM. Before the experiment, all proteins were dialyzed against 20 mM of sodium chloride, pH 7.0. Chitohexaose (Megazyme, Wicklow, Ireland) was freshly dissolved in the dialysis buffer. Data were analyzed using Origin (OriginLab, Northampton, MA, USA) and fitted to a one-binding-site model. Before and after the experiment protein samples were analysed on SDS-PAGE gel and stained with Coomassie Brilliant Blue (S3C Fig).
Dynamic light scattering (DLS) measurements
Mg1LysM, RiSLM and Ecp6 were dialyzed overnight against water, and subsequently incubated with 0.01% Triton X-100 for 4 hours to improve protein solubility. Next, chitohexaose (Megazyme, Wicklow, Ireland) was added in a molar ratio of 1:0, 1:2 and 1:5 (protein:chitin) and incubated overnight. Particle size distribution was measured by a SpectroSize 300 (Xtal Concepts, Hamburg, Germany).
Polymerization assay
Mg1LysM and Ecp6 were adjusted to a concentration of 400 μM and 100 μL of each protein was incubated with 100 μL of 4 mM chitohexaose (Megazyme, Wicklow, Ireland), or 100 μL water as control, at room temperature overnight. Similarly, 100 μL of RiSLM (600 μM) was incubated with 75 μL of 4 mM chitohexaose (Megazyme, Wicklow, Ireland) or water as control in a total volume of 200 μL. The next day, 2 μL of 0.2% methylene blue (Sigma-Aldrich, Missouri, USA) was added and incubated for 10 min after which protein solutions were centrifuged at 20,000 g for 15 min. Photos were taken with a ChemiDoc MP system (Bio-Rad, California, USA) with custom setting for red fluorescent protein (RFP).
Localisation of Mg1LysM and Ecp6
Labelling of effector proteins with BODIPY TMR-X amine reactive probe (Invitrogen, Carlsbad, CA, USA) was performed as described previously [18,19]. For localisation assay of Mg1LysM, conidiospores of Trichoderma viride were harvested from five-day-old potato dextrose agar (PDA) plates and adjusted to 106 conidiospores/mL with half-strength PDB. The conidiospore solution was pipetted into a 96-well microtiter plate in aliquots of 50 μL and the plate was incubated overnight at room temperature for germination. The next day, BODIPY-labeled Mg1LysM was applied at a final concentration of 8 μM and incubated for 4 hrs at room temperature in the dark. Microscopic analysis was performed using a Nikon Eclipse Ti microscope using a 100× Plan apo oil immersion objective (NA 1.4) and a 561 nm laserline. Pictures were processed and analysed with ImageJ (http://rsbweb.nih.gov/ij/).
For localisation study of Ecp6, conidiospores of a V. dahliae transformant were grown in a few micro liters of PDB on a glass slide with coverslip. To prevent the samples from drying out, the slides were mounted on top of moistened tissue in an empty pipette box with water on the bottom. After approximately 6 hr of growth at room temperature, the slides were used for localization studies. Conidiospores of Botrytis cinerea were harvested and germinated overnight in PDB at room temperature. BODIPY-labeled proteins were applied at a concentration of 4 μM and incubated for 2–3 hrs at room temperature in the dark. The localisation studies were performed using a Nikon eclipse 90i UV microscope and NIS-Elements AR 2.3 software (Nikon Instruments Inc., Melville, USA).
Assembly and alignment of Mg1LysM sequences
Illumina whole-genome sequencing data from a global collection of Z. tritici isolates was used to extract Mg1LysM sequences [46]. We used the SPAdes assembler version 3.6.2 [47] to generate de-novo genome assemblies. The SPAdes pipeline includes the BayesHammer read error correction module to build contigs in a stepwise procedure based on increasing k-mer lengths. We defined the k-mer range as 21, 35, 49, 63 and 77. We used the “—careful” option to reduce mismatches and indel errors in the assembly. Polished assemblies were then used to retrieve the contigs containing Mg1LysM orthologs based on blastn [48]. High-confidence sequence matches were extracted with samtools [49] from each draft assembly and aligned using MAFFT version 7.305b [50] using the—auto option and 1,000 iterative refinement cycles. Alignments were processed using JalView [51] and CLC Genomic Workbench 9 (Qiagen, Venlo, Netherlands).
Supporting information
S1 Fig [a]
Protein alignment of Ecp6 and Mg1LysM.
S2 Fig [tif]
2|F|–|F| map.
S4 Fig [he]
Mg1LysM mutants are impaired in protection of against chitinases.
S5 Fig [a]
. LysM effector Ecp6 protein localizes to fungal cell walls.
S1 Table [xlsx]
Mg1LysM protein alignment.
S2 Table [docx]
Primers used in this study.
Zdroje
1. Shibuya N, Kaku H, Kuchitsu K and Maliarik MJ. (1993) Identification of a novel high-affinity binding site for N-acetylchitooligosaccharide elicitor in the membrane fraction from suspension-cultured rice cells. FEBS Lett. 329: 75–78. doi: 10.1016/0014-5793(93)80197-3 8354412
2. Boller T. (2003) Chemoperception of microbial signals in plant cells. Annu. Rev. Plant Physiol. 46: 189–214.
3. Sanchez-Vallet A, Mesters JR and Thomma BPHJ. (2015) The battle for chitin recognition in plant-microbe interactions. FEMS Microbiol Rev. 39: 171–183. doi: 10.1093/femsre/fuu003 25725011
4. Schlumbaum A, Mauch F, Vögeli U and Boller T. (1986) Plant chitinases are potent inhibitors of fungal growth. Nature. 324: 365–367.
5. van Loon LC, Rep M and Pieterse CM. (2006) Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 44: 135–162. doi: 10.1146/annurev.phyto.44.070505.143425 16602946
6. Adrangi S and Faramarzi MA. (2013) From bacteria to human: a journey into the world of chitinases. Biotechnol Adv. 31: 1786–1795. doi: 10.1016/j.biotechadv.2013.09.012 24095741
7. Hamid R, Khan MA, Ahmad M, Ahmad MM, Abdin MZ, Musarrat J, et al. (2013) Chitinases: An update. J Pharm Bioallied Sci. 5: 21–29. doi: 10.4103/0975-7406.106559 23559820
8. Kasprzewska A. (2003) Plant chitinases—regulation and function. Cell Mol Biol Lett. 8: 809–824. 12949620
9. Liu T, Chen L, Ma Q, Shen X and Yang Q. (2014) Structural insights into chitinolytic enzymes and inhibition mechanisms of selective inhibitors. Curr Pharm Des. 20: 754–770. doi: 10.2174/138161282005140214164730 23688083
10. Rovenich H, Zuccaro A and Thomma BP. (2016) Convergent evolution of filamentous microbes towards evasion of glycan-triggered immunity. New Phytol. 212: 896–901. doi: 10.1111/nph.14064 27329426
11. Kaku H, Nishizawa Y, Ishii-Minami N, Akimoto-Tomiyama C, Dohmae N, Takio K, et al. (2006) Plant cells recognize chitin fragments for defense signaling through a plasma membrane receptor. Proc Natl Acad Sci USA. 103: 11086–11091. doi: 10.1073/pnas.0508882103 16829581
12. Miya A, Albert P, Shinya T, Desaki Y, Ichimura K, Shirasu K, et al. (2007) CERK1, a LysM receptor kinase, is essential for chitin elicitor signaling in Arabidopsis. Proc Natl Acad Sci USA. 104: 19613–19618. doi: 10.1073/pnas.0705147104 18042724
13. Liu T, Liu Z, Song C, Hu Y, Han Z, She J, et al. (2012) Chitin-induced dimerization activates a plant immune receptor. Science. 336: 1160–1164. doi: 10.1126/science.1218867 22654057
14. de Jonge R, van Esse HP, Kombrink A, Shinya T, Desaki Y, Bours R, et al. (2010) Conserved fungal LysM effector Ecp6 prevents chitin-triggered immunity in plants. Science. 329: 953–955. doi: 10.1126/science.1190859 20724636
15. Marshall R, Kombrink A, Motteram J, Loza-Reyes E, Lucas J, Hammond-Kosack KE, et al. (2011) Analysis of two in planta expressed LysM effector homologs from the fungus Mycosphaerella graminicola reveals novel functional properties and varying contributions to virulence on wheat. Plant Physiol. 156: 756–769. doi: 10.1104/pp.111.176347 21467214
16. Okmen B, Kemmerich B, Hilbig D, Wemhoner R, Aschenbroich J, Perrar A, et al. (2018) Dual function of a secreted fungalysin metalloprotease in Ustilago maydis. New Phytol. 220: 249–261. doi: 10.1111/nph.15265 29916208
17. de Wit PJ. (2016) Cladosporium fulvum effectors: Weapons in the arms race with tomato. Annu Rev Phytopathol. 54: 1–23. doi: 10.1146/annurev-phyto-011516-040249 27215970
18. van den Burg HA, Harrison SJ, Joosten MH, Vervoort J and de Wit PJ. (2006) Cladosporium fulvum Avr4 protects fungal cell walls against hydrolysis by plant chitinases accumulating during infection. Mol Plant Microbe Interact. 19: 1420–1430. doi: 10.1094/MPMI-19-1420 17153926
19. van Esse HP, Bolton MD, Stergiopoulos I, de Wit PJ and Thomma BPHJ. (2007) The chitin-binding Cladosporium fulvum effector protein Avr4 is a virulence factor. Mol Plant Microbe Interact. 20: 1092–1101. doi: 10.1094/MPMI-20-9-1092 17849712
20. Sanchez-Vallet A, Saleem-Batcha R, Kombrink A, Hansen G, Valkenburg DJ, Thomma BPHJ, et al. (2013) Fungal effector Ecp6 outcompetes host immune receptor for chitin binding through intrachain LysM dimerization. Elife. 2: e00790. doi: 10.7554/eLife.00790 23840930
21. van den Burg HA, Spronk CA, Boeren S, Kennedy MA, Vissers JP, Vuister GW, et al. (2004) Binding of the AVR4 elicitor of Cladosporium fulvum to chitotriose units is facilitated by positive allosteric protein-protein interactions: the chitin-binding site of AVR4 represents a novel binding site on the folding scaffold shared between the invertebrate and the plant chitin-binding domain. J Biol Chem. 279: 16786–16796. doi: 10.1074/jbc.M312594200 14769793
22. Kohler AC, Chen LH, Hurlburt N, Salvucci A, Schwessinger B, Fisher AJ, et al. (2016) Structural analysis of an Avr4 effector ortholog offers insight into chitin binding and recognition by the Cf-4 receptor. Plant Cell. 28: 1945–1965. doi: 10.1105/tpc.15.00893 27401545
23. Hurlburt NK, Chen L-H, Stergiopoulos I and Fisher AJ. (2018) Structure of the Cladosporium fulvum Avr4 effector in complex with (GlcNAc) 6 reveals the ligand-binding mechanism and uncouples its intrinsic function from recognition by the Cf-4 resistance protein. PLoS Pathog. 14: e1007263. doi: 10.1371/journal.ppat.1007263 30148881
24. Stergiopoulos I, van den Burg HA, Okmen B, Beenen HG, van Liere S, Kema GH, et al. (2010) Tomato Cf resistance proteins mediate recognition of cognate homologous effectors from fungi pathogenic on dicots and monocots. Proc Natl Acad Sci USA. 107: 7610–7615. doi: 10.1073/pnas.1002910107 20368413
25. Bolton MD, van Esse HP, Vossen JH, de Jonge R, Stergiopoulos I, Stulemeijer IJ, et al. (2008) The novel Cladosporium fulvum lysin motif effector Ecp6 is a virulence factor with orthologues in other fungal species. Mol Microbiol. 69: 119–136. doi: 10.1111/j.1365-2958.2008.06270.x 18452583
26. de Jonge R and Thomma BPHJ. (2009) Fungal LysM effectors: extinguishers of host immunity? Trends Microbiol. 17: 151–157. doi: 10.1016/j.tim.2009.01.002 19299132
27. Kombrink A and Thomma BPHJ. (2013) LysM effectors: secreted proteins supporting fungal life. PLoS pathog. 9: e1003769. doi: 10.1371/journal.ppat.1003769 24348247
28. Kombrink A, Rovenich H, Shi-Kunne X, Rojas-Padilla E, van den Berg GC, Domazakis E, et al. (2017) Verticillium dahliae LysM effectors differentially contribute to virulence on plant hosts. Mol Plant Pathol. 18: 596–608. doi: 10.1111/mpp.12520 27911046
29. Mentlak TA, Kombrink A, Shinya T, Ryder LS, Otomo I, Saitoh H, et al. (2012) Effector-mediated suppression of chitin-triggered immunity by Magnaporthe oryzae is necessary for rice blast disease. Plant Cell. 24: 322–335. doi: 10.1105/tpc.111.092957 22267486
30. Takahara H, Hacquard S, Kombrink A, Hughes HB, Halder V, Robin GP, et al. (2016) Colletotrichum higginsianum extracellular LysM proteins play dual roles in appressorial function and suppression of chitin-triggered plant immunity. New Phytol. 211: 1323–1337. doi: 10.1111/nph.13994 27174033
31. Bergfors T. (2003) Seeds to crystals. J. Struct. Biol. 142: 66–76. doi: 10.1016/s1047-8477(03)00039-x 12718920
32. Bateman A and Bycroft M. (2000) The structure of a LysM domain from E. coli membrane-bound lytic murein transglycosylase D (MltD). J Mol Biol. 299: 1113–1119. doi: 10.1006/jmbi.2000.3778 10843862
33. Bielnicki J, Devedjiev Y, Derewenda U, Dauter Z, Joachimiak A and Derewenda ZS. (2006) B. subtilis ykuD protein at 2.0 A resolution: insights into the structure and function of a novel, ubiquitous family of bacterial enzymes. Proteins. 62: 144–151. doi: 10.1002/prot.20702 16287140
34. Koharudin LM, Viscomi AR, Montanini B, Kershaw MJ, Talbot NJ, Ottonello S, et al. (2011) Structure-function analysis of a CVNH-LysM lectin expressed during plant infection by the rice blast fungus Magnaporthe oryzae. Structure. 19: 662–674. doi: 10.1016/j.str.2011.03.004 21565701
35. Bozsoki Z, Cheng J, Feng F, Gysel K, Vinther M, Andersen KR, et al. (2017) Receptor-mediated chitin perception in legume roots is functionally separable from Nod factor perception. Proc Natl Acad Sci USA. 114: E8118–E8127. doi: 10.1073/pnas.1706795114 28874587
36. Krissinel E and Henrick K. (2007) Inference of macromolecular assemblies from crystalline state. J Mol Biol. 372: 774–797. doi: 10.1016/j.jmb.2007.05.022 17681537
37. Winn MD, Ballard CC, Cowtan KD, Dodson EJ, Emsley P, Evans PR, et al. (2011) Overview of the CCP4 suite and current developments. Acta Cryst. D. 67: 235–242.
38. Hartmann FE, Sanchez-Vallet A, McDonald BA and Croll D. (2017) A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements. ISME J. 11: 1189–1204. doi: 10.1038/ismej.2016.196 28117833
39. Zeng T, Rodriguez-Moreno L, Mansurkhodzaev A, Wang P, van den Berg W, Gasciolli V, et al. (2020) A lysin motif effector subverts chitin-triggered immunity to facilitate arbuscular mycorrhizal symbiosis. New Phytol. 225: 448–460. doi: 10.1111/nph.16245 31596956
40. Newman J, Egan D, Walter TS, Meged R, Berry I, Ben Jelloul M, et al. (2005) Towards rationalization of crystallization screening for small- to medium-sized academic laboratories: the PACT/JCSG+ strategy. Acta Cryst. D. 61: 1426–1431.
41. Gerlach M, Mueller U and Weiss MS. (2016) The MX Beamlines BL14.1–3 at BESSY II. Journal of large-scale research facilities JLSRF. 2.
42. Adams PD, Afonine PV, Bunkoczi G, Chen VB, Davis IW, Echols N, et al. (2010) PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Cryst. D. 66: 213–221.
43. Murshudov GN, Skubak P, Lebedev AA, Pannu NS, Steiner RA, Nicholls RA, et al. (2011) REFMAC5 for the refinement of macromolecular crystal structures. Acta Cryst. D. 67: 355–367.
44. Emsley P, Lohkamp B, Scott WG and Cowtan K. (2010) Features and development of Coot. Acta Cryst. D. 66: 486–501.
45. Chen VB, Arendall WB 3rd, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, et al. (2010) MolProbity: all-atom structure validation for macromolecular crystallography. Acta Cryst. D. 66: 12–21.
46. Hartmann FE and Croll D. (2017) Distinct trajectories of massive recent gene gains and losses in populations of a microbial eukaryotic pathogen. Mol Biol Evol. 34: 2808–2822. doi: 10.1093/molbev/msx208 28981698
47. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. (2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 19: 455–477. doi: 10.1089/cmb.2012.0021 22506599
48. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. (2009) BLAST+: architecture and applications. BMC Bioinformatics. 10: 421. doi: 10.1186/1471-2105-10-421 20003500
49. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25: 2078–2079. doi: 10.1093/bioinformatics/btp352 19505943
50. Katoh K and Standley DM. (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 30: 772–780. doi: 10.1093/molbev/mst010 23329690
51. Waterhouse AM, Procter JB, Martin DM, Clamp M and Barton GJ. (2009) Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 25: 1189–1191. doi: 10.1093/bioinformatics/btp033 19151095
Článek vyšel v časopise
PLOS Pathogens
2020 Číslo 6
- Antibiotika na nachlazení nezabírají! Jak můžeme zpomalit šíření rezistence?
- Prof. Jan Škrha: Metformin je bezpečný, ale je třeba jej bezpečně užívat a léčbu kontrolovat
- FDA varuje před selfmonitoringem cukru pomocí chytrých hodinek. Jak je to v Česku?
- Jak a kdy u celiakie začíná reakce na lepek? Možnou odpověď poodkryla čerstvá kanadská studie
- Ibuprofen jako alternativa antibiotik při léčbě infekcí močových cest
Nejčtenější v tomto čísle
- Exploring potential of vaginal Lactobacillus isolates from South African women for enhancing treatment for bacterial vaginosis
- Microbiome factors in HPV-driven carcinogenesis and cancers
- Biological sex impacts COVID-19 outcomes
- Bacterial killing by complement requires direct anchoring of membrane attack complex precursor C5b-7