Cloning and characterization of a pectin lyase gene from Colletotrichum lindemuthianumand comparative phylogenetic/structural analyses with genes from phytopathogenic and saprophytic/opportunistic microorganisms
© Lara-Márquez et al; licensee BioMed Central Ltd. 2011
Received: 3 June 2011
Accepted: 9 December 2011
Published: 9 December 2011
Microorganisms produce cell-wall-degrading enzymes as part of their strategies for plant invasion/nutrition. Among these, pectin lyases (PNLs) catalyze the depolymerization of esterified pectin by a β-elimination mechanism. PNLs are grouped together with pectate lyases (PL) in Family 1 of the polysaccharide lyases, as they share a conserved structure in a parallel β-helix. The best-characterized fungal pectin lyases are obtained from saprophytic/opportunistic fungi in the genera Aspergillus and Penicillium and from some pathogens such as Colletotrichum gloeosporioides.
The organism used in the present study, Colletotrichum lindemuthianum, is a phytopathogenic fungus that can be subdivided into different physiological races with different capacities to infect its host, Phaseolus vulgaris. These include the non-pathogenic and pathogenic strains known as races 0 and 1472, respectively.
Here we report the isolation and sequence analysis of the Clpnl2 gene, which encodes the pectin lyase 2 of C. lindemuthianum, and its expression in pathogenic and non-pathogenic races of C. lindemuthianum grown on different carbon sources. In addition, we performed a phylogenetic analysis of the deduced amino acid sequence of Clpnl2 based on reported sequences of PNLs from other sources and compared the three-dimensional structure of Clpnl2, as predicted by homology modeling, with those of other organisms. Both analyses revealed an early separation of bacterial pectin lyases from those found in fungi and oomycetes. Furthermore, two groups could be distinguished among the enzymes from fungi and oomycetes: one comprising enzymes from mostly saprophytic/opportunistic fungi and the other formed mainly by enzymes from pathogenic fungi and oomycetes. Clpnl2 was found in the latter group and was grouped together with the pectin lyase from C. gloeosporioides.
The Clpnl2 gene of C. lindemuthianum shares the characteristic elements of genes coding for pectin lyases. A time-course analysis revealed significant differences between the two fungal races in terms of the expression of Clpnl2 encoding for pectin lyase 2. According to the results, pectin lyases from bacteria and fungi separated early during evolution. Likewise, the enzymes from fungi and oomycetes diverged in accordance with their differing lifestyles. It is possible that the diversity and nature of the assimilatory carbon substrates processed by these organisms played a determinant role in this phenomenon.
Pectin is one of the major components of the primary cell wall of plants and is also found in dividing cells and in the areas of contact between cells that have a secondary cell wall, including xylem and the fibrous cells of woody tissue. Pectin comprises approximately 35% of the primary cell wall of dicots and non-graminaceous monocots. Although its content in secondary walls is greatly reduced, it is believed that pectin plays an important role in the structure and function of both primary and secondary cell walls. The functions of pectin in cell walls are diverse and include plant growth and development, morphogenesis, defense, cell adhesion, cell wall structure, cellular expansion, porosity, ion binding, hydration of seeds, leaf abscission and fruit development, among others [1, 2]. In general, pectin is considered to be a group of polysaccharides that are rich in galacturonic acid (GalA) and present in the form of covalently linked structural domains: homogalacturonan (HG), xylogalacturonan (XGA), rhamnogalacturonan I (RG-I) and rhamnogalacturonan II (RG-II) [1, 2]. The main enzymes involved in the degradation of the HG backbone of pectin are polygalacturonases (PGA, E.C. 126.96.36.199 and XPG, E.C. 188.8.131.52), pectate lyases (PL, E.C. 184.108.40.206 and 220.127.116.11) and pectin lyases (PNL, E.C. 18.104.22.168) .
Pectin lyases (PNLs) catalyze the degradation of pectin through β-elimination; they remove a proton and generate an unsaturated bond between the C-4 and C-5 carbons of the non-reducing end of pectin, which is a neutral form of pectate in which the uronic acid moiety of galacturonic residues has been methyl-esterified. The activity of PNLs is highly dependent on the distribution of the methyl esters over the homogalacturonan backbone. PNLs exhibit pH optima in the range of 6.0-8.5 and, unlike PLs, their activity is independent of Ca2+ ions; it is believed, however, that the residue Arg236 plays a role similar to that of Ca+2 [4, 5]. Pectinase gene expression is regulated at the transcriptional level by the pH of the medium and by carbon sources, as it is induced by pectin and pectic components and repressed by glucose [6–8].
PNLs are grouped into Family 1 of the polysaccharide lyases  and into the pectate lyase superfamily that, in addition to pectin lyases and pectate lyases, also includes plant pollen/style proteins. The three-dimensional structures of five members of the pectate lyase superfamily have been determined. These include Erwinia chrysanthemi pectate lyase C (PELC)  and pectate lyase E (PELE) , Bacillus subtilis pectate lyase  and Aspergillus niger pectin lyase A (PLA)  and pectin lyase B (PLB) . These enzymes fold into a parallel β-helix, which is a topology in which parallel β-strands are wound into a large right-handed coil . Although PLs and PNLs exhibit a similar structural architecture and related catalysis mechanisms, they nonetheless diverge significantly in their carbohydrate binding strategy [4, 13]. Currently, strategies are available for developing functional information from three-dimensional images of enzymes. The growing number of databases on the structure of pectinolytic enzymes has facilitated the analysis of minor structural differences that are responsible for the specific recognition of a unique oligosaccharide sequence in a heterogeneous mixture .
Most of the available information about fungal PNLs and their corresponding encoding genes has been obtained from saprophytic/opportunistic fungi such as Aspergillus niger [16–19], A. orizae [20, 21], A. fumigatus , Penicillium griseoroseum , P. occitanis  and to a lesser extent from the phytopathogenic fungi Glomerella cingulata  and C. gloeosporioides .
The ascomycete C. lindemuthianum is an economically important phytopathogen, and along with its host Phaseolus vulgaris, it provides a convenient model to study the physiological and molecular bases of plant-pathogen interactions . It is an intracellular hemibiotrophic pathogen with physiological races that invade the plant in an interaction consistent to the gene-for-gene model , and monogenic dominant resistance in common bean cultivars leads to the appearance of localized necrotic spots typical of the hypersensitive response (HR) . After penetration of a host epidermal cell in a susceptible cultivar, the pathogenic races of C. lindemuthianum develop an infection vesicle and extend into adjacent cells by producing large primary hyphae, which invaginate without penetrating the host cell membrane and thus persist as a biotrophic interaction. Once a large area of the plant tissue has been colonized, necrotrophic hyphae develop , and this step closely correlates with the production of a number of host cell-wall-degrading enzymes that are characteristic of phytopathogenic fungi [30–32]. Up to know, race 0 is the only strain of C. lindemuthianum unable to infect P. vulgaris, which contrasts with 1472, one of the most virulent races isolated in México . This difference makes the two races an excellent model to investigate the role of pectinolytic enzymes in virulence of C. lindemuthianum. Previous results from this laboratory revealed significant differences between pathogenic (1472) and non-pathogenic (0) races of C. lindemuthianum in terms of growth and production of extracellular PNL activity on different carbon and nitrogen sources in liquid culture. Accordingly, race 1472 grew faster in media containing glucose or polygalacturonic acid, and on 92%-esterified pectin, it produced levels of PNL activity that were approximately 2-fold higher than those produced by race 0. In contrast, cell walls isolated from P. vulgaris hypocotyls and, to a lesser degree, from cellulose sustained the growth of both races but induced PNL only in the pathogenic race .
Here we report the isolation and sequence analysis of the Clpnl2 gene, which encodes pectin lyase 2 of C. lindemuthianum, and its expression in pathogenic and non-pathogenic races of C. lindemuthianum in response to cultivation on different carbon sources. To determine the relationship among the three-dimensional structures of PNLs and the lifestyle of PNL-producing microorganisms, we performed a phylogenetic analysis using protein sequences and deduced amino acid sequences reported for PNLs. A comparative analysis of the three-dimensional structure of the Clpnl2 protein predicted by homology modeling, covering the main body of the protein and the carbohydrate binding site, and the three-dimensional structures of the PNLs used in the phylogenetic analysis was also performed.
Strain and culture conditions
C. lindemuthianum races 0 (non-pathogenic) and 1472 (pathogenic) were kindly provided by Dr. June Simpson (CINVESTAV-IPN, Unidad Irapuato, México) and maintained on potato dextrose agar (PDA, Difco) at 20°C. For DNA extraction, mycelia from C. lindemuthianum race 1472 grown on potato dextrose (PD) for 9 days at 20°C with continuous shaking (150 rpm), was recovered by filtration through Whatman paper No. 1 and stored at -85°C. For induction, 1.6 mg (about 5 cm2) of mycelia from races 0 and 1472 were inoculated in 250 ml-Erlenmeyer flasks containing 50 ml of PD medium and shaken (150 rpm) at 20°C. After 9 days, mycelia was collected by filtration, washed with water and transferred to 250 ml-Erlenmeyer flasks containing 50 ml of modified Mathur's medium (10 mM MgSO4.7H2O, 20 mM KH2PO4, 36 mM L-glutamic acid, distilled water up to 1 L; final pH, 5.5)  supplemented with either 2.5% glucose, 92%-esterified pectin or cell walls from P. vulgaris. Flasks were shaken (150 rpm) at 20°C and after different periods of growth, mycelia was collected by filtration, washed with water and stored at -85°C until use.
Preparation of plant cell walls
Seedlings of P. vulgaris cv. Flor de Mayo were grown for 7 days, and cell walls were extracted and purified from hypocotyls as described elsewhere .
DNA and RNA isolation
Genomic DNA was isolated from C. lindemuthianum mycelia that had been grown for 9 days in PD medium according to standard protocols . Total RNA was purified from mycelia using TRIzol reagent (Invitrogen). RNA samples were treated with DNAse I according to manufacturer's instructions (Invitrogen) to eliminate DNA. The quality and concentration of total RNA were verified using the RNA 6000 Nano LabChip kit (2100 Agilent Bioanalyzer).
Isolation of the homologous DNA Clpnl2 probe from C. lindemuthianum
Genomic DNA from race 1472 was amplified by PCR using the upstream primer pnlD (5'-CAGTACGTCTGGGGTGGTGA-3') and downstream primer pnlR (5'-AAGTAGTTGTTGACGACGTGG-3', which are homologous to sequences between 595 and 614 nt and 891 and 911 nt, respectively, of exon 3 of the Clpnl2 gene from C. gloeosporioides [GenBank: AAD43565]. The PCR incubation mixture was heated at 95°C for 5 min in a thermocycler (Eppendorf Master Cycler Gradient, Brinkmann, Westbury, NY), followed by denaturation for 1 min at 95°C, annealing for 2 min at 48°C and extension for 2 min at 72°C. PCR was then performed for 35 cycles, followed by a final extension for 10 min at 72°C. A PCR product of 383 bp corresponding to pnl2 gene (clpnl2 fragment) was ligated into the pCR 2.1 vector and introduced into E. coli TOP 10 strain from the TOPO TA Cloning kit (Invitrogen).
Genomic DNA library construction and screening
Partial Sau3AI digestion of genomic DNA from race 1472 was used to construct a genomic library in Lambda DASH II/BamHI according to manufacturer's instructions (Stratagene). Screening was performed using 15 × 104 UFP with three rounds of hybridization filters and the homologous Clpnl2 fragment, which was 32P-radiolabeled using the Radprime DNA Labeling System Life Technologies Kit (Tech-Line).
Molecular cloning of the Clpnl2full-length cDNA and expression analyses
The cDNA was amplified by RT-PCR as specified by the manufacturer. SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen) was used to prepare cDNA from total RNA. PCR was performed using the upstream primer Pnl67 (5'-ATGAAGTCTACCATCTTCTCCG-3') and downstream primer Pnl1569 (5'-TTAGATCTTGCGAAACCGGC-3') designed from the DNA Clpnl2 genomic sequence of C. lindemuthianum. The PCR incubation mixture was heated at 94°C for 5 min in a thermocycler (Eppendorf Master Cycler Gradient, Brinkmann, Westbury, NY), followed by 30 cycles of denaturation for 20 sec at 94°C, annealing for 30 sec at 54°C, extension for 1.5 min at 72°C and then by a final extension for 7 min at 72°C. A PCR product of 1,140 bp obtained from total RNA of race 1472 induced with pectin for 4 h and corresponding to the Clpnl2 gene, was ligated into the pCR 2.1 vector (Invitrogen) and three clones were selected and sequenced. The 5' end of cDNA was amplified by 5'RACE as specified by the manufacturer (5'RACE System for Rapid Amplification of cDNA Ends, Invitrogen), with total RNA from race 1472 induced for 4 h with 92%-esterified pectin, using the specific reverse primers Pnl1249 (5'-GTA GTT GTT GAC GAC GTG GAC G-3') and Pnl975 (5'-CGA TGT GCT GGC GGC CG-3'). The amplification products were cloned and five clones were selected and sequenced. For expression analysis, total cDNA (1140 pb) was amplified with specific primers Pnl67 and Pnl1569 in the same conditions described above using total RNA of mycelia from both races induced with 92%-esterified pectin or cell walls from P. vulgaris for 2, 4, 6, 8, 10 and 12 h. For expression analysis, cDNA obtained from cells grown under different conditions was also amplified by PCR using oligonucleotides prepared from ribosomal 18S RNA as a control (5'- TTAGCATGGAATAATRRAATAGGA-3'and 5'-ATTGCAATGCYCTATCCCCA-3) . The PCR incubation mixture was heated at 94°C for 3 min, followed by 35 cycles of denaturation for 1 min at 94°C, annealing for 1 min at 56°C, extension for 1 min at 72°C and then a final extension for 10 min at 72°C. The cDNA of Clpnl2 and the amplified RT-PCR products were analyzed in a Bioanalyzer using the system for quantification and molecular size Agilent DNA 7500 (2100 Agilent Bioanalyzer).
Southern blot hybridization
Genomic DNA of mycelia from race 1472 was digested with selected restriction endonucleases. Digestion products were size-fractionated on a 0.8% agarose gel, transferred to a nylon membrane (Hybond-N+, Amersham Pharmacia Biotec, England), hybridized and detected with a 32P-radiolabeled Clpnl2 probe. Hybridizations were carried out at 60°C in 2X SSC containing 0.5% blocking agent (Roche) and 0.1% SDS. After hybridization, the blot was washed at 60°C for 15 min with 2X SSC containing 1% SDS and then at 60°C for 15 min with 0.2X SSC containing 0.1% SDS.
Sequencing and DNA analysis
The sequences of both strands of DNA of race 1472 and cDNA of both races were determined by the dideoxy-chain termination method using the ABI Prism Dye Cycle Sequencing Ready Reaction Kit in an ABI PRISM 310 DNA sequencer (Applied Biosystems, Foster City, CA). The nucleotide sequences were analyzed using the DNAsis (Hitachi) and 4Peaks v 1.7.2 software (http://mekentosj.com). In silico analyses of putative transcription factor binding sites were performed using the AliBaba2.1 software  and the Transfac 7.0 database ; the regulatory sequences reported for genes of fungal lytic enzymes were also compared. The N-terminal secretion signal sequence was identified with the SignalP 3.0 web server . The protein molecular mass, pI and N-glycosylation sites were calculated on an ExPASy Proteomics Server .
Nucleotide and protein sequences of reported pectin lyases used for phylogenetic analyses.
GenBank: CAD34589, GenBank: AAW03313, GenBank: CAA39305, GenBank: CAA01023, GenBank: ACE00421, GenBank: AAA32701
GenBank: BAB82468, GenBank: BAB82467
Swiss-Prot: BOYCL3, Swiss-Prot: Q4WV10, GenBank: EAL91586, Swiss-Prot: Q4W156
GenBank: EAU31855, GenBank: EAU37973
GenBank: AAA21817, GenBank: AAD43565, GenBank: AAF22244
GenBank: EAW17753, Swiss-Prot: A1CYC2
GenBank: XP_001934252, GenBank: XP_001930850
GenBank: XP_002909420, GenBank: XP_002903922
GenBank: BAA12119, GenBank: AAB84422
Protein homology modeling
The tertiary structure of the deduced amino acid sequence of Clpnl2 was predicted by homology modeling using the Swiss-Model Server  using Pel B from A. niger (PDB: 1qcxA) as template . The prediction of three-dimensional structures of the deduced amino acid sequences used in the phylogenetic analysis was performed in a similar manner. The structural parameters and prediction quality of the modeled structures were evaluated using the program SPDBV v. 4.01 . The energy minimization of the model was performed by GROMOS96 , which was provided by the SPDBV program. MMV 2010.2.0.0 (Molegro ApS) and SPDBV v. 4.01 were used for visualization of molecular structures.
Multiple comparisons of protein structures
The comparison of protein structures was performed using the Voronoi contact method  with the ProCKSI-Server . Calculations were performed using default parameters, and the resultant similarity matrixes (Voronoi-contacts) were standardized and used as the input for clustering of the protein set using the un-weighted pair group method for the arithmetic mean (UPGMA) .
Results and discussion
Isolation and sequence analysis of the Clpnl2gene
Nine positive clones were isolated from the screening of a C. lindemuthianum genomic library using the 32P-radiolabeled fragment of Clpnl2. Southern blot analysis of the clones allowed the identification of a 4.0-kb fragment that hybridized with the PCR probe. The 4.0-kb fragment was subcloned, and 2,159 bp containing the Clpnl2 gene was sequenced [GenBank: JN034038]. The full-length Clpnl2 cDNA of races 1472 [GenBank: JN034039] and 0 [GenBank: JN653459] obtained by RT-PCR from total RNA of mycelium induced for 4 h with 92%-esterified pectin, was sequenced using specific primers designed from the Clpnl2 gene sequence.
The Clpnl2 cDNA contains an ORF of 1140 nucleotides that encodes a putative protein of 379 aa with a N-terminal secretion signal sequence of 19 amino acids, according to the SignalP 3.0 web server . A protein of molecular mass 37.4 kDa and a pI of 9.1 was calculated, and one potential N-glycosylation site was located at position 110 (ExPASy Proteomics Server) . These results are consistent with those reported for the amino acid sequences of Pnl2 and pnlA of C. gloeosporioides [25, 26]. Despite N-glycosylation is common in pectinolytic enzymes and has been reported in several fungal pectin lyases al similar positions, little is known about the function of this posttranslational modification. Although it is believed that it affect enzyme stability and activity [60, 61].
Southern blot analysis
Protein homology modeling
Bayesian analysis allowed the separation of pectin lyases into two groups: one representing bacteria with 100% posterior probability and 100% bootstrap support for MP and NJ analysis, and the other one representing fungi and oomycetes with 100% posterior probability and 98% bootstrap support for NJ. In the group formed by bacteria, sequences from Pectobacterium atrosepticum, P. carotovorum and Bacillus subtilis cluster together with 100% posterior probability. This early separation between amino acid sequences of bacteria and those of oomycetes and fungi can be explained in terms of the evolution of lytic enzymes in these microorganisms for different purposes. Bacteria and some anaerobic fungi produce multi-enzymatic complexes called cellulosomes, which are anchored to the cell surface, allow the microorganisms to bind to lignocellulose substrates and increase the breakdown efficiency of cellulose, hemicellulose and pectin [62, 63]. In contrast, in the majority of fungi and oomycetes, cellulases, pectinases and hemicellulases are not integrated in cellulosome complexes, and the pectin degradation is regulated by a multifunctional control system in which the enzymes act in a synergistic manner and are induced by monosaccharides or small oligosaccharides that are generated as products of the same enzymatic reactions [64, 65].
The inferred tree also showed that the analyzed sequences of saprophytic/opportunistic fungi are clustered into a monophyletic group with 98% posterior probability and 75% and 70% bootstrap support for MP and NJ analyses, respectively. However, phytopathogenic fungi and oomycetes were not clustered together. This may be the result of a reduced representation of sequences in the analysis arising from the few PNL sequences reported for members of these groups.
C. lindemuthianum is found clustered with the amino acid sequences of PnlA and Pnl2 of the fungal pathogen C. gloeosporioides with 100% posterior probability for Bayesian analysis as well as 96% and 99% bootstrap support for MP and NJ analysis, respectively.
Pectin and pectate lyases fold into a parallel β-helix, in which a high structural conservation occurs in regions distant from the active site and particularly in those that contribute to the parallel β-helix architecture. The binding cleft and surroundings constitute the most divergent part of the molecule, which allows variation in substrate specificity [13, 15]. On this background, the results of the phylogenetic analyses and the fact that the classification of the pectin lyases is based both on amino acid sequence similarities as well as their structural features , we believe that a structural comparison would help to strengthen the phylogenetic analysis and to establish a relationship between the genes encoding PNLs with their three-dimensional structures involved in carbohydrate binding.
Multiple comparisons of protein structures
Expression analysis of Clpnl2
The differences between the two races were much more noticeable when 92% esterified pectin was used as the sole carbon source. Transcript expression in the pathogenic race started to increase rapidly, reached the highest levels after 4-6 h and then started to decline, giving a still significant increase at the end of the experimental period (Figure 7B and 7C). The maximum transcript levels on this substrate were clearly higher than those observed on glucose. In contrast, the levels of the Clpnl2 transcript in the non-pathogenic race remained undetectable after 8 h of incubation. Late expression occurred after 10-12 h of cultivation, reaching values comparable to those observed after 4 h of culture on cell walls (Figure 7A and 7C). These results agree with the differences found by Hernández et al. , who analyzed the extracellular activity of pectin lyase in both races of C. lindemuthianum under the same conditions employed in this study. When both races were grown with glucose, extracellular PNL activity was barely detected after 8 (race 1472) and 10 (race 0) days of incubation, as observed in this study. Plant cell walls from P. vulgaris induced a similarly low PNL activity in the two isolates after 7-8 days of incubation. When pectin esterified to 92% was used as the carbon source, the activity in the pathogenic race nearly doubled compared with the activity in the non-pathogenic race. Early transcription of genes encoding lytic enzymes and late detection of the corresponding activities is a well documented phenomenon in different fungi [8, 30, 65, 68]. Apart from the presence of a regulatory system controlling gene expression, the production of active pectinase and probably other lyticases can be modulated by other mechanisms such as postranslational modification and protein transport . These alternatives may help to explain the differences observed in this study.
The pectin lyase of the pathogenic race of C. lindemuthianum is able to degrade highly esterified pectin (92%), unlike that of the non-pathogenic race. Apparently, the differences between the pathogenic and non-pathogenic races of C. lindemuthianum occur as much at the expression level as at the level of enzymatic activity, and it is clear that the non-pathogenic and pathogenic races of C. lindemuthianum respond of different form to the carbon sources (except for glucose, where the mRNA of Clpnl2 and the active enzyme is synthesized at basal levels).
It has been proposed that the basal level of enzymatic activity breaks down the substrate, generating degradation products that further induce enzymatic activity . A similar behavior has been observed in our laboratory for other enzymes that degrade cell walls, such as cellulases and the xylanase and β-xylosidase of C. lindemuthianum (unpublished data).
Several studies have reported that the pectinolytic enzymes play an important role in pathogenesis [70, 71]. These are the first enzymes that act during the infection of the plant, causing extensive degradation of the cell wall and the main symptoms of the disease . However, in addition to enzyme production, the sequence in which the enzymes are produced, the speed of synthesis, concentration and diffusion of enzyme are also fundamental aspects of the pathogenesis process . The non-pathogenic race of C. lindemuthianum used in this work is unable to infect P. vulgaris, and thus its lifestyle is closer to that of a saprophytic fungus. Therefore, it is possible that the differences found between the non-pathogenic (0) and pathogenic (1472) races of C. lindemuthianum are related to the speed of activation of the lytic enzyme genes during the interaction with the host.
The number of pectin lyase sequences corresponding to different species of saprophytic/opportunistic fungi used in our analysis surpassed those of pathogenic oomycetes and fungi. This may be because more species of saprophytic/opportunistic have been studied and their degradation systems are better known. Alternatively, the enzymatic diversity may be the evolutionary effect of the heterogeneity of substrates that were encountered during interactions with an extended variety of hosts. For pectate lyases, it has been proposed that differences in the degree of pectin methylation can explain the existence of isozymes .
Pathogenic fungi and those who have close relationships with their host have developmental strategies that allow them to avoid the plant defenses and penetrate cell walls through the use of lytic enzymes. Plants also rely on strategies that allow them to detect and to defend against the attack of pathogens by producing inhibitors of these enzymes [70, 73, 74]. It is therefore possible that the evolution of unique enzymes was induced in pathogenic fungi and that a greater variability of these enzymes was induced in those fungi with a saprophytic lifestyle, which would explain the presence of amino acid sequences and tertiary structures corresponding to enzymes of saprophytic/opportunistic fungi located between the sequences of pathogenic fungi and oomycetes in the phylogenetic analysis and comparison of structures.
There is evidence that supports a relationship between lytic enzyme production and the lifestyles of fungi and oomycetes. For instance, the genome of the oomycete Hyaloperonospora arabidopsidis has lost several of its hydrolytic enzymes compared with Phytophthora sp., which is likely its ancestor [75, 76]. According to an analysis of the hydrolytic profiles of saprophytic/opportunistic and pathogenic fungi using diverse substrates, the species of phytopathogenic fungi are more active than the non-pathogenic fungi on six of eight tested substrates . It has also been observed that pathogenic fungi of monocotyledonous plants are better adapted to degrade the cell walls of monocotyledonous plants, and pathogens of dicotyledonous plants are better able to degrade the cell walls of dicotyledonous plants, reflecting the host preference .
The Clpnl2 gene, which was cloned from a genomic library of C. lindemuthianum, is a unique copy and contains the characteristic elements of a pectin lyase of Family 1 of polysaccharide lyases.
Phylogenetic analyses showed an early separation between the enzymes of bacteria and those of fungi and oomycetes as well as a tendency of the amino acid sequences of fungi and oomycetes to cluster together according to their lifestyle. These results were confirmed by multiple comparison analysis of structures. According to these results, we believe that is possible that the diversity and nature of the substrates processed by these microorganisms play a determining role in the evolution of their lifestyle. In addition, our results showed that both races of C. lindemuthianum express the Clpnl2 gene, although some differences are observed in the timing and level of expression: the pathogenic race responds faster and at higher levels than the non-pathogenic race. This suggests that there are at least two levels of determination of the lifestyle of the microorganisms: one related to the evolution of the enzymes and one concerning the regulation of the expression of the enzymes. In our model, one race of C. lindemuthianum behaves as a hemibiotrophic pathogen and, according to its inability to infect bean, the other race behaves as a saprophyte. Although this study included the analysis of pectin lyase 2 only, we have observed this behavior with other enzymes of the complex involved in the degradation of the cell wall suggesting that it may be a general phenomenon. The differences at this level can be part of the general response of the fungi to host components. However future studies comparing the enzymatic complex of degradation of more fungi species with different lifestyles are needed to confirm this hypothesis.
Finally, we consider this type of information to be of great importance for the study of the biotechnological potential of these enzymes, as the efficiency of the enzymes could depend on the complexity of the vegetal material to be processed and the lifestyle of organism that is the source of enzymes and/or genes.
The authors thank the financial support provided by the FOMIX CONACYT-Gobierno del Estado de Michoacán (project 2009-05 Clave 116208 to HCC) and CONACYT (scholarship granted to ALM and UCS). We thank Gerardo Vázquez Marrufo by its comments to manuscript.
- Willats WG, McCartney L, Mackie W, Knox JP: Pectin: cell biology and prospects for functional analysis. Plant Mol Biol. 2001, 47: 9-27. 10.1023/A:1010662911148.PubMedView ArticleGoogle Scholar
- Mohnen D: Pectin structure and biosynthesis. Curr Opin Plant Biol. 2008, 11: 266-277. 10.1016/j.pbi.2008.03.006.PubMedView ArticleGoogle Scholar
- de Vries RP, Visser J: Aspergillus enzymes involved in degradation of plant cell wall polysaccharides. Microbiol Mol Biol Rev. 2001, 65: 497-522. 10.1128/MMBR.65.4.497-522.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Herron SR, Benen JA, Scavetta RD, Visser J, Jurnak F: Structure and function of pectic enzymes: virulence factors of plant pathogens. Proc Natl Acad Sci USA. 2000, 97: 8762-8769. 10.1073/pnas.97.16.8762.PubMedPubMed CentralView ArticleGoogle Scholar
- Jayani RS, Saxena S, Gupta R: Microbial pectinolytic enzymes: a review. Process Biochem. 2005, 40: 2931-2944. 10.1016/j.procbio.2005.03.026.View ArticleGoogle Scholar
- Prusky D, McEvoy JL, Leverentz B, Conway WS: Local modulation of host pH by Colletotrichum species as a mechanism to increase virulence. Mol Plant Microbe Interact. 2001, 14: 1105-1113. 10.1094/MPMI.2001.14.9.1105.PubMedView ArticleGoogle Scholar
- Maldonado MC, Strasser de Saad AM, Callieri D: Catabolite repression of the synthesis of inducible polygalacturonase and pectinesterase by Aspergillus niger sp. Curr Microbiol. 1989, 18: 303-306. 10.1007/BF01575945.View ArticleGoogle Scholar
- de Vries RP, Jansen J, Aguilar G, Pařenicová L, Joosten V, Wülfert F, Benen JAE, Visser J: Expression profiling of pectinolytic genes from Aspergillus niger. FEBS lett. 2002, 530: 41-47. 10.1016/S0014-5793(02)03391-4.PubMedView ArticleGoogle Scholar
- Lombard V, Bernard T, Rancurel C: A hierarchical classification of polysaccharide lyases for glycogenomics. Biochem J. 2010, 432: 437-444. 10.1042/BJ20101185.PubMedView ArticleGoogle Scholar
- Yoder MD, Keen NT, Jurnak F: New domain motif: the structure of pectate lyase C, a secreted plant virulence factor. Science. 1993, 260: 1503-1507. 10.1126/science.8502994.PubMedView ArticleGoogle Scholar
- Lietzke SE, Yoder MD, Jurnak F: The three-dimensional structure of pectate lyase E, a plant virulence factor from Erwinia chrysanthemi. Plant Physiol. 1994, 106: 849-862.PubMedPubMed CentralGoogle Scholar
- Pickersgill R, Smith D, Worboys K, Jenkins J: Crystal structure of polygalacturonase from Erwinia carotovora ssp. carotovora. J Biol Chem. 1998, 273: 24660-24664. 10.1074/jbc.273.38.24660.PubMedView ArticleGoogle Scholar
- Mayans O, Scott M, Connerton I, Gravesen T, Benen J, Visser J, Pickersgill R, Jenkins J: Two crystal structures of pectin lyase a from Aspergillus reveal a ph driven conformational change and striking divergence in the substrate-binding clefts of pectin and pectate lyases. Structure. 1997, 5: 677-689. 10.1016/S0969-2126(97)00222-0.PubMedView ArticleGoogle Scholar
- Vitali J, Schick B, Kester HC, Visser J, Jurnak F: The tree-dimensional structure of Aspergillus niger pectin lyase B at 1.7-A resolution. Plant Physiol. 1998, 116: 69-80. 10.1104/pp.116.1.69.PubMedPubMed CentralView ArticleGoogle Scholar
- Herron SR, Jurnak F: Mechanistic lessons from structural studies of the pectate lyases. Advances in pectin and pectinase. Edited by: Voragen F, Schols H, Visser Redited by The Netherlands: Klumer Academic Publishers. 2003, 221-233.View ArticleGoogle Scholar
- Kusters-van Someren MA, Harmsen JAM, Kester HCM, Visser J: Structure of the Aspergillus niger pelA gene and its expression in Aspergillus niger and Aspergillus nidulans. Curr Genet. 1991, 20: 293-299. 10.1007/BF00318518.PubMedView ArticleGoogle Scholar
- Kusters-van Someren M, Flipphi M, de Graaff L, den Broeck van H, Kester H, Hinnen A, Visser J: Characterization of the Aspergillus niger pelB gene: structure and regulation of expression. Mol Gen Genet. 1992, 234: 113-120.PubMedView ArticleGoogle Scholar
- Harmsen JAM, Kusters-van Someren MA, Visser J: Cloning and expression of a second Aspergillus niger pectin lyase gene (pelA): Indications of a pectin lyase gene family in A. niger. Curr Genet. 1990, 18: 161-166. 10.1007/BF00312604.PubMedView ArticleGoogle Scholar
- Gysler C, Harmsen JA, Kester HC, Visser J, Heim J: Isolation and structure of the pectin lyase D-encoding gene from Aspergillus niger. Gene. 1990, 89: 101-108. 10.1016/0378-1119(90)90211-9.PubMedView ArticleGoogle Scholar
- Kitamoto N, Yoshino-Yasuda S, Ohmiya K, Tsukagoshi N: A second pectin lyase gene (pel2) from Aspergillus oryzae KBN616: its sequence analysis and overexpression, and characterization of the gene products. J Biosci Bioeng. 2001, 91: 378-381. 10.1263/jbb.91.378.PubMedView ArticleGoogle Scholar
- Kitamoto N, Yoshino-Yasuda S, Ohmiya K, Tsukagoshi N: Sequence analysis and overexpression of a pectin lyase gene (pel1) from Aspergillus oryzae KBN616. Biosci Biotechnol Biochem. 2001, 65: 209-212. 10.1271/bbb.65.209.PubMedView ArticleGoogle Scholar
- Nierman WC, Pain A, Anderson MJ, Wortman JR, Kim HS, Arroyo J, Berriman M, Abe K, Archer DB, Bermejo C, et al: Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature. 2005, 438: 1151-1156. 10.1038/nature04332.PubMedView ArticleGoogle Scholar
- Bazzolli DS, Ribon AOB, de Queiroz MV, de Araújo EF: Molecular characterization and expression profile of pectin-lyase-encoding genes from Penicillium griseoroseum. Can J Microbiol. 2006, 52: 1070-1077. 10.1139/w06-070.PubMedView ArticleGoogle Scholar
- Trigui-Lahiani H, Gargouri A: Cloning, genomic organisation and mRNA expression of a pectin lyase gene from a mutant strain of Penicillium occitanis. Gene. 2007, 388: 54-60. 10.1016/j.gene.2006.09.022.PubMedView ArticleGoogle Scholar
- Templeton MD, Sharrock KR, Bowen JK, Crowhurst RN, Rikkerink EH: The pectin lyase-encoding gene (pnl) family from Glomerella cingulata: characterization of pnlA and its expression in yeast. Gene. 1994, 142: 141-146. 10.1016/0378-1119(94)90369-7.PubMedView ArticleGoogle Scholar
- Wei Y, Shih J, Li J, Goodwin PH: Two pectin lyase genes, pnl-1 and pnl-2, from Colletotrichum gloeosporioides f. sp. malvae differ in a cellulose-binding domain and in their expression during infection of Malva pusilla. Microbiology. 2002, 148: 2149-2157.PubMedView ArticleGoogle Scholar
- Perfect SE, Hughes HB, O'Connell RJ, Green JR: Colletotrichum: A model genus for studies on pathology and fungal-plant interactions. Fungal Genet Biol. 1999, 27: 186-198. 10.1006/fgbi.1999.1143.PubMedView ArticleGoogle Scholar
- Flor H: Current status of the gene-for-gene concept. Annu Rev Phytopathol. 1971, 9: 275-296. 10.1146/annurev.py.09.090171.001423.View ArticleGoogle Scholar
- O'Connell RJ, Bailey J: Differences in the extent of fungal development, host cell necrosis and symptom expression during race-cultivar interactions between Phaseolus vulgaris and Colletotrichum lindemuthianum. Plant Pathol. 1988, 37: 351-362. 10.1111/j.1365-3059.1988.tb02085.x.View ArticleGoogle Scholar
- Wijesundera R, Bailey J, Byrde R, Fielding A: Cell wall degrading enzymes of Colletotrichum lindemuthianum: their role in the development of bean anthracnose. Physiol Mol Plant Pathol. 1989, 34: 403-413. 10.1016/0885-5765(89)90067-2.View ArticleGoogle Scholar
- Knogge W: Fungal pathogenicity. Curr Opinion Plant Biol. 1998, 1: 324-328. 10.1016/1369-5266(88)80054-2.View ArticleGoogle Scholar
- Dodds PN, Rafiqi M, Gan PHP, Hardham AR, Jones DA, Ellis JG: Effectors of biotrophic fungi and oomycetes: pathogenicity factors and triggers of host resistance. New Phytologist. 2009, 183: 993-1000. 10.1111/j.1469-8137.2009.02922.x.PubMedView ArticleGoogle Scholar
- Rodríguez-Guerra R, Acosta-Gallegos J, González-Chavira M, Simpson J: Patotipos de Colletotrichum lindemuthianum y su implicación en la generación de cultivares resistentes de frijol. Agricultura técnica en México. 2006, 32: 101-114.Google Scholar
- Hernández-Silva L, Piñón-Escobedo C, Cano-Camacho H, Zavala-Páramo MG, Acosta-Rodríguez I, López-Romero E: Comparison of fungal growth and production of extracellular pectin lyase activity by pathogenic and non-pathogenic races of Colletotrichum lindemuthianum cultivated under different conditions. Physiol Mol Plant Pathol. 2007, 70: 88-95. 10.1016/j.pmpp.2007.07.005.View ArticleGoogle Scholar
- Tu J: An improved Mathur's medium for growth, sporulation, and germination of spores of Colletotrichum lindemuthianum. Microbios. 1985, 40: 87-96.Google Scholar
- Fry SC: Wall polymers: chemical characterization. The growing plant cell: chemical and metabolic analysis. 2000, New Jersey: The Blackburn Press, 102-115.Google Scholar
- Ausubel F, Brent R, Kingston R, Moore D, Seidman J, Smith J, Struhl K: Short Protocols in Molecular Biology. 1999, Wiley John & Sons Inc, 1104-4Google Scholar
- Borneman J, Hartin RJ: PCR Primers That Amplify Fungal rRNA Genes from Environmental Samples. App Env Microbiol. 2000, 66: 4356-10.1128/AEM.66.10.4356-4360.2000.View ArticleGoogle Scholar
- Grabe N: AliBaba2: context specific identification of transcription factor binding sites. In Silico Biol. 2002, 2: 1-15.Google Scholar
- Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, et al: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, 34: D108-D110. 10.1093/nar/gkj143.PubMedPubMed CentralView ArticleGoogle Scholar
- Nielsen H, Engelbrecht J, Brunak S, Heijne von G: Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 1997, 10: 1-6. 10.1093/protein/10.1.1.PubMedView ArticleGoogle Scholar
- Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein identification and analysis tools on the ExPASy server. The Proteomic Protocols Handbook. Edited by: Walker JM, Totowa. 2005, NJ: Humana Press Inc, 571-607.View ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal × version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
- Philippe H, Delsuc F, Brinkmann H, Lartillot N: PHYLOGENOMICS. Annu Rev Ecol Evol Syst. 2005, 36: 541-562. 10.1146/annurev.ecolsys.35.112202.130205.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Arnold K, Bordoil L, Kopp J, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2005, 22: 195-201.PubMedView ArticleGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 1997, 18: 2714-2723. 10.1002/elps.1150181505.PubMedView ArticleGoogle Scholar
- Van Gunsteren WF, Billeter WF, Eising AA, Hunenberger PH, Krüger P, Mark AE: Biomolecular simulation: The GROMOS96 manual und user guide. vdf Hochs-chulverlag AG an der ETH Zurich and BIOMOS b v. 1996, Zurich, GroningerGoogle Scholar
- Birzele F, Gewehr JE, Csaba G, Zimmer R: Vorolign-fast structural alignment using Voronoi contacts. Bioinformatics. 2007, 23: 205-211.View ArticleGoogle Scholar
- Barthel D, Hirst JD, Błażewicz J, Burke EK, Krasnogor N: ProCKSI: a decision support system for Protein (Structure) Comparison, Knowledge, Similarity and Information. BMC Bioinformatics. 2007, 8: 416-10.1186/1471-2105-8-416.PubMedPubMed CentralView ArticleGoogle Scholar
- Sokal RR, Michener CD: A statistical method for evaluating systematic relationships. Univ Kans Sci Bull. 1958, 38: 1409-1438.Google Scholar
- Kozak M: Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell. 1986, 44: 283-292. 10.1016/0092-8674(86)90762-2.PubMedView ArticleGoogle Scholar
- Kurtz S, Shore D: RAP1 protein activates and silences transcription of mating-type genes in yeast. Genes Dev. 1991, 5: 616-628. 10.1101/gad.5.4.616.PubMedView ArticleGoogle Scholar
- Andrianopoulos A, Timberlake EE: The Aspergillus nidulans abaA gene encodes a transcriptional activator that acts as a genetic switch to control development. Mol Cell Biol. 1994, 14: 2503-2515.PubMedPubMed CentralView ArticleGoogle Scholar
- Borneman AR, Hynes MJ, Andrianopoulos A: The abaA homologue of Penicillium marneffei participates in two developmental programmes: conidiation and dimorphic growth. Mol Microbiol. 2000, 28: 1034-1047.Google Scholar
- Benoist C, O'Hare K, Breathnach R, Chambon P: The ovalbumin gene-sequence of putative control regions. Nucleic Acids Res. 1980, 8: 127-142. 10.1093/nar/8.1.127.PubMedPubMed CentralView ArticleGoogle Scholar
- Zaret KS, Sherman F: DNA sequence required for efficient transcription termination in yeast. Cell. 1982, 28: 563-573. 10.1016/0092-8674(82)90211-2.PubMedView ArticleGoogle Scholar
- Nevalainen KM, Teo VSJ, Bergquist PL: Heterologous protein expression in filamentous fungi. Trends Biotechnol. 2005, 23: 468-474. 10.1016/j.tibtech.2005.06.002.PubMedView ArticleGoogle Scholar
- Jeoh T, Michener W, Himmel ME, Decker SR, Adney WS: Implications of cellobiohydrolase glycosylation for use in biomass conversion. Biotechnol Biofuels. 2008, 1: 1-10. 10.1186/1754-6834-1-1.View ArticleGoogle Scholar
- Bayer EA, Belaich J-P, Shoham Y, Lamed R: The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu Rev Microbiol. 2004, 58: 521-554. 10.1146/annurev.micro.57.030502.091022.PubMedView ArticleGoogle Scholar
- van Dyk JS, Sakka M, Sakka K, Pletschke BI: Identification of endoglucanases, xylanases, pectinases and mannanases in the multi-enzyme complex of Bacillus licheniformis SVD1. Enzym Microb Tech. 2010, 47: 112-118. 10.1016/j.enzmictec.2010.05.004.View ArticleGoogle Scholar
- de Vries RP: Regulation of Aspergillus genes encoding plant cell wall polysaccharide-degrading enzymes; relevance for industrial production. Appl Microbiol Biotechnol. 2003, 61: 10-20.PubMedView ArticleGoogle Scholar
- Martens-Uzunova ES, Schaap PJ: Assessment of the pectin degrading enzyme network of Aspergillus niger by functional genomics. Fungal Genet Biol. 2009, 46: 170-179. 10.1016/j.fgb.2008.07.021.View ArticleGoogle Scholar
- Zhu J, Weng Z: FAST: A novel protein structure alignment algorithm. Proteins. 2004, 58: 618-627. 10.1002/prot.20331.View ArticleGoogle Scholar
- Yona G, Kedem K: The URMS-RMS hybrid algorithm for fast and sensitive local protein structure alignment. J Comput Biol. 2005, 12: 12-32. 10.1089/cmb.2005.12.12.PubMedView ArticleGoogle Scholar
- Acosta-Rodríguez I, Piñón-Escobedo C, Zavala-Páramo MG, López-Romero E, Cano-Camacho H: Degradation of cellulose by the bean-pathogenic fungus Colletotrichum lindemuthianum. Production of extracellular cellulolytic enzymes by cellulose induction. Antonie Van Leeuwenhoek. 2005, 87: 301-310. 10.1007/s10482-004-6422-6.PubMedView ArticleGoogle Scholar
- Herbert C, Boudart G, Borel C, Jacquet C: Regulation and role of pectinases in phytopathogenic fungi. Advances in pectin and pectinase research. Edited by: Voragen F, Schols H, Visser Redited by Netherlands: Kluwer academic publishers. 2003, 201-201208.View ArticleGoogle Scholar
- Have Ten A, Tenberge K, Benen J: The contribution of cell wall degrading enzymes to pathogenesis of fungal plant pathogens. The Mycota XI. Edited by: KempkenF edited by Berlin, Germany: Springer Verlag. 2002, 341-358.Google Scholar
- Lagaert S, Belien T, Volckaert G: Plant cell walls: Protecting the barrier from degradation by microbial enzymes. Semin Cell Dev Biol. 2009, 20: 1064-1073. 10.1016/j.semcdb.2009.05.008.PubMedView ArticleGoogle Scholar
- Alghisi P, Favaron F: Pectin-degrading enzymes and plant-parasite interactions. Eur J Plant Pathol. 1995, 101: 365-375. 10.1007/BF01874850.View ArticleGoogle Scholar
- Maulik A, Ghosh H, Basu S: Comparative study of protein-protein interaction observed in Polygalacturonase-inhibiting proteins from Phaseolus vulgaris and Glycine max and Polygalacturonase from Fusarium moniliforme. BMC Genomics. 2009, 10: S19-PubMedPubMed CentralView ArticleGoogle Scholar
- King BC, Waxman KD, Nenni NV, Walker LP, Bergstrom GC, Gibson DM: Arsenal of plant cell wall degrading enzymes reflects host preference among plant pathogenic fungi. Biotechnol Biofuels. 2011, 4: 4-10.1186/1754-6834-4-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Dodds PN: Genome Evolution in Plant Pathogens. Science. 2010, 330: 1486-1487. 10.1126/science.1200245.PubMedPubMed CentralView ArticleGoogle Scholar
- Baxter L, Tripathy S, Ishaque N, Boot N, Cabral A, Kemen E, Thines M, Ah-Fong A, Anderson R, Badejoko W, et al: Signatures of adaptation to obligate biotrophy in the Hyaloperonospora arabidopsidis genome. Science. 2010, 330: 1549-1551. 10.1126/science.1195203.PubMedPubMed CentralView ArticleGoogle Scholar
- Huson D, Richter D, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.