Acetyl-CoA carboxylase in higher plants: most plants other than gramineae have both the prokaryotic and the eukaryotic forms of this enzyme

Phylogenetic analysis and structural studies of heteromeric acetyl-CoA carboxylase from the oleaginous Amazonian microalgae Ankistrodesmus sp.: Insights into the BC and BCCP subunits

Acetyl-CoA carboxylase (ACC) is an essential enzyme in fatty acid biosynthesis that catalyzes the formation of malonyl-CoA from acetyl-CoA. While structural studies on ACC components have largely focused on prokaryotes and higher plants, the assembly and molecular adaptations of ACC in microalgae remain underexplored. This study aimed to fill this gap by providing the first structural and evolutionary characterization of both biotin carboxylase (BC) and biotin carboxyl carrier protein (BCCP) from a microalga (Ankistrodesmus sp.). Phylogenetic analysis revealed distinct evolutionary trajectories for BC and BCCP, with BC forming a chlorophyte-specific clade closely related to other oleaginous species, while BCCP displayed two distinct isoforms within green algae, resulting from gene duplication. The crystallographic structure of BC was solved in its apo (1.75 Å) and ADP-Mg2+-bound (1.90 Å) states, revealing conserved conformational changes associated with cofactor binding. BCCP from Ankistrodesmus sp. displayed a unique QLGTF/H motif instead of the canonical AMKXM biotinylation motif, suggesting loss of biotinylation capacity. However, the presence of three additional lysines in the protruding thumb loop, with Lys95 as a candidate for biotin attachment, indicates potential compensatory adaptations. SEC-MALS and pull-down assays confirmed the formation of a stable 1:1 BC-BCCP complex, and circular dichroism showed increased thermal stability of the complex, supporting its structural stability. This study highlights unique structural adaptations in Ankistrodesmus sp. ACC, emphasizing the evolutionary plasticity of BC and BCCP. These insights provide a foundation for future investigations into ACC regulation in photosynthetic organisms and offer potential biotechnological applications for optimizing lipid production in microalgae.

Solution Structure and NMR Chemical Shift Perturbations of the Arabidopsis BCCP1 Identify Intersubunit Interactions Potentially Involved in the Assembly of the Heteromeric Acetyl-CoA Carboxylase

Biotin carboxyl carrier protein (BCCP) is a subunit of the heteromeric acetyl-CoA carboxylase (htACCase), and it chemically links the two half-reactions that constitute the formation of malonyl-CoA from acetyl-CoA, a critical reaction in fatty acid biosynthesis. Because plants are a major source of edible fats and oils, it is important to understand the structural organization of the plant htACCase, relative to its potential to regulate fatty acid biosynthesis in plant plastids. Moreover, unique to the plant htACCase, noncatalytic subunits called biotin attachment domain-containing (BADC) proteins are important in the assembly of the holoenzyme, and they specifically interact with the bcCP and the biotin carboxylase (BC) subunits. We report herein NMR structural studies of the Arabidopsis BCCP isozymes (bcCP1 and BCCP2). We calculated the structure of C-terminal domain of BCCP1 (K200-P280) and explored structural changes in the BCCP1 protein upon its interactions with bc and BADC. The chemical shift perturbation experiments identified potential surface residues on the BCCP1 protein that may facilitate physical interactions between BC and BADC proteins. These studies indicate that the BADC protein interacts with a "thumb"-like protrusion, which is a common structural feature of the bacterial and plant bcCPs, and thereby acts as a potential "cap" to facilitate the assembly of a BC-BCCP-BADC complex.

Validation of assay for measuring acetyl-coenzyme a carboxylase activity in grasses using malachite green

Acetyl-CoA carboxylase (ACCase) is one of the most important enzymes as a herbicide target in gramineous plant species, however, assay methods for the enzyme are primarily limited to those using radioisotopes (RI). Typically, the measurement method that uses RI necessitates specialized facilities and equipment, and involves complex procedures throughout the experiment. As another method for detecting ACCase activity, the colorimetric method using malachite green (MG) is also known. However, reports on this method are limited, and information regarding the simplicity of the procedure and the scope of its application remains unclear. To better understand the method using MG and to develop a simpler assay method, crude enzymes extracted from various target-site resistant (TSR) biotypes of blackgrass (Alopecurus myosuroides) were examined in enzyme inhibition tests. As a result, this method was able to accurately detect the relationship between the chemical classes of ACC inhibitors and cross-resistance to specific TSRs. Moreover, the ACCase activity of other grass species was also examined using this method. By using crude enzymes from various species and a commercially available phosphatase kit containing MG, ACCase activity was detectable with good accuracy. In addition, enzyme inhibition studies using ACCase inhibiting herbicides revealed that this method reproduced results similar to those obtained with the RI method. The Z'-factor, an index of high-throughput screening (HTS), was around 0.7, indicating that it is an excellent screening system. These results suggest that the assay method using MG is very simple, labor-saving, and accurate with a throughput much higher than that of the existing RI method. Therefore, it is strongly suggested that the method could replace the RI method in most cases. These results indicate that it is applicable to HTS for ACCase inhibitors.

Metproxybicyclone, a Novel Carbocyclic Aryl-dione Acetyl-CoA Carboxylase-Inhibiting Herbicide for the Management of Sensitive and Resistant Grass Weeds

Postemergence control of grass weeds has become problematic due to the evolution of resistance to 5-enolpyruvylshikimate-3-phosphate synthase, acetyl-CoA carboxylase (ACCase), and acetolactate synthase-inhibiting herbicides. Herein we describe the invention and synthesis journey toward metproxybicyclone, the first commercial carbocyclic aryl-dione ACCase-inhibiting herbicide for the cost-effective management of grass weeds in dicotyledonous crops and in preplant burndown applications. Glasshouse and field experiments have shown that metproxybicyclone is safe for use on soybean, cotton, and sugar beet, among other crops. It is effective on a variety of key grass weeds including Eleusine indica, Digitaria insularis, Sorghum halepense, and Echinochloa crus-galli. Importantly, metproxybicyclone was more efficacious at killing resistant grass weed populations than current ACCase herbicides. Metproxybicyclone controlled the main ACCase target-site and nontarget site resistant mechanisms in characterized Lolium multiflorum and E. indica populations under glasshouse conditions. Excellent control of a broad resistance-causing D2078G target-site mutant E. indica population was also observed under field conditions.

Bioinformatics Identification and Expression Analysis of Acetyl-CoA Carboxylase Reveal Its Role in Isoflavone Accumulation during Soybean Seed Development

Isoflavones belong to the class of flavonoid compounds, which are important secondary metabolites that play a crucial role in plant development and defense. Acetyl-CoA carboxylase (ACCase) is a biotin-dependent enzyme that catalyzes the conversion of Acetyl-CoA into Malonyl-CoA in plants. It is a key enzyme in fatty acid synthesis and also catalyzes the production of various secondary metabolites. However, information on the ACC gene family in the soybean (Glycine max L. Merr.) genome and the specific members involved in isoflavone biosynthesis is still lacking. In this study, we identified 20 ACC family genes (GmACCs) from the soybean genome and further characterized their evolutionary relationships and expression patterns. Phylogenetic analysis showed that the GmACCs could be divided into five groups, and the gene structures within the same groups were highly conserved, indicating that they had similar functions. The GmACCs were randomly distributed across 12 chromosomes, and collinearity analysis suggested that many GmACCs originated from tandem and segmental duplications, with these genes being under purifying selection. In addition, gene expression pattern analysis indicated that there was functional divergence among GmACCs in different tissues. The GmACCs reached their peak expression levels during the early or middle stages of seed development. Based on the transcriptome and isoflavone content data, a weighted gene co-expression network was constructed, and three candidate genes (Glyma.06G105900, Glyma.13G363500, and Glyma.13G057400) that may positively regulate isoflavone content were identified. These results provide valuable information for the further functional characterization and application of GmACCs in isoflavone biosynthesis in soybean.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

The Cys-2088-Arg mutation in the ACCase gene and enhanced metabolism confer cyhalofop-butyl resistance in Chinese sprangletop (Leptochloa chinensis)

Acetyl-CoA carboxylase (ACCase)-inhibiting herbicides are among the most commonly used herbicides to control grassy weeds, especially Leptochloa chinensis, in rice fields across China. Herein, we collected a suspected resistant (R) population of L. chinensis (HFLJ16) from Lujiang county in Anhui Province. Whole plant dose response tests showed that, compared with the susceptible (S) population, the R population showed high resistance to cyhalofop-butyl (22-fold) and displayed cross-resistance to metamifop (9.7-fold), fenoxaprop-P-ethyl (18.7-fold), quizalofop-P-ethyl (7.6-fold), clodinafop-propargyl (12-fold) and clethodim (8.4-fold). We detected an amino acid substitution (Cys-2088-Arg) in the ACCase of resistant L. chinensis. However, ACCase gene expression levels were not significantly different (P > 0.05) between R plants and S plants, without or with cyhalofop-butyl treatment. Furthermore, pretreatment with piperonyl butoxide (PBO, a cytochrome P450 monooxygenase (CYP450) inhibitor) or 4-chloro-7-nitrobenzoxadiazole (NBD-Cl, a glutathione-S-transferase (GST) inhibitor), inhibited the resistance of the R population to cyhalofop-butyl significantly (by approximately 60% and 26%, respectively). Liquid chromatography tandem mass spectrometry analysis showed that R plants metabolized cyhalofop-butyl and cyhalofop acid (its metabolite) significantly faster than S plants. Three CYP450 genes, one GST gene, and two ABC transporter genes were induced by cyhalofop-butyl and were overexpressed in the R population. Overall, GST-associated detoxification, CYP450 enhancement, and target-site gene mutation are responsible for the resistance of L. chinensis to cyhalofop-butyl.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Overexpression of genes involved in fatty acid biosynthesis increases lipid content in the NaHCO3-tolerant Chlorella sp. JB6

IMPORTANCE

Fatty acid (FA) contents can be altered in Chlorella JB6 in the presence of sodium bicarbonate (NaHCO3). Overexpression of the FA de novo synthesis genes inhibited the growth of JB6 cells and decreased their resistance to NaHCO3, but these transgenic JB6 strains could grow in a medium containing as high as 300 mM NaHCO3. In JB6, ectopic expression of the FA de novo synthesis genes increased the synthesis of very long-chain saturated FA (> 20C).

The RNA landscape of Dunaliella salina in response to short-term salt stress

Using the halotolerant green microalgae Dunaliella salina as a model organism has special merits, such as a wide range of salt tolerance, unicellular organism, and simple life cycle and growth conditions. These unique characteristics make it suitable for salt stress study. In order to provide an overview of the response of Dunaliella salina to salt stress and hopefully to reveal evolutionarily conserved mechanisms of photosynthetic organisms in response to salt stress, the transcriptomes and the genome of the algae were sequenced by the second and the third-generation sequencing technologies, then the transcriptomes under salt stress were compared to the transcriptomes under non-salt stress with the newly sequenced genome as the reference genome. The major cellular biological processes that being regulated in response to salt stress, include transcription, protein synthesis, protein degradation, protein folding, protein modification, protein transport, cellular component organization, cell redox homeostasis, DNA repair, glycerol synthesis, energy metabolism, lipid metabolism, and ion homeostasis. This study gives a comprehensive overview of how Dunaliella salina responses to salt stress at transcriptomic level, especially characterized by the nearly ubiquitous up-regulation of the genes involving in protein folding, DNA repair, and cell redox homeostasis, which may confer the algae important mechanisms to survive under salt stress. The three fundamental biological processes, which face huge challenges under salt stress, are ignored by most scientists and are worth further deep study to provide useful information for breeding economic important plants competent in tolerating salt stress, other than only depending on the commonly acknowledged osmotic balance and ion homeostasis.

The Sapria himalayana genome provides new insights into the lifestyle of endoparasitic plants

BACKGROUND

Sapria himalayana (Rafflesiaceae) is an endoparasitic plant characterized by a greatly reduced vegetative body and giant flowers; however, the mechanisms underlying its special lifestyle and greatly altered plant form remain unknown. To illustrate the evolution and adaptation of S. himalayasna, we report its de novo assembled genome and key insights into the molecular basis of its floral development, flowering time, fatty acid biosynthesis, and defense responses.

RESULTS

The genome of S. himalayana is ~ 1.92 Gb with 13,670 protein-coding genes, indicating remarkable gene loss (~ 54%), especially genes involved in photosynthesis, plant body, nutrients, and defense response. Genes specifying floral organ identity and controlling organ size were identified in S. himalayana and Rafflesia cantleyi, and showed analogous spatiotemporal expression patterns in both plant species. Although the plastid genome had been lost, plastids likely biosynthesize essential fatty acids and amino acids (aromatic amino acids and lysine). A set of credible and functional horizontal gene transfer (HGT) events (involving genes and mRNAs) were identified in the nuclear and mitochondrial genomes of S. himalayana, most of which were under purifying selection. Convergent HGTs in Cuscuta, Orobanchaceae, and S. himalayana were mainly expressed at the parasite-host interface. Together, these results suggest that HGTs act as a bridge between the parasite and host, assisting the parasite in acquiring nutrients from the host.

CONCLUSIONS

Our results provide new insights into the flower development process and endoparasitic lifestyle of Rafflesiaceae plants. The amount of gene loss in S. himalayana is consistent with the degree of reduction in its body plan. HGT events are common among endoparasites and play an important role in their lifestyle adaptation.

Study on molecular mechanism of volatiles variation during Bupleurum scorzonerifolium root development based on metabolome and transcriptome analysis

Bupleurum scorzonerifolium Willd. is a medicinal herb. Its root has a high content of volatile oil (BSVO), which shows a variety of biological activities. Currently, BSVO in the injectable form is used for treating fever in humans and livestock. The yield and quality of volatile oils depends on the developmental stages of plants. However, the changes in BSVO yield and quality during root development in Bupleurum scorzonerifolium and the underlying molecular regulatory mechanisms remain unclear. This knowledge gap is limiting the improvement in the quality of BSVO. In the present study, B. scorzonerifolium root was collected at germinative, vegetative, florescence, fruiting and defoliating stages. The yield of BSVO, metabolic profile of volatile components and transcriptome of root samples at various developmental stages were comprehensively determined and compared. BSVO continuously accumulated from the germinative to fruiting stages, and its level slightly decreased from the fruiting to defoliating stages. A total of 82 volatile components were detected from B. scorzonerifolium root, of which 22 volatiles were identified as differentially accumulated metabolites (DAMs) during the root development. Of these volatiles, fatty acids and their derivatives accounted for the largest proportion. The contents of most major volatiles were highest at the fruiting stage. A large number of differentially expressed genes (DEGs) were detected during B. scorzonerifolium root development, of which 65 DEGs encoded various enzymes and transcription factors regulating the biosynthesis of fatty acids and their derivatives. In further analysis, 42 DEGs were identified to be significantly correlated with DAMs, and these DEGs may be the key genes for the biosynthesis of volatiles. To the best of our knowledge, this is the first study to comprehensively report the changes in the composition and content of volatiles and underlying mechanism during B. scorzonerifolium root development. This study provided important reference for future studies to determine the harvest time of B. scorzonerifolium roots and improve the quality of BSVO.

The PRK/Rubisco shunt strongly influences Arabidopsis seed metabolism and oil accumulation, affecting more than carbon recycling

The carbon efficiency of storage lipid biosynthesis from imported sucrose in green Brassicaceae seeds is proposed to be enhanced by the PRK/Rubisco shunt, in which ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) acts outside the context of the Calvin-Benson-Bassham cycle to recycle CO2 molecules released during fatty acid synthesis. This pathway utilizes metabolites generated by the nonoxidative steps of the pentose phosphate pathway. Photosynthesis provides energy for reactions such as the phosphorylation of ribulose 5-phosphate by phosphoribulokinase (PRK). Here, we show that loss of PRK in Arabidopsis thaliana (Arabidopsis) blocks photoautotrophic growth and is seedling-lethal. However, seeds containing prk embryos develop normally, allowing us to use genetics to assess the importance of the PRK/Rubisco shunt. Compared with nonmutant siblings, prk embryos produce one-third less lipids-a greater reduction than expected from simply blocking the proposed PRK/Rubisco shunt. However, developing prk seeds are also chlorotic and have elevated starch contents compared with their siblings, indicative of secondary effects. Overexpressing PRK did not increase embryo lipid content, but metabolite profiling suggested that Rubisco activity becomes limiting. Overall, our findings show that the PRK/Rubisco shunt is tightly integrated into the carbon metabolism of green Arabidopsis seeds, and that its manipulation affects seed glycolysis, starch metabolism, and photosynthesis.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Genome-wide identification of biotin carboxyl carrier subunits of acetyl-CoA carboxylase in Brassica and their role in stress tolerance in oilseed Brassica napus

Background

Biotin carboxyl carrier protein (BCCP) is a subunit of Acetyl CoA-carboxylase (ACCase) which catalyzes the conversion of acetyl-CoA to malonyl-CoA in a committed step during the de novo biosynthesis of fatty acids. Lipids, lipid metabolites, lipid-metabolizing and -modifying enzymes are known to play a role in biotic and abiotic stress tolerance in plants. In this regard, an understanding of the Brassica napus BCCP genes will aid in the improvement of biotic and abiotic stress tolerance in canola.

Results

In this study, we identified 43 BCCP genes in five Brassica species based on published genome data. Among them, Brassica rapa, Brassica oleracea, Brassica nigra, Brassica napus and Brassica juncea had six, seven, seven, 10 and 13 BCCP homologs, respectively. Phylogenetic analysis categorized them into five classes, each with unique conserved domains. The promoter regions of all BCCP genes contained stress-related cis-acting elements as determined by cis-element analysis. We identified four and three duplicated gene pairs (segmental) in B. napus and B. juncea respectively, indicating the role of segmental duplication in the expansion of this gene family. The Ka/Ks ratios of orthologous gene pairs between Arabidopsis thaliana and five Brassica species were mostly less than 1.0, implying that purifying selection, i.e., selective removal of deleterious alleles, played a role during the evolution of Brassica genomes. Analysis of 10 BnaBCCP genes using qRT-PCR showed a different pattern of expression because of exposure of the plants to biotic stresses, such as clubroot and sclerotinia diseases, and abiotic stresses such as drought, low temperature and salinity stresses.

Conclusions

The identification and functional analysis of the Brassica BCCPs demonstrated that some of these genes might play important roles in biotic and abiotic stress responses. Results from this study could lay the foundation for a better understanding of these genes for the improvement of Brassica crops for stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08920-y.

Overexpression of pea α-carboxyltransferase in Arabidopsis and camelina increases fatty acid synthesis leading to improved seed oil content

Heteromeric acetyl‐CoA carboxylase (htACCase) catalyzes the committed step of de novo fatty acid biosynthesis in most plant plastids. Plant htACCase is comprised of four subunits: α‐ and β‐carboxyltransferase (α‐ and β‐CT), biotin carboxylase, and biotin carboxyl carrier protein. Based on in vivo absolute quantification of htACCase subunits, α‐CT is 3‐ to 10‐fold less abundant than its partner subunit β‐CT in developing Arabidopsis seeds [Wilson and Thelen, J. Proteome Res., 2018, 17 (5)]. To test the hypothesis that low expression of α‐CT limits htACCase activity and flux through fatty acid synthesis in planta, we overexpressed Pisum sativum α‐CT, either with or without its C‐terminal non‐catalytic domain, in both Arabidopsis thaliana and Camelina sativa. First‐generation Arabidopsis seed of 35S::Ps α‐CT (n = 25) and 35S::Ps α‐CTΔ406‐875 (n = 47) were on average 14% higher in oil content (% dry weight) than wild type co‐cultivated in a growth chamber. First‐generation camelina seed showed an average 8% increase compared to co‐cultivated wild type. Biochemical analyses confirmed the accumulation of Ps α‐CT and Ps α‐CTΔ406‐875 protein and higher htACCase activity in overexpression lines during early seed development. Overexpressed Ps α‐CT co‐migrated with native At β‐CT during anion exchange chromatography, indicating co‐association. By successfully increasing seed oil content upon heterologous overexpression of α‐CT, we demonstrate how absolute quantitation of in vivo protein complex stoichiometry can be used to guide rational metabolic engineering.

Advances in Understanding the Genetic Basis of Fatty Acids Biosynthesis in Perilla: An Update

Perilla, also termed as purple mint, Chinese basil, or Perilla mint, is a flavoring herb widely used in East Asia. Both crude oil and essential oil are employed for consumption as well as industrial purposes. Fatty acids (FAs) biosynthesis and oil body assemblies in Perilla have been extensively investigated over the last three decades. Recent advances have been made in order to reveal the enzymes involved in the fatty acid biosynthesis in Perilla. Among those fatty acids, alpha-linolenic acid retained the attention of scientists mainly due to its medicinal and nutraceutical properties. Lipids synthesis in Perilla exhibited similarities with Arabidopsis thaliana lipids’ pathway. The homologous coding genes for polyunsaturated fatty acid desaturases, transcription factors, and major acyl-related enzymes have been found in Perilla via de novo transcriptome profiling, genome-wide association study, and in silico whole-genome screening. The identified genes covered de novo fatty acid synthesis, acyl-CoA dependent Kennedy pathway, acyl-CoA independent pathway, Triacylglycerols (TAGs) assembly, and acyl editing of phosphatidylcholine. In addition to the enzymes, transcription factors including WRINKLED, FUSCA3, LEAFY COTYLEDON1, and ABSCISIC ACID INSENSITIVE3 have been suggested. Meanwhile, the epigenome aspect impacting the transcriptional regulation of FAs is still unclear and might require more attention from the scientific community. This review mainly outlines the identification of the key gene master players involved in Perilla FAs biosynthesis and TAGs assembly that have been identified in recent years. With the recent advances in genomics resources regarding this orphan crop, we provided an updated overview of the recent contributions into the comprehension of the genetic background of fatty acid biosynthesis. The provided resources can be useful for further usage in oil-bioengineering and the design of alpha-linolenic acid-boosted Perilla genotypes in the future.

Gene loss, genome rearrangement, and accelerated substitution rates in plastid genome of Hypericum ascyron (Hypericaceae)

BACKGROUND

Comparative genomic analysis exhibits dynamic evolution of plastid genome (plastome) in the clusioid clade of Malpighiales, which comprise five families, including multiple inversions and gene losses. Little is known about the plastome evolution in Hypericaceae, a large family in the clade. Only the plastome of one species, Cratoxylum cochinchinense, has been published.

RESULTS

We generated a complete plastome sequence for Hypericum ascyron, providing the first complete plastome from the tribe Hypericeae (Hypericaceae). The H. ascyron plastome exhibits dynamic changes in gene and intron content, structure, and sequence divergence compared to the C. cochinchinense plastome from the tribe Cratoxyleae (Hypericaceae). Transcriptome data determined the evolutionary fate of the missing plastid genes infA, rps7, rps16, rpl23, and rpl32 in H. ascyron. Putative functional transfers of infA, rps7, and rpl32 were detected to the nucleus, whereas rps16 and rpl23 were substituted by nuclear-encoded homologs. The plastid rpl32 was integrated into the nuclear-encoded SODcp gene. Our findings suggested that the transferred rpl32 had undergone subfunctionalization by duplication rather than alternative splicing. The H. ascyron plastome rearrangements involved seven inversions, at least three inverted repeat (IR) boundary shifts, which generated gene relocations and duplications. Accelerated substitution rates of plastid genes were observed in the H. ascyron plastome compared with that of C. cochinchinense plastid genes. The higher substitution rates in the accD and clpP were correlated with structural change, including a large insertion of amino acids and losses of two introns, respectively. In addition, we found evidence of positive selection of the clpP, matK, and rps3 genes in the three branches related to H. ascyron. In particular, the matK gene was repeatedly under selection within the family Hypericaceae. Selective pressure in the H. ascyron matK gene was associated with the loss of trnK-UUU and relocation into the IR region.

CONCLUSIONS

The Hypericum ascyron plastome sequence provides valuable information for improving the understanding of plastome evolution among the clusioid of the Malpighiales. Evidence for intracellular gene transfer from the plastid to the nucleus was detected in the nuclear transcriptome, providing insight into the evolutionary fate of plastid genes in Hypericaceae.

Integrated transcriptome and proteome analysis of developing embryo reveals the mechanisms underlying the high levels of oil accumulation in Carya cathayensis Sarg

Hickory (Carya cathayensis Sarg.) is an extraordinary nut-bearing deciduous arbor with high content of oil in its embryo. However, the molecular mechanism underlying high oil accumulation is mostly unknown. Here, we reported that the lipid droplets and oil accumulation gradually increased with the embryo development and the oil content was up to ~76% at maturity. Furthermore, transcriptome and proteome analysis of developing hickory embryo identified 32,907 genes and 9857 proteins. Time-series analysis of gene expressions showed that these genes were divided into 12 clusters and lipid metabolism-related genes were enriched in Cluster 3, with the highest expression levels at 95 days after pollination (S2). Differentially expressed genes and proteins indicated high correlation, and both were enriched in the lipid metabolism. Notably, the genes involved in biosynthesis, transport of fatty acid/lipid and lipid droplets formation had high expression levels at S2, while the expression levels of other genes required for suberin/wax/cutin biosynthesis and lipid degradation were very low at all the sampling time points, ultimately promoting the accumulation of oil. Quantitative reverse-transcription PCR analysis also verified the results of RNA-seq. The co-regulatory networks of lipid metabolism were further constructed and WRINKLED1 (WRI1) was a core transcriptional factor located in the nucleus. Of note, CcWRI1A/B could directly activate the expression of some genes (CcBCCP2A, CcBCCP2B, CcFATA and CcFAD3) required for fatty acid synthesis. These results provided in-depth evidence for revealing the molecular mechanism of high oil accumulation in hickory embryo.

The plastome of the arctic Oxytropis arctobia (Fabaceae) has large differences compared with that of O. splendens and those of related species

Anatomical and physiological specializations for plant adaptation to harsh climates are the results of molecular mechanisms that can be nuclear or organellar encoded. In this study, the complete plastomes of an arctic species, Oxytropis arctobia Bunge (Fabaceae,) and a closely related temperate species, O. splendens Douglas ex Hook., were assembled, annotated and analyzed to search for differences that might help explain their adaptation to different environments. Consistently with the previously sequenced O. bicolor DC. and O. glabra plastomes, the O. arctobia and O. splendens plastomes both have the common features of the inverted repeat-lacking clade (IRLC), as well as the atpF intron loss, which is unique to the genus. However, significant differences distinguishes the O. arctobia from O. splendens and other closely related plastomes (Oxytropis spp. and Astragalus spp.), including a 3 kb inversion, two large insertions (>1 kb), significant modifications of the accD gene, and an overall larger size.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

Gene duplication and rate variation in the evolution of plastid ACCase and Clp genes in angiosperms

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome.

Biochemical and structural characterization of quizalofop-resistant wheat acetyl-CoA carboxylase

A novel nucleotide mutation in ACC1 resulting in an alanine to valine amino acid substitution in acetyl-CoA carboxylase (ACCase) at position 2004 of the Alopecurus myosuroides reference sequence (A2004V) imparts quizalofop resistance in wheat. Genotypes endowed with the homozygous mutation in one or two ACC1 homoeologs are seven- and 68-fold more resistant to quizalofop than a wildtype winter wheat in greenhouse experiments, respectively. In vitro ACCase activities in soluble protein extracts from these varieties are 3.8- and 39.4-fold more resistant to quizalofop with the homozygous mutation in either one or two genomes, relative to the wildtype. The A2004V mutation does not alter the specific activity of wheat ACCase, suggesting that this resistance trait does not affect the catalytic functions of ACCase. Modeling of wildtype and quizalofop-resistant wheat ACCase demonstrates that the A2004V amino acid substitution causes a reduction in the volume of the binding pocket that hinders quizalofop’s interaction with ACCase. Docking studies confirm that the mutation reduces the binding affinity of quizalofop. Interestingly, the models suggest that the A2004V mutation does not affect haloxyfop binding. Follow up in vivo and in vitro experiments reveal that the mutation, in fact, imparts negative cross-resistance to haloxyfop, with quizalofop-resistant varieties exhibiting higher sensitivity to haloxyfop than the wildtype winter wheat line.

Lineage-Specific Variation in IR Boundary Shift Events, Inversions, and Substitution Rates among Caprifoliaceae s.l. (Dipsacales) Plastomes

Caprifoliaceae s.l. plastid genomes (plastomes) show that one inversion and two inverted repeat boundary shifts occurred in the common ancestor of this family, after which the plastomes are generally conserved. This study reports plastome sequences of five additional species, Fedia cornucopiae, Valeriana fauriei, and Valerianella locusta from the subfamily Valerianoideae, as well as Dipsacus japonicus and Scabiosa comosa from the subfamily Dipsacoideae. Combined with the published plastomes, these plastomes provide new insights into the structural evolution of plastomes within the family. Moreover, the three plastomes from the subfamily Valerianoideae exhibited accelerated nucleotide substitution rates, particularly at synonymous sites, across the family. The patterns of accD sequence divergence in the family are dynamic with structural changes, including interruption of the conserved domain and increases in nonsynonymous substitution rates. In particular, the Valeriana accD gene harbors a large insertion of amino acid repeat (AAR) motifs, and intraspecific polymorphism with a variable number of AARs in the Valeriana accD gene was detected. We found a correlation between intron losses and increased ratios of nonsynonymous to synonymous substitution rates in the clpP gene with intensified positive selection. In addition, two Dipsacoideae plastomes revealed the loss of the plastid-encoded rps15, and a potential functional gene transfer to the nucleus was confirmed.

Characterization and Dynamics of Intracellular Gene Transfer in Plastid Genomes of Viola (Violaceae) and Order Malpighiales

Functional gene transfer from organelles to the nucleus, known as intracellular gene transfer (IGT), is an ongoing process in flowering plants. The complete plastid genomes (plastomes) of two Ulleung island endemic violets, Viola ulleungdoensis and V. woosanensis, were characterized, revealing a lack of the plastid-encoded infA, rpl32, and rps16 genes. In addition, functional replacement of the three plastid-encoded genes in the nucleus was confirmed within the genus Viola and the order Malpighiales. Three strategies for the acquisition of a novel transit peptide for successful IGT were identified in the genus Viola. Nuclear INFA acquired a novel transit peptide with very low identity between these proteins, whereas the nuclear RPL32 gene acquired an existing transit peptide via fusion with the nuclear-encoded plastid-targeted SOD gene (Cu-Zn superoxide dismutase superfamily) as one exon, and translated both proteins in the cytosol using alternative mRNA splicing. Nuclear RPS16 contains an internal transit peptide without an N-terminal extension. Gene loss or pseudogenization of the plastid-borne rpl32 and rps16 loci was inferred to occur in the common ancestor of the genus Viola based on an infrageneric phylogenetic framework in Korea. Although infA was lost in the common ancestor of the order Malpighiales, the rpl32 and rps16 genes were lost multiple times independently within the order. Our current study sheds additional light on plastid genome composition and IGT mechanisms in the violet genus and in the order Malpighiales.

Rampant Nuclear Transfer and Substitutions of Plastid Genes in Passiflora

Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora.

Transcriptomic analysis of α-linolenic acid content and biosynthesis in Paeonia ostii fruits and seeds

BACKGROUND

Paeonia ostii is a potentially important oilseed crop because its seed yield is high, and the seeds are rich in α-linolenic acid (ALA). However, the molecular mechanisms underlying ALA biosynthesis during seed kernel, seed testa, and fruit pericarp development in this plant are unclear. We used transcriptome data to address this knowledge gap.

RESULTS

Gas chromatograph-mass spectrometry indicated that ALA content was highest in the kernel, moderate in the testa, and lowest in the pericarp. Therefore, we used RNA-sequencing to compare ALA synthesis among these three tissues. We identified 227,837 unigenes, with an average length of 755 bp. Of these, 1371 unigenes were associated with lipid metabolism. The fatty acid (FA) biosynthesis and metabolism pathways were significantly enriched during the early stages of oil accumulation in the kernel. ALA biosynthesis was significantly enriched in parallel with increasing ALA content in the testa, but these metabolic pathways were not significantly enriched during pericarp development. By comparing unigene transcription profiles with patterns of ALA accumulation, specific unigenes encoding crucial enzymes and transcription factors (TFs) involved in de novo FA biosynthesis and oil accumulation were identified. Specifically, the bell-shaped expression patterns of genes encoding SAD, FAD2, FAD3, PDCT, PDAT, OLE, CLE, and SLE in the kernel were similar to the patterns of ALA accumulation in this tissue. Genes encoding BCCP, BC, KAS I- III, and FATA were also upregulated during the early stages of oil accumulation in the kernel. In the testa, the upregulation of the genes encoding SAD, FAD2, and FAD3 was followed by a sharp increase in the concentrations of ALA. In contrast, these genes were minimally expressed (and ALA content was low) throughout pericarp development.

CONCLUSIONS

We used three tissues with high, moderate, and low ALA concentrations as an exemplar system in which to investigate tissue-specific ALA accumulation mechanisms in P. ostii. The genes and TFs identified herein might be useful targets for future studies of ALA accumulation in the tree peony. This study also provides a framework for future studies of FA biosynthesis in other oilseed plants.

Explore the gene network regulating the composition of fatty acids in cottonseed

BACKGROUND

Cottonseed is one of the major sources of vegetable oil. Analysis of the dynamic changes of fatty acid components and the genes regulating the composition of fatty acids of cottonseed oil is of great significance for understanding the biological processes underlying biosynthesis of fatty acids and for genetic improving the oil nutritional qualities.

RESULTS

In this study, we investigated the dynamic relationship of 13 fatty acid components at 12 developmental time points of cottonseed (Gossypium hirsutum L.) and generated cottonseed transcriptome of the 12 time points. At 5-15 day post anthesis (DPA), the contents of polyunsaturated linolenic acid (C18:3n-3) and saturated stearic acid (C18:0) were higher, while linoleic acid (C18:2n-6) was mainly synthesized after 15 DPA. Using 5 DPA as a reference, 15,647 non-redundant differentially expressed genes were identified in 10-60 DPA cottonseed. Co-expression gene network analysis identified six modules containing 3275 genes significantly associated with middle-late seed developmental stages and enriched with genes related to the linoleic acid metabolic pathway and α-linolenic acid metabolism. Genes (Gh_D03G0588 and Gh_A02G1788) encoding stearoyl-ACP desaturase were identified as hub genes and significantly up-regulated at 25 DPA. They seemed to play a decisive role in determining the ratio of saturated fatty acids to unsaturated fatty acids. FAD2 genes (Gh_A13G1850 and Gh_D13G2238) were highly expressed at 25-50 DPA, eventually leading to the high content of C18:2n-6 in cottonseed. The content of C18:3n-3 was significantly decreased from 5 DPA (7.44%) to 25 DPA (0.11%) and correlated with the expression characteristics of Gh_A09G0848 and Gh_D09G0870.

CONCLUSIONS

These results contribute to our understanding on the relationship between the accumulation pattern of fatty acid components and the expression characteristics of key genes involved in fatty acid biosynthesis during the entire period of cottonseed development.

The evolutionary fate of rpl32 and rps16 losses in the Euphorbia schimperi (Euphorbiaceae) plastome

Gene transfers from mitochondria and plastids to the nucleus are an important process in the evolution of the eukaryotic cell. Plastid (pt) gene losses have been documented in multiple angiosperm lineages and are often associated with functional transfers to the nucleus or substitutions by duplicated nuclear genes targeted to both the plastid and mitochondrion. The plastid genome sequence of Euphorbia schimperi was assembled and three major genomic changes were detected, the complete loss of rpl32 and pseudogenization of rps16 and infA. The nuclear transcriptome of E. schimperi was sequenced to investigate the transfer/substitution of the rpl32 and rps16 genes to the nucleus. Transfer of plastid-encoded rpl32 to the nucleus was identified previously in three families of Malpighiales, Rhizophoraceae, Salicaceae and Passifloraceae. An E. schimperi transcript of pt SOD-1-RPL32 confirmed that the transfer in Euphorbiaceae is similar to other Malpighiales indicating that it occurred early in the divergence of the order. Ribosomal protein S16 (rps16) is encoded in the plastome in most angiosperms but not in Salicaceae and Passifloraceae. Substitution of the E. schimperi pt rps16 was likely due to a duplication of nuclear-encoded mitochondrial-targeted rps16 resulting in copies dually targeted to the mitochondrion and plastid. Sequences of RPS16-1 and RPS16-2 in the three families of Malpighiales (Salicaceae, Passifloraceae and Euphorbiaceae) have high sequence identity suggesting that the substitution event dates to the early divergence within Malpighiales.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Studies of an insecticidal inhibitor of acetyl-CoA carboxylase in the nematode C. elegans

We have studied the mode of action of the insecticide spirotetramat in the nematode Caenorhabditis elegans. A combination of symptomology, forward genetics and genome editing show that spirotetramat acts on acetyl-CoA carboxylase (ACC) in C. elegans, as it does in insects. We found C. elegans embryos exposed to spirotetramat show a cell division defect which closely resembles the phenotype of loss-of-function mutations in the gene pod-2, which encodes ACC. We then identified two mutations in the carboxyl transferase domain of pod-2 (ACC) which confer resistance and were confirmed using CRISPR/Cas9. One of these mutations substitutes an invertebrate-specific amino acid with one ubiquitous in other taxa; this residue may, therefore, be a determinant of the selectivity of spirotetramat for invertebrates. Such a mutation may also be the target of selection for resistance in the field. Our study is a further demonstration of the utility of C. elegans in studying bioactive chemicals.

Phosphorus supply via a fed-batch strategy improves lipid heterotrophic production of Chlorella regularis

Intracellular phosphorus (P) accumulation can improve microalgal growth and lipid synthesis. However, large excess of P causes cell poisoning. This study utilized a P-fed-batch strategy to investigate its potential to improve the utilization of the excessive P, while avoiding toxic side effects. This strategy contributed to a more complete utilization of the intracellularly stored P, which enhanced the microalgae biomass by 10-15% by upregulating the brassinosteroid growth hormone gene at a P-fed-batch frequency of 2-8. Furthermore, the lipid content increased by 4-16% via upregulation of lipid synthesis-related genes. As a result, the P-fed-batch strategy significantly increased the lipid production by 13-19%. The content of saturated fatty acid increased by ~ 100%, implying improved combustibility and oxidative stability. This is the first study of this P-fed-batch strategy and provides a new concept for the complete utilization of excessive P.

Molecular characteristics of the first case of haloxyfop-resistant Poa annua

Haloxyfop is one of two acetyl-coenzyme A carboxylase (ACCase) inhibitors that is recommended for controlling Poa annua. We have characterised a population of P. annua that had developed resistance to haloxyfop. This resistant population was found to be almost 20 times less sensitive to haloxyfop than a susceptible population based on percentage survival of individuals in two dose-response experiments. However, the haloxyfop-resistant population was still susceptible to clethodim. Pre-treatment of resistant individuals with a cytochrome P450 inhibitor, malathion, did not change the sensitivity level of the resistant plants to haloxyfop, suggesting that a non-target site mechanism of resistance involving enhanced metabolism, was not responsible for this resistance in P. annua. Gene sequencing showed that a target site mutation at position 2041, which replaced isoleucine with threonine in the carboxyltransferase (CT) domain of the ACCase enzyme, was associated with resistance to haloxyfop in the resistant population. An evaluation of the 3-D structure of the CT domain suggested that, unlike Asn-2041, which is the most common mutation at this position reported to date, Thr-2041 does not change the conformational structure of the CT domain. This is the first study investigating the molecular mechanism involved with haloxyfop resistance in P. annua.

Plant Unsaturated Fatty Acids: Biosynthesis and Regulation

In most plants, major unsaturated fatty acids (UFAs) are three C18 species, namely, oleic (18:1), linoleic (18:2), and α-linolenic (18:3) acids. These simple compounds play multiple crucial roles in planta and are also important economic traits of oil crops. The enzymatic steps of C18 UFA biosynthesis have been well established. However, the associated FA/lipid trafficking between the plastid and the endoplasmic reticulum remains largely unclear, as does the regulation of the expression and activities of the involved enzymes. In this review, we will revisit the biosynthesis of C18 UFAs with an emphasis on the trafficking, and present an overview of the key enzymes and their regulation. Of particular interest is the emerging regulatory network composed of transcriptional factors and upstream signaling pathways. The review thereby provides the promise of using physical, biochemical and/or genetic means to manipulate FA composition and increase oil yield in crop improvement.

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

Cross-resistance pattern to ACCase-inhibiting herbicides in a novel Trp1999Leu mutation American sloughgrass (Beckmannia syzigachne) population

The plastid acetyl coenzyme carboxylase (ACCase) Trp₁₉₉₉Leu mutation was identified in a Beckmannia syzigachne population resistant to fenoxaprop-p-ethyl. The pattern of cross-resistance for the Trp₁₉₉₉Leu mutation is still ambiguous. In this paper, mutant homozygote (1999Leu/Leu, RR) and wild type (1999Trp/Trp, SS) B. syzigachne plants with the same genetic background were purified from the JS-26 population using the dCAPS method. The activity of ACCase in RR and SS was determined. Then, the cross-resistance pattern to ACCase inhibiting herbicides of the Trp₁₉₉₉Leu mutation was determined using the whole-plant method. ACCase activity showed that the Trp₁₉₉₉Leu mutation decreased ACCase sensitivity to fenoxaprop-p-ethyl by 2.73-fold. A dose-response experiment indicated that the Trp₁₉₉₉Leu mutation conferred high resistance to quizalofop-p-ethyl (20.29-fold), metamifop (12.22-fold) and pinoxaden (18.60-fold), moderate resistance to fenoxaprop-p-ethyl (8.20-fold) and sethoxydim (6.38-fold), low resistance to cyhalofop-butyl (2.73-fold) and no resistance to clodinafop-propargyl (1.42 fold) and clethodim (1.59-fold). This is the first report of the role of Trp₁₉₉₉Leu in fenoxaprop-p-ethyl resistance and of the patterns of cross-resistance to ACCase-inhibiting herbicides in B. syzigachne.

Understanding evolution in Poales: Insights from Eriocaulaceae plastome

In this study, we report the plastome of Eriocaulon decemflorum (Eriocaulaceae) and make an effort to understand the genome evolution, structural rearrangements and gene content of the order Poales by comparing it with other available plastomes. The size of complete E. decemflorum plastome is 151,671 bp with an LSC (81,477bp), SSC (17,180bp) and a pair of IRs (26,507 bp). The plastome exhibits GC content of 35.8% and 134 protein-coding genes with 19 genes duplicated in the IR region. The Eriocaulaceae plastome is characterized by the presence of accD, ycf1 and ycf2 genes and presence of introns in clpP and rpoC1 genes which have been lost in the Graminid plastomes. Phylogenomic analysis based on 81 protein-coding genes placed Eriocaulaceae sister to Mayacaceae. The present study enhances our understanding of the evolution of Poales by analyzing the plastome data from the order.

Derived Polymorphic Amplified Cleaved Sequence (dPACS): A Novel PCR-RFLP Procedure for Detecting Known Single Nucleotide and Deletion-Insertion Polymorphisms

Most methods developed for detecting known single nucleotide polymorphisms (SNP) and deletion-insertion polymorphisms (DIP) are dependent on sequence conservation around the SNP/DIP and are therefore not suitable for application to heterogeneous organisms. Here we describe a novel, versatile and simple PCR-RFLP procedure baptised 'derived Polymorphic Amplified Cleaved Sequence' (dPACS) for genotyping individual samples. The notable advantage of the method is that it employs a pair of primers that cover the entire fragment to be amplified except for one or few diagnostic bases around the SNP/DIP being investigated. As such, it provides greater opportunities to introduce mismatches in one or both of the 35-55 bp primers for creating a restriction site that unambiguously differentiates wild from mutant sequences following PCR-RFLP and horizontal MetaPhorTM gel electrophoresis. Selection of effective restriction enzymes and primers is aided by the newly developed dPACS 1.0 software. The highly transferable dPACS procedure is exemplified here with the positive detection (in up to 24 grass and broadleaf species tested) of wild type proline106 of 5-enolpyruvylshikimate-3-phosphate synthase and its serine, threonine and alanine variants that confer resistance to glyphosate, and serine264 and isoleucine2041 which are key target-site determinants for weed sensitivities to some photosystem II and acetyl-CoA carboxylase inhibiting herbicides, respectively.

Role of CYP81A cytochrome P450s in clomazone metabolism in Echinochloa phyllopogon

Clomazone is a herbicide used in the cultivation of numerous crops due to its unique site of action and effectiveness on weeds. The differences in clomazone susceptibility among plants have been attributed to the differences in their complex clomazone metabolic pathways that are not fully understood. We previously identified two CYP81A cytochrome P450 monooxygenases that metabolize five chemically unrelated herbicides in multiple-herbicide resistant Echinochloa phyllopogon. Since the resistant E. phyllopogon have decreased clomazone susceptibility, involvement of these P450s in clomazone resistance was suggested. In this study, we revealed that each P450 gene endowed Arabidopsis thaliana (Arabidopsis) with clomazone resistance. Consistent with this, clomazone resistance co-segregated with resistance to other herbicides in F₆ progenies of crosses between susceptible and resistant E. phyllopogon, suggesting that the P450s are involved in differential clomazone susceptibility in E. phyllopogon. Arabidopsis transformations of the other seven CYP81As of E. phyllopogon found that two more genes, CYP81A15 and CYP81A24, decreased Arabidopsis susceptibility to clomazone. Differences in substrate preference between clomazone and a herbicide that inhibits acetolactate synthase were suggested among the four CYP81A P450s. This study provides insights into clomazone metabolism in plants.

Plastidial and Mitochondrial Malonyl CoA-ACP Malonyltransferase is Essential for Cell Division and Its Overexpression Increases Storage Oil Content

Malonyl-acyl carrier protein (ACP) is a key building block for the synthesis of fatty acids, which are important components of cell membranes, storage oils and lipid-signaling molecules. Malonyl CoA-ACP malonyltransferase (MCAMT) catalyzes the production of malonyl-ACP and CoA from malonyl-CoA and ACP. Here, we report that MCAMT plays a critical role in cell division and has the potential to increase the storage oil content in Arabidopsis. The quantitative real-time PCR and MCAMT promoter:GUS analyses showed that MCAMT is predominantly expressed in shoot and root apical meristems, leaf hydathodes and developing embryos. The fluorescent signals of MCAMT:eYFP were observed in both chloroplasts and mitochondria of tobacco leaf protoplasts. In particular, the N-terminal region (amino acid residues 1-30) of MCAMT was required for mitochondrial targeting. The Arabidopsis mcamt-1 and -2 mutants exhibited an embryo-lethal phenotype because of the arrest of embryo development at the globular stage. The transgenic Arabidopsis expressing antisense MCAMT RNA showed growth retardation caused by the defects in cell division. The overexpression of MCAMT driven by the promoter of the senescence-associated 1 (SEN1) gene, which is predominantly expressed in developing seeds, increased the seed yield and storage oil content of Arabidopsis. Taken together, the plastidial and mitochondrial MCAMT is essential for Arabidopsis cell division and is a novel genetic resource useful for enhancing storage oil content in oilseed crops.

Highly accelerated rates of genomic rearrangements and nucleotide substitutions in plastid genomes of Passiflora subgenus Decaloba

Plastid genomes (plastomes) of photosynthetic angiosperms are for the most part highly conserved in their organization, mode of inheritance and rates of nucleotide substitution. A small number of distantly related lineages share a syndrome of features that deviate from this general pattern, including extensive genomic rearrangements, accelerated rates of nucleotide substitution, biparental inheritance and plastome-genome incompatibility. Previous studies of plastomes in Passiflora with limited taxon sampling suggested that the genus exhibits this syndrome. To examine this phenomenon further, 15 new plastomes from Passiflora were sequenced and combined with previously published data to examine the phylogenetic relationships, genome organization and evolutionary rates across all five subgenera and the sister genus Adenia. Phylogenomic analyses using 68 protein-coding genes shared by Passiflora generated a fully resolved and strongly supported tree that is congruent with previous phylogenies based on a few plastid and nuclear loci. This phylogeny was used to examine the distribution of plastome rearrangements across Passiflora. Multiple gene and intron losses and inversions were identified in Passiflora with some occurring in parallel and others that extended across the Passifloraceae. Furthermore, extensive expansions and contractions of the inverted repeat (IR) were uncovered and in some cases this resulted in exclusion of all ribosomal RNA genes from the IR. The most highly rearranged lineage was subgenus Decaloba, which experienced extensive IR expansion that incorporated up to 25 protein-coding genes usually located in large single copy region. Nucleotide substitution rate analyses of 68 protein-coding genes across the genus showed lineage- and locus-specific acceleration. Significant increase in dS, dN and dN/dS was detected for clpP across the genus and for ycf4 in certain lineages. Significant increases in dN and dN/dS for ribosomal subunits and plastid-encoded RNA polymerase genes were detected in the branch leading to the expanded IR-clade in subgenus Decaloba. This subgenus displays the syndrome of unusual features, making it an ideal system to investigate the dynamic evolution of angiosperm plastomes.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.

Transcriptome Analysis Reveals Dynamic Fat Accumulation in the Walnut Kernel

Walnut (Juglans regia L.) is an important woody oilseed species cultivated throughout the world. In this study, comparative transcript profiling was performed using high-throughput RNA sequencing technology at the following three stages of walnut fat synthesis in the “Lvling” walnut cultivar: the initial developmental stage (L1), the fast developing stage (L2), and the last developing stage (L3). A total of 68.18 GB of data were obtained on the three developmental stages, and 92% to 94% of clean data were able to be located to the reference genome. Further comparisons of the transcripts in the three libraries revealed that 724, 2027, and 4817 genes were differentially expressed between the L2 and L1 (L2vsL1), L3 and L2 (L3vsL2), and L3 and L1 (L3vsL1) samples, respectively. Through the GO gene enrichment analysis, differentially expressed genes (DEGs) in L2vsL1, L3vsL2, and L3vsL1 were enriched into 3, 0, and 2 functional categories, respectively. According to the KEGG enrichment analysis, DEGs in L2vsL1, L3vsL2, and L3vsL1 were annotated into 77, 110, and 3717 taxonomic metabolic pathways in the KEGG database, respectively. Next, we analyzed expression levels of genes related to fat synthesis. Our results indicated that ACCase, LACS, and FAD7 were the key genes related to fat synthesis. The high-throughput transcriptome sequencing of walnut in different developmental stages has greatly enriched the current genomic available resources. The comparison of DEGs under different developmental stages identified a wealth of candidate genes involved in fat synthesis, which will facilitate further genetic improvement and molecular studies of the walnut.

The Repeat Sequences and Elevated Substitution Rates of the Chloroplast accD Gene in Cupressophytes

The plastid accD gene encodes a subunit of the acetyl-CoA carboxylase (ACCase) enzyme. The length of accD gene has been supposed to expand in Cryptomeria japonica, Taiwania cryptomerioides, Cephalotaxus, Taxus chinensis, and Podocarpus lambertii, and the main reason for this phenomenon was the existence of tandemly repeated sequences. However, it is still unknown whether the accD gene length in other cupressophytes has expanded. Here, in order to investigate how widespread this phenomenon was, 18 accD sequences and its surrounding regions of cupressophyte were sequenced and analyzed. Together with 39 GenBank sequence data, our taxon sampling covered all the extant gymnosperm orders. The repetitive elements and substitution rates of accD among 57 gymnosperm species were analyzed, the results show: (1) Reading frame length of accD gene in 18 cupressophytes species has also expanded. (2) Many repetitive elements were identified in accD gene of cupressophyte lineages. (3) The synonymous and non-synonymous substitution rates of accD were accelerated in cupressophytes. (4) accD was located in rearrangement endpoints. These results suggested that repetitive elements may mediate the chloroplast genome rearrangement and accelerated the substitution rates.

Transcriptomic analysis of Perilla frutescens seed to insight into the biosynthesis and metabolic of unsaturated fatty acids

Background

Perilla frutescens is well known for its high α-linolenic acid (ALA) accumulation in seeds and medicinal values as well as a source of edible and general-purpose oils. However, the regulatory mechanisms of the biosynthesis of fatty acid in its seeds remain poorly understood due to the lacking of sequenced genome. For better understanding the regulation of lipid metabolism and further increase its oil content or modify oil composition, time-course transcriptome and lipid composition analyses were performed.

Results

Analysis of fatty acid content and composition showed that the α-linolenic acid and oleic acid accumulated rapidly from 5 DAF to 15 DAF and then kept relatively stable. However, the amount of palmitic acid and linoleic acid decreased quickly from 5 DAF to 15DAF. No significant variation of stearic acid content was observed from 5 DAF to 25DAF. Our transcriptome data analyses revealed that 110,176 unigenes were generated from six seed libraries at 5, 10, 20 DAF. Of these, 53 (31 up, 22 down) and 653 (259 up, 394 down) genes showed temporal and differentially expression during the seed development in 5 DAF vs 10 DAF, 20 vs 10 DAF, respectively. The differentially expressed genes were annotated and found to be involved in distinct functional categories and metabolic pathways. Deep mining of transcriptome data led to the identification of key genes involved in fatty acid and triacylglycerol biosynthesis and metabolism. Thirty seven members of transcription factor family AP2, B3 and NFYB putatively involved in oil synthesis and deposition were differentially expressed during seed development. The results of qRT-PCR for selected genes showed a strong positive correlation with the expression abundance measured in RNA-seq analysis.

Conclusions

The present study provides valuable genomic resources for characterizing Perilla seed gene expression at the transcriptional level and will extend our understanding of the complex molecular and cellular events of oil biosynthesis and accumulation in oilseed crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4595-z) contains supplementary material, which is available to authorized users.

Molecular Evolution of Chloroplast Genomes of Orchid Species: Insights into Phylogenetic Relationship and Adaptive Evolution

Orchidaceae is the 3rd largest family of angiosperms, an evolved young branch of monocotyledons. This family contains a number of economically-important horticulture and flowering plants. However, the limited availability of genomic information largely hindered the study of molecular evolution and phylogeny of Orchidaceae. In this study, we determined the evolutionary characteristics of whole chloroplast (cp) genomes and the phylogenetic relationships of the family Orchidaceae. We firstly characterized the cp genomes of four orchid species: Cremastra appendiculata, Calanthe davidii, Epipactis mairei, and Platanthera japonica. The size of the chloroplast genome ranged from 153,629 bp (C. davidi) to 160,427 bp (E. mairei). The gene order, GC content, and gene compositions are similar to those of other previously-reported angiosperms. We identified that the genes of ndhC, ndhI, and ndhK were lost in C. appendiculata, in that the ndh I gene was lost in P. japonica and E. mairei. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in orchid species. E. mairei had the highest number of repeats (81), while C. davidii had the lowest number (57). The total number of Simple Sequence Repeats is at least 50 in C. davidii, and, at most, 78 in P. japonica. Interestingly, we identified 16 genes with positive selection sites (the psbH, petD, petL, rpl22, rpl32, rpoC1, rpoC2, rps12, rps15, rps16, accD, ccsA, rbcL, ycf1, ycf2, and ycf4 genes), which might play an important role in the orchid species’ adaptation to diverse environments. Additionally, 11 mutational hotspot regions were determined, including five non-coding regions (ndhB intron, ccsA-ndhD, rpl33-rps18, ndhE-ndhG, and ndhF-rpl32) and six coding regions (rps16, ndhC, rpl32, ndhI, ndhK, and ndhF). The phylogenetic analysis based on whole cp genomes showed that C. appendiculata was closely related to C. striata var. vreelandii, while C. davidii and C. triplicate formed a small monophyletic evolutionary clade with a high bootstrap support. In addition, five subfamilies of Orchidaceae, Apostasioideae, Cypripedioideae, Epidendroideae, Orchidoideae, and Vanilloideae, formed a nested evolutionary relationship in the phylogenetic tree. These results provide important insights into the adaptive evolution and phylogeny of Orchidaceae.

The complete chloroplast genome sequence of the CAM epiphyte Spanish moss (Tillandsia usneoides, Bromeliaceae) and its comparative analysis

Spanish moss (Tillandsia usneoides) is an epiphytic bromeliad widely distributed throughout tropical and warm temperate America. This plant is highly adapted to extreme environmental conditions. Striking features of this species include specialized trichomes (scales) covering the surface of its shoots aiding the absorption of water and nutrients directly from the atmosphere and a specific photosynthesis using crassulacean acid metabolism (CAM). Here we report the plastid genome of Spanish moss and present the comparison of genome organization and sequence evolution within Poales. The plastome of Spanish moss has a quadripartite structure consisting of a large single copy (LSC, 87,439 bp), two inverted regions (IRa and IRb, 26,803 bp) and short single copy (SSC, 18,612 bp) region. The plastid genome had 37.2% GC content and 134 genes with 88 being unique protein-coding genes and 20 of these are duplicated in the IR, similar to other reported bromeliads. Our study shows that early diverging lineages of Poales do not have high substitution rates as compared to grasses, and plastid genomes of bromeliads show structural features considered to be ancestral in graminids. These include the loss of the introns in the clpP and rpoC1 genes and the complete loss or partial degradation of accD and ycf genes in the Graminid clade. Further structural rearrangements appeared in the graminids lacking in Spanish moss, which include a 28-kb inversion between the trnG-UCC-rps14 region and 6-kb in the trnG-UCC-psbD, followed by a third <1kb inversion in the trnT sequence.

Acetyl-CoA carboxylase overexpression in herbicide-resistant large crabgrass (Digitaria sanguinalis)

BACKGROUND

The occurrence of herbicide-resistant weed biotypes is increasing and this report of an acetyl-CoA carboxylase (ACCase) inhibitor-resistant Digitaria sanguinalis L. Scop. from southwestern Ontario is another example. The identified weed escaped control in an onion and carrot rotation in which graminicides were used for several consecutive years. Our goal was to characterize the level and mechanism of resistance of the biotype.

RESULTS

The biotype was resistant to all five ACCase inhibitor herbicides tested. Gene-expression profiling was performed because none of the mutations known to confer resistance in the ACCase gene were detected. RNASeq and quantitative reverse-transcriptase PCR (qRT-PCR) results indicated that transcription of ACCase was 3.4-9.3 times higher in the resistant biotype than the susceptible biotype. ACCase gene copy number was determined by qPCR to be five to seven times higher in the resistant compared with the susceptible biotype. ACCase gene overexpression was directly related to the increase of the ACCase gene copy number.

CONCLUSION

Our results are consistent with the hypothesis that overexpression of the herbicide target gene ACCase confers resistance to the herbicide. This is the first reported case of target gene duplication conferring resistance to a herbicide other than glyphosate. © 2017 Society of Chemical Industry See related Article.

Fluazifop-P-butyl induced ROS generation with IAA (indole-3-acetic acid) oxidation in Acanthospermum hispidum D.C

Acanthospermum hispidum D.C. was particularly susceptible to fluazifop-P-butyl, an aryloxyphenoxypropionate herbicide, and the primary action site for the herbicide was shoot apical meristem, which is also the main site of indole-3-acetic acid (IAA) biosynthesis and action. Membrane lipid peroxidation caused by increasing levels of reactive oxygen species (ROS) was considered as an action mechanism of fluazifop-P-butyl in A. hispidum. To further clarify the ROS inducing mechanism of fluazifop-P-butyl in the plant, the interactions between fluazifop-P-butyl and auxin compounds IAA or 2,4-dichlorophenoxyacetic acid (2,4-D) were studied. Haloxyfop-P-methyl, an AOPP herbicide which is inactive on A. hispidum, was used for comparison. The results showed that the growth inhibition and malondialdehyde or H₂O₂ increases induced by fluazifop-P-butyl on A. hispidum were reversed by IAA or 2,4-D. The IAA content was decreased but the contents of three IAA oxidation metabolites, indole-3-methanol, indole-3-aldehyde and indole-3-carboxylic acid were increased by fluazifop-P-butyl in A. hispidum, but not by haloxyfop-P-methyl. The growth of A. hispidum was not inhibited by three IAA oxidative compounds. Moreover, the activities of IAA oxidase and peroxidase were increased by fluazifop-P-butyl but not by haloxyfop-P-methyl, and the increase was reversed by IAA or 2,4-D. We suggest that there is an antagonistic effect between fluazifop-P-butyl and IAA or 2,4-D, and the IAA oxidation may be involved in the action mechanism of fluazifop-P-butyl in A. hispidum.