Crystal structure of the Chlamydomonas starch debranching enzyme isoamylase ISA1 reveals insights into the mechanism of branch trimming and complex assembly

The mechanism of carbon and nitrogen metabolism under low temperature and low light stress after heading in late indica rice

Double-season late indica rice frequently experiences low temperature accompanying with low light stress during the grain filling stage in southern China, which alters the carbon and nitrogen metabolism of the rice grains, thereby impacting both grain yield and quality. However, the physiological mechanism is still unclear. A pot experiment using two late indica rice cultivars (high-quality and common-quality rice cultivars) was conducted under control (CK), low temperature (LT) and low temperature and light (LT + LL) to investigate the grain filling, photosynthetic characteristics, carbon and nitrogen metabolic enzymes and related gene expression. The results indicated that both LT and LT + LL treatments primarily reduced grain weight in the two cultivars compared to CK. This effect was particularly pronounced for LT + LL, which showed a significant difference in both superior and inferior grains of the common-quality rice cultivar. This reduction in grain weight was attributed to an average decrease in the photosynthetic rate of rice leaves by 23.5%, a decrease in Fv/Fm by 6.2%, and a decrease in the original grain filling rate by 10.3%. LT and LT + LL treatments decreased the activities of soluble starch synthase (SSS), granules bound starch synthetase (GBSS) and starch branch enzyme (SBE) in the early grain filling stage, and its related gene expression including OsAGPL2, OsAGPS2b, OsGBSSI and OsSSIIIa, and the influence degree was intensified by LT + LL, compared with LT. For nitrogen metabolism, the activity of glutamic oxalic aminotransferase (GOT) also significantly decreased by LT and LT + LL during grain filling in the cultivars, which was relatively lower under LT + LL, and only decreased the glutamic pyruvic transaminase (GPT) in the early grain filling stage but increased in the later period. However, the related gene expressions of OsGOT1B and OsGS1;3 were enhanced significantly by LT and LT + LL treatments, and reached the highest level under LT + LL treatment, implying the disordered process of nitrogen metabolism. The results suggest that low temperature decreased photosynthetic traits to hinder grain filling in the cultivars, mainly derived from the worsened carbon and nitrogen metabolism, and LT + LL combined stress aggravated the damage degree.

Carbohydrate binding modules: Compact yet potent accessories in the specific substrate binding and performance evolution of carbohydrate-active enzymes

Carbohydrate binding modules (CBMs) are independent non-catalytic domains widely found in carbohydrate-active enzymes (CAZymes), and they play an essential role in the substrate binding process of CAZymes by guiding the appended catalytic modules to the target substrates. Owing to their precise recognition and selective affinity for different substrates, CBMs have received increasing research attention over the past few decades. To date, CBMs from different origins have formed a large number of families that show a variety of substrate types, structural features, and ligand recognition mechanisms. Moreover, through the modification of specific sites of CBMs and the fusion of heterologous CBMs with catalytic domains, improved enzymatic properties and catalytic patterns of numerous CAZymes have been achieved. Based on cutting-edge technologies in computational biology, gene editing, and protein engineering, CBMs as auxiliary components have become portable and efficient tools for the evolution and application of CAZymes. With the aim to provide a theoretical reference for the functional research, rational design, and targeted utilization of novel CBMs in the future, we systematically reviewed the function-related characteristics and potentials of CAZyme-derived CBMs in this review, including substrate recognition and binding mechanisms, non-catalytic contributions to enzyme performances, module modifications, and innovative applications in various fields.

Origin and evolution of the main starch biosynthetic enzymes

Starch, a semi-crystalline energy storage form primarily found in plant plastids plays a crucial role in various food or no-food applications. Despite the starch biosynthetic pathway's main enzymes have been characterized, their origin and evolution remained a subject of debate. In this study, we conducted the comprehensive phylogenetic and structural analysis of three types of starch biosynthetic enzymes: starch synthase (SS), starch branching enzyme (SBE) and isoamylase-type debranching enzyme (ISA) from 51,151 annotated genomes. Our findings provide valuable insights into the possible scenario for the origin and evolution of the starch biosynthetic pathway. Initially, the ancestor of SBE can be traced back to an unidentified bacterium that existed before the formation of the last eukaryotic common ancestor (LECA) via horizontal gene transfer (HGT). This transfer event likely provided the eukaryote ancestor with the ability to synthesize glycogen. Furthermore, during the emergence of Archaeplastida, one clade of SS was transferred from Deltaproteobacteria by HGT, while ISA and the other clade of SS originated from Chlamydiae through endosymbiosis gene transfer (EGT). Both these transfer events collectively contributed to the establishment of the original starch biosynthetic pathway. Subsequently, after the divergence of Viridiplantae from Rhodophyta, all three enzymes underwent multiple duplications and N-terminus extension domain modifications, resulting in the formation of functionally specialized isoforms and ultimately leading to the complete starch biosynthetic pathway. By shedding light on the evolutionary origins of key enzymes involved in the starch biosynthetic pathway, this study provides important insights into the evolutionary events of plants.

Distribution of Aromatic Amino Acid Residues in Substrate-Binding Regions Modulates Substrate Specificity of Microbial Debranching Enzymes

Debranching enzymes (DBEs) directly hydrolyze α-1,6-glucosidic linkages in glycogen, starch, and related polysaccharides, making them important in the starch processing industry. However, the ambiguous substrate specificity usually restricts synergistic catalysis with other amylases for improving starch utilization. Herein, a glycogen-debranching enzyme from Saccharolobus solfataricus (SsGDE) and two isoamylases from Pseudomonas amyloderamosa (PaISO) and Chlamydomonas reinhardtii (CrISO) were used to investigate the molecular mechanism of substrate specificity. Along with the structure-based computational analysis, the aromatic residues in the substrate-binding region of DBEs played an important role in binding substrates. The aromatic residues in SsGDE appeared clustered, contributing to a small substrate-binding region. In contrast, the aromatic residues in isoamylase were distributed dispersedly, forming a large active site. The distinct characteristics of substrate-binding regions in SsGDE and isoamylase might explain their substrate preferences for maltodextrin and amylopectin, respectively. By modulating the substrate-binding region of SsGDE, variants Y323F and V375F were obtained with significantly enhanced activities, and the activities of Y323F and V375F increased by 30 and 60% for amylopectin, and 20 and 23% for DE4 maltodextrin, respectively. This study revealed the molecular mechanisms underlying the substrate specificity for SsGDE and isoamylases, providing a route for engineering enzymes to achieve higher catalytic performance.

Structure and genetic regulation of starch formation in sorghum (Sorghum bicolor (L.) Moench) endosperm: A review

This review focuses on the structure and genetic regulation of starch formation in sorghum (Sorghum bicolor (L.) Moench) endosperm. Sorghum is an important cereal crop that is well suited to grow in regions with high temperatures and limited water resources due to its C4 metabolism. The endosperm of sorghum kernels is a rich source of starch, which is composed of two main components: amylose and amylopectin. The synthesis of starch in sorghum endosperm involves multiple enzymatic reactions, which are regulated by complex genetic and environmental factors. Recent research has identified several genes involved in the regulation of starch synthesis in sorghum endosperm. In addition, the structure and properties of sorghum starch can also be influenced by environmental factors such as temperature, water availability, and soil nutrients. A better understanding of the structure and genetic regulation of starch formation in sorghum endosperm can have important implications for the development of sorghum-based products with improved quality and nutritional value. This review provides a comprehensive summary of the current knowledge on the structure and genetic regulation of starch formation in sorghum endosperm and highlights the potential for future research to further improve our understanding of this important process.

FLOURY ENDOSPERM 6 mutations enhance the sugary phenotype caused by the loss of ISOAMYLASE1 in barley

Barley double mutants in two genes involved in starch granule morphology, HvFLO6 and HvISA1, had impaired starch accumulation and higher grain sugar levels than either single mutant. Starch is a biologically and commercially important glucose polymer synthesized by plants as semicrystalline starch granules (SGs). Because SG morphology affects starch properties, mutants with altered SG morphology may be useful in breeding crops with desirable starch properties, including potentially novel properties. In this study, we employed a simple screen for mutants with altered SG morphology in barley (Hordeum vulgare). We isolated mutants that formed compound SGs together with the normal simple SGs in the endosperm and found that they were allelic mutants of the starch biosynthesis genes ISOAMYLASE1 (HvISA1) and FLOURY ENDOSPERM 6 (HvFLO6), encoding starch debranching enzyme and CARBOHYDRATE-BINDING MODULE 48-containing protein, respectively. We generated the hvflo6 hvisa1 double mutant and showed that it had significantly reduced starch biosynthesis and developed shrunken grains. In contrast to starch, soluble α-glucan, phytoglycogen, and sugars accumulated to higher levels in the double mutant than in the single mutants. In addition, the double mutants showed defects in SG morphology in the endosperm and in the pollen. This novel genetic interaction suggests that hvflo6 acts as an enhancer of the sugary phenotype caused by hvisa1 mutation.

An analysis of sugary endosperm in sorghum: Characterization of mutant phenotypes depending on alleles of the corresponding starch debranching enzyme

Sorghum is the fifth most important cereal crop. Here we performed molecular genetic analyses of the 'SUGARY FETERITA' (SUF) variety, which shows typical sugary endosperm traits (e.g., wrinkled seeds, accumulation of soluble sugars, and distorted starch). Positional mapping indicated that the corresponding gene was located on the long arm of chromosome 7. Within the candidate region of 3.4 Mb, a sorghum ortholog for maize Su1 (SbSu) encoding a starch debranching enzyme ISA1 was found. Sequencing analysis of SbSu in SUF uncovered nonsynonymous single nucleotide polymorphisms (SNPs) in the coding region, containing substitutions of highly conserved amino acids. Complementation of the rice sugary-1 (osisa1) mutant line with the SbSu gene recovered the sugary endosperm phenotype. Additionally, analyzing mutants obtained from an EMS-induced mutant panel revealed novel alleles with phenotypes showing less severe wrinkles and higher Brix scores. These results suggested that SbSu was the corresponding gene for the sugary endosperm. Expression profiles of starch synthesis genes during the grain-filling stage demonstrated that a loss-of-function of SbSu affects the expression of most starch synthesis genes and revealed the fine-tuned gene regulation in the starch synthetic pathway in sorghum. Haplotype analysis using 187 diverse accessions from a sorghum panel revealed the haplotype of SUF showing severe phenotype had not been used among the landraces and modern varieties. Thus, weak alleles (showing sweet and less severe wrinkles), such as in the abovementioned EMS-induced mutants, are more valuable for grain sorghum breeding. Our study suggests that more moderate alleles (e.g. produced by genome editing) should be beneficial for improving grain sorghum.

Crystal structure of a novel homodimeric D-allulose 3-epimerase from a Clostridia bacterium

D-Allulose, a low-calorie rare sugar with various physiological functions, is mainly produced through the isomerization of D-fructose by ketose 3-epimerases (KEases), which exhibit various substrate specificities. A novel KEase from a Clostridia bacterium (CDAE) was identified to be a D-allulose 3-epimerase and was further characterized as thermostable and metal-dependent. In order to explore its structure–function relationship, the crystal structure of CDAE was determined using X-ray diffraction at 2.10 Å resolution, revealing a homodimeric D-allulose 3-epimerase structure with extensive interactions formed at the dimeric interface that contribute to structure stability. Structural analysis identified the structural features of CDAE, which displays a common (β/α)8-TIM barrel and an ordered Mn2+-binding architecture at the active center, which may explain the positive effects of Mn2+ on the activity and stability of CDAE. Furthermore, comparison of CDAE and other KEase structures revealed several structural differences, highlighting the remarkable differences in enzyme–substrate binding at the O4, O5 and O6 sites of the bound substrate, which are mainly induced by distinct hydrophobic pockets in the active center. The shape and hydrophobicity of this pocket appear to produce the differences in specificity and affinity for substrates among KEase family enzymes. Exploration of the crystal structure of CDAE provides a better understanding of its structure–function relationship, which might provide a basis for molecular modification of CDAE and further provides a reference for other KEases.

Starch synthase II plays a crucial role in starch biosynthesis and the formation of multienzyme complexes in cassava storage roots

Starch is a glucose polymer synthesized by green plants for energy storage and is crucial for plant growth and reproduction. The biosynthesis of starch polysaccharides is mediated by members of the large starch synthase (SS) protein superfamily. Here, we showed that in cassava storage roots, soluble starch synthase II (MeSSII) plays an important role in starch biosynthesis and the formation of protein complexes with other starch biosynthetic enzymes by directly interacting with MeSSI, MeSBEII, and MeISAII. MeSSII-RNAi cassava lines showed increased amylose content and reduced biosynthesis of the intermediate chain of amylopectin (B1 type) in their storage roots, leading to altered starch physicochemical properties. Furthermore, gel permeation chromatography analysis of starch biosynthetic enzymes between wild type and MeSSII-RNAi lines confirmed the key role of MeSSII in the organization of heteromeric starch synthetic protein complexes. The lack of MeSSII in cassava also reduced the capacity of MeSSI, MeSBEII, MeISAI, and MeISAII to bind to starch granules. These findings shed light on the key components of the starch biosynthesis machinery in root crops.

The amino acid on the top of the active groove allosterically modulates product specificity of the 1,4-α-glucan branching enzyme

The 1,4-α-glucan branching enzymes (GBEs, EC 2.4.1.18) catalyse the formation of α-1,6 branching points in starch, presenting several potential applications in modifying starch. Previous study proved that W285 is considered to act as a "switch" to stop extension of substrates in the structure of GBE from Cyanothece sp. (cceBE). In the structure of GBE from Rhodothermus obamensis STB05 (RoGBE), the amino acid 160 site is structurally similar to the W285 in cceBE. In order to explore the role of this site in RoGBE, several engineered mutants individually substituted with Arg, Phe and Ala at G160 were studied in our research. The results show that substitution with Arg and Phe increased branching activity significantly, and the ratio of short glucan chains among all oligosaccharides increased. Finally, we proposed that the G160 is a 'door model' to elucidate introduced mutagenesis that triggers and controls the length of binding glucan chains of starch.

Production and structure prediction of amylases from Chlorella vulgaris

Amylases are enzymes required for starch degradation and are naturally produced by many microorganisms. These enzymes are used in several fields such as food processing, beverage, and medicine as well as in the formulation of enzymatic detergents proving their significance in modern biotechnology. In this study, a three-stage growth mode was applied to enhance starch production and amylase detection from Chlorella vulgaris. Stress conditions applied in the second stage of cultivation led to an accumulation of proteins (75% DW) and starch (21% DW) and a decrease in biomass. Amylase activities were detected and they showed high production levels especially on day 3 (35 U/ml) and day 5 (22.5 U/ml) of the second and third stages, respectively. The bioinformatic tools used to seek amylase protein sequences from TSA database of C. vulgaris revealed 7 putative genes encoding for 4 α-amylases, 2 β-amylases, and 1 isoamylase. An in silico investigation showed that these proteins are different in their lengths as well as in their cellular localizations and oligomeric states though they share common features like CSRs of GH13 family or active site of GH14 family. In brief, this study allowed for the production and in silico characterization of amylases from C. vulgaris.

Allelic variation in sugary1 gene affecting kernel sweetness among diverse-mutant and -wild-type maize inbreds

Sweet corn is popular worldwide as vegetable. Though large numbers of sugary1 (su1)-based sweet corn germplasm are available, allelic diversity in su1 gene encoding SU1 isoamylase among diverse maize inbreds has not been analyzed. Here, we characterized the su1 gene in maize and compared with allied species. The entire su1 gene (11,720 bp) was sequenced among six mutant (su1) and five wild (Su1) maize inbreds. Fifteen InDels of 2-45 bp were selected to develop markers for studying allelic diversity in su1 gene among 19 mutant- (su1) and 29 wild-type (Su1) inbreds. PIC ranged from 0.15 (SU-InDel7) to 0.37 (SU-InDel13). Major allele frequency varied from 0.52 to 0.90, while gene diversity ranged from 0.16 to 0.49. Phylogenetic tree categorized 48 maize inbreds in two clusters each for wild- type (Su1) and mutant (su1) types. 44 haplotypes of su1 were observed, with three haplotypes (Hap6, Hap22 and Hap29) sharing more than one genotype. Further, comparisons were made with 23 orthologues of su1 from 16 grasses and Arabidopsis. Maize possessed 15-19 exons in su1, while it was 11-24 exons among orthologues. Introns among the orthologues were longer (77-2206 bp) than maize (859-1718 bp). SU1 protein of maize and orthologues had conserved α-amylase and CBM_48 domains. The study also provided physicochemical properties and secondary structure of SU1 protein in maize and its orthologues. Phylogenetic analysis showed closer relationship of maize SU1 protein with P. hallii, S. bicolor and E. tef than Triticum sp. and Oryza sp. The study showed that presence of high allelic diversity in su1 gene which can be utilized in the sweet corn breeding program. This is the first report of comprehensive characterization of su1 gene and its allelic forms in diverse maize and related orthologues.

Microbial starch debranching enzymes: Developments and applications

Starch debranching enzymes (SDBEs) hydrolyze the α-1,6 glycosidic bonds in polysaccharides such as starch, amylopectin, pullulan and glycogen. SDBEs are also important enzymes for the preparation of sugar syrup, resistant starch and cyclodextrin. As the synergistic catalysis of SDBEs and other starch-acting hydrolases can effectively improve the raw material utilization and production efficiency during starch processing steps such as saccharification and modification, they have attracted substantial research interest in the past decades. The substrate specificities of the two major members of SDBEs, pullulanases and isoamylases, are quite different. Pullulanases generally require at least two α-1,4 linked glucose units existing on both sugar chains linked by the α-1,6 bond, while isoamylases require at least three units of α-1,4 linked glucose. SDBEs mainly belong to glycoside hydrolase (GH) family 13 and 57. Except for GH57 type II pullulanse, GH13 pullulanases and isoamylases share plenty of similarities in sequence and structure of the core catalytic domains. However, the N-terminal domains, which might be one of the determinants contributing to the substrate binding of SDBEs, are distinct in different enzymes. In order to overcome the current defects of SDBEs in catalytic efficiency, thermostability and expression level, great efforts have been made to develop effective enzyme engineering and fermentation strategies. Herein, the diverse biochemical properties and distinct features in the sequence and structure of pullulanase and isoamylase from different sources are summarized. Up-to-date developments in the enzyme engineering, heterologous production and industrial applications of SDBEs is also reviewed. Finally, research perspective which could help understanding and broadening the applications of SDBEs are provided.

Enhancing carbohydrate repartitioning into lipid and carotenoid by disruption of microalgae starch debranching enzyme

Light/dark cycling is an inherent condition of outdoor microalgae cultivation, but is often unfavorable for lipid accumulation. This study aims to identify promising targets for metabolic engineering of improved lipid accumulation under outdoor conditions. Consequently, the lipid-rich mutant Chlamydomonas sp. KOR1 was developed through light/dark-conditioned screening. During dark periods with depressed CO₂ fixation, KOR1 shows rapid carbohydrate degradation together with increased lipid and carotenoid contents. KOR1 was subsequently characterized with extensive mutation of the ISA1 gene encoding a starch debranching enzyme (DBE). Dynamic time-course profiling and metabolomics reveal dramatic changes in KOR1 metabolism throughout light/dark cycles. During light periods, increased flux from CO₂ through glycolytic intermediates is directly observed to accompany enhanced formation of small starch-like particles, which are then efficiently repartitioned in the next dark cycle. This study demonstrates that disruption of DBE can improve biofuel production under light/dark conditions, through accelerated carbohydrate repartitioning into lipid and carotenoid.

A Review of Starch Biosynthesis in Relation to the Building Block-Backbone Model

Starch is a water-insoluble polymer of glucose synthesized as discrete granules inside the stroma of plastids in plant cells. Starch reserves provide a source of carbohydrate for immediate growth and development, and act as long term carbon stores in endosperms and seed tissues for growth of the next generation, making starch of huge agricultural importance. The starch granule has a highly complex hierarchical structure arising from the combined actions of a large array of enzymes as well as physicochemical self-assembly mechanisms. Understanding the precise nature of granule architecture, and how both biological and abiotic factors determine this structure is of both fundamental and practical importance. This review outlines current knowledge of granule architecture and the starch biosynthesis pathway in relation to the building block-backbone model of starch structure. We highlight the gaps in our knowledge in relation to our understanding of the structure and synthesis of starch, and argue that the building block-backbone model takes accurate account of both structural and biochemical data.

Functional Analysis of Starch Metabolism in Plants

In plants, starch is synthesized in leaves during the day-time from fixed carbon through photosynthesis and is mobilized at night to support continued respiration, sucrose export, and growth in the dark. The main crops where starch is biosynthesized and stored are corn, rice, wheat, and potatoes, and they are mainly used as food resources for humankind. There are many genes that are involved in starch biosynthesis from cytosol to storage organs in plants. ADP-glucose, UDP- glucose, and glucose-6-phosphate are synthesized catalyzed by UDP-invertase, AGPase, hexokinase, and P- hexose-isomerase in cytosol. Starch composed of amylopectin and amylose is synthesized by starch synthase, granule bound starch synthase, starch-branching enzyme, debranching enzyme, and pullulanase, which is primarily responsible for starch production in storage organs. Recently, it has been uncovered that structural genes are controlled by proteins derived from other genes such as transcription factors. To obtain more precise information on starch metabolism, the functions of genes and transcription factors need to be studied to understand their roles and functions in starch biosynthesis in plants. However, the roles of genes related to starch biosynthesis are not yet clearly understood. The papers of this special issue contain reviews and research articles on these topics and will be a useful resource for researchers involved in the quality improvement of starch storage crops.

Theoretical and experimental approaches to understand the biosynthesis of starch granules in a physiological context

Starch, a plant-derived insoluble carbohydrate composed of glucose polymers, is the principal carbohydrate in our diet and a valuable raw material for industry. The properties of starch depend on the arrangement of glucose units within the constituent polymers. However, key aspects of starch structure and the underlying biosynthetic processes are not well understood, limiting progress towards targeted improvement of our starch crops. In particular, the major component of starch, amylopectin, has a complex three-dimensional, branched architecture. This architecture stems from the combined actions of a multitude of enzymes, each having broad specificities that are difficult to capture experimentally. In this review, we reflect on experimental approaches and limitations to decipher the enzymes' specificities and explore possibilities for in silico simulations of these activities. We believe that the synergy between experimentation and simulation is needed for the correct interpretation of experimental data and holds the potential to greatly advance our understanding of the overall starch biosynthetic process. We furthermore propose that the formation of glucan secondary structures, concomitant with its synthesis, is a previously overlooked factor that directly affects amylopectin architecture through its impact on enzyme function.

Structure-function analysis of silkworm sucrose hydrolase uncovers the mechanism of substrate specificity in GH13 subfamily 17 exo-α-glucosidases

The domestic silkworm Bombyx mori expresses two sucrose-hydrolyzing enzymes, BmSUH and BmSUC1, belonging to glycoside hydrolase family 13 subfamily 17 (GH13_17) and GH32, respectively. BmSUH has little activity on maltooligosaccharides, whereas other insect GH13_17 α-glucosidases are active on sucrose and maltooligosaccharides. Little is currently known about the structural mechanisms and substrate specificity of GH13_17 enzymes. In this study, we examined the crystal structures of BmSUH without ligands; in complexes with substrates, products, and inhibitors; and complexed with its covalent intermediate at 1.60-1.85 Å resolutions. These structures revealed that the conformations of amino acid residues around subsite -1 are notably different at each step of the hydrolytic reaction. Such changes have not been previously reported among GH13 enzymes, including exo- and endo-acting hydrolases, such as α-glucosidases and α-amylases. Amino acid residues at subsite +1 are not conserved in BmSUH and other GH13_17 α-glucosidases, but subsite -1 residues are absolutely conserved. Substitutions in three subsite +1 residues, Gln¹⁹¹, Tyr²⁵¹, and Glu⁴⁴⁰, decreased sucrose hydrolysis and increased maltase activity of BmSUH, indicating that these residues are key for determining its substrate specificity. These results provide detailed insights into structure-function relationships in GH13 enzymes and into the molecular evolution of insect GH13_17 α-glucosidases.

Starch-binding domains as CBM families-history, occurrence, structure, function and evolution

The term "starch-binding domain" (SBD) has been applied to a domain within an amylolytic enzyme that gave the enzyme the ability to bind onto raw, i.e. thermally untreated, granular starch. An SBD is a special case of a carbohydrate-binding domain, which in general, is a structurally and functionally independent protein module exhibiting no enzymatic activity but possessing potential to target the catalytic domain to the carbohydrate substrate to accommodate it and process it at the active site. As so-called families, SBDs together with other carbohydrate-binding modules (CBMs) have become an integral part of the CAZy database (http://www.cazy.org/). The first two well-described SBDs, i.e. the C-terminal Aspergillus-type and the N-terminal Rhizopus-type have been assigned the families CBM20 and CBM21, respectively. Currently, among the 85 established CBM families in CAZy, fifteen can be considered as families having SBD functional characteristics: CBM20, 21, 25, 26, 34, 41, 45, 48, 53, 58, 68, 69, 74, 82 and 83. All known SBDs, with the exception of the extra long CBM74, were recognized as a module consisting of approximately 100 residues, adopting a β-sandwich fold and possessing at least one carbohydrate-binding site. The present review aims to deliver and describe: (i) the SBD identification in different amylolytic and related enzymes (e.g., CAZy GH families) as well as in other relevant enzymes and proteins (e.g., laforin, the β-subunit of AMPK, and others); (ii) information on the position in the polypeptide chain and the number of SBD copies and their CBM family affiliation (if appropriate); (iii) structure/function studies of SBDs with a special focus on solved tertiary structures, in particular, as complexes with α-glucan ligands; and (iv) the evolutionary relationships of SBDs in a tree common to all SBD CBM families (except for the extra long CBM74). Finally, some special cases and novel potential SBDs are also introduced.

Structural basis of glycogen metabolism in bacteria

The evolution of metabolic pathways is a major force behind natural selection. In the spotlight of such process lies the structural evolution of the enzymatic machinery responsible for the central energy metabolism. Specifically, glycogen metabolism has emerged to allow organisms to save available environmental surplus of carbon and energy, using dedicated glucose polymers as a storage compartment that can be mobilized at future demand. The origins of such adaptive advantage rely on the acquisition of an enzymatic system for the biosynthesis and degradation of glycogen, along with mechanisms to balance the assembly and disassembly rate of this polysaccharide, in order to store and recover glucose according to cell energy needs. The first step in the classical bacterial glycogen biosynthetic pathway is carried out by the adenosine 5'-diphosphate (ADP)-glucose pyrophosphorylase. This allosteric enzyme synthesizes ADP-glucose and acts as a point of regulation. The second step is carried out by the glycogen synthase, an enzyme that generates linear α-(1→4)-linked glucose chains, whereas the third step catalyzed by the branching enzyme produces α-(1→6)-linked glucan branches in the polymer. Two enzymes facilitate glycogen degradation: glycogen phosphorylase, which functions as an α-(1→4)-depolymerizing enzyme, and the debranching enzyme that catalyzes the removal of α-(1→6)-linked ramifications. In this work, we rationalize the structural basis of glycogen metabolism in bacteria to the light of the current knowledge. We describe and discuss the remarkable progress made in the understanding of the molecular mechanisms of substrate recognition and product release, allosteric regulation and catalysis of all those enzymes.

Functional Roles of Starch Binding Domains and Surface Binding Sites in Enzymes Involved in Starch Biosynthesis

Biosynthesis of starch is catalyzed by a cascade of enzymes. The activity of a large number of these enzymes depends on interaction with polymeric substrates via carbohydrate binding sites, which are situated outside of the catalytic site and its immediate surroundings including the substrate-binding crevice. Such secondary binding sites can belong to distinct starch binding domains (SBDs), classified as carbohydrate binding modules (CBMs), or be surface binding sites (SBSs) exposed on the surface of catalytic domains. Currently in the Carbohydrate-Active enZYmes (CAZy) database SBDs are found in 13 CBM families. Four of these families; CBM20, CBM45, CBM48, and CBM53 are represented in enzymes involved in starch biosynthesis, namely starch synthases, branching enzymes, isoamylases, glucan, water dikinases, and α-glucan phosphatases. A critical role of the SBD in activity has not been demonstrated for any of these enzymes. Among the well-characterized SBDs important for starch biosynthesis are three CBM53s of Arabidopsis thaliana starch synthase III, which have modest affinity. SBSs, which are overall less widespread than SBDs, have been reported in some branching enzymes, isoamylases, synthases, phosphatases, and phosphorylases active in starch biosynthesis. SBSs appear to exert roles similar to CBMs. SBSs, however, have also been shown to modulate specificity for example by discriminating the length of chains transferred by branching enzymes. Notably, the difference in rate of occurrence between SBDs and SBSs may be due to lack of awareness of SBSs. Thus, SBSs as opposed to CBMs are not recognized at the protein sequence level, which hampers their identification. Moreover, only a few SBSs in enzymes involved in starch biosynthesis have been functionally characterized, typically by structure-guided site-directed mutagenesis. The glucan phosphatase Like SEX4 2 from A. thaliana has two SBSs with weak affinity for β-cyclodextrin, amylose and amylopectin, which were indicated by mutational analysis to be more important than the active site for initial substrate recognition. The present review provides an update on occurrence of functional SBDs and SBSs in enzymes involved in starch biosynthesis.

Starch formation inside plastids of higher plants

Starch is a water-insoluble polyglucan synthesized inside the plastid stroma within plant cells, serving a crucial role in the carbon budget of the whole plant by acting as a short-term and long-term store of energy. The highly complex, hierarchical structure of the starch granule arises from the actions of a large suite of enzyme activities, in addition to physicochemical self-assembly mechanisms. This review outlines current knowledge of the starch biosynthetic pathway operating in plant cells in relation to the micro- and macro-structures of the starch granule. We highlight the gaps in our knowledge, in particular, the relationship between enzyme function and operation at the molecular level and the formation of the final, macroscopic architecture of the granule.

Cloning of the full-length isoamylase3 gene from cassava Manihot esculenta Crantz 'KU50' and its heterologous expression in E. coli

Isoamylase (EC.3.2.1.68), an essential enzyme in starch metabolism, catalyses the cleavage of α-1,6 glucosidic linkages of branched α-polyglucans such as beta-limit dextrin and amylopectin, but not pullulan. Three different isoamylase isoforms have been reported in plants and algae. We herein report on the first success in preparation of full-length isoamylase3 gene (MeISA3) of cassava Manihot esculenta Crantz 'KU50' from 5' Rapid Amplification of cDNA Ends (5' RACE). The MeISA3 was cloned to pET21b and expressed in E. coli. The Histrap^TM-purified rMeISA3 appeared as a single band protein with approximate molecular size of 75 kDa on SDS-PAGE and Western blot, while 80 kDa was shown by gel filtration chromatography. This indicated the existence of a monomeric enzyme. Biochemical characterisation of rMeISA3 showed that the enzyme was specific towards beta-limit dextrin, with optimal activity at 37 °C pH 6.0. Activity of rMeISA3 could be significantly promoted by Mg²⁺ and Co²⁺. rMeISA3 debranched glucan chains of amylopectin were confirmed by HPAEC-PAD analysis.

Elucidation of the mechanism of interaction between Klebsiella pneumoniae pullulanase and cyclodextrin

Crystal structures of Klebsiella pneumoniae pullulanase (KPP) in complex with α-cyclodextrin (α-CD), β-cyclodextrin (β-CD) and γ-cyclodextrin (γ-CD) were refined at around 1.98–2.59 Å resolution from data collected at SPring-8. In the structures of the complexes obtained with 1 mM α-CD or γ-CD, one molecule of CD was found at carbohydrate-binding module 41 only (CBM41). In the structures of the complexes obtained with 1 mM β-CD or with 10 mM α-CD or γ-CD, two molecules of CD were found at CBM41 and in the active-site cleft, where the hydrophobic residue of Phe746 occupies the inside cavity of the CD rings. In contrast to α-CD and γ-CD, one β-CD molecule was found at the active site only in the presence of 0.1 mM β-CD. These results were coincident with the solution experiments, which showed that β-CD inhibits this enzyme more than a thousand times more potently than α-CD and γ-CD. The strong inhibition of β-CD is caused by the optimized interaction between β-CD and the side chain of Phe746. The increased K i values of the F746A mutant for β-CD supported the importance of Phe746 in the strong interaction of pullulanase with β-CD.

Evolutionary, structural and expression analysis of core genes involved in starch synthesis

Starch is the main storage carbohydrate in plants and an important natural resource for food, feed and industrial raw materials. However, the details regarding the pathway for starch biosynthesis and the diversity of biosynthetic enzymes involved in this process are poorly understood. This study uses a comprehensive phylogenetic analysis of 74 sequenced plant genomes to revisit the evolutionary history of the genes encoding ADP-glucose pyrophosphorylase (AGPase), starch synthase (SS), starch branching enzyme (SBE) and starch de-branching enzyme (DBE). Additionally, the protein structures and expression patterns of these four core genes in starch biosynthesis were studied to determine their functional differences. The results showed that AGPase, SS, SBE and DBE have undergone complicated evolutionary processes in plants and that gene/genome duplications are responsible for the observed differences in isoform numbers. A structure analysis of these proteins suggested that the deletion/mutation of amino acids in some active sites resulted in not only structural variation but also sub-functionalization or neo-functionalization. Expression profiling indicated that AGPase-, SS-, SBE- and DBE-encoding genes exhibit spatio-temporally divergent expression patterns related to the composition of functional complexes in starch biosynthesis. This study provides a comprehensive atlas of the starch biosynthetic pathway, and these data should support future studies aimed at increasing understanding of starch biosynthesis and the functional evolutionary divergence of AGPase, SS, SBE, and DBE in plants.

Heterologous co-expression in E. coli of isoamylase genes from cassava Manihot esculenta Crantz 'KU50' achieves enzyme-active heteromeric complex formation

Cloning of two isoamylase genes, MeISA1 and MeISA2, from cassava (Manihot esculenta Crantz) tubers, accompanied by their co-expression in E. coli demonstrates a requirement for heteromeric complex formation to achieve debranching activity. Starch debranching enzyme (DBE) or isoamylase (ISA) (EC.3.2.1.68), an important enzyme in starch metabolism, catalyses the hydrolysis of α-1,6 glycosidic linkages of amylopectin. Isoforms of ISAs have been reported in higher plants and algae (Fujita et al. in Planta 208:283-293, 1999; Hussain et al. in Plant Cell 15:133-149, 2003; Ishizaki et al. in Agric Biol Chem 47:771-779, 1983; Mouille et al. in Plant Cell 8:1353-1366, 1996). In the current work, cassava ISA genes were isolated from cDNA generated from total RNA from tubers of Manihot esculanta Crantz cultivar KU50. MeISA1 and MeISA2 were successfully amplified and cloned into a pETDuet1 vector. The putative MeISA1 and MeISA2 proteins comprised 763 and 882 amino acids, with substantial similarity to StISA1 and StISA2 from potato (84.4% and 68.9%, respectively). Recombinant MeISA1 and MeISA2 were co-expressed in Escherichia coli SoluBL21 (DE3). Histrap^TM-Purified rMeISA1 and rMeISA2 showed approximate molecular weights of 87 and 99 kDa, respectively, by SDS-PAGE. Debranching activity was only detectable in the column fractions where both recombinant ISA isoforms were present. The heteromeric DBE from crude extracts of 4-5 h induced cultures analysed by gel filtration chromatography and western blot showed combinations of rMeISA1 and rMeISA2 at ratios of 1:1 to 4:1. Pooled fractions with DBE activity were used for enzyme characterisation, which showed that the enzyme was specific for amylopectin, with optimum activity at 37 °C and pH 7.0. Enzyme activity was enhanced by Co²⁺, Mg²⁺ and Ca²⁺, but was strongly inhibited by Cu²⁺. Debranched amylopectin products showed chain length distributions typical of plant DBE.

Starch Biosynthesis in the Developing Endosperms of Grasses and Cereals

The starch-rich endosperms of the Poaceae, which includes wild grasses and their domesticated descendents the cereals, have provided humankind and their livestock with the bulk of their daily calories since the dawn of civilization up to the present day. There are currently unprecedented pressures on global food supplies, largely resulting from population growth, loss of agricultural land that is linked to increased urbanization, and climate change. Since cereal yields essentially underpin world food and feed supply, it is critical that we understand the biological factors contributing to crop yields. In particular, it is important to understand the biochemical pathway that is involved in starch biosynthesis, since this pathway is the major yield determinant in the seeds of six out of the top seven crops grown worldwide. This review outlines the critical stages of growth and development of the endosperm tissue in the Poaceae, including discussion of carbon provision to the growing sink tissue. The main body of the review presents a current view of our understanding of storage starch biosynthesis, which occurs inside the amyloplasts of developing endosperms.

Functional and structural characterization of plastidic starch phosphorylase during barley endosperm development

The production of starch is essential for human nutrition and represents a major metabolic flux in the biosphere. The biosynthesis of starch in storage organs like barley endosperm operates via two main pathways using different substrates: starch synthases use ADP-glucose to produce amylose and amylopectin, the two major components of starch, whereas starch phosphorylase (Pho1) uses glucose-1-phosphate (G1P), a precursor for ADP-glucose production, to produce α-1,4 glucans. The significance of the Pho1 pathway in starch biosynthesis has remained unclear. To elucidate the importance of barley Pho1 (HvPho1) for starch biosynthesis in barley endosperm, we analyzed HvPho1 protein production and enzyme activity levels throughout barley endosperm development and characterized structure-function relationships of HvPho1. The molecular mechanisms underlying the initiation of starch granule biosynthesis, that is, the enzymes and substrates involved in the initial transition from simple sugars to polysaccharides, remain unclear. We found that HvPho1 is present as an active protein at the onset of barley endosperm development. Notably, purified recombinant protein can catalyze the de novo production of α-1,4-glucans using HvPho1 from G1P as the sole substrate. The structural properties of HvPho1 provide insights into the low affinity of HvPho1 for large polysaccharides like starch or amylopectin. Our results suggest that HvPho1 may play a role during the initiation of starch biosynthesis in barley.

Arabidopsis thaliana FAR-RED ELONGATED HYPOCOTYLS3 (FHY3) and FAR-RED-IMPAIRED RESPONSE1 (FAR1) modulate starch synthesis in response to light and sugar

In living organisms, daily light/dark cycles profoundly affect cellular processes. In plants, optimal growth and development, and adaptation to daily light-dark cycles, require starch synthesis and turnover. However, the underlying molecular mechanisms coordinating daily starch metabolism remain poorly understood. To explore the roles of Arabidopsis thaliana light signal transduction proteins FAR-RED ELONGATED HYPOCOTYLS3 (FHY3) and FAR-RED-IMPAIRED RESPONSE1 (FAR1) in starch metabolism, the contents of starch and water-soluble polysaccharides, and the structure of starch granules were investigated in fhy3, far1 and fhy3 far1 mutant plants. Disruption of FHY3 or FAR1 reduced starch accumulation and altered starch granule structure in the fhy3-4, far1-2, and fhy3-4 far1-2 mutant plants. Furthermore, molecular and genetic evidence revealed that the gene encoding the starch-debranching enzyme ISOAMYLASE2 (ISA2) is a direct target of FHY3 and FAR1, and functions in light-induced starch synthesis. Our data establish the first molecular link between light signal transduction and starch synthesis, suggesting that the light-signaling proteins FHY3 and FAR1 influence starch synthesis and starch granule formation through transcriptional activation of ISA2.

Bound Substrate in the Structure of Cyanobacterial Branching Enzyme Supports a New Mechanistic Model

Branching enzyme (BE) catalyzes the formation of α-1,6-glucosidic linkages in amylopectin and glycogen. The reaction products are variable, depending on the organism sources, and the mechanistic basis for these different outcomes is unclear. Although most cyanobacteria have only one BE isoform belonging to glycoside hydrolase family 13, Cyanothece sp. ATCC 51142 has three isoforms (BE1, BE2, and BE3) with distinct enzymatic properties, suggesting that investigations of these enzymes might provide unique insights into this system. Here, we report the crystal structure of ligand-free wild-type BE1 (residues 5-759 of 1-773) at 1.85 Å resolution. The enzyme consists of four domains, including domain N, carbohydrate-binding module family 48 (CBM48), domain A containing the catalytic site, and domain C. The central domain A displays a (β/α)₈-barrel fold, whereas the other domains adopt β-sandwich folds. Domain N was found in a new location at the back of the protein, forming hydrogen bonds and hydrophobic interactions with CBM48 and domain A. Site-directed mutational analysis identified a mutant (W610N) that bound maltoheptaose with sufficient affinity to enable structure determination at 2.30 Å resolution. In this structure, maltoheptaose was bound in the active site cleft, allowing us to assign subsites -7 to -1. Moreover, seven oligosaccharide-binding sites were identified on the protein surface, and we postulated that two of these in domain A served as the entrance and exit of the donor/acceptor glucan chains, respectively. Based on these structures, we propose a substrate binding model explaining the mechanism of glycosylation/deglycosylation reactions catalyzed by BE.

Structure and function of α-glucan debranching enzymes

α-Glucan debranching enzymes hydrolyse α-1,6-linkages in starch/glycogen, thereby, playing a central role in energy metabolism in all living organisms. They belong to glycoside hydrolase families GH13 and GH57 and several of these enzymes are industrially important. Nine GH13 subfamilies include α-glucan debranching enzymes; isoamylase and glycogen debranching enzymes (GH13_11); pullulanase type I/limit dextrinase (GH13_12-14); pullulan hydrolase (GH13_20); bifunctional glycogen debranching enzyme (GH13_25); oligo-1 and glucan-1,6-α-glucosidases (GH13_31); pullulanase type II (GH13_39); and α-amylase domains (GH13_41) in two-domain amylase-pullulanases. GH57 harbours type II pullulanases. Specificity differences, domain organisation, carbohydrate binding modules, sequence motifs, three-dimensional structures and specificity determinants are discussed. The phylogenetic analysis indicated that GH13_39 enzymes could represent a "missing link" between the strictly α-1,6-specific debranching enzymes and the enzymes with dual specificity and α-1,4-linkage preference.

Formation of starch in plant cells

Starch-rich crops form the basis of our nutrition, but plants have still to yield all their secrets as to how they make this vital substance. Great progress has been made by studying both crop and model systems, and we approach the point of knowing the enzymatic machinery responsible for creating the massive, insoluble starch granules found in plant tissues. Here, we summarize our current understanding of these biosynthetic enzymes, highlighting recent progress in elucidating their specific functions. Yet, in many ways we have only scratched the surface: much uncertainty remains about how these components function together and are controlled. We flag-up recent observations suggesting a significant degree of flexibility during the synthesis of starch and that previously unsuspected non-enzymatic proteins may have a role. We conclude that starch research is not yet a mature subject and that novel experimental and theoretical approaches will be important to advance the field.

Comparison of Chain-Length Preferences and Glucan Specificities of Isoamylase-Type α-Glucan Debranching Enzymes from Rice, Cyanobacteria, and Bacteria

It has been believed that isoamylase (ISA)-type α-glucan debranching enzymes (DBEs) play crucial roles not only in α-glucan degradation but also in the biosynthesis by affecting the structure of glucans, although molecular basis on distinct roles of the individual DBEs has not fully understood. In an attempt to relate the roles of DBEs to their chain-length specificities, we analyzed the chain-length distribution of DBE enzymatic reaction products by using purified DBEs from various sources including rice, cyanobacteria, and bacteria. When DBEs were incubated with phytoglycogen, their chain-length specificities were divided into three groups. First, rice endosperm ISA3 (OsISA3) and Eschericia coli GlgX (EcoGlgX) almost exclusively debranched chains having degree of polymerization (DP) of 3 and 4. Second, OsISA1, Pseudomonas amyloderamosa ISA (PsaISA), and rice pullulanase (OsPUL) could debranch a wide range of chains of DP≧3. Third, both cyanobacteria ISAs, Cyanothece ATCC 51142 ISA (CytISA) and Synechococcus elongatus PCC7942 ISA (ScoISA), showed the intermediate chain-length preference, because they removed chains of mainly DP3-4 and DP3-6, respectively, while they could also react to chains of DP5-10 and 7–13 to some extent, respectively. In contrast, all these ISAs were reactive to various chains when incubated with amylopectin. In addition to a great variation in chain-length preferences among various ISAs, their activities greatly differed depending on a variety of glucans. Most strikingly, cyannobacteria ISAs could attack branch points of pullulan to a lesser extent although no such activity was found in OsISA1, OsISA3, EcoGlgX, and PsaISA. Thus, the present study shows the high possibility that varied chain-length specificities of ISA-type DBEs among sources and isozymes are responsible for their distinct functions in glucan metabolism.

Phylogeny-structured carbohydrate metabolism across microbiomes collected from different units in wastewater treatment process

BACKGROUND

With respect to global priority for bioenergy production from plant biomass, understanding the fundamental genetic associations underlying carbohydrate metabolisms is crucial for the development of effective biorefinery process. Compared with gut microbiome of ruminal animals and wood-feed insects, knowledge on carbohydrate metabolisms of engineered biosystems is limited.

RESULTS

In this study, comparative metagenomics coupled with metabolic network analysis was carried out to study the inter-species cooperation and competition among carbohydrate-active microbes in typical units of wastewater treatment process including activated sludge and anaerobic digestion. For the first time, sludge metagenomes demonstrated rather diverse pool of carbohydrate-active genes (CAGs) comparable to that of rumen microbiota. Overall, the CAG composition correlated strongly with the microbial phylogenetic structure across sludge types. Gene-centric clustering analysis showed the carbohydrate pathways of sludge systems were shaped by different environmental factors, including dissolved oxygen and salinity, and the latter showed more determinative influence of phylogenetic composition. Eventually, the highly clustered co-occurrence network of CAGs and saccharolytic phenotypes, revealed three metabolic modules in which the prevalent populations of Actinomycetales, Clostridiales and Thermotogales, respectively, play significant roles as interaction hubs, while broad negative co-exclusion correlations observed between anaerobic and aerobic microbes, probably implicated roles of niche separation by dissolved oxygen in determining the microbial assembly.

CONCLUSIONS

Sludge microbiomes encoding diverse pool of CAGs was another potential source for effective lignocellulosic biomass breakdown. But unlike gut microbiomes in which Clostridiales, Lactobacillales and Bacteroidales play a vital role, the carbohydrate metabolism of sludge systems is built on the inter-species cooperation and competition among Actinomycetales, Clostridiales and Thermotogales.

Amylopectin biosynthetic enzymes from developing rice seed form enzymatically active protein complexes

Amylopectin is a highly branched, organized cluster of glucose polymers, and the major component of rice starch. Synthesis of amylopectin requires fine co-ordination between elongation of glucose polymers by soluble starch synthases (SSs), generation of branches by branching enzymes (BEs), and removal of misplaced branches by debranching enzymes (DBEs). Among the various isozymes having a role in amylopectin biosynthesis, limited numbers of SS and BE isozymes have been demonstrated to interact via protein-protein interactions in maize and wheat amyloplasts. This study investigated whether protein-protein interactions are also found in rice endosperm, as well as exploring differences between species. Gel permeation chromatography of developing rice endosperm extracts revealed that all 10 starch biosynthetic enzymes analysed were present at larger molecular weights than their respective monomeric sizes. SSIIa, SSIIIa, SSIVb, BEI, BEIIb, and PUL co-eluted at mass sizes >700kDa, and SSI, SSIIa, BEIIb, ISA1, PUL, and Pho1 co-eluted at 200-400kDa. Zymogram analyses showed that SSI, SSIIIa, BEI, BEIIa, BEIIb, ISA1, PUL, and Pho1 eluted in high molecular weight fractions were active. Comprehensive co-immunoprecipitation analyses revealed associations of SSs-BEs, and, among BE isozymes, BEIIa-Pho1, and pullulanase-type DBE-BEI interactions. Blue-native-PAGE zymogram analyses confirmed the glucan-synthesizing activity of protein complexes. These results suggest that some rice starch biosynthetic isozymes are physically associated with each other and form active protein complexes. Detailed analyses of these complexes will shed light on the mechanisms controlling the unique branch and cluster structure of amylopectin, and the physicochemical properties of starch.

Sequence variation, differential expression, and divergent evolution in starch-related genes among accessions of Arabidopsis thaliana

Transitory starch metabolism is a nonlinear and highly regulated process. It originated very early in the evolution of chloroplast-containing cells and is largely based on a mosaic of genes derived from either the eukaryotic host cell or the prokaryotic endosymbiont. Initially located in the cytoplasm, starch metabolism was rewired into plastids in Chloroplastida. Relocation was accompanied by gene duplications that occurred in most starch-related gene families and resulted in subfunctionalization of the respective gene products. Starch-related isozymes were then evolutionary conserved by constraints such as internal starch structure, posttranslational protein import into plastids and interactions with other starch-related proteins. 25 starch-related genes in 26 accessions of Arabidopsis thaliana were sequenced to assess intraspecific diversity, phylogenetic relationships, and modes of selection. Furthermore, sequences derived from additional 80 accessions that are publicly available were analyzed. Diversity varies significantly among the starch-related genes. Starch synthases and phosphorylases exhibit highest nucleotide diversities, while pyrophosphatases and debranching enzymes are most conserved. The gene trees are most compatible with a scenario of extensive recombination, perhaps in a Pleistocene refugium. Most genes are under purifying selection, but disruptive selection was inferred for a few genes/substitutiones. To study transcript levels, leaves were harvested throughout the light period. By quantifying the transcript levels and by analyzing the sequence of the respective accessions, we were able to estimate whether transcript levels are mainly determined by genetic (i.e., accession dependent) or physiological (i.e., time dependent) parameters. We also identified polymorphic sites that putatively affect pattern or the level of transcripts.