Functional Classification and Characterization of the Fungal Glycoside Hydrolase 28 Protein Family.

Pectin is a major constituent of the plant cell wall, comprising compounds with important industrial applications such as homogalacturonan, rhamnogalacturonan and xylogalacturonan. A large array of enzymes is involved in the degradation of this amorphous substrate. The Glycoside Hydrolase 28 (GH28) family includes polygalacturonases (PG), rhamnogalacturonases (RG) and xylogalacturonases (XG) that share a structure of three to four pleated ß-sheets that form a rod with the catalytic site amidst a long, narrow groove. Although these enzymes have been studied for many years, there has been no systematic analysis. We have collected a comprehensive set of GH28 encoding sequences to study their evolution in fungi, directed at obtaining a functional classification, as well as at the identification of substrate specificity as functional constraint. Computational tools such as Alphafold, Consurf and MEME were used to identify the subfamilies' characteristics. A hierarchic classification defines the major classes of endoPG, endoRG and endoXG as well as three exoPG classes. Ascomycete endoPGs are further classified in two subclasses whereas we identify four exoRG subclasses. Diversification towards exomode is explained by loops that appear inserted in a number of turns. Substrate-driven diversification can be identified by various specificity determining positions that appear to surround the binding groove.

Structure-function analysis of Sedolisins: evolution of tripeptidyl peptidase and endopeptidase subfamilies in fungi.

BACKGROUND: Sedolisins are acid proteases that are related to the basic subtilisins. They have been identified in all three superkingdoms but are not ubiquitous, although fungi that secrete acids as part of their lifestyle can have up to six paralogs. Both TriPeptidyl Peptidase (TPP) and endopeptidase activity have been identified and it has been suggested that these correspond to separate subfamilies. RESULTS: We studied eukaryotic sedolisins by computational analysis. A maximum likelihood tree shows one major clade containing non-fungal sequences only and two major as well as two minor clades containing only fungal sequences. One of the major fungal clades contains all known TPPs whereas the other contains characterized endosedolisins. We identified four Cluster Specific Inserts (CSIs) in endosedolisins, of which CSIs 1, 3 and 4 appear as solvent exposed according to structure modeling. Part of CSI2 is exposed but a short stretch forms a novel and partially buried α-helix that induces a conformational change near the binding pocket. We also identified a total of 15 specificity determining positions (SDPs) of which five, identified in two independent analyses, form highly connected SDP sub-networks. Modeling of virtual mutants suggests a key role for the W307A or F307A substitution. The remaining four key SDPs physically interact at the interface of the catalytic domain and the enzyme's prosegment. Modeling of virtual mutants suggests these SDPs are indeed required to compensate the conformational change induced by CSI2 and the A307. One of the two small fungal clades concerns a subfamily that contains 213 sequences, is mostly similar to the major TPP subfamily but differs, interestingly, in position 307, showing mostly isoleucine and threonine. CONCLUSIONS: Analysis confirms there are at least two sedolisin subfamilies in fungi: TPPs and endopeptidases, and suggests a third subfamily with unknown characteristics. Sequence and functional diversification was centered around buried SDP307 and resulted in a conformational change of the pocket. Mutual Information network analysis forms a useful instrument in the corroboration of predicted SDPs.

Molecular dynamics and structure function analysis show that substrate binding and specificity are major forces in the functional diversification of Eqolisins.

BACKGROUND: Eqolisins are rare acid proteases found in archaea, bacteria and fungi. Certain fungi secrete acids as part of their lifestyle and interestingly these also have many eqolisin paralogs, up to nine paralogs have been recorded. This suggests a process of functional redundancy and diversification has occurred, which was the subject of the research we performed and describe here. RESULTS: We identified eqolisin homologs by means of iterative HMMER analysis of the NR database. The identified sequences were scrutinized for which new hallmarks were identified by molecular dynamics simulations of mutants in highly conserved positions, using the structure of an eqolisin that was crystallized in the presence of a transition state inhibitor. Four conserved glycines were shown to be important for functionality. A substitution of W67F is shown to be accompanied by the L105W substitution. Molecular dynamics shows that the W67 binds to the substrate via a π-π stacking and a salt bridge, the latter being stronger in a virtual W67F/L105W double mutant of the resolved structure of Scytalido-carboxyl peptidase-B (PDB ID: 2IFW). Additional problematic mutations are discussed. Upon sequence scrutiny we obtained a set of 233 sequences that was used to reconstruct a Bayesian phylogenetic tree. We identified 14 putative specificity determining positions (SDPs) of which four are explained by mere structural explanations and nine seem to correspond to functional diversification related with substrate binding and specificity. A first sub-network of SDPs is related to substrate specificity whereas the second sub-network seems to affect the dynamics of three loops that are involved in substrate binding. CONCLUSION: The eqolisins form a small superfamily of acid proteases with nevertheless many paralogs in acidic fungi. Functional redundancy has resulted in diversification related to substrate specificity and substrate binding.

HMMER Cut-off Threshold Tool (HMMERCTTER): Supervised classification of superfamily protein sequences with a reliable cut-off threshold.

BACKGROUND: Protein superfamilies can be divided into subfamilies of proteins with different functional characteristics. Their sequences can be classified hierarchically, which is part of sequence function assignation. Typically, there are no clear subfamily hallmarks that would allow pattern-based function assignation by which this task is mostly achieved based on the similarity principle. This is hampered by the lack of a score cut-off that is both sensitive and specific. RESULTS: HMMER Cut-off Threshold Tool (HMMERCTTER) adds a reliable cut-off threshold to the popular HMMER. Using a high quality superfamily phylogeny, it clusters a set of training sequences such that the cluster-specific HMMER profiles show cluster or subfamily member detection with 100% precision and recall (P&R), thereby generating a specific threshold as inclusion cut-off. Profiles and thresholds are then used as classifiers to screen a target dataset. Iterative inclusion of novel sequences to groups and the corresponding HMMER profiles results in high sensitivity while specificity is maintained by imposing 100% P&R self detection. In three presented case studies of protein superfamilies, classification of large datasets with 100% precision was achieved with over 95% recall. Limits and caveats are presented and explained. CONCLUSIONS: HMMERCTTER is a promising protein superfamily sequence classifier provided high quality training datasets are used. It provides a decision support system that aids in the difficult task of sequence function assignation in the twilight zone of sequence similarity. All relevant data and source codes are available from the Github repository at the following URL: https://github.com/BBCMdP/HMMERCTTER.

Chlorogenic acid, anthocyanin and flavan-3-ol biosynthesis in flesh and skin of Andean potato tubers (Solanum tuberosum subsp. andigena).

Natural variation of Andean potato was used to study the biosynthesis of phenolic compounds. Levels of phenolic compounds and corresponding structural gene transcripts were examined in flesh and skin of tubers. Phenolic acids, mainly chlorogenic acid (CGA), represent the major compounds, followed by anthocyanins and flavan-3-ols. High-anthocyanin varieties have high levels of CGA. Both metabolite and transcript levels were higher in skin than in flesh and showed a good correspondence. Two hydroxycinnamoyl-CoA transferases (HCT/HQT) have been involved in CGA production, of which HCT reflects CGA levels. Catechin was found in pigmented tissues whereas epicatechin was restricted to tuber skin. Transcripts of leucoanthocyanidin reductase (LCR), which generates catechin, could not be detected. Anthocyanidin reductase (ANR) transcripts, the enzyme responsible for epicatechin production, showed similar levels among samples. These data suggest that the biosynthesis of flavan-3-ols in potato tuber would require ANR but not LCR and that an epimerization process is involved.

Evolutionary and Functional Relationships in the Truncated Hemoglobin Family.

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.

Phospholipase D δ knock-out mutants are tolerant to severe drought stress.

Phospholipase D (PLD) is involved in different plant processes, ranging from responses to abiotic and biotic stress to plant development. Phospholipase Dδ (PLDδ) is activated in dehydration and salt stress, producing the lipid second messenger phosphatidic acid. In this work we show that pldδ Arabidopsis mutants were more tolerant to severe drought than wild-type plants. PLDδ has been shown to be required for ABA regulation of stomatal closure of isolated epidermal peels. However, there was no significant difference in stomatal conductance at the whole plant level between wild-type and pldδ mutants. Since PLD hydrolyses structural phospholipids, then we looked at membrane integrity. Ion leakage measurements showed that during dehydration of leaf discs pldδ mutant has less membrane degradation compared to the wild-type. We further analyzed the mutants and showed that pldδ have higher mRNA levels of RAB18 and RD29A compared to wild-type plants under normal growth conditions. Transient expression of AtPLDδ in Nicotiana benthamiana plants induced a wilting phenotype. These findings suggest that, in wt plants PLDδ disrupt membranes in severe drought stress and, in the absence of the protein (PLDδ knock-out) might drought-prime the plants, making them more tolerant to severe drought stress. The results are discussed in relation to PLDδ role in guard cell signaling and drought tolerance.

Chlorogenic Acid Biosynthesis Appears Linked with Suberin Production in Potato Tuber (Solanum tuberosum).

Potato (Solanum tuberosum L.) is a good source of dietary antioxidants. Chlorogenic acid (CGA) and caffeic acid (CA) are the most abundant phenolic acid antioxidants in potato and are formed by the phenylpropanoid pathway. A number of CGA biosynthetic routes that involve hydroxycinnamoyl-CoA quinate hydroxycinnamoyl transferase (HQT) and/or hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) have been proposed, but little is known about their path in potato. CA production requires a caffeoyl shikimate esterase (CSE), and CA serves as a substrate of lignin precursor ferulic acid via the action of caffeic/5-hydroxyferulic acid O-methyltransferase (COMT I). CGA is precursor of caffeoyl-CoA and, via caffeoyl-CoA O-methyltransferase (CCoAOMT), of feruloyl-CoA. Feruloyl-CoA is required for lignin and suberin biosynthesis, crucial for tuber development. Here, metabolite and transcript levels of the mentioned and related enzymes, such as cinnamate 4-hydroxylase (C4H), were determined in the flesh and skin of fresh and stored tubers. Metabolite and transcript levels were higher in skin than in flesh, irrespective of storage. CGA and CA production appear to occur via p-coumaroyl-CoA, using HQT and CSE, respectively. HCT is likely involved in CGA remobilization toward suberin. The strong correlation between CGA and CA, the correspondence with C4H, HQT, CCoAOMT2, and CSE, and the negative correlation of HCT and COMT I in potato tubers suggest a major flux toward suberin.

The diversity of algal phospholipase D homologs revealed by biocomputational analysis.

Phospholipase D (PLD) participates in the formation of phosphatidic acid, a precursor in glycerolipid biosynthesis and a second messenger. PLDs are part of a superfamily of proteins that hydrolyze phosphodiesters and share a catalytic motif, HxKxxxxD, and hence a mechanism of action. Although HKD-PLDs have been thoroughly characterized in plants, animals and bacteria, very little is known about these enzymes in algae. To fill this gap in knowledge, we performed a biocomputational analysis by means of HMMER iterative profiling, using most eukaryotic algae genomes available. Phylogenetic analysis revealed that algae exhibit very few eukaryotic-type PLDs but possess, instead, many bacteria-like PLDs. Among algae eukaryotic-type PLDs, we identified C2-PLDs and PXPH-like PLDs. In addition, the dinoflagellate Alexandrium tamarense features several proteins phylogenetically related to oomycete PLDs. Our phylogenetic analysis also showed that algae bacteria-like PLDs (proteins with putative PLD activity) fall into five clades, three of which are novel lineages in eukaryotes, composed almost entirely of algae. Specifically, Clade II is almost exclusive to diatoms, whereas Clade I and IV are mainly represented by proteins from prasinophytes. The other two clades are composed of mitochondrial PLDs (Clade V or Mito-PLDs), previously found in mammals, and a subfamily of potentially secreted proteins (Clade III or SP-PLDs), which includes a homolog formerly characterized in rice. In addition, our phylogenetic analysis shows that algae have non-PLD members within the bacteria-like HKD superfamily with putative cardiolipin synthase and phosphatidylserine/phosphatidylglycerophosphate synthase activities. Altogether, our results show that eukaryotic algae possess a moderate number of PLDs that belong to very diverse phylogenetic groups.

The tomato phosphatidylinositol-phospholipase C2 (SlPLC2) is required for defense gene induction by the fungal elicitor xylanase.

The tomato [Solanum lycopersicum (Sl)] phosphatidylinositol-phospholipase C (PI-PLC) gene family is composed of six members, named SlPLC1 to SlPLC6, differentially regulated upon pathogen attack. We have previously shown that the fungal elicitor xylanase rapidly induces nitric oxide (NO), which is required for PI-PLCs activity and downstream defense responses in tomato cell suspensions. Here, we show that all six SlPLC genes are expressed in tomato cell suspensions. Treatment of the cells with xylanase induces an early increase in SlPLC5 transcript levels, followed by a raise of the amount of SlPLC2 transcripts. The production of NO is required to augment SlPLC5 transcript levels in xylanase-treated tomato cells. Xylanase also induces SlPLC2 and SlPLC5 transcript levels in planta. We knocked-down the expression of SlPLC2 and SlPLC5 by virus-induced gene silencing. We found that SlPLC2 is required for xylanase-induced expression of the defense-related genes PR1 and HSR203J.

Extensive expansion of A1 family aspartic proteinases in fungi revealed by evolutionary analyses of 107 complete eukaryotic proteomes.

The A1 family of eukaryotic aspartic proteinases (APs) forms one of the 16 AP families. Although one of the best characterized families, the recent increase in genome sequence data has revealed many fungal AP homologs with novel sequence characteristics. This study was performed to explore the fungal AP sequence space and to obtain an in-depth understanding of fungal AP evolution. Using a comprehensive phylogeny of approximately 700 AP sequences from the complete proteomes of 87 fungi and 20 nonfungal eukaryotes, 11 major clades of APs were defined of which clade I largely corresponds to the A1A subfamily of pepsin-archetype APs. Clade II largely corresponds to the A1B subfamily of nepenthesin-archetype APs. Remarkably, the nine other clades contain only fungal APs, thus indicating that fungal APs have undergone a large sequence diversification. The topology of the tree indicates that fungal APs have been subject to both "birth and death" evolution and "functional redundancy and diversification." This is substantiated by coclustering of certain functional sequence characteristics. A meta-analysis toward the identification of Cluster Determining Positions (CDPs) was performed in order to investigate the structural and biochemical basis for diversification. Seven CDPs contribute to the secondary structure of the enzyme. Three other CDPs are found in the vicinity of the substrate binding cleft. Tree topology, the large sequence variation among fungal APs, and the apparent functional diversification suggest that an amendment to update the current A1 AP classification based on a comprehensive phylogenetic clustering might contribute to refinement of the classification in the MEROPS peptidase database.

Evolution and functional diversification of the small heat shock protein/α-crystallin family in higher plants.

Small heat shock proteins (sHSPs) are chaperones that play an important role in stress tolerance. They consist of an alpha-crystallin domain (ACD) flanked by N- and C-terminal regions. However, not all proteins that contain an ACD, hereafter referred to as ACD proteins, are sHSPs because certain ACD proteins are known to have different functions. Furthermore, since not all ACD proteins have been identified yet, current classifications are incomplete. A total of 17 complete plant proteomes were screened for the presence of ACD proteins by HMMER profiling and the identified ACD protein sequences were classified by maximum likelihood phylogeny. Differences among and within groups were analysed, and levels of functional constraint were determined. There are 29 different classes of ACD proteins, eight of which contain classical sHSPs and five likely chaperones. The other classes contain proteins with uncharacterised or poorly characterised functions. N- and C-terminal sequences are conserved within the phylogenetic classes. Phylogenetics suggests a single duplication of the CI sHSP ancestor that occurred prior to the speciation of mono- and dicotyledons. This was followed by a number of more recent duplications that resulted in the presence of many paralogues. The results suggest that N- and C-terminal sequences of sHSPs play a role in class-specific functionality and that non-sHSP ACD proteins have conserved but unexplored functions, which are mainly determined by subsequences other than that of the ACD.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea.

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

The aspartic proteinase family of three Phytophthora species.

BACKGROUND: Phytophthora species are oomycete plant pathogens with such major social and economic impact that genome sequences have been determined for Phytophthora infestans, P. sojae and P. ramorum. Pepsin-like aspartic proteinases (APs) are produced in a wide variety of species (from bacteria to humans) and contain conserved motifs and landmark residues. APs fulfil critical roles in infectious organisms and their host cells. Annotation of Phytophthora APs would provide invaluable information for studies into their roles in the physiology of Phytophthora species and interactions with their hosts. RESULTS: Genomes of Phytophthora infestans, P. sojae and P. ramorum contain 11-12 genes encoding APs. Nine of the original gene models in the P. infestans database and several in P. sojae and P. ramorum (three and four, respectively) were erroneous. Gene models were corrected on the basis of EST data, consistent positioning of introns between orthologues and conservation of hallmark motifs. Phylogenetic analysis resolved the Phytophthora APs into 5 clades. Of the 12 sub-families, several contained an unconventional architecture, as they either lacked a signal peptide or a propart region. Remarkably, almost all APs are predicted to be membrane-bound. CONCLUSIONS: One of the twelve Phytophthora APs is an unprecedented fusion protein with a putative G-protein coupled receptor as the C-terminal partner. The others appear to be related to well-documented enzymes from other species, including a vacuolar enzyme that is encoded in every fungal genome sequenced to date. Unexpectedly, however, the oomycetes were found to have both active and probably-inactive forms of an AP similar to vertebrate BACE, the enzyme responsible for initiating the processing cascade that generates the Aß peptide central to Alzheimer's Disease. The oomycetes also encode enzymes similar to plasmepsin V, a membrane-bound AP that cleaves effector proteins of the malaria parasite Plasmodium falciparum during their translocation into the host red blood cell. Since the translocation of Phytophthora effector proteins is currently a topic of intense research activity, the identification in Phytophthora of potential functional homologues of plasmepsin V would appear worthy of investigation. Indeed, elucidation of the physiological roles of the APs identified here offers areas for future study. The significant revision of gene models and detailed annotation presented here should significantly facilitate experimental design.

Phosphatidic acid production in chitosan-elicited tomato cells, via both phospholipase D and phospholipase C/diacylglycerol kinase, requires nitric oxide.

Nitric oxide (NO) and the lipid second messenger phosphatidic acid (PA) are involved in plant defense responses during plant-pathogen interactions. NO has been shown to be involved in the induction of PA production in response to the pathogen associated molecular pattern (PAMP) xylanase in tomato cells. It was shown that NO is critical for PA production induced via phospholipase C (PLC) in concerted action with diacylglycerol kinase (DGK) but not for the xylanase-induced PA via phospholipase D (PLD). In order to study whether this is a general phenomenon during PAMP perception or if it is particular for xylanase, we studied the effect of the PAMP chitosan in tomato cell suspensions. We observed a rapid NO production in tomato cells treated with chitosan. Chitosan induced the formation of PA by activating both PLD and PLC/DGK. The activation of either phospholipase-mediated signaling pathway was inhibited in cells treated with the NO scavenger cPTIO. This indicates that NO is required for PA generation via both the PLD and PLC/DGK pathway during plant defense response in chitosan elicited cells. Responses downstream PA were studied. PLC inhibitors neomycin and U73122 inhibited chitosan-induced ROS production. Differences between xylanase and chitosan-induced phospholipid signaling pathways are discussed.

Identification of tomato phosphatidylinositol-specific phospholipase-C (PI-PLC) family members and the role of PLC4 and PLC6 in HR and disease resistance.

The perception of pathogen-derived elicitors by plants has been suggested to involve phosphatidylinositol-specific phospholipase-C (PI-PLC) signalling. Here we show that PLC isoforms are required for the hypersensitive response (HR) and disease resistance. We characterised the tomato [Solanum lycopersicum (Sl)] PLC gene family. Six Sl PLC-encoding cDNAs were isolated and their expression in response to infection with the pathogenic fungus Cladosporium fulvum was studied. We found significant regulation at the transcriptional level of the various SlPLCs, and SlPLC4 and SlPLC6 showed distinct expression patterns in C. fulvum-resistant Cf-4 tomato. We produced the encoded proteins in Escherichia coli and found that both genes encode catalytically active PI-PLCs. To test the requirement of these Sl PLCs for full Cf-4-mediated recognition of the effector Avr4, we knocked down the expression of the encoding genes by virus-induced gene silencing. Silencing of SlPLC4 impaired the Avr4/Cf-4-induced HR and resulted in increased colonisation of Cf-4 plants by C. fulvum expressing Avr4. Furthermore, expression of the gene in Nicotiana benthamiana enhanced the Avr4/Cf-4-induced HR. Silencing of SlPLC6 did not affect HR, whereas it caused increased colonisation of Cf-4 plants by the fungus. Interestingly, Sl PLC6, but not Sl PLC4, was also required for resistance to Verticillium dahliae, mediated by the transmembrane Ve1 resistance protein, and to Pseudomonas syringae, mediated by the intracellular Pto/Prf resistance protein couple. We conclude that there is a differential requirement of PLC isoforms for the plant immune response and that Sl PLC4 is specifically required for Cf-4 function, while Sl PLC6 may be a more general component of resistance protein signalling.

The Botrytis cinerea aspartic proteinase family.

The ascomycete plant pathogen Botrytis cinerea secretes aspartic proteinase (AP) activity. Functional analysis was carried out on five aspartic proteinase genes (Bcap1-5) reported previously. Single and double mutants lacking these five genes showed neither a reduced secreted proteolytic activity, nor a reduction in virulence and they showed no alteration in sensitivity to antifungal proteins purified from grape juice. Scrutiny of the B. cinerea genome revealed the presence of nine additional Bcap genes, denoted Bcap6-14. The product of the Bcap8 gene was found to constitute up to 23% of the total protein secreted by B. cinerea. Bcap8-deficient mutants secreted approximately 70% less AP activity but were just as virulent as the wild-type strain. Phylogenetic analysis showed that Bcap8 has orthologs in many basidiomycetes but only few ascomycetes including the biocontrol fungus Trichoderma harzanium. Potential functions of the 14 APs in B. cinerea are discussed based on their sequence characteristics, phylogeny and predicted localization.

Three QTLs for Botrytis cinerea resistance in tomato.

Tomato (Solanum lycopersicum) is susceptible to grey mold (Botrytis cinerea). Partial resistance to this fungus was identified in accessions of wild relatives of tomato such as S. habrochaites LYC4. In order to identify loci involved in quantitative resistance (QTLs) to B. cinerea, a population of 174 F(2) plants was made originating from a cross between S. lycopersicum cv. Moneymaker and S. habrochaites LYC4. The population was genotyped and tested for susceptibility to grey mold using a stem bioassay. Rbcq1, a QTL reducing lesion growth (LG) and Rbcq2, a QTL reducing disease incidence (DI) were identified. Rbcq1 is located on Chromosome 1 and explained 12% of the total phenotypic variation while Rbcq2 is located on Chromosome 2 and explained 15% of the total phenotypic variation. Both QTL effects were confirmed by assessing disease resistance in two BC(2)S(1) progenies segregating for either of the two QTLs. One additional QTL, Rbcq4 on Chromosome 4 reducing DI, was identified in one of the BC(2)S(1) progenies. F(2) individuals, homozygous for the Rbcq2 and Rbcq4 alleles of S. habrochaites showed a reduction of DI by 48%. QTLs from S. habrochaites LYC4 offer good perspectives for breeding B. cinerea resistant tomato cultivars.

An aspartic proteinase gene family in the filamentous fungus Botrytis cinerea contains members with novel features.

Botrytis cinerea, an important fungal plant pathogen, secretes aspartic proteinase (AP) activity in axenic cultures. No cysteine, serine or metalloproteinase activity could be detected. Proteinase activity was higher in culture medium containing BSA or wheat germ extract, as compared to minimal medium. A proportion of the enzyme activity remained in the extracellular glucan sheath. AP was also the only type of proteinase activity in fluid obtained from B. cinerea-infected tissue of apple, pepper, tomato and zucchini. Five B. cinerea genes encoding an AP were cloned and denoted Bcap1-5. Features of the encoded proteins are discussed. BcAP1, especially, has novel characteristics. A phylogenetic analysis was performed comprising sequences originating from different kingdoms. BcAP1 and BcAP5 did not cluster in a bootstrap-supported clade. BcAP2 clusters with vacuolar APs. BcAP3 and BcAP4 cluster with secreted APs in a clade that also contains glycosylphosphatidylinositol-anchored proteinases from Saccharomyces cerevisiae and Candida albicans. All five Bcap genes are expressed in liquid cultures. Transcript levels of Bcap1, Bcap2, Bcap3 and Bcap4 are subject to glucose and peptone repression. Transcripts from all five Bcap genes were detected in infected plant tissue, indicating that at least part of the AP activity in planta originates from the pathogen.

A tomato metacaspase gene is upregulated during programmed cell death in Botrytis cinerea-infected leaves.

Programmed cell death (PCD) in plant cells is often accompanied by biochemical and morphological hallmarks similar to those of animal apoptosis. However, orthologs of animal caspases, cysteinyl aspartate-specific proteases that constitute the core component of animal apoptosis, have not yet been identified in plants. Recent studies have revealed the presence of a family of genes encoding proteins with distant homology to mammalian caspases, designated metacaspases, in the Arabidopsis thaliana genome. Here, we describe the isolation of LeMCA1, a type-II metacaspase cDNA clone from tomato (Lycopersicon esculentum Mill.). BLAST analysis demonstrated that the LeMCA1 gene is located in close vicinity of several genes that have been linked with PCD. Southern analysis indicated the existence of at least one more metacaspase in the tomato genome. LeMCA1 mRNA levels rapidly increased upon infection of tomato leaves with Botrytis cinerea, a fungal pathogen that induces cell death in several plant species. LeMCA1 was not upregulated during chemical-induced PCD in suspension-cultured tomato cells.