Two long-chain acyl-CoA synthetases from Arabidopsis thaliana involved in peroxisomal fatty acid beta-oxidation

Evolutionary origin and gradual accumulation with plant evolution of the LACS family

BACKGROUND

LACS (long-chain acyl-CoA synthetase) genes are widespread in organisms and have multiple functions in plants, especially in lipid metabolism. However, the origin and evolutionary dynamics of the LACS gene family remain largely unknown.

RESULTS

Here, we identified 1785 LACS genes in the genomes of 166 diverse plant species and identified the clades (I, II, III, IV, V, VI) of six clades for the LACS gene family of green plants through phylogenetic analysis. Based on the evolutionary history of plant lineages, we found differences in the origins of different clades, with Clade IV originating from chlorophytes and representing the origin of LACS genes in green plants. The structural characteristics of different clades indicate that clade IV is relatively independent, while the relationships between clades (I, II, III) and clades (V, VI) are closer. Dispersed duplication (DSD) and transposed duplication (TRD) are the main forces driving the evolution of plant LACS genes. Network clustering analysis further grouped all LACS genes into six main clusters, with genes within each cluster showing significant co-linearity. Ka/Ks results suggest that LACS family genes underwent purifying selection during evolution. We analyzed the phylogenetic relationships and characteristics of six clades of the LACS gene family to explain the origin, evolutionary history, and phylogenetic relationships of different clades and proposed a hypothetical evolutionary model for the LACS family of genes in plants.

CONCLUSIONS

Our research provides genome-wide insights into the evolutionary history of the LACS gene family in green plants. These insights lay an important foundation for comprehensive functional characterization in future research.

The acyl-CoA synthetase TgACS1 allows neutral lipid metabolism and extracellular motility in Toxoplasma gondii through relocation via its peroxisomal targeting sequence (PTS) under low nutrient conditions

Apicomplexa parasites cause major diseases such as toxoplasmosis and malaria that have major health and economic burdens. These unicellular pathogens are obligate intracellular parasites that heavily depend on lipid metabolism for the survival within their hosts. Their lipid synthesis relies on an essential combination of fatty acids (FAs) obtained from both de novo synthesis and scavenging from the host. The constant flux of scavenged FA needs to be channeled toward parasite lipid storage, and these FA storages are timely mobilized during parasite division. In eukaryotes, the utilization of FA relies on their obligate metabolic activation mediated by acyl-co-enzyme A (CoA) synthases (ACSs), which catalyze the thioesterification of FA to a CoA. Besides the essential functions of FA for parasite survival, the presence and roles of ACS are yet to be determined in Apicomplexa. Here, we identified TgACS1 as a Toxoplasma gondii cytosolic ACS that is involved in FA mobilization in the parasite specifically during low host nutrient conditions, especially in extracellular stages where it adopts a different localization. Heterologous complementation of yeast ACS mutants confirmed TgACS1 as being an Acyl-CoA synthetase of the bubble gum family that is most likely involved in β-oxidation processes. We further demonstrate that TgACS1 is critical for gliding motility of extracellular parasite facing low nutrient conditions, by relocating to peroxisomal-like area.IMPORTANCEToxoplasma gondii, causing human toxoplasmosis, is an Apicomplexa parasite and model within this phylum that hosts major infectious agents, such as Plasmodium spp., responsible for malaria. The diseases caused by apicomplexans are responsible for major social and economic burdens affecting hundreds of millions of people, like toxoplasmosis chronically present in about one-third of the world's population. Lack of efficient vaccines, rapid emergence of resistance to existing treatments, and toxic side effects of current treatments all argue for the urgent need to develop new therapeutic tools to combat these diseases. Understanding the key metabolic pathways sustaining host-intracellular parasite interactions is pivotal to develop new efficient ways to kill these parasites. Current consensus supports parasite lipid synthesis and trafficking as pertinent target for novel treatments. Many processes of this essential lipid metabolism in the parasite are not fully understood. The capacity for the parasites to sense and metabolically adapt to the host physiological conditions has only recently been unraveled. Our results clearly indicate the role of acyl-co-enzyme A (CoA) synthetases for the essential metabolic activation of fatty acid (FA) used to maintain parasite propagation and survival. The significance of our research is (i) the identification of seven of these enzymes that localize at different cellular areas in T. gondii parasites; (ii) using lipidomic approaches, we show that TgACS1 mobilizes FA under low host nutrient content; (iii) yeast complementation showed that acyl-CoA synthase 1 (ACS1) is an ACS that is likely involved in peroxisomal β-oxidation; (iv) the importance of the peroxisomal targeting sequence for correct localization of TgACS1 to a peroxisomal-like compartment in extracellular parasites; and lastly, (v) that TgACS1 has a crucial role in energy production and extracellular parasite motility.

Carotenoid retention during post-harvest storage of Capsicum annuum: the role of the fruit surface structure

In this study, a chilli pepper (Capsicum annuum) panel for post-harvest carotenoid retention was studied to elucidate underlying mechanisms associated with this commercial trait of interest. Following drying and storage, some lines within the panel had an increase in carotenoids approaching 50% compared with the initial content at the fresh fruit stage. Other lines displayed a 25% loss of carotenoids. The quantitative determination of carotenoid pigments with concurrent cellular analysis indicated that in most cases, pepper fruit with thicker (up to 4-fold) lipid exocarp layers and smooth surfaces exhibit improved carotenoid retention properties. Total cutin monomer content increased in medium/high carotenoid retention fruits and subepidermal cutin deposits were responsible for the difference in exocarp thickness. Cutin biosynthesis and cuticle precursor transport genes were differentially expressed between medium/high and low carotenoid retention genotypes, and this supports the hypothesis that the fruit cuticle can contribute to carotenoid retention. Enzymatic degradation of the cuticle and cell wall suggests that in Capsicum the carotenoids (capsanthin and its esters) are embedded in the lipidic exocarp layer. This was not the case in tomato. Collectively, the data suggest that the fruit cuticle could provide an exploitable resource for the enhancement of fruit quality.

Genome-wide association study and pathway analysis to decipher loci associated with Fusarium ear rot resistance in tropical maize germplasm

Breeding for host resistance is the most efficient and environmentally safe method to curb the spread of fusarium ear rot (FER). However, conventional breeding for resistance to FER is hampered by the complex polygenic nature of this trait, which is highly influenced by environmental conditions. This study aimed to identify genomic regions, single nucleotide polymorphisms (SNPs), and putative candidate genes associated with FER resistance as well as candidate metabolic pathways and pathway genes involved in it. A panel of 151 tropical inbred maize lines were used to assess the genetic architecture of FER resistance over two seasons. During the study period, seven SNPs associated with FER resistance were identified on chromosomes 1, 2, 4, 5, and 9, accounting for 4-11% of the phenotypic variance. These significant markers were annotated into four genes. Seven significant metabolic pathways involved in FER resistance were identified using the Pathway Association Study Tool, the most significant being the superpathway of the glyoxylate cycle. Overall, this study confirmed that resistance to FER is indeed a complex mechanism controlled by several small to medium-effect loci. Our findings may contribute to fast-tracking the efforts to develop disease-resistant maize lines through marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10722-023-01793-4.

Metabolomics-mediated elucidation of rice responses to salt stress

MAIN CONCLUSION

Role of salinity responsive metabolites of rice and its wild species has been discussed. Salinity stress is one of the important environmental stresses that severely affects rice productivity. Although, several vital physio-biochemical and molecular responses have been activated in rice under salinity stress which were well described in literatures, the mechanistic role of salt stress and microbes-induced metabolites to overcome salt stress in rice are less studied. Nevertheless, over the years, metabolomic studies have allowed a comprehensive analyses of rice salt stress responses. Hence, we review the salt stress-triggered alterations of various metabolites in rice and discuss their significant roles toward salinity tolerance. Some of the metabolites such as serotonin, salicylic acid, ferulic acid and gentisic acid may act as signaling molecules to activate different downstream salt-tolerance mechanisms; whereas, the other compounds such as amino acids, sugars and organic acids directly act as protective agents to maintain osmotic balance and scavenger of reactive oxygen species during the salinity stress. The quantity, type, tissues specificity and time of accumulation of metabolites induced by salinity stress vary between salt-sensitive and tolerant rice genotypes and thus, contribute to their different degrees of salt tolerance. Moreover, few tolerance metabolites such as allantoin, serotonin and melatonin induce unique pathways for activation of defence mechanisms in salt-tolerant varieties of rice, suggesting their potential roles as the universal biomarkers for salt tolerance. Therefore, these metabolites can be applied exogenously to the sensitive genotypes of rice to enhance their performance under salt stress. Furthermore, the microbes of rhizosphere also participated in rice salt tolerance either directly or indirectly by regulating their metabolic pathways. Thus, this review for the first time offers valuable and comprehensive insights into salt-induced spatio-temporal and genotype-specific metabolites in different genotypes of rice which provide a reference point to analyze stress-gene-metabolite relationships for the biomarker designing in rice. Further, it can also help to decipher several metabolic systems associated with salt tolerance in rice which will be useful in developing salt-tolerance cultivars by conventional breeding/genetic engineering/exogenous application of metabolites.

Mapping the castor bean endosperm proteome revealed a metabolic interaction between plastid, mitochondria, and peroxisomes to optimize seedling growth

In this work, we studied castor-oil plant Ricinus communis as a classical system for endosperm reserve breakdown. The seeds of castor beans consist of a centrally located embryo with the two thin cotyledons surrounded by the endosperm. The endosperm functions as major storage tissue and is packed with nutritional reserves, such as oil, proteins, and starch. Upon germination, mobilization of the storage reserves requires inter-organellar interplay of plastids, mitochondria, and peroxisomes to optimize growth for the developing seedling. To understand their metabolic interactions, we performed a large-scale organellar proteomic study on castor bean endosperm. Organelles from endosperm of etiolated seedlings were isolated and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS). Computer-assisted deconvolution algorithms were applied to reliably assign the identified proteins to their correct subcellular localization and to determine the abundance of the different organelles in the heterogeneous protein samples. The data obtained were used to build a comprehensive metabolic model for plastids, mitochondria, and peroxisomes during storage reserve mobilization in castor bean endosperm.

Genome-Wide Identification, Characterization, and Expression Analysis of Long-Chain Acyl-CoA Synthetases in Carya illinoinensis under Different Treatments

As crucial enzymes in the lipid metabolic network, long-chain acyl-CoA synthases (LACSs) are members of the acyl-activated enzyme superfamily and play a crucial role in epidermal wax synthesis, plant lipid anabolic metabolism, and stress tolerance. In this study, 11 pecan LACS genes were identified and categorized into five groups and located on nine chromosomes. The significant degree of conservation in the AtLACS and CiLACS protein sequences was demonstrated by multiple sequence alignment and conserved motif analysis. Cis-acting element analysis identified numerous stress-responsive and hormone-inducible elements in the promoter regions of CiLACS genes. The expression levels of CiLACS9 and CiLACS9-1 were considerably up-regulated under salt and drought stress, according to the qRT-RCR study. Treatment with ABA also led to increased expression levels of CiLACS1, CiLACS1-1, CiLACS2, and CiLACS9-1. Notably, CiLACS4, CiLACS4-1, CiLACS9, and CiLACS9-1 exhibited peak expression levels at 135 days after anthesis and are likely to have been crucial in the accumulation of seed kernel oil. Moreover, the CiLACS9 gene was shown to be located in the cytoplasm. These findings offer a theoretical framework for clarifying the roles of LACS genes in the processes of pecan kernel oil synthesis and response to abiotic stressors.

Genome-Wide Analysis and Functional Characterization of LACS Gene Family Associated with Lipid Synthesis in Cotton (Gossypium spp.)

Cotton (Gossypium spp.) is the fifth largest oil crop in the world, and cottonseed provides abundant vegetable oil resources and industrial bioenergy fuels for people; therefore, it is of practical significance to increase the oil content of cotton seeds for improving the oil yield and economic benefits of planting cotton. Long-chain acyl-coenzyme A (CoA) synthetase (LACS) capable of catalyzing the formation of acyl-CoAs from free fatty acids has been proven to significantly participate in lipid metabolism, of which whole-genome identification and functional characterization of the gene family have not yet been comprehensively analyzed in cotton. In this study, a total of sixty-five LACS genes were confirmed in two diploid and two tetraploid Gossypium species, which were divided into six subgroups based on phylogenetic relationships with twenty-one other plants. An analysis of protein motif and genomic organizations displayed structural and functional conservation within the same group but diverged among the different group. Gene duplication relationship analysis illustrates the LACS gene family in large scale expansion through WGDs/segmental duplications. The overall Ka/Ks ratio indicated the intense purifying selection of LACS genes in four cotton species during evolution. The LACS genes promoter elements contain numerous light response cis-elements associated with fatty acids synthesis and catabolism. In addition, the expression of almost all GhLACS genes in high seed oil were higher compared to those in low seed oil. We proposed LACS gene models and shed light on their functional roles in lipid metabolism, demonstrating their engineering potential for modulating TAG synthesis in cotton, and the genetic engineering of cottonseed oil provides a theoretical basis.

Peroxisomal KAT2 (3-ketoacyl-CoA thiolase 2) gene has a key role in gingerol biosynthesis in ginger (Zingiber officinale Rosc.)

Ginger is an important spice crop with medicinal values and gingerols are the most abundant pungent polyphenols present in ginger, responsible for most of its pharmacological properties. The present study focuses on the molecular mechanism of gingerol biosynthesis in ginger using transcriptome analysis. Suppression Subtractive Hybridization (SSH) was done in leaf and rhizome tissues using high gingerol-producing ginger somaclone B3 as the tester and parent cultivar Maran as the driver and generated high-quality leaf and rhizome Expressed Sequence Tags (ESTs). The Blast2GO annotations of the ESTs revealed the involvement of leaf ESTs in secondary metabolite production, identifying the peroxisomal KAT2 gene (Leaf EST 9) for the high gingerol production in ginger. Rhizome ESTs mostly coded for DNA metabolic processes and differential genes for high gingerol production were not observed in rhizomes. In the qRT-PCR analysis, somaclone B3 had shown high chalcone synthase (CHS: rate-limiting gene in gingerol biosynthetic pathway) activity (0.54 fold) in the leaves of rhizome sprouts. The presence of a high gingerol gene in leaf ESTs and high expression of CHS in leaves presumed that the site of synthesis of gingerols in ginger is the leaves. A modified pathway for gingerol/polyketide backbone formation has been constructed explaining the involvement of KAT gene isoforms KAT2 and KAT5 in gingerol/flavonoid biosynthesis, specifically the KAT2 gene which is otherwise thought to be involved mainly in β-oxidation. The results of the present investigations have the potential of utilizing KAT/thiolase superfamily enzymes for protein/metabolic pathway engineering in ginger for large-scale production of gingerols.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13562-022-00825-x.

Overexpression the BnLACS9 could increase the chlorophyll and oil content in Brassica napus

BACKGROUND

Chlorophyll is a very important pigment involved in photosynthesis, while plant acyl-CoA biosynthesis is derived from plastid-localized fatty acids (FAs). Until now, the regulation of the acyl-CoA pathway for chlorophyll biosynthesis is still unknown.

RESULTS

Here, we identified a long-chain acyl-CoA synthetase (LACS) gene BnLACS9 from Brassica napus. BnLACS9 complemented a LACS-deficient yeast strain YB525, which indicated that BnLACS9 has the LACS function. BnLACS9 was localized in the chloroplast envelope membrane, while mainly expressed in young leaves and flowers. Overexpression of BnLACS9 in Nicotiana benthamiana resulted in an increase in total CoA and MGDG content. In B. napus with overexpression of BnLACS9, the number of chloroplast grana lamellae and the chlorophyll content, as well as the MGDG and DGDG contents, increased compared to wild type. The net photosynthetic rate, dry weight of the entire plant and oil content of seeds increased significantly, accompanied by an increase in chlorophyll content. Transcriptome analysis revealed that overexpression of BnLACS9 improved the pathway of acyl-CoA biosynthesis and further improved the enzymes in the glycolipid synthesis pathway, while acyl-CoA was the substrate for glycolipid synthesis. The increased glycolipids, especially MGDG and DGDG, accelerated the formation of the chloroplast grana lamellae, which increased the number of chloroplast thylakoid grana lamella and further lead to increased chlorophyll content.

CONCLUSIONS

In the present study, we demonstrated that BnLACS9 played a crucial role in glycolipids and chlorophyll biosynthesis in B. napus. The results also provide a new direction and theoretical basis for the improvement of the agronomic traits of plants.

Plastidic membrane lipids are oxidized by a lipoxygenase in Lobosphaera incisa

Green microalgae can accumulate neutral lipids, as part of a general lipid remodeling mechanism under stress such as nitrogen starvation. Lobosphaera incisa is of special interest because of its unique TAG acyl chain composition, especially 20:4 (n-6) can reach up to 21% of dry weight after nitrogen starvation. In order to identify factors that may influence the accumulation of polyunsaturated fatty acids (PUFAs), we identified recently a linoleate 13-lipoxygenase (LiLOX). It shares highest identity with plastidic enzymes from vascular plants and is induced upon nitrogen starvation. Here, we confirmed the localization of LiLOX in the stroma of plastids via transient expression in epithelial onion cells. In order to further characterize this enzyme, we focused on the identification of the endogenous substrate of LiLOX. In this regard, an ex vivo enzymatic assay, coupled with non-targeted analysis via mass spectrometry allowed the identification of MGDG, DGDG and PC as three substrate candidates, later confirmed via in vitro assays. Further investigation revealed that LiLOX has preferences towards the lipid class MGDG, which seems in agreement with its localization in the galactolipid rich plastid. Altogether, this study shows the first characterization of plastidic LOX from green algae, showing preference for MGDGs. However, lipidomics analysis did neither reveal an endogenous LiLOX product nor the final end product of MGDG oxidation. Nevertheless, the latter is a key to understanding the role of this enzyme and since its expression is highest during the degradation of the plastidic membrane, it is tempting to assume its involvement in this process.

The mechanism analysis of exogenous melatonin in limiting pear fruit aroma decrease under low temperature storage

Exogenous melatonin (MT) is widely used in fruit preservation, and can increase the storage time and delay the quality deterioration. Firstly, it was found that 150 μM MT was the optimal concentration to treat 'Xinli No.7' under storage at 4 °C for 60 days. MT could significantly improve oxidase activity and inhibit the reduction of physiological indexes, including pulp hardness, weight loss, titratable acid and soluble solid content. MT could also reduce ethylene release and limit the reduction of fruit aroma. The average content of fruit aroma substance increased by 43.53%. A relevant RNA-Seq database was built to further explore the regulation mechanism of MT. A total of 2,761 differentially expressed genes (DEGs) were identified. DEGs were enriched in 64 functional groups and 191 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. DEGs were mainly enriched in alpha-linolenic acid metabolism, fatty acid metabolism and plant hormone signal transduction pathway. The gene pycom09g05270 belonging to long chain acyl-CoA synthetase family and participating in fatty acid metabolism pathway was identified, and its expression level was consistent with fragments per kilobase per million mapped reads (FPKM) values, implying that pycom09g05270 might play a vital role in maintaining quality during the storage process.

Long-Chain Acyl-CoA Synthetases Promote Poplar Resistance to Abiotic Stress by Regulating Long-Chain Fatty Acid Biosynthesis

Long-chain acyl-CoA synthetases (LACSs) catalyze fatty acids (FAs) to form fatty acyl-CoA thioesters, which play essential roles in FA and lipid metabolisms and cuticle wax biosynthesis. Although LACSs from Arabidopsis have been intensively studied, the characterization and function of LACSs from poplar are unexplored. Here, 10 poplar PtLACS genes were identified from the poplar genome and distributed to eight chromosomes. A phylogenetic tree indicated that PtLACSs are sorted into six clades. Collinearity analysis and duplication events demonstrated that PtLACSs expand through segmental replication events and experience purifying selective pressure during the evolutionary process. Expression patterns revealed that PtLACSs have divergent expression changes in response to abiotic stress. Interaction proteins and GO analysis could enhance the understanding of putative interactions among protein and gene regulatory networks related to FA and lipid metabolisms. Cluster networks and long-chain FA (LCFA) and very long-chain FA (VLCFA) content analysis revealed the possible regulatory mechanism in response to drought and salt stresses in poplar. The present study provides valuable information for the functional identification of PtLACSs in response to abiotic stress metabolism in poplar.

Synthesis of C20-38 Fatty Acids in Plant Tissues

Very-long-chain fatty acids (VLCFA) are involved in a number of important plant physiological functions. Disorders in the expression of genes involved in the synthesis of VLCFA lead to a number of phenotypic consequences, ranging from growth retardation to the death of embryos. The elongation of VLCFA in the endoplasmic reticulum (ER) is carried out by multiple elongase complexes with different substrate specificities and adapted to the synthesis of a number of products required for a number of metabolic pathways. The information about the enzymes involved in the synthesis of VLCFA with more than 26 atoms of Carbon is rather poor. Recently, genes encoding enzymes involved in the synthesis of both regular-length fatty acids and VLCFA have been discovered and investigated. Polyunsaturated VLCFA in plants are formed mainly by 20:1 elongation into new monounsaturated acids, which are then imported into chloroplasts, where they are further desaturated. The formation of saturated VLCFA and their further transformation into a number of aliphatic compounds included in cuticular waxes and suberin require the coordinated activity of a large number of different enzymes.

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.

Identification of Arabidopsis Protein Kinases That Harbor Functional Type 1 Peroxisomal Targeting Signals

Peroxisomes are eukaryotic specific organelles that perform diverse metabolic functions including fatty acid β-oxidation, reactive species metabolism, photorespiration, and responses to stress. However, the potential regulation of these functions by post-translational modifications, including protein phosphorylation, has had limited study. Recently, we identified and catalogued a large number of peroxisomal phosphorylated proteins, implicating the presence of protein kinases in this organelle. Here, we employed available prediction models coupled with sequence conservation analysis to identify 31 protein kinases from the Arabidopsis kinome (all protein kinases) that contain a putative, non-canonical peroxisomal targeting signal type 1 (PTS1). From this, twelve C-terminal domain-PTS1s were demonstrated to be functional in vivo, targeting enhanced yellow fluorescent protein to peroxisomes, increasing the list of presumptive peroxisomal protein kinases to nineteen. Of the twelve protein kinases with functional PTS1s, we obtained full length clones for eight and demonstrated that seven target to peroxisomes in vivo. Screening homozygous mutants of the presumptive nineteen protein kinases revealed one candidate (GPK1) that harbors a sugar-dependence phenotype, suggesting it is involved in regulating peroxisomal fatty acid β-oxidation. These results present new opportunities for investigating the regulation of peroxisome functions.

Long-chain acyl-CoA synthetases activate fatty acids for lipid synthesis, remodeling and energy production in Chlamydomonas

Long-chain acyl-CoA synthetases (LACSs) play many roles in mammals, yeasts and plants, but knowledge on their functions in microalgae remains fragmented. Here via genetic, biochemical and physiological analyses, we unraveled the function and roles of LACSs in the model microalga Chlamydomonas reinhardtii. In vitro assays on purified recombinant proteins revealed that CrLACS1, CrLACS2 and CrLACS3 all exhibited bona fide LACS activities toward a broad range of free fatty acids. The Chlamydomonas mutants compromised in CrLACS1, CrLACS2 or CrLACS3 did not show any obvious phenotypes in lipid content or growth under nitrogen (N)-replete condition. But under N-deprivation, CrLACS1 or CrLACS2 suppression resulted in c. 50% less oil, yet with a higher amount of chloroplast lipids. By contrast, CrLACS3 suppression impaired oil remobilization and cell growth severely during N-recovery, supporting its role in fatty acid β-oxidation to provide energy and carbon sources for regrowth. Transcriptomics analysis suggested that the observed lipid phenotypes are likely not due to transcriptional reprogramming but rather a shift in metabolic adjustment. Taken together, this study provided solid experimental evidence for essential roles of the three Chlamydomonas LACS enzymes in lipid synthesis, remodeling and catabolism, and highlighted the importance of lipid homeostasis in cell growth under nutrient fluctuations.

Transcriptome and physiological analyses provide insights into the leaf epicuticular wax accumulation mechanism in yellowhorn

Plantations and production of yellowhorn, one of the most important woody oil and urban greening trees widely cultivated in northern China, have gradually become limited by drought stress. The epicuticular wax layer plays a key role in the protection of yellowhorn trees from drought and other stresses. However, there is no research on the mechanism of wax loading in yellowhorn trees. In this study, we investigated the anatomical and physiological characteristics of leaves from different germplasm resources and different parts of the same tree and compared their cuticle properties. In addition, the different expression patterns of genes involved in wax accumulation were analyzed, and a coexpression network was built based on transcriptome sequencing data. Morphological and physiological comparisons found that the sun leaves from the outer part of the crown had thicker epicuticular wax, which altered the permeability and improved the drought resistance of leaves, than did shade leaves. Based on transcriptome data, a total of 3008 and 1324 differentially expressed genes (DEGs) were identified between the sun leaves and shade leaves in glossy- and non-glossy-type germplasm resources, respectively. We identified 138 DEGs involved in wax biosynthesis and transport, including structural genes (such as LACS8, ECH1, and ns-LTP) and transcription factors (such as MYB, WRKY, and bHLH transcription factor family proteins). The coexpression network showed a strong correlation between these DEGs. The differences in gene expression patterns between G- and NG-type germplasm resources under different light conditions were very clear. These results not only provide a theoretical basis for screening and developing drought-resistant yellowhorn germplasm resources but also provide a data platform to reveal the wax accumulation process of yellowhorn leaves.

Genome-wide comparative analysis of long-chain acyl-CoA synthetases (LACSs) gene family: A focus on identification, evolution and expression profiling related to lipid synthesis

In plants, Long-chain acyl-CoA synthetases (LACSs) play key roles in activating fatty acids to fatty acyl-CoA thioesters, which are then further involved in lipid synthesis and fatty acid catabolism. LACSs have been intensively studied in Arabidopsis, but its evolutionary relationship in green plants is unexplored. In this study, we performed a comprehensive genome-wide analysis of the LACS gene family across green plants followed by phylogenetic clustering analysis, gene structure determination, detection of conserved motifs, gene expression in tissues and subcellular localization. Our results identified LACS genes in 122 plant species including algae, low land plants (i.e., mosses and lycophytes), monocots, and eudicots. In total, 697 sequences were identified, and 629 sequences were selected because of alignment and some duplication errors. The retrieved amino acid sequences ranged from 271 to 1056 residues and diversified in intron/exon patterns in different LACSs. Phylogenetic clustering grouped LACS gene family into six major clades with distinct potential functions. This classification is well supported by examining gene structure and conserved motifs. Also, gene expression analysis and subcellular localization substantiate with clade division in the phylogeny, indicating that the evolutionary pattern is visible in their functionality. Additionally, experimental analysis of lacs2 mutant validated that LACS2 plays key roles in suberin synthesis. Thus, our study not only provides an evolutionary mechanism underlying functional diversification but also lays the foundation for further elucidation of the LACS gene family.

Dynamic transcriptome analysis identifies genes related to fatty acid biosynthesis in the seeds of Prunus pedunculata Pall

BACKGROUND

Prunus pedunculata Pall, the deciduous shrub of Amygdalus subgenus in Rosaceae, is a new kind of desert oil-bearing tree. It has a long story of being planted in the West and North of China for sand fixation and desert control. In addition, the seeds of P. pedunculata are rich of oil, especially the monounsaturated fatty acid and polyunsaturated fatty acid. However, little is known about the molecular mechanisms of oil accumulation during the seed development of P. pedunculata.

RESULTS

The seeds of P. pedunculata from three independent plants at 10, 18, 24, 31, 39, 45, 59 and 73 days after flowering (DAF) were obtained and the oil compositions were evaluated. It showed that oleic acid was the dominant type of oil content in the mature seeds (from 32.724% at 10DAF to 72.06% at 73DAF). Next, transcriptome sequencing for the developing seeds produced 988.795 million high quality reads and TRINITY assembled 326,271 genes for the first transcriptome for P. pedunculata. After the assembled transcriptome was evaluated by BUSCO with 85.9% completeness, we identified 195,342, 109,850 and 121,897 P. pedunculata genes aligned to NR, GO and KEGG pathway databases, respectively. Then, we predicted 23,229 likely proteins from the assembled transcriptome and identified 1917 signal peptides and 5512 transmembrane related proteins. In the developing seeds we detected 91,362 genes (average FPKM > 5) and correlation analysis indicated three possible development stages - early (10 ~ 24DAF), middle (31 ~ 45DAF) and late (59 ~ 73DAF). We next analyzed the differentially expressed genes (DEGs) in the developing seeds. Interestingly, compared to 10DAF the number of DEGs was increased from 4406 in 18DAF to 27,623 in 73DAF. Based on the gene annotation, we identified 753, 33, 8 and 645 DEGs related to the fatty acid biosynthesis, lipid biosynthesis, oil body and transcription factors. Notably, GPAT, DGD1, LACS2, UBC and RINO were highly expressed at the early development stage, ω6-FAD, SAD, ACP, ACCA and AHG1 were highly expressed at the middle development stage, and LACS6, DGD1, ACAT1, AGPAT, WSD1, EGY2 and oleosin genes were highly expressed at the late development stage.

CONCLUSIONS

This is the first time to study the developing seed transcriptome of P. pedunculata and our findings will provide a valuable resource for future studies. More importantly, it will improve our understanding of molecular mechanisms of oil accumulation in P. pedunculata.

Acyl-CoA thioesterase activity of peroxisomal ABC protein ABCD1 is required for the transport of very long-chain acyl-CoA into peroxisomes

The ABCD1 protein, one of the four ATP-binding cassette (ABC) proteins in subfamily D, is located on the peroxisomal membrane and is involved in the transport of very long chain fatty acid (VLCFA)-CoA into peroxisomes. Its mutation causes X-linked adrenoleukodystophy (X-ALD): an inborn error of peroxisomal β-oxidation of VLCFA. Whether ABCD1 transports VLCFA-CoA as a CoA ester or free fatty acid is controversial. Recently, Comatose (CTS), a plant homologue of human ABCD1, has been shown to possess acyl-CoA thioesterase (ACOT) activity, and it is suggested that this activity is required for transport of acyl-CoA into peroxisomes. However, the precise transport mechanism is unknown. Here, we expressed human His-tagged ABCD1 in methylotrophic yeast, and characterized its ACOT activity and transport mechanism. The expressed ABCD1 possessed both ATPase and ACOT activities. The ACOT activity of ABCD1 was inhibited by p-chloromercuribenzoic acid (pCMB), a cysteine-reactive compound. Furthermore, we performed a transport assay with ABCD1-containing liposomes using 7-nitro-2–1,3-benzoxadiazol-4-yl (NBD)-labeled acyl-CoA as the substrate. The results showed that the fatty acid produced from VLCFA-CoA by ABCD1 is transported into liposomes and that ACOT activity is essential during this transport process. We propose a detailed mechanism of VLCFA-CoA transport by ABCD1.

Functional Role of Long-Chain Acyl-CoA Synthetases in Plant Development and Stress Responses

Fatty acids (FAs) play vital roles in plants as components of lipid membranes that demarcate cells and organelles, as sources of stored energy in the form of neutral lipids, and as signaling molecules that elicit plant responses to adverse conditions. The activation of FAs through the formation of acyl-CoA intermediates by acyl-CoA synthetase (ACS) family enzymes is required for their synthesis and degradation. Long-chain ACSs (LACSs) represent a small subgroup of ACS enzymes that specifically convert long-chain or very-long-chain FAs into corresponding thioesters for multiple lipid-associated processes. Alteration of LACS activity often results in pleiotropic phenotypes such as male sterility, organ fusion, aberrant cuticular structure, delayed seed germination, altered seed oil content, and plant capacity to respond to various environmental stresses. This review provides a comprehensive analysis of LACS family enzymes including substrate specificity, tissue-specific expression patterns, and distinct subcellular localization highlighting their specific roles in lipid synthesis and degradation, the effects of altered LACS activity on plant development, the relationship between LACS activity and stress resistance, and the regulation of LACS activity. Finally, we pose several major questions to be addressed, which would advance our current understanding of LACS function in plants.

Evolutionary Maintenance of the PTS2 Protein Import Pathway in the Stramenopile Alga Nannochloropsis

The stramenopile alga Nannochloropsis evolved by secondary endosymbiosis of a red alga by a heterotrophic host cell and emerged as a promising organism for biotechnological applications, such as the production of polyunsaturated fatty acids and biodiesel. Peroxisomes play major roles in fatty acid metabolism but experimental analyses of peroxisome biogenesis and metabolism in Nannochloropsis are not reported yet. In fungi, animals, and land plants, soluble proteins of peroxisomes are targeted to the matrix by one of two peroxisome targeting signals (type 1, PTS1, or type 2, PTS2), which are generally conserved across kingdoms and allow the prediction of peroxisomal matrix proteins from nuclear genome sequences. Because diatoms lost the PTS2 pathway secondarily, we investigated its presence in the stramenopile sister group of diatoms, the Eustigmatophyceae, represented by Nannochloropsis. We detected a full-length gene of a putative PEX7 ortholog coding for the cytosolic receptor of PTS2 proteins and demonstrated its expression in Nannochloropsis gaditana. The search for predicted PTS2 cargo proteins in N. gaditana yielded several candidates. In vivo subcellular targeting analyses of representative fusion proteins in different plant expression systems demonstrated that two predicted PTS2 domains were indeed functional and sufficient to direct a reporter protein to peroxisomes. Peroxisome targeting of the predicted PTS2 cargo proteins was further confirmed in Nannochloropsis oceanica by confocal and transmission electron microscopy. Taken together, the results demonstrate for the first time that one group of stramenopile algae maintained the import pathway for PTS2 cargo proteins. To comprehensively map and model the metabolic capabilities of Nannochloropsis peroxisomes, in silico predictions needs to encompass both the PTS1 and the PTS2 matrix proteome.

A combination of genome-wide association study and transcriptome analysis in leaf epidermis identifies candidate genes involved in cuticular wax biosynthesis in Brassica napus

BACKGROUND

Brassica napus L. is one of the most important oil crops in the world. However, climate-change-induced environmental stresses negatively impact on its yield and quality. Cuticular waxes are known to protect plants from various abiotic/biotic stresses. Dissecting the genetic and biochemical basis underlying cuticular waxes is important to breed cultivars with improved stress tolerance.

RESULTS

Here a genome-wide association study (GWAS) of 192 B. napus cultivars and inbred lines was used to identify single-nucleotide polymorphisms (SNPs) associated with leaf waxes. A total of 202 SNPs was found to be significantly associated with 31 wax traits including total wax coverage and the amounts of wax classes and wax compounds. Next, epidermal peels from leaves of both high-wax load (HW) and low-wax load (LW) lines were isolated and used to analyze transcript profiles of all GWAS-identified genes. Consequently, 147 SNPs were revealed to have differential expressions between HW and LW lines, among which 344 SNP corresponding genes exhibited up-regulated while 448 exhibited down-regulated expressions in LW when compared to those in HW. According to the gene annotation information, some differentially expressed genes were classified into plant acyl lipid metabolism, including fatty acid-related pathways, wax and cutin biosynthesis pathway and wax secretion. Some genes involved in cell wall formation and stress responses have also been identified.

CONCLUSIONS

Combination of GWAS with transcriptomic analysis revealed a number of directly or indirectly wax-related genes and their associated SNPs. These results could provide clues for further validation of SNPs for marker-assisted breeding and provide new insights into the genetic control of wax biosynthesis and improving stress tolerance of B. napus.

Lipid Composition and Associated Gene Expression Patterns during Pollen Germination and Pollen Tube Growth in Olive (Olea europaea L.)

Pollen lipids are essential for sexual reproduction, but our current knowledge regarding lipid dynamics in growing pollen tubes is still very scarce. Here, we report unique lipid composition and associated gene expression patterns during olive pollen germination. Up to 376 genes involved in the biosynthesis of all lipid classes, except suberin, cutin and lipopolysaccharides, are expressed in olive pollen. The fatty acid profile of olive pollen is markedly different compared with other plant organs. Triacylglycerol (TAG), containing mostly C12-C16 saturated fatty acids, constitutes the bulk of olive pollen lipids. These compounds are partially mobilized, and the released fatty acids enter the β-oxidation pathway to yield acetyl-CoA, which is converted into sugars through the glyoxylate cycle during the course of pollen germination. Our data suggest that fatty acids are synthesized de novo and incorporated into glycerolipids by the 'eukaryotic pathway' in elongating pollen tubes. Phosphatidic acid is synthesized de novo in the endomembrane system during pollen germination and seems to have a central role in pollen tube lipid metabolism. The coordinated action of fatty acid desaturases FAD2-3 and FAD3B might explain the increase in linoleic and alpha-linolenic acids observed in germinating pollen. Continuous synthesis of TAG by the action of diacylglycerol acyltransferase 1 (DGAT1) enzyme, but not phosphoplipid:diacylglycerol acyltransferase (PDAT), also seems plausible. All these data allow for a better understanding of lipid metabolism during the olive reproductive process, which can impact, in the future, on the increase in olive fruit yield and, therefore, olive oil production.

Multi-targeted trehalose-6-phosphate phosphatase I harbors a novel peroxisomal targeting signal 1 and is essential for flowering and development

This work reveals information about new peroxisomal targeting signals type 1 and identifies trehalose-6-phosphate phosphatase I as multitargeted and is implicated in plant development, reproduction, and stress response. A putative, non-canonical peroxisomal targeting signal type 1 (PTS1) Pro-Arg-Met > was identified in the extreme C-terminus of trehalose-6-phosphate phosphatase (TPP)I. TPP catalyzes the final step of trehalose synthesis, and the enzyme was previously characterized to be nuclear only (Krasensky et al. in Antioxid Redox Signal 21(9):1289-1304, 2014). Here we show that the TPPI C-terminal decapeptide ending with Pro-Arg-Met > or Pro-Lys-Met > can indeed function as a PTS1. Upon transient expression in two plant expression systems, the free C- or N-terminal end led to the full-length TPPI targeting to peroxisomes and plastids, respectively. The nucleus and nucleolus targeting of the full-length TPPI was observed in both cases. The homozygous T-DNA insertion line of TPPI showed a pleiotropic phenotype including smaller leaves, shorter roots, delayed flowering, hypersensitivity to salt, and a sucrose dependent seedling development. Our results identify novel PTS1s, and TPPI as a protein multi-targeted to peroxisomes, plastids, nucleus, and nucleolus. Altogether our findings implicate an essential role for TPPI in development, reproduction, and cell signaling.

Functional Analysis of Rice Long-Chain Acyl-CoA Synthetase 9 (OsLACS9) in the Chloroplast Envelope Membrane

The long-chain acyl-CoA synthetases (LACSs) are involved in lipid synthesis, fatty acid catabolism, and the transport of fatty acids between subcellular compartments. These enzymes catalyze the critical reaction of fatty acyl chains to fatty acyl-CoAs for the triacylglycerol biosynthesis used as carbon and energy reserves. In Arabidopsis, LACSs are encoded by a family of nine genes, with LACS9 being the only member located in the chloroplast envelope membrane. However, the comprehensive role of LACS9 and its contribution to plant metabolism have not been explored thoroughly. In this study, we report on the identification and characterization of LACS9 mutants in rice plants. Our results indicate that the loss-of-function mutations in OsLACS9 affect the architecture of internodes resulting in dwarf plants with large starch granules in the chloroplast, showing the suppression of starch degradation. Moreover, the plastid localization of α-amylase I-1 (AmyI-1)-a key enzyme involved in starch breakdown in plastids-was suppressed in the lacs9 mutant line. Immunological and confocal laser scanning microscopy analyses showed that OsLACS9-GFP is located in the chloroplast envelope in green tissue. Microscopic analysis showed that OsLACS9s interact with each other in the plastid envelope membrane. Furthermore, OsLACS9 is also one of the proteins transported to plastids without a transit peptide or involvement of the Toc/Tic complex system. To identify the plastid-targeting signal of OsLACS9, the transient expression and localization of a series of N-terminal truncated OsLACS9-green fluorescent protein (GFP) fusion proteins were examined. Truncation analyses identified the N-terminal 30 amino acid residues to be required for OsLACS9 plastid localization. Overall, the data in this study provide an advanced understanding of the function of OsLACS9 and its role in starch degradation and plant growth.

Systems biology approach identifies functional modules and regulatory hubs related to secondary metabolites accumulation after transition from autotrophic to heterotrophic growth condition in microalgae

Heterotrophic growth mode is among the most promising strategies put forth to overcome the low biomass and secondary metabolites productivity challenge. To shedding light on the underlying molecular mechanisms, transcriptome meta-analysis was integrated with weighted gene co-expression network analysis (WGCNA), connectivity analysis, functional enrichment, and hubs identification. Meta-analysis and Functional enrichment analysis demonstrated that most of the biological processes are up-regulated at heterotrophic growth condition, which leads to change of genetic architectures and phenotypic outcomes. WGNCA analysis of meta-genes also resulted four significant functional modules across logarithmic (LG), transition (TR), and production peak (PR) phases. The expression pattern and connectivity characteristics of the brown module as a non-preserved module vary across LG, TR, and PR phases. Functional analysis identified Carotenoid biosynthesis, Fatty acid metabolism and Methane metabolism as enriched pathways in the non-preserved module. Our integrated approach was applied here, identified some hubs, such as a serine hydroxymethyltransferase (SHMT1), which is the best candidate for development of metabolites accumulating strains in microalgae. Current study provided a new insight into underlying metabolite accumulation mechanisms and opens new avenue for the future applied studies in the microalgae field.

Fatty Acyl Synthetases and Thioesterases in Plant Lipid Metabolism: Diverse Functions and Biotechnological Applications

Plants use fatty acids to synthesize acyl lipids for many different cellular, physiological, and defensive roles. These roles include the synthesis of essential membrane, storage, or surface lipids, as well as the production of various fatty acid-derived metabolites used for signaling or defense. Fatty acids are activated for metabolic processing via a thioester linkage to either coenzyme A or acyl carrier protein. Acyl synthetases metabolically activate fatty acids to their thioester forms, and acyl thioesterases deactivate fatty acyl thioesters to free fatty acids by hydrolysis. These two enzyme classes therefore play critical roles in lipid metabolism. This review highlights the surprisingly complex and varying roles of fatty acyl synthetases in plant lipid metabolism, including roles in the intracellular trafficking of fatty acids. This review also surveys the many specialized fatty acyl thioesterases characterized to date in plants, which produce a great diversity of fatty acid products in a tissue-specific manner. While some acyl thioesterases produce fatty acids that clearly play roles in plant-insect or plant-microbial interactions, most plant acyl thioesterases have yet to be fully characterized both in terms of their substrate specificities and their functions. The biotechnological applications of plant acyl thioesterases and synthetases are also discussed, as there is significant interest in these enzymes as catalysts for the sustainable production of fatty acids and their derivatives for industrial uses.

Long-Chain acyl-CoA Synthetase LACS2 Contributes to Submergence Tolerance by Modulating Cuticle Permeability in Arabidopsis

In Arabidopsis thaliana, LONG-CHAIN ACYL-COA SYNTHETASEs (LACSs) catalyze the synthesis of long-chain acyl-CoAs and function in diverse biological processes. We have recently revealed that LACS2 is primarily involved in the production of polyunsaturated linolenoyl-CoA, essential for the activation of ethylene response transcription factors-mediated hypoxia signaling. Here, we further reported the dual role of LACS2 in the regulation of submergence tolerance by modulating cuticle permeability in Arabidopsis cells. LACS2-overexpressors (LACS2-OEs) showed improved tolerance to submergence, with higher accumulation of cuticular wax and cutin in their rosettes. In contrast, knockout of LACS2 in the lacs2-3 mutant resulted in hypersensitivity to submergence with reduced wax crystals and thinner cutin layer. By analyses of plant surface permeability, we observed that the hypoxic sensitivities in the LACS2-OEs and lacs2-3 mutant were physiologically correlated with chlorophyll leaching, water loss rates, ionic leakage, and gas exchange. Thus, our findings suggest the role of LACS2 in plant response to submergence by modulating cuticle permeability in plant cells.

Using Transcriptome Analysis to Screen for Key Genes and Pathways Related to Cytoplasmic Male Sterility in Cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is one of the most important cash crops worldwide. Cytoplasmic male sterility (CMS) is an excellent breeding system for exploitation of heterosis, which has great potential to increase crop yields. To understand the molecular mechanism of CMS in cotton, we compared transcriptome, cytomorphological, physiological and bioinformatics data between the CMS line C2P5A and its maintainer line C2P5B. By using high-throughput sequencing technology, 178,166 transcripts were assembled and 2013 differentially expression genes (DEGs) were identified at three different stages of C2P5A anther development. In this study, we identified DEGs associated with reactive oxygen species (ROS), peroxisomes, aldehyde dehydrogenases (ALDH), cytochrome oxidase subunit VI, and cytochrome P450, and DEGs associated with tapetum development, Jojoba acyl-CoA reductase-related male sterility protein, basic helix-loop-helix (bHLH) and MYB transcription factors. The abnormal expression of one of these genes may be responsible for the CMS C2P5A line. In gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, DEGs were mainly related to carbohydrate metabolism, amino acid metabolism, transport and catabolism, and signal transduction. Carbohydrate metabolism provides energy for anther development, starch and sucrose metabolism, fatty acid biosynthesis and metabolism and ascorbate and aldarate metabolism. These results showed that numerous genes and multiple complex metabolic pathways regulate cotton anther development. Weighted correlation network analysis (WGCNA) indicated that three modules, ‘turquoise,’ ‘blue,’ and ‘green,’ were specific for the CMS C2P5A line. The ‘turquoise’ and ‘blue’ modules were mainly related to carbohydrate metabolism, amino acid metabolism, energy metabolism, peroxisomes, pyruvate metabolism as well as fatty acid degradation. The ‘green’ module was mainly related to energy metabolism, carbon metabolism, translation, and lipid metabolism. RNA-sequencing and WGCNA polymerization modules were screened for key genes and pathways related to CMS in cotton. This study presents a new perspective for further research into the metabolic pathways of pollen abortion in the CMS C2P5A line and also provides a theoretical basis for its breeding and production.

A systems genetics analysis in Eucalyptus reveals coordination of metabolic pathways associated with xylan modification in wood-forming tissues

Acetyl- and methylglucuronic acid decorations of xylan, the dominant hemicellulose in secondary cell walls (SCWs) of woody dicots, affect its interaction with cellulose and lignin to determine SCW structure and extractability. Genes and pathways involved in these modifications may be targets for genetic engineering; however, little is known about the regulation of xylan modifications in woody plants. To address this, we assessed genetic and gene expression variation associated with xylan modification in developing xylem of Eucalyptus grandis × Eucalyptus urophylla interspecific hybrids. Expression quantitative trait locus (eQTL) mapping identified potential regulatory polymorphisms affecting gene expression modules associated with xylan modification. We identified 14 putative xylan modification genes that are members of five expression modules sharing seven trans-eQTL hotspots. The xylan modification genes are prevalent in two expression modules. The first comprises nucleotide sugar interconversion pathways supplying the essential precursors for cellulose and xylan biosynthesis. The second contains genes responsible for phenylalanine biosynthesis and S-adenosylmethionine biosynthesis required for glucuronic acid and monolignol methylation. Co-expression and co-regulation analyses also identified four metabolic sources of acetyl coenxyme A that appear to be transcriptionally coordinated with xylan modification. Our systems genetics analysis may provide new avenues for metabolic engineering to alter wood SCW biology for enhanced biomass processability.

New guidelines for fluorophore application in peroxisome targeting analyses in transient plant expression systems

Peroxisome research has been revolutionized by proteome studies combined with in vivo subcellular targeting analyses. Yellow and cyan fluorescent protein (YFP and CFP) are the classical fluorophores of plant peroxisome research. In the new transient expression system of Arabidopsis seedlings co-cultivated with Agrobacterium we detected the YFP fusion of one candidate protein in peroxisomes, but only upon co-transformation with the peroxisome marker, CFP-PTS1. The data suggested that the YFP fusion was directed to peroxisomes due to its weak heterodimerization ability with CFP-PTS1, allowing piggy-back import into peroxisomes. Indeed, if co-expressed with monomeric Cerulean-PTS1 (mCer-PTS1), the YFP fusion was no longer matrix localized. We systematically investigated the occurrence and extent of dimerization-based piggy-back import for different fluorophore combinations in five major transient plant expression systems. In Arabidopsis seedlings and tobacco leaves both untagged YFP and monomeric Venus were imported into peroxisomes if co-expressed with CFP-PTS1 but not with mCer-PTS1. By contrast, piggy-back import of cytosolic proteins was not observed in Arabidopsis and tobacco protoplasts or in onion epidermal cells for any fluorophore combination at any time point. Based on these important results we formulate new guidelines for fluorophore usage and experimental design to guarantee reliable identification of novel plant peroxisomal proteins.

Molecular and biochemical analysis of the castor caruncle reveals a set of unique genes involved in oil accumulation in non-seed tissues

BACKGROUND

With the increasing demand for vegetative oil and the approach of peak seed oil production, it is important to develop new oil production platforms from non-seed tissues. Castor bean (Ricinus communis) is one of the crops for vegetable oil for industrial applications with yield around 1.4 ton oil per hectare produced in seed. The castor caruncle is a non-seed tissue attached to seed.

RESULTS

Caruncle accumulates up to 40% oil by weight in the form of triacylglycerol (TAG), with a highly contrasting fatty acid composition when compared to the seed oil. Biochemical analysis indicated that the caruncle synthesizes TAGs independent of the seed. Such non-seed tissue has provided an excellent resource for understanding the mechanism of oil accumulation in tissues other than seeds. Transcriptome analysis revealed the key members of gene families involved in fatty acid synthesis and TAG assembly in the caruncle. A transient expression assay of these selected genes resulted in a 20-fold increased TAG accumulation in leaves.

CONCLUSIONS

Castor caruncle utilizes an independent system to synthesize TAGs. Results provide the possibility of exploiting caruncle gene set to engineer oil production in non-seed tissues or microbes.

Functional Overlap of Long-Chain Acyl-CoA Synthetases in Arabidopsis

Long-chain acyl-CoA synthetases (LACSs) play diverse and essential roles in lipid metabolism. The genomes of model eukaryotic organisms encode multiple LACS genes, and the substrate specificities of LACS homologs often overlap substantially. Homologous LACSs tend to differ in their expression patterns, localizations, and, by extension, the metabolic pathways to which they contribute. The Arabidopsis genome encodes a family of nine LACS genes, which have been characterized largely by reverse genetic analysis of mutant phenotypes. Because of redundancy, distinguishing the contributions of some Arabidopsis LACS genes has been challenging. Here, we have attempted to clarify the functions of LACSs that functionally overlap by synopsizing the results of previous work, isolating a suite of higher-order mutants that were previously lacking, and analyzing oil, wax, cutin, cuticle permeability, fertility and growth phenotypes. LACS1, LACS2, LACS4, LACS8 and LACS9 all affect cuticular lipid metabolism, but have different precise roles. Seed set, seed weight and storage oil amounts of higher-order lacs1, lacs2, lacs4, lacs8 and lacs9 mutants vary greatly, with these traits subject to different effects of fertility and oil synthesis defects. LACS4, LACS8 and LACS9 have partially redundant roles in development, as lacs4 lacs8 and lacs4 lacs9 double mutants are dwarf. lacs4 lacs8 lacs9 triple mutants were not recovered, and are assumed to be non-viable. Together, these results sketch a complex network of functions and functional interactions within the Arabidopsis LACS gene family.

A Framework to Investigate Peroxisomal Protein Phosphorylation in Arabidopsis

Peroxisomes perform essential roles in a range of cellular processes, highlighted by lipid metabolism, reactive species detoxification, and response to a variety of stimuli. The ability of peroxisomes to grow, divide, respond to changing cellular needs, interact with other organelles, and adjust their proteome as required, suggest that, like other organelles, their specialized roles are highly regulated. Similar to most other cellular processes, there is an emerging role for protein phosphorylation to regulate these events. In this review, we establish a knowledge framework of key players that control protein phosphorylation events in the plant peroxisome (i.e., the protein kinases and phosphatases), and highlight a vastly expanded set of (phospho)substrates.

Engineering Arabidopsis long-chain acyl-CoA synthetase 9 variants with enhanced enzyme activity

Long-chain acyl-CoA synthetase (LACS, EC 6.2.1.3) catalyzes the ATP-dependent activation of free fatty acid to form acyl-CoA, which, in turn, serves as the major acyl donor for various lipid metabolic pathways. Increasing the size of acyl-CoA pool by enhancing LACS activity appears to be a useful approach to improve the production and modify the composition of fatty acid-derived compounds, such as triacylglycerol. In the present study, we aimed to improve the enzyme activity of Arabidopsis thaliana LACS9 (AtLACS9) by introducing random mutations into its cDNA using error-prone PCR. Two AtLACS9 variants containing multiple amino acid residue substitutions were identified with enhanced enzyme activity. To explore the effect of each amino acid residue substitution, single-site mutants were generated and the amino acid substitutions C207F and D238E were found to be primarily responsible for the increased activity of the two variants. Furthermore, evolutionary analysis revealed that the beneficial amino acid site C207 is conserved among LACS9 from plant eudicots, whereas the other beneficial amino acid site D238 might be under positive selection. Together, our results provide valuable information for the production of LACS variants for applications in the metabolic engineering of lipid biosynthesis in oleaginous organisms.

Characterization of Fatty Acid Exporters involved in fatty acid transport for oil accumulation in the green alga Chlamydomonas reinhardtii

BACKGROUND

In the past few decades, microalgae biofuel has become one of the most interesting sources of renewable energy. However, the higher cost of microalgae biofuel compared to that of petroleum prevented microalgae biofuel production. Therefore, the research on increasing lipid productivity from microalgae becomes more important. The lipid production source, triacylglycerol biosynthesis in microalgae requires short chain fatty acids as substrates, which are synthesized in chloroplasts. However, the transport mechanism of fatty acids from microalgae chloroplasts to cytosol remains unknown.

RESULTS

cDNAs from two homologs of the Arabidopsis fatty acid exporter 1 (FAX1) were cloned from Chlamydomonas reinhardtii and were named crfax1 and crfax2. Both CrFAXs were involved in fatty acid transport, and their substrates were mainly C16 and C18 fatty acids. Overexpression of both CrFAXs increased the accumulation of the total lipid content in algae cells, and the fatty acid compositions were changed under normal TAP or nitrogen deprivation conditions. Overexpression of both CrFAXs also increased the chlorophyll content. The MGDG content was decreased but the TAG, DAG, DGDG and other lipid contents were increased in CrFAXs overexpression strains.

CONCLUSION

These results reveal that CrFAX1 and CrFAX2 were involved in mediating fatty acid export for lipids biosynthesis in C. reinhardtii. In addition, overexpression of both CrFAXs obviously increased the intracellular lipid content, especially the triacylglycerol content in microalgae, which provides a potential technology for the production of more biofuels using microalgae.

Genome-wide identification and characterization of apple long-chain Acyl-CoA synthetases and expression analysis under different stresses

Long-chain acyl-CoA synthetases (LACSs) are members of the acyl-activating enzyme superfamily that have important roles in lipid synthesis and storage, fatty acid catabolism, vectorial acylation, and synthesis of cutin and wax. Here, 11 apple MdLACS genes were identified based on the Malus × domestica reference genome, clustered into six groups and mapped to ten chromosomes. Multiple sequence alignment and conserved motifs analyses showed that the sequences of the AtLACS and MdLACS proteins were highly conserved. A cis-element analysis in the promoter regions of the MdLACS genes revealed various elements related to stress responsiveness and plant hormones. Subsequently, expression analysis demonstrated that the MdLACS genes had different expression profiles in different tissues in response to various abiotic stresses. To further study the function of MdLACS genes in apple, MdLACS1 was isolated to identify its basic function, which the function of MdLACS1 in response to apple abiotic stress resistance was determined by the transgenic method. The results showed the MdLACS1 enhanced tolerance to polyethylene glycol, salt, and abscisic acid in the apple callus, suggesting that MdLACS1 is an important regulator in response to abiotic stresses. Finally, the functional interoperability network among the MdLACS proteins was predicted and analyzed, which could the understanding of the possible interactions among proteins and genes regulatory networks concerned with wax biosynthesis and regulatory mechanisms in response to abiotic stresses in apple.

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty-one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome-wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome-wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well-known tocopherol pathway genes would greatly improve tocopherol content in maize.

Substrate preferences of long-chain acyl-CoA synthetase and diacylglycerol acyltransferase contribute to enrichment of flax seed oil with α-linolenic acid

Seed oil from flax (Linum usitatissimum) is enriched in α-linolenic acid (ALA; 18:3Δ^{9cis,12cis,15cis} ), but the biochemical processes underlying the enrichment of flax seed oil with this polyunsaturated fatty acid are not fully elucidated. Here, a potential process involving the catalytic actions of long-chain acyl-CoA synthetase (LACS) and diacylglycerol acyltransferase (DGAT) is proposed for ALA enrichment in triacylglycerol (TAG). LACS catalyzes the ATP-dependent activation of free fatty acid to form acyl-CoA, which in turn may serve as an acyl-donor in the DGAT-catalyzed reaction leading to TAG. To test this hypothesis, flax LACS and DGAT cDNAs were functionally expressed in Saccharomyces cerevisiae strains to probe their possible involvement in the enrichment of TAG with ALA. Among the identified flax LACSs, LuLACS8A exhibited significantly enhanced specificity for ALA over oleic acid (18:1Δ^9cis ) or linoleic acid (18:2Δ^9cis,12cis ). Enhanced α-linolenoyl-CoA specificity was also observed in the enzymatic assay of flax DGAT2 (LuDGAT2-3), which displayed ∼20 times increased preference toward α-linolenoyl-CoA over oleoyl-CoA. Moreover, when LuLACS8A and LuDGAT2-3 were co-expressed in yeast, both in vitro and in vivo experiments indicated that the ALA-containing TAG enrichment process was operative between LuLACS8A- and LuDGAT2-3-catalyzed reactions. Overall, the results support the hypothesis that the cooperation between the reactions catalyzed by LACS8 and DGAT2 may represent a route to enrich ALA production in the flax seed oil.

Cellular and Molecular Responses of Dunaliella tertiolecta by Expression of a Plant Medium Chain Length Fatty Acid Specific Acyl-ACP Thioesterase

Metabolic engineering of microalgae to accumulate high levels of medium chain length fatty acids (MCFAs) has met with limited success. Traditional approaches employ single introduction of MCFA specific acyl-ACP thioesterases (TEs), but our current research in transgenic Dunaliella tertiolecta line has highlighted that, there is no single rate-limiting approach that can effectively increase MCFA levels. Here, we explore the accumulation of MCFAs in D. tertiolecta after transgenic expression of myristic acid biased TE (C14TE). We observe that the MCFA levels were negatively correlated to the fatty acid (FA) synthesis genes, ketoacyl-ACP synthase II (KASII), stearoyl-CoA-9-desaturase (Δ9D), and oleoyl-CoA-12-desaturase (Δ12D). To further examine the molecular mechanism of MCFA accumulation in microalgae, we investigate the transcriptomic dynamics of the MCFA producing strain of D. tertiolecta. At the transcript level, enhanced MCFA accumulation primarily involved up-regulation of photosynthetic genes and down-regulation of genes from central carbon metabolic processes, resulting in an overall decrease in carbon precursors for FA synthesis. We additionally observe that MCFA specific peroxisomal β-oxidation gene (ACX3) was greatly enhanced to prevent excessive build-up of unusual MCFA levels. Besides, long chain acyl-CoA synthetase gene (LACS) was down-regulated, likely in attempt to control fatty acyl supply flux to FA synthesis cycle. This article provides a spatial regulation model of unusual FA accumulation in microalgae and a platform for additional metabolic engineering targeting pathways from FA synthesis, FA transport, and peroxisomal β-oxidation to achieve microalgae oils with higher levels of MCFAs.

Transcriptomic Analyses of Scrippsiella trochoidea Reveals Processes Regulating Encystment and Dormancy in the Life Cycle of a Dinoflagellate, with a Particular Attention to the Role of Abscisic Acid

Due to the vital importance of resting cysts in the biology and ecology of many dinoflagellates, a transcriptomic investigation on Scrippsiella trochoidea was conducted with the aim to reveal the molecular processes and relevant functional genes regulating encystment and dormancy in dinoflagellates. We identified via RNA-seq 3,874 (out of 166,575) differentially expressed genes (DEGs) between resting cysts and vegetative cells; a pause of photosynthesis (confirmed via direct measurement of photosynthetic efficiency); an active catabolism including β-oxidation, glycolysis, glyoxylate pathway, and TCA in resting cysts (tested via measurements of respiration rate); 12 DEGs encoding meiotic recombination proteins and members of MEI2-like family potentially involved in sexual reproduction and encystment; elevated expressions in genes encoding enzymes responding to pathogens (chitin deacetylase) and ROS stress in cysts; and 134 unigenes specifically expressed in cysts. We paid particular attention to genes pertaining to phytohormone signaling and identified 4 key genes regulating abscisic acid (ABA) biosynthesis and catabolism, with further characterization based on their full-length cDNA obtained via RACE-PCR. The qPCR results demonstrated elevated biosynthesis and repressed catabolism of ABA during the courses of encystment and cyst dormancy, which was significantly enhanced by lower temperature (4 ± 1°C) and darkness. Direct measurements of ABA using UHPLC-MS/MS and ELISA in vegetative cells and cysts both fully supported qPCR results. These results collectively suggest a vital role of ABA in regulating encystment and maintenance of dormancy, akin to its function in seed dormancy of higher plants. Our results provided a critical advancement in understanding molecular processes in resting cysts of dinoflagellates.

PLATINUM SENSITIVE 2 LIKE impacts growth, root morphology, seed set, and stress responses

Eukaryotic protein phosphatase 4 (PP4) is a PP2A-type protein phosphatase that is part of a conserved complex with regulatory factors PSY2 and PP4R2. Various lines of Arabidopsis thaliana with mutated PP4 subunit genes were constructed to study the so far completely unknown functions of PP4 in plants. Mutants with knocked out putative functional homolog of the PSY2 LIKE (PSY2L) gene were dwarf and bushy, while plants with knocked out PP4R2 LIKE (PP4R2L) looked very similar to WT. The psy2l seedlings had short roots with disorganized morphology and impaired meristem. Seedling growth was sensitive to the genotoxin cisplatin. Global transcript analysis (RNA-seq) of seedlings and rosette leaves revealed several groups of genes, shared between both types of tissues, strongly influenced by knocked out PSY2L. Receptor kinases, CRINKLY3 and WAG1, important for growth and development, were down-regulated 3–7 times. EUKARYOTIC ELONGATION FACTOR5A1 was down-regulated 4–6 fold. Analysis of hormone sensitive genes indicated that abscisic acid levels were high, while auxin, cytokinin and gibberellic acid levels were low in psy2l. Expression of specific transcription factors involved in regulation of anthocyanin synthesis were strongly elevated, e.g. the master regulator PAP1, and intriguingly TT8, which is otherwise mainly expressed in seeds. The psy2l mutants accumulated anthocyanins under conditions where WT did not, pointing to PSY2L as a possible upstream negative regulator of PAP1 and TT8. Expression of the sugar-phosphate transporter GPT2, important for cellular sugar and phosphate homeostasis, was enhanced 7–8 times. Several DNA damage response genes, including the cell cycle inhibitor gene WEE1, were up-regulated in psy2l. The activation of DNA repair signaling genes, in combination with phenotypic traits showing aberrant root meristem and sensitivity to the genotoxic cisplatin, substantiate the involvement of Arabidopsis PSY2L in maintenance of genome integrity.

Modeling the Metabolism of Arabidopsis thaliana: Application of Network Decomposition and Network Reduction in the Context of Petri Nets

Motivation:Arabidopsis thaliana is a well-established model system for the analysis of the basic physiological and metabolic pathways of plants. Nevertheless, the system is not yet fully understood, although many mechanisms are described, and information for many processes exists. However, the combination and interpretation of the large amount of biological data remain a big challenge, not only because data sets for metabolic paths are still incomplete. Moreover, they are often inconsistent, because they are coming from different experiments of various scales, regarding, for example, accuracy and/or significance. Here, theoretical modeling is powerful to formulate hypotheses for pathways and the dynamics of the metabolism, even if the biological data are incomplete. To develop reliable mathematical models they have to be proven for consistency. This is still a challenging task because many verification techniques fail already for middle-sized models. Consequently, new methods, like decomposition methods or reduction approaches, are developed to circumvent this problem. Methods: We present a new semi-quantitative mathematical model of the metabolism of Arabidopsis thaliana. We used the Petri net formalism to express the complex reaction system in a mathematically unique manner. To verify the model for correctness and consistency we applied concepts of network decomposition and network reduction such as transition invariants, common transition pairs, and invariant transition pairs. Results: We formulated the core metabolism of Arabidopsis thaliana based on recent knowledge from literature, including the Calvin cycle, glycolysis and citric acid cycle, glyoxylate cycle, urea cycle, sucrose synthesis, and the starch metabolism. By applying network decomposition and reduction techniques at steady-state conditions, we suggest a straightforward mathematical modeling process. We demonstrate that potential steady-state pathways exist, which provide the fixed carbon to nearly all parts of the network, especially to the citric acid cycle. There is a close cooperation of important metabolic pathways, e.g., the de novo synthesis of uridine-5-monophosphate, the γ-aminobutyric acid shunt, and the urea cycle. The presented approach extends the established methods for a feasible interpretation of biological network models, in particular of large and complex models.

Tobacco pollen tubes - a fast and easy tool for studying lipid droplet association of plant proteins

In recent years, lipid droplets have emerged as dynamic organelles rather than inactive storage sites for triacylglycerol. The number of proteins known to be associated with lipid droplets has increased, but remains small in comparison with those found with other organelles. Also the mechanisms of how lipid droplets are recognized and bound by proteins need deeper investigation. Here, we present a fast, simple and inexpensive approach to assay proteins for their association with lipid droplets in vivo that can help to screen protein candidates or mutated variants of proteins for their association in an efficient manner. For this, a system to transiently transform Nicotiana tabacum pollen grains was used because these naturally contain lipid droplets. We designed vectors for fast cloning of genes as fusions with either mVenus or mCherry. This allowed us to assay colocalization with lipid droplets stained with Nile Red and Bodipy 505/515, respectively. We successfully tested our system not only for proteins from Arabidopsis thaliana, but also for proteins from the moss Physcomitrella patens and the alga Chlamydomonas reinhardtii. The small size of the vector used allows easy exchange of codons by site-directed mutagenesis. We used this to show that two proline residues in the proline knot of a caleosin are not essential for the binding of lipid droplets. We also demonstrated that peroxisomes are not associated with the lipid droplets in tobacco pollen tubes, which reduces the risk of false interpretation of microscopic data in our system.

High light stress triggers distinct proteomic responses in the marine diatom Thalassiosira pseudonana

Background

Diatoms are able to acclimate to frequent and large light fluctuations in the surface ocean waters. However, the molecular mechanisms underlying these acclimation responses of diaotms remain elusive.

Results

In this study, we investigated the mechanism of high light protection in marine diatom Thalassiosira pseudonana using comparative proteomics in combination with biochemical analyses. Cells treated under high light (800 μmol photons m⁻²s⁻¹) for 10 h were subjected to proteomic analysis. We observed that 143 proteins were differentially expressed under high light treatment. Light-harvesting complex proteins, ROS scavenging systems, photorespiration, lipid metabolism and some specific proteins might be involved in light protection and acclimation of diatoms. Non-photochemical quenching (NPQ) and relative electron transport rate could respond rapidly to varying light intensities. High-light treatment also resulted in increased diadinoxanthin + diatoxanthin content, decreased Fv/Fm, increased triacylglycerol and altered fatty acid composition. Under HL stress, levels of C14:0 and C16:0 increased while C20:5ω3 decreased.

Conclusions

We demonstrate that T. pseudonana has efficient photoprotective mechanisms to deal with HL stress. De novo synthesis of Ddx/Dtx and lipid accumulation contribute to utilization of the excess energy. Our data will provide new clues for in-depth study of photoprotective mechanisms in diatoms.

Electronic supplementary material

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

Oil Biosynthesis in Underground Oil-Rich Storage Vegetative Tissue: Comparison of Cyperus esculentus Tuber with Oil Seeds and Fruits

Cyperus esculentus is unique in that it can accumulate rich oil in its tubers. However, the underlying mechanism of tuber oil biosynthesis is still unclear. Our transcriptional analyses of the pathways from pyruvate production up to triacylglycerol (TAG) accumulation in tubers revealed many distinct species-specific lipid expression patterns from oil seeds and fruits, indicating that in C. esculentus tuber: (i) carbon flux from sucrose toward plastid pyruvate could be produced mostly through the cytosolic glycolytic pathway; (ii) acetyl-CoA synthetase might be an important contributor to acetyl-CoA formation for plastid fatty acid biosynthesis; (iii) the expression pattern for stearoyl-ACP desaturase was associated with high oleic acid composition; (iv) it was most likely that endoplasmic reticulum (ER)-associated acyl-CoA synthetase played a significant role in the export of fatty acids between the plastid and ER; (v) lipid phosphate phosphatase (LPP)-δ was most probably related to the formation of the diacylglycerol (DAG) pool in the Kennedy pathway; and (vi) diacylglyceroltransacylase 2 (DGAT2) and phospholipid:diacylglycerolacyltransferase 1 (PDAT1) might play crucial roles in tuber oil biosynthesis. In contrast to oil-rich fruits, there existed many oleosins, caleosins and steroleosins with very high transcripts in tubers. Surprisingly, only a single ortholog of WRINKLED1 (WRI1)-like transcription factor was identified and it was poorly expressed during tuber development. Our study not only provides insights into lipid metabolism in tuber tissues, but also broadens our understanding of TAG synthesis in oil plants. Such knowledge is of significance in exploiting this oil-rich species and manipulating other non-seed tissues to enhance storage oil production.

The Roles of β-Oxidation and Cofactor Homeostasis in Peroxisome Distribution and Function in Arabidopsis thaliana

Key steps of essential metabolic pathways are housed in plant peroxisomes. We conducted a microscopy-based screen for anomalous distribution of peroxisomally targeted fluorescence in Arabidopsis thaliana This screen uncovered 34 novel alleles in 15 genes affecting oil body mobilization, fatty acid β-oxidation, the glyoxylate cycle, peroxisome fission, and pexophagy. Partial loss-of-function of lipid-mobilization enzymes conferred peroxisomes clustered around retained oil bodies without other notable defects, suggesting that this microscopy-based approach was sensitive to minor perturbations, and that fatty acid β-oxidation rates in wild type are higher than required for normal growth. We recovered three mutants defective in PECTIN METHYLESTERASE31, revealing an unanticipated role in lipid mobilization for this cytosolic enzyme. Whereas mutations reducing fatty acid import had peroxisomes of wild-type size, mutations impairing fatty acid β-oxidation displayed enlarged peroxisomes, possibly caused by excess fatty acid β-oxidation intermediates in the peroxisome. Several fatty acid β-oxidation mutants also displayed defects in peroxisomal matrix protein import. Impairing fatty acid import reduced the large size of peroxisomes in a mutant defective in the PEROXISOMAL NAD⁺ TRANSPORTER (PXN), supporting the hypothesis that fatty acid accumulation causes pxn peroxisome enlargement. The diverse mutants isolated in this screen will aid future investigations of the roles of β-oxidation and peroxisomal cofactor homeostasis in plant development.