BACKGROUND:Soybean is an economically important crop which flowers predominantly in response to photoperiod. Several major loci controlling the quantitative trait for reproductive timing have been identified, of which allelic combinations at three of these loci, E1, E2, and E3, are the dominant factors driving time to flower and reproductive period. However, functional genomics studies have identified additional loci which affect reproductive timing, many of which are less understood. A better characterization of these genes will enable fine-tuning of adaptation to various production environments. Two such genes, E1La and E1Lb, have been implicated in flowering by previous studies, but their effects have yet to be assessed under natural photoperiod regimes.RESULTS:Natural and induced variants of E1La and E1Lb were identified and introgressed into lines harboring either E1 or its early flowering variant, e1-as. Lines were evaluated for days to flower and maturity in a Maturity Group (MG) III production environment. These results revealed that variation in E1La and E1Lb promoted earlier flowering and maturity, with stronger effects in e1-as background than in an E1 background. The geographic distribution of E1La alleles among wild and cultivated soybean revealed that natural variation in E1La likely contributed to northern expansion of wild soybean, while breeding programs in North America exploited e1-as to develop cultivars adapted to northern latitudes.CONCLUSION:This research identified novel alleles of the E1 paralogues, E1La and E1Lb, which promote flowering and maturity under natural photoperiods. These loci represent sources of genetic variation which have been under-utilized in North American breeding programs to control reproductive timing, and which can be valuable additions to a breeder's molecular toolbox.

Article Figures and data Abstract eLife digest Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity. eLife digest Plants produce a vast variety of different molecules known as secondary or specialized metabolites to attract pollinating insects, such as bees, or protect themselves against herbivores and pests. The secondary metabolites are made from simple building blocks that are readily available in plants, including amino acids, fatty acids and sugars. Different species of plant, and even different parts of the same plant, produce their own sets of secondary metabolites. For example, the hairs on the surface of tomatoes and other members of the nightshade family of plants make metabolites known as acylsugars. These chemicals deter herbivores and pests from damaging the plants. To make acylsugars, the plants attach long chains known as fatty acyl groups to molecules of sugar, such as sucrose. Some members of the nightshade family produce acylsugars with longer chains than others. In particular, acylsugars with long chains are only found in tomatoes and other closely-related species. It remained unclear how the nightshade family evolved to produce acylsugars with chains of different lengths. To address this question, Fan et al. used genetic and biochemical approaches to study tomato plants and other members of the nightshade family. The experiments identified two genes known as AACS and AECH in tomatoes that produce acylsugars with long chains. These two genes originated from the genes of older enzymes that metabolize fatty acids – the building blocks of fats – in plant cells. Unlike the older genes, AACS and AECH were only active at the tips of the hairs on the plant’s surface. Fan et al. then investigated the evolutionary relationship between 11 members of the nightshade family and two other plant species. This revealed that AACS and AECH emerged in the nightshade family around the same time that longer chains of acylsugars started appearing. These findings provide insights into how plants evolved to be able to produce a variety of secondary metabolites that may protect them from a broader range of pests. The gene cluster identified in this work could be used to engineer other species of crop plants to start producing acylsugars as natural pesticides. Introduction Despite the enormous structural diversity of plant specialized metabolites, they are derived from a relatively small number of primary metabolites, such as sugars, amino acids, nucleotides, and fatty acids (Maeda, 2019). These lineage-, tissue- or cell- type specific specialized metabolites mediate environmental interactions, such as herbivore and pathogen deterrence or pollinator and symbiont attraction (Mithöfer and Boland, 2012; Pichersky and Lewinsohn, 2011). Specialized metabolism evolution is primarily driven by gene duplication (Moghe and Last, 2015; Panchy et al., 2016), and relaxed selection of the resulting gene pairs allows modification of cell- and tissue-specific gene expression and changes in enzymatic activity. This results in expanded substrate recognition and/or diversified product formation (Khersonsky and Tawfik, 2010; Leong and Last, 2017). The neofunctionalized enzymes can prime the origin and diversification of specialized metabolic pathways (Schenck and Last, 2020; Weng et al., 2012; Weng, 2014). There are many examples of mechanisms that lead to novel enzymatic activities in specialized cell- or tissue-types, however, the principles that govern assembly of multi-enzyme specialized metabolic pathways are less well established. One appealing hypothesis involves the stepwise recruitment of pathway enzymes (Noda-Garcia et al., 2018). In rare cases, non-homologous specialized metabolic enzyme genes occur in proximity to each other in a genomic region, forming a biosynthetic gene cluster (Nützmann et al., 2016; Nützmann and Osbourn, 2014; Rokas et al., 2018). In recent years, an increasing number of specialized metabolic gene clusters (SMGCs) were experimentally identified or bioinformatically predicted in plants (Boutanaev et al., 2015; Castillo et al., 2013; Schläpfer et al., 2017). However, although most experimentally characterized plant SMGCs are co-expressed, the majority of the bioinformatically predicted ones do not show coexpression under global network analysis (Wisecaver et al., 2017). While examples of SMGCs are still relatively rare in plants, experimentally validated cases were reported for a surprisingly diverse group of pathways. These include terpenes (Chae et al., 2014; Prisic et al., 2004; Qi et al., 2004; Wilderman et al., 2004), cyclic hydroxamic acids (Frey et al., 1997), biosynthetically unrelated alkaloids (Itkin et al., 2013; Winzer et al., 2012), polyketides (Schneider et al., 2016), cyanogenic glucosides (Takos et al., 2011), and modified fatty acids (Jeon et al., 2020). However, whereas each cluster encodes multiple non-homologous enzymes of a biosynthetic pathway, evolution of their assembly is not well understood. Acylsugars are a group of insecticidal (Leckie et al., 2016) and anti-inflammatory (Herrera-Salgado et al., 2005) chemicals mainly observed in glandular trichomes of Solanaceae species (Fan et al., 2019; Schuurink and Tissier, 2020). These specialized metabolites are sugar aliphatic esters with three levels of structural diversity across the Solanaceae family: acyl chain length, acylation position, and sugar core (Fan et al., 2019). The primary metabolites sucrose and aliphatic acyl-CoAs are the biosynthetic precursors of acylsucroses in plants as evolutionarily divergent as the cultivated tomato Solanum lycopersicum (Fan et al., 2016; Figure 1), Petunia axillaris (Nadakuduti et al., 2017) and Salpiglossis sinuata (Moghe et al., 2017). The core tomato acylsucrose biosynthetic pathway involves four BAHD [BEAT, AHCT, HCBT, DAT (D'Auria, 2006) family acylsucrose acyltransferases (Sl-ASAT1 through Sl-ASAT4), which are specifically expressed in the type I/IV trichome tip cells (Fan et al., 2016; Schilmiller et al., 2015; Schilmiller et al., 2012). These enzymes catalyze consecutive reactions utilizing sucrose and acyl-CoA substrates to produce the full set of cultivated tomato acylsucroses in vitro (Fan et al., 2016). Figure 1 Download asset Open asset Primary metabolites are biosynthetic precursors of tomato trichome acylsugars. In cultivated tomatoes, the trichome acylsucroses are synthesized by four Sl-ASATs using the primary metabolites – sucrose and different types of acyl-CoAs – as substrates. In this study we provide evidence that medium chain fatty acids are converted to acyl-CoAs by an acyl-CoA synthetase for medium chain acylsugar biosynthesis. Co-option of primary metabolic enzymes contributed to the evolution of acylsugar biosynthesis and led to interspecific structural diversification across the Solanum tomato clade. One example is an invertase-like enzyme originating from carbohydrate metabolism that generates acylglucoses in the wild tomato S. pennellii through cleavage of the acylsucrose glycosidic bond (Leong et al., 2019). In another case, allelic variation of a truncated isopropylmalate synthase-like enzyme (IPMS3) – from branched chain amino acid metabolism – leads to acylsugar iC4/iC5 (2-methylpropanoic/3-methylbutanoic acid) acyl chain diversity in S. pennellii and S. lycopersicum (Ning et al., 2015). Acylsugar structural diversity is even more striking across the family. Previous studies revealed variation in acyl chain length (Ghosh et al., 2014; Liu et al., 2017; Moghe et al., 2017): Nicotiana, Petunia and Salpiglossis species were reported to accumulate acylsugars containing only short acyl chains (carbon number, C ≤ 8). In contrast, some species in Solanum and other closely related genera produce acylsugars with medium acyl chains (C ≥ 10). These results are consistent with the hypothesis that the capability to produce medium chain acylsugars varies across the Solanaceae family. In this study, we identify a metabolic gene cluster on tomato chromosome 7 containing two non-homologous genes – acylsugar acyl-CoA synthetase (AACS) and acylsugar enoyl-CoA hydratase (AECH) – affecting medium chain acylsugar biosynthesis. Genetic and biochemical results show that the trichome enriched AACS and AECH are involved in generating medium chain acyl-CoAs, which are donor substrates for acylsugar biosynthesis. Genomic analysis revealed a syntenic region on chromosome 12, where the acylsucrose biosynthetic Sl-ASAT1 is located (Fan et al., 2016). Phylogenetic analysis of the syntenic regions in Solanaceae and beyond led to evolutionary reconstruction of the origin of the acylsugar gene cluster. We infer that sequential gene insertion facilitated emergence of this gene cluster in tomato. These results provide insights into specialized metabolic evolution through emergence of cell-type specific gene expression, the formation of metabolic gene clusters and illuminates additional examples of primary metabolic enzymes being co-opted into specialized metabolism. Results Identification of a metabolic gene cluster that affects tomato trichome medium chain acylsugar biosynthesis S. pennellii natural accessions (Mandal et al., 2020), as well as the S. lycopersicum M82 ×S. pennellii LA0716 chromosomal substitution introgression lines (ILs) (Eshed and Zamir, 1995), offer convenient resources to investigate interspecific genetic variation that affects acylsugar metabolic diversity (Mandal et al., 2020; Schilmiller et al., 2010). In a rescreen of ILs for S. pennellii genetic regions that alter trichome acylsugar profiles (Schilmiller et al., 2010), IL7-4 was found to accumulate increased C10 medium chain containing acylsugars compared with M82 (Figure 2, A and B). The genetic locus that contributes to the acylsugar phenotype was narrowed down to a 685 kb region through screening selected backcross inbred lines (BILs) (Ofner et al., 2016) that have recombination breakpoints on chromosome 7 (Figure 2C). Because tomato acylsucrose biosynthesis occurs in trichomes, candidate genes in this region were filtered based on their trichome-specific expression patterns. This analysis identified a locus containing multiple tandemly duplicated genes of three families – an acyl-CoA synthetase (ACS), enoyl-CoA hydratase (ECH), and BAHD acyltransferase. Our analysis (Moore et al., 2020) revealed co-expression of four Sl-ASATs (Fan et al., 2019) and three genes at the locus – Solyc07g043630, Solyc07g043660, and Solyc07g043680 (Supplementary file 1 and Figure 2—figure supplement 1). Expression of these three genes was trichome enriched (Figure 2D), and thus they were selected for further analysis. Figure 2 with 1 supplement see all Download asset Open asset Mapping of a genetic locus related to acylsugar variations in tomato interspecific introgression lines. (A) Electrospray ionization negative (ESI-) mode, base-peak intensity (BPI) LC/MS chromatogram of trichome metabolites from cultivated tomato S. lycopersicum M82 and introgression line IL7-4. The orange bars highlight two acylsugars that have higher abundance in IL7-4 than in M82. For the acylsucrose nomenclature, ‘S’ refers to a sucrose backbone, ‘3:22’ means three acyl chains with twenty-two carbons in total. The length of each acyl chain is shown in the parentheses. (B) Peak area percentage of seven major trichome acylsugars in M82 and IL7-4. The sum of the peak area percentage of each acylsugar is equal to 100% in each sample. The data is shown for three plants ± SEM. **p<0.01, Welch two-sample t test. Figure 2—source data 1 includes values for the analysis. (C) Mapping the genetic locus contributing to the IL7-4 acylsugar phenotype using selected backcross inbred lines (BILs) that have recombination break points within the introgression region of IL7-4. (D) Narrowing down candidate genes in the locus using trichome/stem RNA-seq datasets generated from previous study (Ning et al., 2015). A region with duplicated genes of three types – acyl-CoA synthetase (ACS), BAHD acyltransferase, and enoyl-CoA hydratase (ECH) – is shown. The red-blue color gradient provides a visual marker to rank the expression levels represented by Fragments Per Kilobase of transcript per Million mapped reads (FPKM). Coexpression analysis of tomato ACS, ECH, and BAHD acyltransferase family genes is shown in Figure 2—figure supplement 1. Figure 2—source data 1 Data used to make Figure 2B. Peak area percentage of seven major trichome acylsugars in M82 and IL7-4. https://cdn.elifesciences.org/articles/56717/elife-56717-fig2-data1-v2.xlsx Download elife-56717-fig2-data1-v2.xlsx The three candidate genes were tested for involvement in tomato acylsugar biosynthesis by making loss of function mutations using the CRISPR-Cas9 gene editing system. Two guide RNAs (gRNAs) were designed to target one or two exons of each gene to assist site-specific DNA cleavage by hCas (Brooks et al., 2014; Figure 3—figure supplement 1, A–C). In the self-crossed T1 progeny of stably transformed M82 plants, at least two homozygous mutants were obtained in Solyc07g043630, Solyc07g043660, and Solyc07g043680 (Figure 3—figure supplement 1, A–C), and these were analyzed for leaf trichome acylsugar changes. Altered acylsugar profiles were observed in the ACS-annotated Solyc07g043630 or ECH-annotated Solyc07g043680 mutants (Figure 3, A and B), but not in the ACS-annotated Solyc07g043660 mutant (Figure 3—figure supplement 1D). Despite carrying mutations in distinctly annotated genes (ACS or ECH), the two mutants exhibited the same phenotype – no detectable medium acyl chain (C10 or C12) containing acylsugars (Figure 3, A and B). We renamed Solyc07g043630 as acylsugar acyl-CoA synthetase 1 (Sl-AACS1) and Solyc07g043680 as acylsugar enoyl-CoA hydratase 1 (Sl-AECH1) based on this analysis. Figure 3 with 2 supplements see all Download asset Open asset CRISPR/Cas9-mediated gene knockout of tomato Sl-AACS1 or Sl-AECH1 eliminates detectable medium chain containing acylsugars. (A) Combined LC/MS extracted ion chromatograms of trichome metabolites from CRISPR mutants sl-aacs1 and sl-aech1. The medium chain acylsugars that are not detected in the two mutants are denoted by pairs of vertical dotted lines. Figure 3—figure supplement 1 describes the design of the gRNAs and details of the gene edits. (B) Quantification of seven major trichome acylsugars in sl-aacs1 and sl-aech1 mutants. Two independent T2 generation transgenic lines for each mutant were used for analysis. The peak area/internal standard (IS) normalized by leaf dry weight (DW) is shown from six plants ± SEM. Figure 3—source data 1 includes values for the analysis. (C) Confocal fluorescence images showing that GFP fluorescence driven by Sl-AACS1 or Sl-AECH1 is located in the tip cells of type I/IV trichomes. Their tissue specific expressions are similar to Sl-ASAT1 (Fan et al., 2016), which locates in a chromosome 12 region that is syntenic to the locus containing Sl-AACS1 and Sl-AECH1. Figure 3—figure supplement 2 provides the detailed information of the syntenic region. Sl-AACS1, Sl-AECH1, and Sl-ASAT1 are the only gene models with demonstrated functions in acylsugar biosynthesis. Figure 3—source data 1 Data used to make Figure 3B. Quantification of seven major trichome acylsugars in the CRISPR mutants sl-aacs1 and sl-aech1, as well as the parent M82. https://cdn.elifesciences.org/articles/56717/elife-56717-fig3-data1-v2.xlsx Download elife-56717-fig3-data1-v2.xlsx Further genomic analysis revealed that Sl-AACS1 and Sl-AECH1 belong to a syntenic region shared with a locus on chromosome 12, where Sl-ASAT1 is located (Figure 3C and Figure 3—figure supplement 2). Sl-ASAT1 is specifically expressed in trichome tip cells and encodes the enzyme catalyzing the first step of tomato acylsucrose biosynthesis (Fan et al., 2016). This led us to test the cell-type expression pattern of Sl-AACS1 and Sl-AECH1. Like Sl-ASAT1, the promoters of both genes drove GFP expression in the trichome tip cells of stably transformed M82 plants (Figure 3C). This supports our hypothesis that Sl-AACS1 and Sl-AECH1 are involved in tomato trichome acylsugar biosynthesis. Taken together, we identified a metabolic gene cluster involved in medium chain acylsugar biosynthesis, which is composed of two cell-type specific genes. In vitro analysis of Sl-AACS1 and Sl-AECH1 implicates their roles in medium chain acyl-CoA metabolism ACS and ECH are established to function in multiple cell compartments for the metabolism of acyl-CoA (Buchanan et al., 2015), the acyl donor substrates for ASAT enzymes. We sought to understand the organelle targeting of Sl-AACS1 and Sl-AECH1, to advance our knowledge of acylsugar machinery at the subcellular level. We constructed expression cassettes of Sl-AACS1, Sl-AECH1 and Solyc07g043660 with C-terminal cyan fluorescent protein (CFP), hypothesizing that the targeting peptides reside at the N-terminus of precursor proteins. When co-expressed in tobacco leaf epidermal cells, three CFP-tagged recombinant proteins co-localized with the mitochondrial marker MT-RFP (Nelson et al., 2007; Figure 4A and Figure 4—figure supplement 1A). To rule out the possibility of peroxisomal localization, we fused Sl-AACS1, Sl-AECH1, or Solyc07g043660 with N-terminus fused yellow fluorescent protein (YFP), considering that potential peroxisomal targeting peptides are usually located on the C-terminus (Brocard and Hartig, 2006). The expressed YFP-recombinant proteins were not co-localized with the peroxisomal marker RFP-PTS (Nelson et al., 2007; Figure 4—figure supplement 1B). Instead, they appeared distributed in the cytosol (Figure 4—figure supplement 1B), presumably because the N-terminal YFP blocked the mitochondria targeting signal. Taken together, protein expression and co-localization analyses suggest that Sl-AACS1, Sl-AECH1, and Solyc07g043660 encode enzymes targeted to mitochondria. Figure 4 with 2 supplements see all Download asset Open asset Functional analysis of Sl-AACS1 and Sl-AECH1 in N. benthamiana and recombinant Sl-AACS1 enzyme analysis. (A) Confocal images of co-expression analysis in tobacco leaf epidermal cells using C-terminal CFP-tagged either Sl-AACS1 or Sl-AECH1 and the mitochondrial marker MT-RFP. Arrowheads point to mitochondria that are indicated by MT-RFP fluorescent signals. Scale bar equals 10 μm. Figure 4—figure supplement 1B describes that the expressed YFP-recombinant proteins were not co-localized with the peroxisomal marker RFP-PTS (B) Aliphatic fatty acids of different chain lengths were used as the substrates to test Sl-AACS1 acyl-CoA synthetase activity. Mean amount of acyl-CoAs generated (nmol min−1 mg−1 proteins) was used to represent enzyme activities. The results are from three measurements ± SEM. Figure 4—source data 1 includes values for the measurements. (C) Enzyme activity of Sl-AACS1 for six fatty acid substrates. (D) Identification of membrane lipid phosphatidylcholine (PC), which contains medium acyl chains, following transient expression of Sl-AECH1 in N. benthamiana leaves. The results from expressing Sl-AACS1 and co-expressing both Sl-AECH1 and Sl-AACS1 are also shown. Mole percentage (Mol %) of the acyl chains from membrane lipids with carbon number 12, 14, 16, and 18 are shown for three biological replicates ± SEM. *p<0.05, **p<0.01. Welch two-sample t test was performed comparing with the empty vector control. Figure 4—source data 2 includes values for the lipid analysis. Acyl groups of the same chain lengths with saturated and unsaturated bonds were combined in the calculation. Figure 4—figure supplement 2 shows that the putative Sl-AECH1 orthologs from S. pennellii and S. quitoense generated medium chain lipids in the infiltrated leaves. Figure 4—source data 1 Data used to make Figure 4B. Aliphatic fatty acids of different chain lengths were used as the substrates to test Sl-AACS1 acyl-CoA synthetase activity. https://cdn.elifesciences.org/articles/56717/elife-56717-fig4-data1-v2.xlsx Download elife-56717-fig4-data1-v2.xlsx Figure 4—source data 2 Data used to make Figure 4D and Figure 4—figure supplement 1C. N. benthamiana leaf membrane lipid acyl chain composition. https://cdn.elifesciences.org/articles/56717/elife-56717-fig4-data2-v2.xlsx Download elife-56717-fig4-data2-v2.xlsx Sl-AACS1 belongs to a group of enzymes that activate diverse carboxylic acid substrates to produce acyl-CoAs. We hypothesized that Sl-AACS1 uses medium chain fatty acids as substrates, because ablation of Sl-AACS1 eliminated acylsugars with medium acyl chains. To characterize the in vitro activity of Sl-AACS1, we purified recombinant His-tagged proteins from Escherichia coli. Enzyme assays were performed by supplying fatty acid substrates with even carbon numbers from C2 through C18 (Figure 4B). The results showed that Sl-AACS1 utilized fatty acid substrates with lengths ranging from C6 to C12, including those with a terminal branched carbon (iC10:0) or an unsaturated bond (trans-2-decenoic acid, C10:1) (Figure 4, B and C). However, no activity was observed with the 3-hydroxylated C12 and C14 fatty acids as substrates (Figure 4B). These results support our hypothesis that Sl-AACS1 produces medium chain acyl-CoAs, which are in vivo substrates for acylsugar biosynthesis. To test whether Sl-AACS1 and Sl-AECH1 can produce medium chain acyl-CoAs in planta, we transiently expressed these genes in Nicotiana benthamiana leaves using Agrobacterium-mediated infiltration (Sainsbury et al., 2009). It is challenging to directly measure plant acyl-CoAs, due to their low concentration and separate organellar pools. We used an alternative approach and characterized membrane lipids, which are produced from acyl-CoA intermediates. We took advantage of the observation that N. benthamiana membrane lipids do not accumulate detectable acyl chains of 12 carbons or shorter. N. benthamiana leaves were infiltrated with constructs containing Sl-AACS1 or Sl-AECH1 individually, or together (Figure 4D). In contrast to the empty vector control, infiltration of Sl-AECH1 led to detectable levels of C12 acyl chains in the leaf membrane lipid phosphatidylcholine (PC) (Figure 4D). We also observed increased C14 acyl chains in PC, phosphatidylglycerol (PG), sulfoquinovosyl diacylglycerol (SQDG), and digalactosyldiacylglycerol (DGDG) in Sl-AECH1 infiltrated plants (Figure 4D and Figure 4—figure supplement 1C). These results suggest that Sl-AECH1 participates in generation of medium chain acyl-CoAs in planta, which are channeled into lipid biosynthesis. No medium chain acylsugars were detected, presumably due to the lack of core acylsugar biosynthetic machinery in N. benthamiana mesophyll cells. We asked whether the closest known homologs of Sl-AECH1 from Solanum species can generate medium chain lipids when transiently expressed in N. benthamiana. Two SQDGs with C12 chains were monitored by LC/MS as peaks diagnostic of lipids containing medium chain fatty acids (Figure 4—figure supplement 2, A and B). The results showed that only the putative Sl-AECH1 orthologs Sopen07g023250 (Sp-AECH1) and Sq_c37194 (Sq-AECH1) – from S. pennellii and S. quitoense respectively – generated medium chain lipids in the infiltrated leaves (Figure 4—figure supplement 2C). This confirms that not all ECHs can produce medium chain lipids and suggests that the function of Sl-AECH1 evolved recently, presumably as a result of neofunctionalization after gene duplication (Figure 4—figure supplement 2C). AACS1 and AECH1 are evolutionarily conserved in the Solanum Medium chain acylsugars were documented in Solanum species besides cultivated tomato, including S. pennellii (Leong et al., 2019), S. nigrum (Moghe et al., 2017), as well as the more distantly related S. quitoense (Leong et al., 2020; Hurney, 2018) and S. lanceolatum (Herrera-Salgado et al., 2005). We hypothesized that evolution of AACS1 and AECH1 contributed to medium chain acylsugar biosynthesis in Solanum. As a test, we analyzed the genomes of Solanum species other than cultivated tomato. Indeed, the acylsugar related synteny containing ACS and ECH was found in both S. pennellii and S. melongena (eggplant), suggesting that the cluster assembly evolved before divergence of the tomato and eggplant lineage (Figure 5A). Figure 5 with 2 supplements see all Download asset Open asset AACS1 and AECH1 are evolutionarily conserved in Solanum plants. (A) A conserved syntenic genomic region containing AACS1 and AECH1 was found in three selected Solanum species. Nodes representing estimated dates since the last common ancestors (Särkinen et al., 2013) shown on the left. The closest homologs of AACS1 and AECH1 in Solanum quitoense are shown without genomic context because the genes were identified from RNA-seq and genome sequences are not available. The lines connect genes representing putative orthologs across the four species. The trichome/stem RNA-seq data of two biological S. pennellii replicates are summarized (Supplementary file 2) for genes in the syntenic region. The red-blue color gradient provides a visual marker to rank the expression levels in FPKM. Structures of representative medium chain acylsugars from S. quitoense (acylinositol, I4:26) (Hurney, 2018) and S. pennellii (acylglucose, G3:19) (Leong et al., 2019) are on the right. Figure 5—figure supplement 1 shows that stable Sp-AACS1 transformation of the M82 CRISPR mutant sl-aacs1 restores C12 containing acylsugars (B) CRISPR/Cas9-mediated gene knockout of Sp-AACS1 or Sp-AECH1 in S. pennellii produce no detectable medium chain containing acylsugars. The ESI+ mode LC/MS extracted ion chromatograms of trichome metabolites are shown for each mutant. The m/z 127.01 (left panel) corresponds to the glucopyranose ring fragment that both acylsucroses and acylglucoses generate under high collision energy positive-ion mode. The m/z 155.14 (center panel) and 183.17 (right panel) correspond to the acylium ions from acylsugars with chain length of C10 and C12, respectively. Figure 5—figure supplement 2A–C describes the design of the gRNAs and the detailed information of gene edits. (C) Silencing Sq-AACS1 (Sq-c34025) or Sq-AECH1 (Sq-c37194) in S. quitoense using VIGS leads to reduction of total acylsugars. The peak area/internal standard (IS) normalized by leaf dry weight was shown from sixteen plants ± SEM. ***p<0.001, Welch two-sample t test. Figure 5—figure supplement 2E and F describes the VIGS experimental design and the representative LC/MS extracted ion chromatograms of S. quitoense major acylsugars. (D) Reduced gene expression of Sq-AACS1 or Sq-AECH1 correlates with decreased acylsugar levels in S. quitoense. The qRT-PCR gene expression data are plotted with acylsugar levels of the same leaf as described in Figure 5—figure supplement 2E. Figure 5—source data 1 includes raw data for the S. quitoense VIGS experiments. Figure 5—source data 1 Data used to make Figure 5C and D. Silencing Sq-AACS1 or Sq-AECH1 in S. quitoense using VIGS leads to reduction of total acylsugars. https://cdn.elifesciences.org/articles/56717/elife-56717-fig5-data1-v2.xlsx Download elife-56717-fig5-data1-v2.xlsx We applied gene expression and genetic approaches to test the in vivo functions of ACS and ECH in selected Solanum species. To explore the expression pattern of S. pennellii ACS and ECH cluster genes, we performed RNA-seq analysis on trichomes and shaved stems to identify acylsugar biosynthetic candidates (Supplementary file 2). The expression pattern of S. pennellii cluster genes is strikingly similar to S. lycopersicum: one ECH and two ACS genes are highly enriched in trichomes, including the orthologs of Sl-AACS1 and Sl-AECH1. Sp-AACS1 function (Sopen07g023200) was first tested by asking whether it can reverse the cultivated tomato sl-aacs1 mutant acylsugar phenotype. Indeed, Sp-AACS1 restored C12 containing acylsugars in the stably transformed sl-aacs1 plants (Figure 5—figure supplement 1). To directly test Sp-AACS1 and Sp-AECH1 function, we used CRISPR-Cas9 to make single mutants in S. pennellii LA0716. No medium chain acylsugars were detected in T0 generation mutants with edits for each gene (Figure 5B and Figure 5—figure supplement 2, A and C). Similar to the ACS-annotated Solyc07g043660 cultivated tomato mutant (Figure 3—figure supplement 1D), deletion of S. pennellii ortholog Sopen07g023220 has no observed effects on S. pennellii trichome acylsugars (Figure 5—figure supplement 2D). The medium chain acylsugar producer S. quitoense (Hurney, 2018) was used for AACS1 and AECH1 fun

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Soybean plays an important role in seed oil production for foods and industrial products in the USA. Chemical hydrogenation of commodity soybean oil increased functionality but unavoidably created trans fat which has been linked to many health issues in humans. An alternative to using hydrogenation of the oil to enhance oxidative stability is to genetically increase the level of oleic acid in the seed oil. Our goal was to develop a soybean germplasm as a source for more functional soybean oil with a high oleic acid, increased stearic acid, and low linolenic acid profile utilizing new and existing variant alleles of key fatty acid desaturase genes. Using forward genetics, novel alleles of FAD2-1A were discovered and characterized from a mutant soybean population already containing elevated seed stearic acid content. One of the two new FAD2-1A alleles, when combined with existing mutant alleles of FAD2-1B created a high oleic acid seed oil phenotype in soybean germplasm lines. The other new FAD2-1A allele produced lower and more variable oleic acid levels when combined with existing mutant alleles of FAD2-1B. Seed stearic acid increased to ~ 10–11% in lines containing combinations of FAD2-1A and FAD2-1B mutant alleles plus the SACPD-C missense mutant alleles; however, the increase in stearic acid was associated with a decrease in the oleic acid content and did not meet the target of at least 20% stearic acid in the seed oil. In addition, existing FAD3 mutant alleles were incorporated into the novel high oleic acid and high oleic/elevated stearic acid soybean lines with four or five targeted alleles combined (FAD2-1A, FAD2-1B, FAD3A, FAD3C, and SACPD-C) to create soybean germplasm with more functional soybean oil. This study is beneficial for improving the quality of soybean oil based on nutritional value and oxidative stability.