Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut (Arachis hypogaea L.)

The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils. Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect-increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes. With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanut, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold (α = 1%). We detected 7,682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (p-value < 1.0E-16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Discovery and Engineering of Pathways for Production of α-Branched Organic Acids

Cell-based synthesis offers many opportunities for preparing small molecules from simple renewable carbon sources by telescoping multiple reactions into a single fermentation step. One challenge in this area is the development of enzymatic carbon-carbon bond forming cycles that enable a modular disconnection of a target structure into cellular building blocks. In this regard, synthetic pathways based on thiolase enzymes to catalyze the initial carbon-carbon bond forming step between acyl coenzyme A (CoA) substrates offer a versatile route for biological synthesis, but the substrate diversity of such pathways is currently limited. In this report, we describe the identification and biochemical characterization of a thiolase-ketoreductase pair involved in production of branched acids in the roundworm, Ascaris suum, that demonstrates selectivity for forming products with an α-methyl branch using a propionyl-CoA extender unit. Engineering synthetic pathways for production of α-methyl acids in Escherichia coli using these enzymes allows the construction of microbial strains that produce either chiral 2-methyl-3-hydroxy acids (1.1 ± 0.2 g L^-1) or branched enoic acids (1.12 ± 0.06 g L^-1) in the presence of a dehydratase at 44% and 87% yield of fed propionate, respectively. In vitro characterization along with in vivo analysis indicates that the ketoreductase is the key driver for selectivity, forming predominantly α-branched products even when paired with a thiolase that highly prefers unbranched linear products. Our results expand the utility of thiolase-based pathways and provide biosynthetic access to α-branched compounds as precursors for polymers and other chemicals.

Fatty acid synthesis enzyme clans

Ketoacyl reductases (KRs), hydroxyacyl dehydratases (HDs), and enoyl reductases (ERs) are part of the fatty acid/polyketide synthesis cycle. They are known as acyl dehydrogenases, enoyl hydratases, and hydroxyacyl dehydrogenases, respectively, when catalyzing their reverse reactions. Earlier, we classified these enzymes into four KR, eight HD, and five ER families by statistical criteria. Members of all four KR families and three ER families have Rossmann folds, while five HD family members have HotDog folds. This suggests that those proteins with the same folds in different families may be distantly related, and therefore in clans, even though their amino acid sequences may not be homologous. We have now defined two clans containing three of the four KR families and two of the eight HD families, using manual and statistical tests. One of the ER families is related to the KR clan.

Evolution of acyl-ACP-thioesterases and β-ketoacyl-ACP-synthases revealed by protein-protein interactions

The fatty acid synthase (FAS) is a conserved primary metabolic enzyme complex capable of tolerating cross-species engineering of domains for the development of modified and overproduced fatty acids. In eukaryotes, acyl-acyl carrier protein thioesterases (TEs) off-load mature cargo from the acyl carrier protein (ACP), and plants have developed TEs for short/medium-chain fatty acids. We showed that engineering plant TEs into the green microalga Chlamydomonas reinhardtii does not result in the predicted shift in fatty acid profile. Since fatty acid biosynthesis relies on substrate recognition and protein-protein interactions between the ACP and its partner enzymes, we hypothesized that plant TEs and algal ACP do not functionally interact. Phylogenetic analysis revealed major evolutionary differences between FAS enzymes, including TEs and ketoacyl synthases (KSs), in which the former is present only in some species, whereas the latter is present in all, and has a common ancestor. In line with these results, TEs appeared to be selective towards their ACP partners whereas KSs showed promiscuous behavior across bacterial, plant and algal species. Based on phylogenetic analyses, in silico docking, in vitro mechanistic crosslinking and in vivo algal engineering, we propose that phylogeny can predict effective interactions between ACPs and partner enzymes.

Sequence analysis and structure prediction of enoyl-CoA hydratase from Avicennia marina: implication of various amino acid residues on substrate-enzyme interactions

Enoyl-CoA hydratase catalyzes the hydration of 2-trans-enoyl-CoA into 3-hydroxyacyl-CoA. The present study focuses on the correlation between the functional and structural aspects of enoyl-CoA hydratase from Avicennia marina. We have used bioinformatics tools to construct and analyze 3D homology models of A. marina enoyl-CoA hydratase (AMECH) bound to different substrates and inhibitors and studied the residues involved in the ligand-enzyme interaction. Structural information obtained from the models was compared with those of the reported crystal structures. We observed that the overall folds were similar; however, AMECH showed few distinct structural changes which include structural variation in the mobile loop, formation and loss of certain interactions between the active site residues and substrates. Some changes were also observed within specific regions of the enzyme. Glu106 is almost completely conserved in sequences of the isomerases/hydratases including AMECH while Glu86 which is the other catalytic residue in most of the isomerases/hydratases is replaced by Gly and shows no interaction with the substrate. Asp114 is located within 4Å distance of the catalytic water which makes it a probable candidate for the second catalytic residue in AMECH. Another prominent feature of AMECH is the presence of structurally distinct mobile loop having a completely different coordination with the hydrophobic binding pocket of acyl portion of the substrate.