Conserved catalytic machinery and the prediction of a common fold for several families of glycosyl hydrolases

Insight into the functional roles of Glu175 in the hyperthermostable xylanase XYL10C-ΔN through structural analysis and site-saturation mutagenesis

BACKGROUND

Improving the hydrolytic performance of hemicellulases to degrade lignocellulosic biomass is of considerable importance for second-generation biorefinery. Xylanase, as the crucial hemicellulase, must be thermostable and have high activity for its potential use in the bioethanol industry. To obtain excellent xylanase candidates, it is necessary to understand the structure-function relationships to provide a meaningful reference to improve the enzyme properties. This study aimed to investigate the catalytic mechanism of a highly active hyperthermophilic xylanase variant, XYL10C-ΔN, for hemicellulose degradation.

RESULTS

By removing the N-terminal 66 amino acids, the variant XYL10C-ΔN showed a 1.8-fold improvement in catalytic efficiency and could hydrolyze corn stover more efficiently in hydrolysis of corn stover; however, it retained similar thermostability to the wild-type XYL10C. Based on the crystal structures of XYL10C-ΔN and its complex with xylobiose, Glu175 located on loop 3 was found to be specific to GH10 xylanases and probably accounts for the excellent enzyme properties by interacting with Lys135 and Met137 on loop 2. Site-saturation mutagenesis confirmed that XYL10C-ΔN with glutamate acid at position 175 had the highest catalytic efficiency, specific activity, and the broadest pH-activity profile. The functional roles of Glu175 were also verified in the mutants of another two GH10 xylanases, XylE and XynE2, which showed increased catalytic efficiencies and wider pH-activity profiles.

CONCLUSIONS

XYL10C-ΔN, with excellent thermostability, high catalytic efficiency, and great lignocellulose-degrading capability, is a valuable candidate xylanase for the biofuel industry. The mechanism underlying improved activity of XYN10C-ΔN was thus investigated through structural analysis and functional verification, and Glu175 was identified to play the key role in the improved catalytic efficiency. This study revealed the importance of a key residue (Glu175) in XYN10C-ΔN and provides a reference to modify GH10 xylanases for improved catalytic performance.

Gene cloning and expression of fungal lignocellulolytic enzymes from the rumen of gayal (Bos frontalis)

A total of 6,219 positive clones were obtained by constructing a BAC library of uncultured ruminal fungi of gayal, and two clones (xynF1 and eglF2) with lignocellulolytic enzyme activity were selected. The sequencing results showed that xynF1 and eglF2 had 903-bp, and 1,995-bp, open reading frames likely to encode β-xylanase (XynF1) and β-glucosidase (EglF2), respectively. The amino acid sequence of XynF1 had 99% coverage and 95% homology to the endo-β-1,4-xylanase encoded by the cellulase gene of Orpinomyces sp. LT-3 (GenBank accession No. AEO51791.1). The amino acid sequence of EglF2 had 99% coverage and 93% homology to the β-glucosidase encoded by the cellulase gene of Piromyces sp. E2 (GenBank accession No. CAC34952.1). Analysis using the SMART software showed that XynF1 contains a glycoside hydrolase family 11 functional module and a carbohydrate-binding module, while EglF2 contains a glycoside hydrolase family 1 functional module. XynF1 showed the highest relative enzymatic activity, up to 95%, at 45°C and pH 4.2, while EglF2 showed the highest relative enzymatic activity, up to 95%, at 55°C and pH 6.2. In this study, we achieved efficient expression of the xynF1 and eglF2 genes in Pichia pastoris, which laid a foundation for the practical application of the lignocellulolytic enzymes.

Overexpression of a Cellobiose-Glucose-Halotolerant Endoglucanase from Scytalidium thermophilum

Enzyme reaction products and by-products from pretreatment steps can inhibit endoglucanases and are major factors limiting the efficiency of enzymatic lignocellulosic biomass hydrolysis. The gene encoding the endoglucanase from Scytalidium thermophilum (egst) was cloned and expressed as a soluble protein in Pichia pastoris GS115. The recombinant enzyme (Egst) was monomeric (66 kDa) and showed an estimated carbohydrate content of 53.3% (w/w). The optimum temperature and pH of catalysis were 60-70 °C and pH of 5.5, respectively. The enzyme was highly stable at pH 3.0-8.0 with a half-life in water of 100 min at 65 °C. The Egst presented good halotolerance, retaining 84.1 and 71.4% of the control activity in the presence of 0.5 and 2.0 mol L^-1 NaCl, respectively. Hydrolysis of medium viscosity carboxymethylcellulose (CMC) by Egst was stimulated 1.77-, 1.84-, 1.64-, and 1.8-fold by dithiothreitol, β-mercaptoethanol, cysteine, and manganese at 10, 10, 10, and 5 mmol L^-1 concentration, respectively. The enzyme hydrolyzed CMC with maximal velocity and an apparent affinity constant of 432.10 ± 16.76 and 10.5 ± 2.53 mg mL^-1, respectively. Furthermore, the Egst was tolerant to reaction products and able to act on pretreated fractions sugarcane bagasse demonstrating excellent properties for application in the hydrolysis of lignocellulosic biomass.

An Integrated Computational Framework to Assess the Mutational Landscape of α-L-Iduronidase IDUA Gene

Mucopolysaccharidosis type I is a lysosomal genetic disorder caused due to the deficiency of the α-L-iduronidase enzyme (IDUA). Mutations associated with IDUA lead to mild to severe forms of diseases characterized by different clinical features. In the present study, we first performed a comprehensive analysis using various in silico prediction tools to screen and prioritize the missense mutations or nonsynonymous SNPs (nsSNPs) associated with IDUA. Subsequently, statistical analysis was empowered to examine the predictive ability and accuracy of the in silico prediction tool results supporting the disease phenotype ranging from mild to severe. Till date, no study has been carried out in IDUA in analyzing the impact of the nsSNPs at the structural level. In this context with the aid of pathogenic and stability prediction in silico tools, we identified nsSNPs R89Q, R89W, and P533R to be most deleterious and disease-causing having impact on the function of the protein. Extensive molecular dynamics analysis was performed using Gromacs to understand the deleterious nature of the mutants. Variations observed between the trajectory files of native and mutants R89Q, R89W, and P533R using Gromacs utilities enabled us to measure the adverse effects on the protein and could be the underlying reasons for the disease pathogenesis. These findings may be helpful in understanding the genotype-phenotype relationship and molecular basis of the disease to design drugs for better treatment. J. Cell. Biochem. 119: 555-565, 2018. © 2017 Wiley Periodicals, Inc.

Structural basis of the substrate specificity and instability in solution of a glycosidase from Lactobacillus plantarum

Statistics from structural genomics initiatives reveal that around 50-55% of the expressed, non-membrane proteins cannot be purified and therefore structurally characterized due to solubility problems, which emphasized protein solubility as one of the most serious concerns in structural biology projects. Lactobacillus plantarum CECT 748^T produces an aggregation-prone glycosidase (LpBgl) that we crystallized previously. However, this result could not be reproduced due to protein instability and therefore further high-resolution structural analyses of LpBgl were impeded. The obtained crystals of LpBgl diffracted up to 2.48Å resolution and permitted to solve the structure of the enzyme. Analysis of the active site revealed a pocket for phosphate-binding with an uncommon architecture, where a phosphate molecule is tightly bound suggesting the recognition of 6-phosphoryl sugars. In agreement with this observation, we showed that LpBgl exhibited 6-phospho-β-glucosidase activity. Combination of structural and mass spectrometry results revealed the formation of dimethyl arsenic adducts on the solvent exposed cysteine residues Cys211 and Cys292. Remarkably, the double mutant Cys211Ser/Cys292Ser resulted stable in solution at high concentrations indicating that the marginal solubility of LpBgl can be ascribed specifically to these two cysteine residues. The 2.30Å crystal structure of this double mutant showed no disorder around the newly incorporated serine residues and also loop rearrangements within the phosphate-binding site. Notably, LpBgl could be prepared at high yield by proteolytic digestion of the fusion protein LSLt-LpBgl, which raises important questions about potential hysteretic processes upon its initial production as an enzyme fused to a solubility enhancer.

Gene Expression and Molecular Characterization of a Xylanase from Chicken Cecum Metagenome

A xylanase gene xynA_MG1 with a 1,116-bp open reading frame, encoding an endo-β-1,4-xylanase, was cloned from a chicken cecum metagenome. The translated XynA_MG1 protein consisted of 372 amino acids including a putative signal peptide of 23 amino acids. The calculated molecular mass of the mature XynA_MG1 was 40,013 Da, with a theoretical pI value of 5.76. The amino acid sequence of XynA_MG1 showed 59% identity to endo-β-1,4-xylanase from Prevotella bryantii and Prevotella ruminicola and 58% identity to that from Prevotella copri. XynA_MG1 has two conserved motifs, DVVNE and TEXD, containing two active site glutamates and an invariant asparagine, characteristic of GH10 family xylanase. The xynA_MG1 gene without signal peptide sequence was cloned and fused with thioredoxin protein (Trx.Tag) in pET-32a plasmid and overexpressed in Escherichia coli Tuner™(DE3)pLysS. The purified mature XynA_MG1 was highly salt-tolerant and stable and displayed higher than 96% of its catalytic activity in the reaction containing 1 to 4 M NaCl. It was only slightly affected by common organic solvents added in aqueous solution to up to 5 M. This chicken cecum metagenome-derived xylanase has potential applications in animal feed additives and industrial enzymatic processes requiring exposure to high concentrations of salt and organic solvents.

Dual role of imidazole as activator/inhibitor of sweet almond (Prunus dulcis) β-glucosidase

The activity of Prunus dulcis (sweet almond) β-glucosidase at the expense of p-nitrophenyl-β-d-glucopyranoside at pH 6 was determined, both under steady-state and pre-steady-state conditions. Using crude enzyme preparations, competitive inhibition by 1-5 mM imidazole was observed under both kinetic conditions tested. However, when imidazole was added to reaction mixtures at 0.125-0.250 mM, we detected a significant enzyme activation. To further inspect this effect exerted by imidazole, β-glucosidase was purified to homogeneity. Two enzyme isoforms were isolated, i.e. a full-length monomer, and a dimer containing a full-length and a truncated subunit. Dimeric β-glucosidase was found to perform much better than the monomeric enzyme, independently of the kinetic conditions used to assay enzyme activity. In addition, the sensitivity towards imidazole was found to differ between the two isoforms. While monomeric enzyme was indeed found to be relatively insensitive to imidazole, dimeric β-glucosidase was observed to be significantly activated by 0.125-0.250 mM imidazole under pre-steady-state conditions. Further, steady-state assays revealed that the addition of 0.125 mM imidazole to reaction mixtures increases the Km of dimeric enzyme from 2.3 to 6.7 mM. The activation of β-glucosidase dimer by imidazole is proposed to be exerted via a conformational transition poising the enzyme towards proficient catalysis.

Isolation and characterization of a novel endo-β-1,4-glucanase from a metagenomic library of the black-goat rumen

The various types of lignocellulosic biomass found in plants comprise the most abundant renewable bioresources on Earth. In this study, the ruminal microbial ecosystem of black goats was explored because of their strong ability to digest lignocellulosic forage. A metagenomic fosmid library containing 115,200 clones was prepared from the black-goat rumen and screened for a novel cellulolytic enzyme. The KG35 gene, containing a novel glycosyl hydrolase family 5 cellulase domain, was isolated and functionally characterized. The novel glycosyl hydrolase family 5 cellulase gene is composed of a 963-bp open reading frame encoding a protein of 320 amino acid residues (35.1 kDa). The deduced amino acid sequence showed the highest sequence identity (58%) for sequences from the glycosyl hydrolase family 5 cellulases. The novel glycosyl hydrolase family 5 cellulase gene was overexpressed in Escherichia coli. Substrate specificity analysis revealed that this recombinant glycosyl hydrolase family 5 cellulase functions as an endo-β-1,4-glucanase. The recombinant KG35 endo-β-1,4-glucanase showed optimal activity within the range of 30–50 °C at a pH of 6–7. The thermostability was retained and the pH was stable in the range of 30–50 °C at a pH of 5–7.

A Bacteroidetes locus dedicated to fungal 1,6-β-glucan degradation: Unique substrate conformation drives specificity of the key endo-1,6-β-glucanase

Glycans are major nutrients available to the human gut microbiota. The Bacteroides are generalist glycan degraders, and this function is mediated largely by polysaccharide utilization loci (PULs). The genomes of several Bacteroides species contain a PUL, PUL1,6-β-glucan, that was predicted to target mixed linked plant 1,3;1,4-β-glucans. To test this hypothesis we characterized the proteins encoded by this locus in Bacteroides thetaiotaomicron, a member of the human gut microbiota. We show here that PUL1,6-β-glucan does not orchestrate the degradation of a plant polysaccharide but targets a fungal cell wall glycan, 1,6-β-glucan, which is a growth substrate for the bacterium. The locus is up-regulated by 1,6-β-glucan and encodes two enzymes, a surface endo-1,6-β-glucanase, BT3312, and a periplasmic β-glucosidase that targets primarily 1,6-β-glucans. The non-catalytic proteins encoded by PUL1,6-β-glucan target 1,6-β-glucans and comprise a surface glycan-binding protein and a SusD homologue that delivers glycans to the outer membrane transporter. We identified the central role of the endo-1,6-β-glucanase in 1,6-β-glucan depolymerization by deleting bt3312, which prevented the growth of B. thetaiotaomicron on 1,6-β-glucan. The crystal structure of BT3312 in complex with β-glucosyl-1,6-deoxynojirimycin revealed a TIM barrel catalytic domain that contains a deep substrate-binding cleft tailored to accommodate the hook-like structure adopted by 1,6-β-glucan. Specificity is driven by the complementarity of the enzyme active site cleft and the conformation of the substrate. We also noted that PUL1,6-β-glucan is syntenic to many PULs from other Bacteroidetes, suggesting that utilization of yeast and fungal cell wall 1,6-β-glucans is a widespread adaptation within the human microbiota.

Transcriptome analysis of the whitefly, Bemisia tabaci MEAM1 during feeding on tomato infected with the crinivirus, Tomato chlorosis virus, identifies a temporal shift in gene expression and differential regulation of novel orphan genes

BACKGROUND

Whiteflies threaten agricultural crop production worldwide, are polyphagous in nature, and transmit hundreds of plant viruses. Little is known how whitefly gene expression is altered due to feeding on plants infected with a semipersistently transmitted virus. Tomato chlorosis virus (ToCV; genus Crinivirus, family Closteroviridae) is transmitted by the whitefly (Bemisia tabaci) in a semipersistent manner and infects several globally important agricultural and ornamental crops, including tomato.

RESULTS

To determine changes in global gene regulation in whiteflies after feeding on tomato plants infected with a crinivirus (ToCV), comparative transcriptomic analysis was performed using RNA-Seq on whitefly (Bemisia tabaci MEAM1) populations after 24, 48, and 72 h acquisition access periods on either ToCV-infected or uninfected tomatoes. Significant differences in gene expression were detected between whiteflies fed on ToCV-infected tomato and those fed on uninfected tomato among the three feeding time periods: 447 up-regulated and 542 down-regulated at 24 h, 4 up-regulated and 7 down-regulated at 48 h, and 50 up-regulated and 160 down-regulated at 72 h. Analysis revealed differential regulation of genes associated with metabolic pathways, signal transduction, transport and catabolism, receptors, glucose transporters, α-glucosidases, and the uric acid pathway in whiteflies fed on ToCV-infected tomatoes, as well as an abundance of differentially regulated novel orphan genes. Results demonstrate for the first time, a specific and temporally regulated response by the whitefly to feeding on a host plant infected with a semipersistently transmitted virus, and advance the understanding of the whitefly vector-virus interactions that facilitate virus transmission.

CONCLUSION

Whitefly transmission of semipersistent viruses is believed to require specific interactions between the virus and its vector that allow binding of virus particles to factors within whitefly mouthparts. Results provide a broader understanding of the potential mechanism of crinivirus transmission by whitefly, aid in discerning genes or loci in whitefly that influence virus interactions or transmission, and subsequently facilitate development of novel, genetics-based control methods against whitefly and whitefly-transmitted viruses.

Effect of light and prey availability on gene expression of the mixotrophic chrysophyte, Ochromonas sp

Background

Ochromonas is a genus of mixotrophic chrysophytes that is found ubiquitously in many aquatic environments. Species in this genus can be important consumers of bacteria but vary in their ability to perform photosynthesis. We studied the effect of light and bacteria on growth and gene expression of a predominantly phagotrophic Ochromonas species. Axenic cultures of Ochromonas sp. were fed with heat-killed bacteria (HKB) and grown in constant light or darkness. RNA was extracted from cultures in the light or in the dark with HKB present (Light + HKB; Dark + HKB), and in the light after HKB were depleted (Light + depleted HKB).

Results

There were no significant differences in the growth or bacterial ingestion rates between algae grown in light or dark conditions. The availability of light led to a differential expression of only 8% of genes in the transcriptome. A number of genes associated with photosynthesis, phagotrophy, and tetrapyrrole synthesis was upregulated in the Light + HKB treatment compared to Dark + HKB. Conversely, the comparison between the Light + HKB and Light + depleted HKB treatments revealed that the presence of HKB led to differential expression of 59% of genes, including the majority of genes involved in major carbon and nitrogen metabolic pathways. Genes coding for unidirectional enzymes for the utilization of glucose were upregulated in the presence of HKB, implying increased glycolytic activities during phagotrophy. Algae without HKB upregulated their expression of genes coding for ammonium transporters, implying uptake of inorganic nitrogen from the culture medium when prey were unavailable.

Conclusions

Transcriptomic results agreed with previous observations that light had minimal effect on the population growth of Ochromonas sp. However, light led to the upregulation of a number of phototrophy- and phagotrophy-related genes, while the availability of bacterial prey led to prominent changes in major carbon and nitrogen metabolic pathways. Our study demonstrated the potential of transcriptomic approaches to improve our understanding of the trophic physiologies of complex mixotrophs, and revealed responses in Ochromonas sp. not apparent from traditional culture studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3549-1) contains supplementary material, which is available to authorized users.

The Structure of an Archaeal β-Glucosaminidase Provides Insight into Glycoside Hydrolase Evolution

The archaeal exo-β-d-glucosaminidase (GlmA) is a dimeric enzyme that hydrolyzes chitosan oligosaccharides into monomer glucosamines. GlmA is a member of the glycosidase hydrolase (GH)-A superfamily-subfamily 35 and is a novel enzyme in terms of its primary structure. Here, we present the crystal structure of GlmA in complex with glucosamine at 1.27 Å resolution. The structure reveals that a monomeric form of GlmA shares structural homology with GH42 β-galactosidases, whereas most of the spatial positions of the active site residues are identical to those of GH35 β-galactosidases. We found that upon dimerization, the active site of GlmA changes shape, enhancing its ability to hydrolyze the smaller substrate in a manner similar to that of homotrimeric GH42 β-galactosidase. However, GlmA can differentiate glucosamine from galactose based on one charged residue while using the "evolutionary heritage residue" it shares with GH35 β-galactosidase. Our study suggests that GH35 and GH42 β-galactosidases evolved by exploiting the structural features of GlmA.

Computational determination of the effects of virulent Escherichia coli and salmonella bacteriophages on human gut

BACKGROUND AND OBJECTIVE

Salmonella and Escherichia coli are different types of bacteria that cause food poisoning in humans. In the elderly, infants and people with chronic conditions, it is very dangerous if Salmonella or E. coli gets into the bloodstream and then they must be treated by phage therapy. Treating Salmonella and E. coli by phage therapy affects the gut flora. This research paper presents a system for detecting the effects of virulent E. coli and Salmonella bacteriophages on human gut.

METHODS

A method based on Domain-Domain Interactions (DDIs) model is implemented in the proposed system to determine the interactions between the proteins of human gut bacteria and the proteins of bacteriophages that infect virulent E. coli and Salmonella. The system helps gastroenterologists to realize the effect of injecting bacteriophages that infect virulent E. coli and Salmonella on the human gut.

RESULTS

By testing the system over Enterobacteria phage 933W, Enterobacteria phage VT2-Sa and Enterobacteria phage P22, it resulted in four interactions between the proteins of the bacteriophages that infect E. coli O157:H7, E. coli O104:H4 and Salmonella typhimurium and the proteins of human gut bacterium strains.

CONCLUSION

Several effects were detected such as: antibacterial activity against a number of bacterial species in human gut, regulation of cellular differentiation and organogenesis during gut, lung, and heart development, ammonia assimilation in bacteria, yeasts, and plants, energizing defense system and its function in the detoxification of lipopolysaccharide, and in the prevention of bacterial translocation in human gut.

Immobilization of Glycoside Hydrolase Families GH1, GH13, and GH70: State of the Art and Perspectives

Glycoside hydrolases (GH) are enzymes capable to hydrolyze the glycosidic bond between two carbohydrates or even between a carbohydrate and a non-carbohydrate moiety. Because of the increasing interest for industrial applications of these enzymes, the immobilization of GH has become an important development in order to improve its activity, stability, as well as the possibility of its reuse in batch reactions and in continuous processes. In this review, we focus on the broad aspects of immobilization of enzymes from the specific GH families. A brief introduction on methods of enzyme immobilization is presented, discussing some advantages and drawbacks of this technology. We then review the state of the art of enzyme immobilization of families GH1, GH13, and GH70, with special attention on the enzymes β-glucosidase, α-amylase, cyclodextrin glycosyltransferase, and dextransucrase. In each case, the immobilization protocols are evaluated considering their positive and negative aspects. Finally, the perspectives on new immobilization methods are briefly presented.

Periplasmic Cytophaga hutchinsonii Endoglucanases Are Required for Use of Crystalline Cellulose as the Sole Source of Carbon and Energy

The soil bacterium Cytophaga hutchinsonii actively digests crystalline cellulose by a poorly understood mechanism. Genome analyses identified nine genes predicted to encode endoglucanases with roles in this process. No predicted cellobiohydrolases, which are usually involved in the utilization of crystalline cellulose, were identified. Chromosomal deletions were performed in eight of the endoglucanase-encoding genes: cel5A, cel5B, cel5C, cel9A, cel9B, cel9C, cel9E, and cel9F Each mutant retained the ability to digest crystalline cellulose, although the deletion of cel9C caused a modest decrease in cellulose utilization. Strains with multiple deletions were constructed to identify the critical cellulases. Cells of a mutant lacking both cel5B and cel9C were completely deficient in growth on cellulose. Cell fractionation and biochemical analyses indicate that Cel5B and Cel9C are periplasmic nonprocessive endoglucanases. The requirement of periplasmic endoglucanases for cellulose utilization suggests that cellodextrins are transported across the outer membrane during this process. Bioinformatic analyses predict that Cel5A, Cel9A, Cel9B, Cel9D, and Cel9E are secreted across the outer membrane by the type IX secretion system, which has been linked to cellulose utilization. These secreted endoglucanases may perform the initial digestion within amorphous regions on the cellulose fibers, releasing oligomers that are transported into the periplasm for further digestion by Cel5B and Cel9C. The results suggest that both cell surface and periplasmic endoglucanases are required for the growth of C. hutchinsonii on cellulose and that novel cell surface proteins may solubilize and transport cellodextrins across the outer membrane.

IMPORTANCE

The bacterium Cytophaga hutchinsonii digests crystalline cellulose by an unknown mechanism. It lacks processive cellobiohydrolases that are often involved in cellulose digestion. Critical cellulolytic enzymes were identified by genetic analyses. Intracellular (periplasmic) nonprocessive endoglucanases performed an important role in cellulose utilization. The results suggest a model involving partial digestion at the cell surface, solubilization and uptake of cellodextrins across the outer membrane by an unknown mechanism, and further digestion within the periplasm. The ability to sequester cellodextrins and digest them intracellularly may limit losses of soluble cellobiose to other organisms. C. hutchinsonii uses an unusual approach to digest cellulose and is a potential source of novel proteins to increase the efficiency of conversion of cellulose into soluble sugars and biofuels.

iTRAQ-Based Quantitative Proteomics Analysis of Black Rice Grain Development Reveals Metabolic Pathways Associated with Anthocyanin Biosynthesis

BACKGROUND

Black rice (Oryza sativa L.), whose pericarp is rich in anthocyanins (ACNs), is considered as a healthier alternative to white rice. Molecular species of ACNs in black rice have been well documented in previous studies; however, information about the metabolic mechanisms underlying ACN biosynthesis during black rice grain development is unclear.

RESULTS

The aim of the present study was to determine changes in the metabolic pathways that are involved in the dynamic grain proteome during the development of black rice indica cultivar, (Oryza sativa L. indica var. SSP). Isobaric tags for relative and absolute quantification (iTRAQ) MS/MS were employed to identify statistically significant alterations in the grain proteome. Approximately 928 proteins were detected, of which 230 were differentially expressed throughout 5 successive developmental stages, starting from 3 to 20 days after flowering (DAF). The greatest number of differentially expressed proteins was observed on 7 and 10 DAF, including 76 proteins that were upregulated and 39 that were downregulated. The biological process analysis of gene ontology revealed that the 230 differentially expressed proteins could be sorted into 14 functional groups. Proteins in the largest group were related to metabolic process, which could be integrated into multiple biochemical pathways. Specifically, proteins with a role in ACN biosynthesis, sugar synthesis, and the regulation of gene expression were upregulated, particularly from the onset of black rice grain development and during development. In contrast, the expression of proteins related to signal transduction, redox homeostasis, photosynthesis and N-metabolism decreased during grain maturation. Finally, 8 representative genes encoding different metabolic proteins were verified via quantitative real-time polymerase chain reaction (qRT-PCR) analysis, these genes had differed in transcriptional and translational expression during grain development.

CONCLUSIONS

Expression analyses of metabolism-related protein groups belonging to different functional categories and subcategories indicated that significantly upregulated proteins were related to flavonoid and starch synthesis. On the other hand, the downregulated proteins were determined to be related to nitrogen metabolism, as well as other functional categories and subcategories, including photosynthesis, redox homeostasis, tocopherol biosynthetic, and signal transduction. The results provide valuable new insights into the characterization and understanding of ACN pigment production in black rice.

Prediction of Certain Well-Characterized Domains of Known Functions within the PE and PPE Proteins of Mycobacteria

The PE and PPE protein family are unique to mycobacteria. Though the complete genome sequences for over 500 M. tuberculosis strains and mycobacterial species are available, few PE and PPE proteins have been structurally and functionally characterized. We have therefore used bioinformatics tools to characterize the structure and function of these proteins. We selected representative members of the PE and PPE protein family by phylogeny analysis and using structure-based sequence annotation identified ten well-characterized protein domains of known function. Some of these domains were observed to be common to all mycobacterial species and some were species specific.

Carbohydrate-active enzymes: sequences, shapes, contortions and cells

The enzyme-catalysed degradation of oligo and polysaccharides is of considerable interest in many fields ranging from the fundamental–understanding the intrinsic chemical beauty–through to the applied, including diverse practical applications in medicine and biotechnology. Carbohydrates are the most stereochemically-complex biopolymer, and myriad different natural polysaccharides have led to evolution of multifaceted enzyme consortia for their degradation. The glycosidic bonds that link sugar monomers are among the most chemically-stable, yet enzymatically-labile, bonds in the biosphere. That glycoside hydrolases can achieve a rate enhancement (kcat/kuncat) >1017-fold provides testament to their remarkable proficiency and the sophistication of their catalysis reaction mechanisms. The last two decades have seen significant advances in the discovery of new glycosidase sequences, sequence-based classification into families and clans, 3D structures and reaction mechanisms, providing new insights into enzymatic catalysis. New impetus to these studies has been provided by the challenges inherent in plant and microbial polysaccharide degradation, both in the context of environmentally-sustainable routes to foods and biofuels, and increasingly in human nutrition. Study of the reaction mechanism of glycoside hydrolases has also inspired the development of enzyme inhibitors, both as mechanistic probes and increasingly as therapeutic agents. We are on the cusp of a new era where we are learning how to dovetail powerful computational techniques with structural and kinetic data to provide an unprecedented view of conformational details of enzyme action.

Antisense down-regulation of the strawberry β-galactosidase gene FaβGal4 increases cell wall galactose levels and reduces fruit softening

Strawberry softening is characterized by an increase in the solubilization and depolymerization of pectins from cell walls. Galactose release from pectin side chains by β-galactosidase enzymes has been proposed as one reason for the increase in soluble pectins. A putative β-galactosidase gene, FaβGal4, has been identified using a custom-made oligonucleotide-based strawberry microarray platform. FaβGal4 was expressed mainly in the receptacle during fruit ripening, and was positively regulated by abscisic acid and negatively regulated by auxins. To ascertain the role of FaβGal4 in strawberry softening, transgenic plants containing an antisense sequence of this gene under the control of the CaMV35S promoter were generated. Phenotypic analyses were carried out in transgenic plants during three consecutive growing seasons, using non-transformed plants as control. Two out of nine independent transgenic lines yielded fruits that were 30% firmer than control at the ripe stage. FaβGal4 mRNA levels were reduced by 70% in ripe fruits from these selected transgenic lines, but they also showed significant silencing of FaβGal1, although the genes did not share significant similarity. These two transgenic lines also showed an increase in pectin covalently bound to the cell wall, extracted using Na2CO3. The amount of galactose in cell walls from transgenic fruits was 30% higher than in control; notably, the galactose increase was larger in the 1 M KOH fraction, which is enriched in hemicellulose. These results suggest that FaβGal4 participates in the solubilization of covalently bound pectins during ripening, reducing strawberry fruit firmness.

Role of the non-conserved amino acid asparagine 285 in the glycone-binding pocket of Neosartorya fischeri β-glucosidase

Neosartorya fischeriβ-glucosidase (NfBGL595) is distinguished from other BGLs by its high turnover forp-nitrophenyl β-d-glucopyranoside (pNPG) and flavones.

NMR analysis of the binding mode of two fungal endo-β-1,4-mannanases from GH5 and GH26 families

A different mode of action among two endo-β-1,4 mannanases from Podospora anserina has been revealed through an accurate NMR analysis.

An Aspergillus nidulans GH26 endo-β-mannanase with a novel degradation pattern on highly substituted galactomannans

The activity and substrate degradation pattern of a novel Aspergillus nidulans GH26 endo-β-mannanase (AnMan26A) was investigated using two galactomannan substrates with varying amounts of galactopyranosyl residues. The AnMan26A was characterized in parallel with the GH26 endomannanase from Podospora anserina (PaMan26A) and three GH5 endomannanases from A. nidulans and Trichoderma reesei (AnMan5A, AnMan5C and TrMan5A). The initial rates and the maximal degree of enzymatically catalyzed conversion of locust bean gum and guar gum galactomannans were determined. The hydrolysis product profile at maximal degree of conversion was determined using DNA sequencer-Assisted Saccharide analysis in High throughput (DASH). This is the first reported use of this method for analyzing galactomannooligosaccharides. AnMan26A and PaMan26A were found to have a novel substrate degradation pattern on the two galactomannan substrates. On the highly substituted guar gum AnMan26A and PaMan26A reached 35-40% as their maximal degree of conversion whereas the three tested GH5 endomannanases only reached 8-10% as their maximal degree of conversion. α-Galactosyl-mannose was identified as the dominant degradation product resulting from AnMan26A and PaMan26A action on guar gum, strongly indicating that these two enzymes can accommodate galactopyranosyl residues in the -1 and in the +1 subsite. The degradation of α-6(4)-6(3)-di-galactosyl-mannopentaose by AnMan26A revealed accommodation of galactopyranosyl residues in the -2, -1 and +1 subsite of the enzyme. Accommodation of galactopyranosyl residues in subsites -2 and +1 has not been observed for other characterized endomannanases to date. Docking analysis of galactomannooligosaccharides in available crystal structures and homology models supported the conclusions drawn from the experimental results. This newly discovered diversity of substrate degradation patterns demonstrates an expanded functionality of fungal endomannanases, than hitherto reported.

Role of BGT-1 and BGT-2, two predicted GPI-anchored glycoside hydrolases/glycosyltransferases, in cell wall remodeling in Neurospora crassa

Neurospora crassa BGT-1 (NCU06381) and BGT-2 (NCU09175) are two putative glycoside hydrolases (GHs) with additional predicted glycosyltransferase activity and binding sites for a glycosyl phosphatidyl inositol (GPI) anchor that would facilitate their attachment to the plasma membrane (PM). To discern their role in key morphogenetic events during vegetative development of N. crassa, BGT-1 and BGT-2 were labeled with the green fluorescent protein (GFP). The gfp was inserted immediately after the signal peptide sequence, within the bgt-1 encoding sequence, or directly before the GPI-binding site in the case of bgt-2. Both BGT-1-GFP and BGT-2-GFP were observed at the PM of the hyphal apical dome, excluding the foremost apical region and the Spitzenkörper (Spk), where chitin and β-1,3-glucan synthases have been previously found. These and previous studies suggest a division of labor of the cell wall synthesizing machinery at the hyphal dome: at the very tip, glucans are synthesized by enzymes that accumulate at the Spk, before getting incorporated into the PM, whereas at the subtending zone below the apex, glucans are presumably hydrolyzed, producing amenable ends for further branching and crosslinking with other cell wall polymers. Additionally, BGT-1-GFP and BGT-2-GFP were observed at the leading edge of new developing septa, at unreleased interconidial junctions, at conidial poles, at germling and hyphal fusion sites, and at sites of branch emergence, all of them processes that seemingly involve cell wall remodeling. Even though single and double mutant strains for the corresponding genes did not show a drastic reduction of growth rate, bgt-2Δ and bgt-1Δ::bgt-2Δ strains exhibited an increased resistance to the cell wall stressors calcofluor white (CW) and congo red (CR) than the reference strain, which suggests they present significant architectural changes in their cell wall. Furthermore, the conidiation defects observed in the mutants indicate a significant role of BGT-1 and BGT-2 on the re-arrangement of glucans needed at the conidiophore cell wall to allow conidial separation.

Hooked on α-d-galactosidases: from biomedicine to enzymatic synthesis

α-d-Galactosidases (EC 3.2.1.22) are enzymes employed in a number of useful bio-based applications. We have depicted a comprehensive general survey of α-d-galactosidases from different origin with special emphasis on marine example(s). The structures of natural α-galactosyl containing compounds are described. In addition to 3D structures and mechanisms of action of α-d-galactosidases, different sources, natural function and genetic regulation are also covered. Finally, hydrolytic and synthetic exploitations as free or immobilized biocatalysts are reviewed. Interest in the synthetic aspects during the next years is anticipated for access to important small molecules by green technology with an emphasis on alternative selectivity of this class of enzymes from different sources.

Fucoidanases

In recent years, the research of fucoidans has steadily increased. The interest in these substances is due to their various biological activities. Despite a wide range of biological activity and the lack of oral toxicity, fucoidans remain relatively unexploited as a source of medicines because of their heterogeneity. Enzymes that degrade polyanionic polysaccharides are widely used for establishing their structures and structure-activity relationships. Sometimes, to obtain preparations of polysaccharides with standard characteristics, for example, medicines and food supplements, enzymatic treatment can be also applied. Only a few sources of enzymes with fucoidanase activity have been described, and only a few studies regarding the isolation and characterization of fucoidanases have been performed. The data on the specificity of fucoidanases: the type of cleaved glycoside bond, the relation between catalytic activity and the degree of substrate sulfation are scarce. The review summarizes achievements in the research of fucoidanases, mechanisms of enzymatic degradation of fucoidans, as well as of structures of sulfated fucooligosaccharides obtained under the action of fucoidanases.

Cloning and Genomic Organization of a Rhamnogalacturonase Gene from Locally Isolated Strain of Aspergillus niger

The rhg gene encoding a rhamnogalacturonase was isolated from the novel strain A1 of Aspergillus niger. It consists of an ORF of 1.505 kb encoding a putative protein of 446 amino acids with a predicted molecular mass of 47 kDa, belonging to the family 28 of glycosyl hydrolases. The nature and position of amino acids comprising the active site as well as the three-dimensional structure were well conserved between the A. niger CTM10548 and fungal rhamnogalacturonases. The coding region of the rhg gene is interrupted by three short introns of 56 (introns 1 and 3) and 52 (intron 2) bp in length. The comparison of the peptide sequence with A. niger rhg sequences revealed that the A1 rhg should be an endo-rhamnogalacturonases, more homologous to rhg A than rhg B A. niger known enzymes. The comparison of rhg nucleotide sequence from A. niger A1 with rhg A from A. niger shows several base changes. Most of these changes (59 %) are located at the third base of codons suggesting maintaining the same enzyme function. We used the rhamnogalacturonase A from Aspergillus aculeatus as a template to build a structural model of rhg A1 that adopted a right-handed parallel β-helix.

Modeling of cooked starch digestion process using recombinant human pancreatic α-amylase and maltase-glucoamylase for in vitro evaluation of α-glucosidase inhibitors

In human, digestion of cooked starch mainly involves breaking down of α-amylase to α-limit dextrins and small linear malto-oligosaccharides, which are in turn hydrolyzed to glucose by the gut mucosal maltase-glucoamylase (MGAM). Human pancreatic α-amylase (HPA), amino- and carboxyl-terminal portions of MGAM (ntMGAM and ctMGAM) catalyze the hydrolysis of α-D-(1,4) glycosidic linkages in starch, playing a crucial role in the production of glucose in the human lumen. Accordingly, these enzymes are effective drug targets for the treatments of type 2 diabetes and obesity. In this study, a Plackett-Burman based statistical screening procedure was adopted to determine the most critical factors affecting cooked starch digestion by the combination of HPA, ctMGAM and ntMGAM. Six factors were tested and experimental results showed that pH and temperature were the major influencing factors, with optimal pH and temperature at 6.0 and 50 °C, respectively. Surprisingly, ntMGAM had no significant contribution to the glucose production from starch digestion compared to the HPA and ctMGAM. The optimal proportion of HPA and ctMGAM in a starch digestion system was further determined by response surface methodology. Results showed a maximum starch digestion (88.05%) within 0.5 h when used HPA:ctMGAM=1:9 (U). The inhibitory effects of various inhibitors on the cooked starch digestion by HPA1/ctMGAM9 were evaluated by determining their half maximal inhibitory concentration (IC50) values. Acarviostatin II03 showed the highest inhibitory activity, with 67 times higher potency than acarbose. Moreover, acarviostatin II03 could significantly depress postprandial blood glucose levels in mice, better than that by acarbose. These findings suggest that our in vitro enzymatic system can simulate in vivo starch digestion process, and thus can be used to screen and evaluate α-glucosidase inhibitors.

Towards New Drug Targets? Function Prediction of Putative Proteins of Neisseria meningitidis MC58 and Their Virulence Characterization

Neisseria meningitidis is a Gram-negative aerobic diplococcus, responsible for a variety of meningococcal diseases. The genome of N. meningitidis MC58 is comprised of 2114 genes that are translated into 1953 proteins. The 698 genes (∼35%) encode hypothetical proteins (HPs), because no experimental evidence of their biological functions are available. Analyses of these proteins are important to understand their functions in the metabolic networks and may lead to the discovery of novel drug targets against the infections caused by N. meningitidis. This study aimed at the identification and categorization of each HP present in the genome of N. meningitidis MC58 using computational tools. Functions of 363 proteins were predicted with high accuracy among the annotated set of HPs investigated. The reliably predicted 363 HPs were further grouped into 41 different classes of proteins, based on their possible roles in cellular processes such as metabolism, transport, and replication. Our studies revealed that 22 HPs may be involved in the pathogenesis caused by this microorganism. The top two HPs with highest virulence scores were subjected to molecular dynamics (MD) simulations to better understand their conformational behavior in a water environment. We also compared the MD simulation results with other virulent proteins present in N. meningitidis. This study broadens our understanding of the mechanistic pathways of pathogenesis, drug resistance, tolerance, and adaptability for host immune responses to N. meningitidis.

Identification and Characterization of Two Endogenous β-Glucosidases from the Termite Coptotermes formosanus

Coptotermes formosanus is a well-known wood-feeding termite that can degrade lignocellulose polysaccharides efficiently with its unique multi-enzyme catalysis system. β-glucosidase (BG) is one of the important cellulases in its enzyme system. However, there may present multiple endogenous BGs in termite digestive system for various properties and functions. This study aims to characterize two BG homologs and reveal their potential coordinative effect. In this study, two endogenous BG homologs (CfGlu1B and CfGlu1C) from C. formosanus showed different substrate specificity. CfGlu1B favors cellobiose while CfGlu1C favors sucrose. Besides, CfGlu1C exhibited higher alkali resistance than CfGlu1B. Kinetic analysis revealed that CfGlu1B enzyme's activity toward p-NP-β-D-glucopyranoside (p-NPG) was higher than that of CfGlu1C, and the difference mainly attributes to the turnover number (K cat). In addition, the activity assay showed significant synergistic action of CfGlu1B and CfGlu1C in degrading D-lactosum. For effect of metals, Cu(2+) inhibited both enzymes and Ca(2+) increased the activity of CfGlu1C but not CfGlu1B. Site-directed mutagenesis analysis indicated that both enzymes lost activities when residues E190 of CfGlu1B and E168 of CfGlu1C were mutated to alanine, respectively, which were essential active centers of the GHF1 enzymes. Moreover, mutation H252N increased the activity of enzyme CfGlu1C by 2.1-fold. This study implies interesting possibilities for better practical biotechnological use in green energy production.

Exchange of active site residues alters substrate specificity in extremely thermostable β-glycosidase from Thermococcus kodakarensis KOD1

β-Glycosidase from Thermococcus kodakarensis KOD1 is a hyperthermophilic enzyme with β-glucosidase, β-mannosidase, β-fucosidase and β-galactosidase activities. Sequence alignment with other β-glycosidases from hyperthermophilic archaea showed two unique active site residues, Gln77 and Asp206. These residues were represented by Arg and Asp in all other hyperthermophilic β-glycosidases. The two active site residues were mutated to Q77R, D206N and D206Q, to study the role of these unique active site residues in catalytic activity and to alter the substrate specificity to enhance its β-glucosidase activity. The secondary structure analysis of all the mutants showed no change in their structure and exhibited in similar conformation like wild-type as they all existed in dimer form in an SDS-PAGE under non-reducing conditions. Q77R and D206Q affected the catalytic activity of the enzyme whereas the D206N altered the catalytic turn-over rate for glucosidase and mannosidase activities with fucosidase activity remain unchanged. Gln77 is reported to interact with catalytic nucleophile and Asp206 with axial C2-hydroxyl group of substrates. Q77R might have made some changes in three dimensional structure due to its electrostatic effect and lost its catalytic activity. The extended side chains of D206Q is predicted to affect the substrate binding during catalysis. The high-catalytic turn-over rate by D206N for β-glucosidase activity makes it a useful enzyme in cellulose degradation at high temperatures.

Differential expression of α-L-arabinofuranosidases during maize (Zea mays L.) root elongation

Specific α- l -arabinofuranosidases are involved in the realisation of elongation growth process in cells with type II cell walls. Elongation growth in a plant cell is largely based on modification of the cell wall. In type II cell walls, the Ara/Xyl ratio is known to decrease during elongation due to the partial removal of Ara residues from glucuronoarabinoxylan. We searched within the maize genome for the genes of all predicted α-L-arabinofuranosidases that may be responsible for such a process and related their expression to the activity of the enzyme and the amount of free arabinose measured in six zones of a growing maize root. Eight genes of the GH51 family (ZmaABFs) and one gene of the GH3 family (ZmaARA-I) were identified. The abundance of ZmaABF1 and 3-6 transcripts was highly correlated with the measured enzymatic activity and free arabinose content that significantly increased during elongation. The transcript abundances also coincided with the pattern of changes in the Ara/Xyl ratio of the xylanase-extractable glucuronoarabinoxylan described in previous studies. The expression of ZmaABF3, 5 and 6 was especially up-regulated during elongation although corresponding proteins are devoid of the catalytic glutamate at the proper position. ZmaABF2 transcripts were specifically enriched in the root cap and meristem. A single ZmaARA-I gene was not expressed as a whole gene but instead as splice variants that encode the C-terminal end of the protein. Changes in the ZmaARA-I transcript level were rather moderate and had no significant correlation with free arabinose content. Thus, elongation growth of cells with type II cell walls is accompanied by the up-regulation of specific and predicted α-L-arabinofuranosidase genes, and the corresponding activity is indeed pronounced and is important for the modification of glucuronoarabinoxylan, which plays a key role in the modification of the cell wall supramolecular organisation.

Family 46 Carbohydrate-binding Modules Contribute to the Enzymatic Hydrolysis of Xyloglucan and β-1,3-1,4-Glucans through Distinct Mechanisms

Structural carbohydrates comprise an extraordinary source of energy that remains poorly utilized by the biofuel sector as enzymes have restricted access to their substrates within the intricacy of plant cell walls. Carbohydrate active enzymes (CAZYmes) that target recalcitrant polysaccharides are modular enzymes containing noncatalytic carbohydrate-binding modules (CBMs) that direct enzymes to their cognate substrate, thus potentiating catalysis. In general, CBMs are functionally and structurally autonomous from their associated catalytic domains from which they are separated through flexible linker sequences. Here, we show that a C-terminal CBM46 derived from BhCel5B, a Bacillus halodurans endoglucanase, does not interact with β-glucans independently but, uniquely, acts cooperatively with the catalytic domain of the enzyme in substrate recognition. The structure of BhCBM46 revealed a β-sandwich fold that abuts onto the region of the substrate binding cleft upstream of the active site. BhCBM46 as a discrete entity is unable to bind to β-glucans. Removal of BhCBM46 from BhCel5B, however, abrogates binding to β-1,3-1,4-glucans while substantially decreasing the affinity for decorated β-1,4-glucan homopolymers such as xyloglucan. The CBM46 was shown to contribute to xyloglucan hydrolysis only in the context of intact plant cell walls, but it potentiates enzymatic activity against purified β-1,3-1,4-glucans in solution or within the cell wall. This report reveals the mechanism by which a CBM can promote enzyme activity through direct interaction with the substrate or by targeting regions of the plant cell wall where the target glucan is abundant.

Glandular β-glucosidases in juvenile Chrysomelina leaf beetles support the evolution of a host-plant-dependent chemical defense

Plant-feeding insects are spread across the entire plant kingdom. Because they chew externally on leaves, leaf beetle of the subtribe Chrysomelina sensu stricto are constantly exposed to life-threatening predators and parasitoids. To counter these pressures, the juveniles repel their enemies by displaying glandular secretions that contain defensive compounds. These repellents can be produced either de novo (iridoids) or by using plant-derived precursors. The autonomous production of iridoids pre-dates the evolution of phytochemical-based defense strategies. Both strategies include hydrolysis of the secreted non-toxic glycosides in the defensive exudates. By combining in vitro as well as in vivo experiments, we show that iridoid de novo producing as well as sequestering species rely on secreted β-glucosidases to cleave the pre-toxins. Our phylogenetic analyses support a common origin of chrysomeline β-glucosidases. The kinetic parameters of these β-glucosidases demonstrated substrate selectivity which reflects the adaptation of Chrysomelina sensu stricto to the chemistry of their hosts during the course of evolution. However, the functional studies also showed that the broad substrate selectivity allows building a chemical defense, which is dependent on the host plant, but does not lead to an "evolutionary dead end".

An update on Acanthamoeba keratitis: diagnosis, pathogenesis and treatment

Free-living amoebae of the genus Acanthamoeba are causal agents of a severe sight-threatening infection of the cornea known as Acanthamoeba keratitis. Moreover, the number of reported cases worldwide is increasing year after year, mostly in contact lens wearers, although cases have also been reported in non-contact lens wearers. Interestingly, Acanthamoeba keratitis has remained significant, despite our advances in antimicrobial chemotherapy and supportive care. In part, this is due to an incomplete understanding of the pathogenesis and pathophysiology of the disease, diagnostic delays and problems associated with chemotherapeutic interventions. In view of the devastating nature of this disease, here we present our current understanding of Acanthamoeba keratitis and molecular mechanisms associated with the disease, as well as virulence traits of Acanthamoeba that may be potential targets for improved diagnosis, therapeutic interventions and/or for the development of preventative measures. Novel molecular approaches such as proteomics, RNAi and a consensus in the diagnostic approaches for a suspected case of Acanthamoeba keratitis are proposed and reviewed based on data which have been compiled after years of working on this amoebic organism using many different techniques and listening to many experts in this field at conferences, workshops and international meetings. Altogether, this review may serve as the milestone for developing an effective solution for the prevention, control and treatment of Acanthamoeba infections.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.