Skip to main content

Main menu

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae

User menu

  • My alerts
  • Log in

Search

  • Advanced search
Plant Cell
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae
  • My alerts
  • Log in
Plant Cell

Advanced Search

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Follow PlantCell on Twitter
  • Visit PlantCell on Facebook
  • Visit Plantae
Research ArticleLARGE-SCALE BIOLOGY ARTICLES
Open Access

Candidate Gene Networks for Acylsugar Metabolism and Plant Defense in Wild Tomato Solanum pennellii

Sabyasachi Mandal, Wangming Ji, Thomas D. McKnight
Sabyasachi Mandal
Department of Biology, Texas A&M University, College Station, Texas 77843
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sabyasachi Mandal
Wangming Ji
Department of Biology, Texas A&M University, College Station, Texas 77843
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Wangming Ji
Thomas D. McKnight
Department of Biology, Texas A&M University, College Station, Texas 77843
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Thomas D. McKnight
  • For correspondence: mcknight@bio.tamu.edu

Published January 2020. DOI: https://doi.org/10.1105/tpc.19.00552

  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading
  • © 2020 American Society of Plant Biologists. All rights reserved.

Abstract

Many solanaceous plants secrete acylsugars, which are branched-chain and straight-chain fatty acids esterified to Glu or Suc. These compounds have important roles in plant defense and potential commercial applications. However, several acylsugar metabolic genes remain unidentified, and little is known about regulation of this pathway. Comparative transcriptomics between low- and high-acylsugar-producing accessions of Solanum pennellii revealed that expression levels of known and novel candidate genes (putatively encoding beta-ketoacyl-(acyl-carrier-protein) synthases, peroxisomal acyl-activating enzymes, ATP binding cassette (ABC) transporters, and central carbon metabolic proteins) were positively correlated with acylsugar accumulation, except two genes previously reported to be involved in acylglucose biosynthesis. Genes putatively encoding oxylipin metabolic proteins, subtilisin-like proteases, and other antimicrobial defense proteins were upregulated in low-acylsugar–producing accessions. Transcriptome analysis after biochemical inhibition of biosynthesis of branched-chain amino acids (precursors to branched-chain fatty acids) by imazapyr showed concentration-dependent downregulation of known and most acylsugar candidate genes, but not defense genes. Weighted gene correlation network analysis identified separate coexpressed gene networks for acylsugar metabolism (including six transcription factor genes and flavonoid metabolic genes) and plant defense (including genes putatively encoding NB-ARC and leucine-rich repeat sequences, protein kinases and defense signaling proteins, and previously mentioned defense proteins). Additionally, virus-induced gene silencing of two trichomes preferentially expressed candidate genes for straight-chain fatty acid biosynthesis confirmed their role in acylsugar metabolism.

INTRODUCTION

Plants synthesize a wide variety of specialized metabolites that provide selective advantages in specific environmental conditions. Different classes of specialized metabolites are found in specific taxonomic groups, and many of them are valuable phytochemicals. Acylsugars are nonvolatile and viscous specialized metabolites secreted through trichomes of solanaceous plants (Fobes et al., 1985; Kroumova et al., 2016; Moghe et al., 2017), and their roles in plant defense have been studied extensively. For example, acylsugars from Solanum pennellii act as feeding and oviposition deterrents for insect pests and exert toxic effects on different herbivores (Hawthorne et al., 1992; Juvik et al., 1994; Liedl et al., 1995). Acylsugars from other genera also have important roles in providing protection against herbivores and plant pathogens (Chortyk et al., 1997; Hare, 2005; Luu et al., 2017). Additionally, amphipathic acylsugar molecules are believed to reduce the surface tension of adsorbed dew, thereby providing an additional source of water for the plants (Fobes et al., 1985). Acylsugars also have potential applications as pesticides (Puterka et al., 2003), food additives (Hill and Rhode, 1999), cosmetic and personal care products (Hill and Rhode, 1999), antibiotics (Chortyk et al., 1993), and anti-inflammatory compounds (Pérez-Castorena et al., 2010). Due to these important biological roles and potential commercial applications of acylsugars, a detailed understanding of the genetic network involved in acylsugar metabolism is required for successful crop breeding programs and metabolic engineering of acylsugar production.

Figure1
  • Download figure
  • Open in new tab
  • Download powerpoint

In S. pennellii, acylsugars are mostly 2,3,4-tri-O-acylated Glu esters, and some Suc esters, with C4 to C12 fatty acids (Walters and Steffens, 1990; Shapiro et al., 1994; Schilmiller et al., 2015). Predominant branched-chain fatty acids (BCFAs) include 2-methylpropanoic acid, 3-methylbutanoic acid, 2-methylbutanoic acid, and 8-methylnonanoic acid, whereas predominant straight-chain fatty acids (SCFAs) include n-decanoic acid and n-dodecanoic acid. BCFAs are derived from branched-chain amino acids (BCAAs), whereas SCFAs are hypothesized to be derived by a fatty acid synthase (FAS)-mediated de novo fatty acid biosynthetic process (Walters and Steffens, 1990).

Biosynthesis of acylsugars can be divided into two phases: (1) synthesis of fatty acyl chains, and (2) esterification of these acyl molecules to Glu or Suc (Figure 1A). Three acylsugar acyltransferases (ASATs) are involved in phase 2 of acylsucrose biosynthesis (Schilmiller et al., 2015; Fan et al., 2016, 2017). During the preparation of this article, an invertase (acylsucrose fructofuranosidase1; ASFF1) was reported to be capable of producing acylglucose from acylsucrose (Leong et al., 2019). Previous biochemical studies reported a UDP-glucose:fatty acid glucosyltransferase (UDP-Glc:FA GT; Ghangas and Steffens, 1993; Kuai et al., 1997) and a serine carboxypeptidase-like glucose acyltransferase (SCPL GAT; Li et al., 1999; Li and Steffens, 2000) that initiate acylglucose biosynthesis in S. pennellii (Supplemental Figure 1). However, to our knowledge, no genetic evidence is available to support this model.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Current Model of Acylsugar Production in S. pennellii and Workflow to Identify Candidate Genes for Acylsugar Metabolism.

(A) BCFAs are derived from BCAAs, and SCFAs are assumed to be produced by FAS. ASATs use fatty acyl-CoAs for acylsucrose biosynthesis. Recently, an invertase-like enzyme (ASFF1) was reported, which produces acylglucose from acylsucrose. Enzymes are in red font. Unidentified enzymes are indicated with a question mark. Double arrows indicate more than one enzymatic step. C2, 3, 4 indicates esterification at respective positions on Glu. P2, 3, 4 indicates esterification at respective positions on the Suc pyranose ring.

(B) Workflow to identify acylsugar candidate genes. Comparative transcriptomics identified 1,679 DEGs between three low-acylsugar–producing accessions (“LOW”) and three high-acylsugar–producing accessions (“HIGH”). An independent comparison with higher sequencing coverage identified 3,524 DEGs between another low-acylsugar–producing accession, LA1920, and the high-acylsugar–producing accession LA0716. In total, 1,087 DEGs were common to both comparisons. Expression levels of selected candidate genes were determined in LA1302, which accumulates an intermediate level of acylsugars (“MEDIUM”). Analysis of gene expression in response to inhibition of BCAA biosynthesis narrowed the list of candidates, and WGCNA (Langfelder and Horvath, 2008) was performed to identify the coexpressed gene network associated with acylsugar accumulation. Additionally, VIGS led to functional validation of two candidate genes.

Acylsugars are exuded through type-IV trichomes in S. pennellii (Fobes et al., 1985; Slocombe et al., 2008), and many acylsucrose biosynthetic genes are expressed in trichome tip cells (Schilmiller et al., 2012, 2015; Fan et al., 2016). However, analysis of periclinal chimeras indicates that acylglucose biosynthesis in S. pennellii is not trichome-specific (Kuai et al., 1997). Therefore, we used leaf samples for comparative transcriptomics studies reported here. Although this approach underrepresents actual expression differences of trichome-specific genes, it ensures that a complete snapshot of candidate gene networks could be captured.

Different accessions of S. pennellii produce different amounts of total acylsugars (Shapiro et al., 1994). For example, LA0716 produces acylsugars up to 20% of its leaf dry weight (Fobes et al., 1985), whereas in many low-acylsugar–producing accessions, acylsugars make up <1% of leaf dry weight (Shapiro et al., 1994). Here, we exploited this natural genetic variation among S. pennellii accessions to identify candidate biosynthetic and regulatory genes for this trait. Comparative transcriptomics after biochemical inhibition of BCAA biosynthesis further refined the list of candidate genes and illuminated possible regulatory mechanisms of acylsugar production. Functional analysis through virus-induced gene silencing (VIGS) of two candidate genes for SCFA biosynthesis corroborated their involvement in acylsugar metabolism and validated our approach to identifying genes involved in this pathway. Additionally, gene coexpression network analysis revealed genetic networks involved in acylsugar metabolism and plant defense pathways.

RESULTS

Comparative Transcriptomics between Low- and High-Acylsugar–Producing Accessions

We performed differential gene expression analysis between a group of three low-acylsugar–producing accessions (LA1911, LA1912, and LA1926; collectively referred to as “LOW” in this article) and a group of three high-acylsugar–producing accessions (LA1941, LA1946, and LA0716; collectively referred to as “HIGH” in this article) of S. pennellii (Figure 1B; Supplemental Table 1; Shapiro et al., 1994). Three individual plants for each of the six accessions, with an average of >25 million processed paired-end reads for each plant, were used for this analysis (Supplemental Table 2). A total of 19,379 genes passed our filtering criteria for minimum expression levels, and we obtained 1,679 differentially expressed genes (DEGs; see “Methods”; Supplemental Data Set 1: Sheet 1). Of the 1,679 DEGs, 931 were upregulated and 748 were downregulated in the “HIGH” group.

We also compared transcriptomes between another low-acylsugar–producing accession, LA1920, and the high-acylsugar–producing accession LA0716, with each accession having four biological replicates and slightly higher sequencing coverage, as an independent examination of differential gene expression. In this analysis, 19,967 genes passed our filtering criteria for minimum expression levels (see “Methods”), and we identified 3,524 DEGs (1,465 upregulated and 2,059 downregulated genes in LA0716; Supplemental Data Set 1: Sheet 2). To identify candidate genes most likely to be involved in acylsugar metabolism, we compared “LOW” versus “HIGH” DEGs with LA1920 versus LA0716 DEGs (Figure 1B). A total of 1,087 DEGs were common to both comparisons (586 upregulated and 501 downregulated genes in high-acylsugar–producing accessions; Supplemental Data Set 1: Sheet 3). Enrichment analysis revealed significant overrepresentation of gene ontology (GO) terms such as “secondary metabolite biosynthetic process” (GO:0044550), “fatty acid biosynthetic process” (GO:0006633), “branched-chain amino acid biosynthetic process” (GO:0009082), and “defense response” (GO:0006952) in the list of common DEGs (Supplemental Figure 2).

We also examined the leaf transcriptome of accession LA1302 (three biological replicates), which accumulates an intermediate level of acylsugars (referred to as “MEDIUM” in this article; Supplemental Table 1; Shapiro et al., 1994), to determine whether expression levels of candidate genes were consistent with the amount of acylsugar production (Figure 1B). Multidimensional scaling plots showed a clear distinction in gene expression profiles among groups of accessions that accumulate different amounts of acylsugars (Supplemental Figures 3A and 3B). Volcano plots showed high log2(fold-change) values associated with low false discovery rate (FDR) values, indicating robust differential expression (Supplemental Figures 3C and 3D). On the other hand, 10 “housekeeping genes” showed similar expression levels across all biological groups (Supplemental Figure 4).

BCAA Metabolic Genes Are Upregulated in High-Acylsugar–Producing Accessions

Because BCAAs are precursors to acylsugar BCFAs (Walters and Steffens, 1990), we examined genes involved in BCAA metabolism. Many genes involved in both biosynthesis of BCAAs in plastids and conversion of BCAAs to branched-chain acyl molecules in mitochondria (Binder, 2010; Maloney et al., 2010) were upregulated in high-acylsugar–producing accessions (Table 1; Supplemental Figure 5). Most of these BCAA/BCFA metabolic DEGs showed intermediate levels of expression in the “MEDIUM” accession LA1302 (Figure 2A). Correlation analysis (Spearman’s rank correlation coefficient, SRCC; denoted by ρ) using gene expression profiles in 29 samples revealed a strong positive correlation among most of these genes (Figure 3).

View this table:
  • View inline
  • View popup
Table 1. DEGs Involved in Phase 1 of Acylsugar Production (Acyl Chain Synthesis)
Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Heatmaps Showing Expression Levels of Genes with Known and Putative Functions in Acylsugar Metabolism. Genes are designated by their gene identifier numbers (Sopen IDs).

(A) Acylsugar phase-1–related genes (Table 1). Expression levels of many candidate genes in LA1302 (“MEDIUM” accession) were intermediate between low- and high-acylsugar–producing accessions.

(B) Acylsugar phase-2–related genes (Table 2). Sopen04g001210 (marked with an asterisk) was identified as a carboxylesterase gene related to ASH (Schilmiller et al., 2016). Sopen03g040490 (marked with a double-asterisk) has been recently reported as a trichome-expressed invertase gene that is capable of producing acylglucose from acylsucrose (Leong et al., 2019). AG, acylglucose phase 2 (previous model).

(C) Other genes related to acylsugar metabolism (Table 3). Two genes (Sopen04g023150 and Sopen01g047950) were removed from the heatmap due to their very high expression levels. Sopen04g023150 had FPKM values of 570, 624, 455, 179, and 254 in “LOW,” LA1920, LA1302, LA0716, and “HIGH” groups, respectively. Sopen01g047950 had FPKM values of 61, 203, 79, 29, and 19 in “LOW,” LA1920, LA1302, LA0716, and “HIGH” groups, respectively. AAE1, acyl-activating enzyme 1; ABC, ABC transporters; AG, acylglucose phase 2 (previous model); ASAT, acylsucrose acyltransferase; ASH, acylsugar hydrolase; BCAA, branched-chain amino acid; CCM, central carbon metabolism; FAS, fatty acid synthase components; FPKM, fragments per kilobase of transcript per million mapped reads.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Correlation Among Expression Profiles of Selected Genes.

FPKM values for 29 phase-1- and phase-2–related genes (from Tables 1 and 2, respectively) in our 29 samples were used to determine pairwise SRCC. A [29 × 29] matrix was created to visualize correlation. Red, white, and blue colors indicate positive, zero, and negative correlations, respectively. Darker colors indicate stronger correlations. SRCC values and associated P values are given in Supplemental Data Set 2. Sopen04g001210 (marked with an asterisk) is a carboxylesterase gene related to ASH. Sopen03g040490 (marked with a double-asterisk) is the recently identified invertase gene. Abbreviations: AG, acylglucose phase 2 (previous model).

Genes Putatively Encoding a KAS IV/KAS II-like Enzyme, Other FAS Components, and AAEs Are Upregulated in High-Acylsugar–Producing Accessions

Sopen12g004230, predicted to encode a beta-ketoacyl-(acyl-carrier-protein) synthase II (KAS II)-like enzyme, had 39-fold higher expression in the “HIGH” group (FDR = 4.40E−18). The SCFA profile of S. pennellii acylsugars (predominantly C10 and C12, with a trace amount of C8 in LA1302) is strikingly similar to that of the seed oil in many Cuphea species (Shapiro et al., 1994; Dehesh et al., 1998; Slabaugh et al., 1998; Schütt et al., 2002). A specialized version of KAS II, referred to as KAS IV, is an important determinant of chain length in Cuphea medium-chain (C8 to C12) fatty acids (Dehesh et al., 1998; Slabaugh et al., 1998; Schütt et al., 2002). We compared amino acid sequences of five KAS IV/KAS II-like enzymes from four Cuphea species with the Sopen12g004230 sequence, but no similarity was found. However, sequence analysis of the adjacent Sopen12g004240 coding region (1,785-amino acids; annotated as reverse transcriptase) revealed a KAS domain at the C-terminal end (1,374 to 1,784), and BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) analysis of Cuphea KAS IV/KAS II-like enzymes showed high sequence similarity with Sopen12g004240 (77 to 83% identity, 87 to 92% similarity, 82 to 90% coverage, e-value 0.0; Supplemental Table 3). Apparently, a transposon insertion into an intron of Sopen12g004230-Sopen12g004240 led to misannotation of this locus as two separate genes (Supplemental Figure 6A). We confirmed that this locus produces a single transcript with reverse transcription (RT)-PCR analysis, and Sanger sequencing confirmed that this single transcript can produce a functional KAS protein (Supplemental Figures 6B and 6C). An intraspecific F2 population derived from the cross between accessions LA0716 and LA1912 identified a quantitative trait locus on the upper arm of chromosome 12 that is associated with C10 SCFA production (Blauth et al., 1999). The locus containing Sopen12g004230-Sopen12g004240 is in this region.

DEGs putatively encoding a KAS III, other components of the FAS complex, and acyl-activating enzymes (AAEs; ATP/AMP binding proteins that form fatty acyl-CoA molecules from free fatty acids, ATP and CoA; Shockey et al., 2003) were also upregulated in high-acylsugar–producing accessions (Table 1). Expression profiles of these genes showed a very strong positive correlation among themselves, as well as with BCAA/BCFA metabolic DEGs (Figure 3; Supplemental Data Set 2).

Acylsugar Phase-2 Metabolic Genes Are Upregulated in High-Acylsugar–Producing Accessions and Show a Strong Positive Correlation with Phase-1 DEGs

Genes encoding ASATs and the recently identified invertase (SpASFF1), which are involved in phase 2 of acylsugar biosynthesis (Schilmiller et al., 2015; Fan et al., 2016, 2017; Leong et al., 2019), were upregulated in high-acylsugar–producing accessions (Table 2). ASATs use acyl-CoA molecules as their substrates (Schilmiller et al., 2012, 2015; Fan et al., 2016, 2017), and free SCFAs produced by FAS-catalyzed de novo fatty acid biosynthesis must be activated to their acyl-CoA derivatives (SCFA-CoA) before they can be used by ASATs (Figure 1A). All three ASAT genes exhibited very strong positive correlation with both acyl-activating enzyme1 (AAE1) genes and FAS components (Figure 3; Supplemental Data Set 2).

View this table:
  • View inline
  • View popup
Table 2. Genes Involved in Phase 2 of Acylsugar Production

Recently, three genes encoding acylsugar hydrolases (ASHs; carboxylesterases that remove acyl groups from acylsucrose molecules) and a related carboxylesterase gene (Sopen04g001210) were reported in cultivated and wild tomato (Schilmiller et al., 2016). Expression of these genes, except Sopen05g030120 (Sp-ASH1) showed strong positive correlation with genes putatively involved in phase-1 acyl chain synthesis (Figure 3; Supplemental Data Set 2). The expression profile of ASH1 in different tissues of Solanum lycopersicum led the authors to propose additional non-trichome–localized functions for Sl-ASH1 (Schilmiller et al., 2016); our results further support the idea that ASH1 may not have a critical role in trichome acylsucrose metabolism.

Expression Levels of Two Previously Reported Acylglucose Biosynthetic Genes Are Not Consistent with Acylsugar Levels

Expression profiles of two genes encoding UDP-Glc:FA GT and SCPL GAT, respectively (previous model of acylglucose biosynthesis; Supplemental Figure 1; Ghangas and Steffens, 1993; Kuai et al., 1997; Li et al., 1999; Li and Steffens, 2000) were not positively correlated with acylsugar levels (Table 2; Figure 2) and other acylsugar metabolic genes (Figure 3; Supplemental Data Set 2).

Genes Putatively Encoding ATP Binding Cassette Transporters and Central Carbon Metabolic Proteins Are Upregulated in High-Acylsugar–Producing Accessions

Five genes belonging to the ATP-binding cassette (ABC) transporter family showed differential expression, and three of them (Sopen12g034820, Sopen04g005380, and Sopen03g001870) were upregulated in high-acylsugar–producing accessions (Table 3). Recently, the role of central carbon metabolism (CCM) in supporting specialized metabolism in glandular trichomes of S.lycopersicum and Solanum habrochaites was investigated, and it was concluded that glandular trichomes import most of their fixed carbon from underlying leaf tissues, despite being able to perform some photosynthetic activities (Balcke et al., 2017). Special roles of chloroplast malic enzyme in supplying pyruvate and glutathione s-transferase (GST) in alleviating oxidative stress were also reported. We did find genes putatively encoding these central carbon metabolic (CCM) proteins in the list of 1,087 DEGs (Table 3; Figure 2C).

View this table:
  • View inline
  • View popup
Table 3. degs Involved in Transport Process, ccm, and Regulation of Gene Expression

Trichome-Enriched Expression of Acylsugar Candidate Genes

Genes encoding ASATs and the recently identified invertase (SpASFF1) are expressed in apical cells of type-I/IV glandular trichomes (Schilmiller et al., 2012, 2015; Fan et al., 2016; Leong et al., 2019). Using data of Ning et al. (2015), who published expression profiles of S. lycopersicum genes in isolated stem trichomes versus shaved stems, we investigated if S. lycopersicum homologs of candidate genes show a trichome-enriched expression pattern. We used reciprocal BLAST to identify putative S. lycopersicum orthologs of candidate genes (Supplemental Data Set 3), and they showed trichome-enriched expression pattern (26- to 562-fold for AAE1 members, 143- to 470-fold for ABC transporter genes, and 16- to 287-fold for CCM genes; Supplemental Table 4).

Genes Putatively Encoding Oxylipin Metabolic Proteins and Other Defense Proteins Are Upregulated in Low-Acylsugar–Producing Accessions

Defense response genes were present in the list of 1,087 DEGs (GO:0006952; Supplemental Figure 2), and we found that genes putatively involved in oxylipin biosynthesis (fatty acid desaturase [FAD] and lipoxygenase [LOX]), metabolism (oxo-phytodienoate reductase and GST), and signaling pathways (TGA transcription factors; Mueller et al., 2008), as well as those putatively encoding subtilisin-like proteases (SLPs) and additional defense proteins were upregulated in low-acylsugar–producing accessions (Table 4; Figure 4). Based on expression profiles, these putative defense genes showed moderate to strong correlation with each other (Supplemental Figure 7), indicating an active network of innate immunity in low-acylsugar–producing accessions. According to data of Ning et al. (2015), putative S. lycopersicum orthologs of these defense DEGs do not show trichome-enriched expression (except two of 28 DEGs; Sopen09g006360 and Sopen10g019020; Supplemental Table 4).

View this table:
  • View inline
  • View popup
Table 4. degs Putatively Involved in Oxylipin-Related and Amine Oxidase-Related Defense Systems
Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Heatmap Showing Expression Levels of Genes Putatively Involved in Plant Defense (Table 4 Genes).

Expression levels of many genes in LA1302 were similar to low-acylsugar–producing accessions but different from high-acylsugar–producing accessions. Some oxylipin-related defense genes, which were upregulated in low-acylsugar–producing accessions, had higher expression levels than amine oxidase (AOX)-related defense genes, which were upregulated in high-acylsugar–producing accessions. OLM, oxylipin metabolism.

Amine Oxidase-Related Defense Response Genes Are Upregulated in High-Acylsugar–Producing Accessions

Genes putatively encoding amine oxidase domains (FLOWERING LOCUS D-like proteins; required for systemic acquired resistance in Arabidopsis [Arabidopsis thaliana]; Singh et al., 2013) and a corresponding aldehyde oxidase (Cona et al., 2006) were upregulated in high-acylsugar–producing accessions (Table 4). H2O2 generated by amine oxidase- and aldehyde oxidase-catalyzed reactions acts as a signal molecule and plays an important role during pathogen invasion by acting on extensin proteins (Niebel et al., 1993; Jackson et al., 2001; Cona et al., 2006). Two genes predicted to encode extensin-like glycoproteins were also upregulated in high-acylsugar–producing accessions (Table 4), and one of them (Sopen05g003410) showed strong positive correlation with amine oxidase- and aldehyde oxidase-encoding genes (Supplemental Figure 7). Although amine oxidase-related defense genes were upregulated in high-acylsugar–producing accessions, it is important to note that their expression levels were much lower than some oxylipin-related defense genes in low-acylsugar–producing accessions (Figure 4). Among these amine oxidase-related DEGs, we identified a putative S. lycopersicum ortholog (Solyc01g005850) for only one gene (Sopen01g001850; extension-like protein), and according to the data of Ning et al. (2015), Solyc01g005850 does not show trichome-enriched expression (Supplemental Table 4).

Biochemical Inhibition of BCAA Biosynthesis Changes Transcript Levels of Many Acylsugar Candidate DEGs, But Not Defense DEGs

Acetolactate synthase (ALS) catalyzes the first common step in the biosynthesis of BCAAs (Supplemental Figure 5A), and biochemical inhibition of this enzyme significantly lowers acylsugar production in S. pennellii (Walters and Steffens, 1990). To determine whether inhibition of this enzyme changed expression of ALS and other candidate genes involved in acylsugar metabolism, we treated leaves of S. pennellii (the high-acylsugar–producing accession LA0716) with the ALS inhibitor imazapyr at 0.1 mM and 1 mM concentrations.

Expression of 171 of the 1,087 DEGs previously identified by comparative transcriptomics was affected by imazapyr treatment at both concentrations, including many of the DEGs involved in BCAA/BCFA metabolism (Table 5; Supplemental Data Set 4: Sheet 1). It should be noted that at 1 mM imazapyr, expression of Sopen04g026270 and Sopen01g028100 (branched-chain keto acid dehydrogenase complex E1 and E2 subunits, respectively) was significantly increased, possibly because the amount of 2-ketobutyrate, a substrate of ALS that also can be used directly by the dehydrogenase complex to make propionyl-CoA, builds up (Walters and Steffens, 1990).

View this table:
  • View inline
  • View popup
Table 5. Response of Selected Genes to Imazapyr Treatment

DEGs putatively involved in SCFA metabolism (FAS componentscand AAE1 members) also showed significant downregulation at both 0.1 mM and 1 mM imazypyr in a concentration-dependent manner. AAE1 proteins have a peroxisomal location in Arabidopsis (Reumann et al., 2009); one gene (Sopen11g007710) predicted to encode a PMP22/Mpv17 family peroxisomal membrane protein had 7.7-fold (FDR = 1.17e−08) and 165-fold (FDR = 3.62E−24) lower expression at 0.1 mM and 1 mM imazapyr treatment, respectively (Supplemental Data Set 4: Sheet 1). Sopen11g007710 was upregulated 29-fold in the “HIGH” group (FDR = 2.80E–11; Supplemental Data Set 1), and its putative ortholog in S. lycopersicum (Solyc11g012990; Supplemental Data Set 3) has 1,420-fold higher expression in isolated trichomes than shaved stems, consistent with trichome-enriched expression of AAE1 members (Ning et al., 2015). PMP22 is presumably involved in controlling permeability of peroxisomal membrane (Brosius et al., 2002). These results suggest a role of peroxisome in acylsugar metabolism.

DEGs encoding ASATs and the invertase showed significant decrease in gene expression level in a concentration-dependent manner (by as much as 94-fold, FDR = 1.65E–14), as did ASHs and the related carboxylesterase Sopen04g001210 (by as much as 70-fold, FDR = 8.82E–50; Table 5). On the other hand, Sopen01g049990 and Sopen10g020280 (encoding UDP-Glc:FA GT and SCPL GAT, respectively; previous model of acylglucose biosynthesis) showed slightly higher and similar expression levels, respectively, in response to imazapyr treatment.

DEGs putatively encoding CCM proteins and three ABC transporters were significantly downregulated (by as much as 45-fold, FDR = 4.99E−10 for RUBISCO small subunit; Table 5). Imazapyr treatment in Arabidopsis resulted in significant induction of nine genes (out of the 10 genes reported) encoding ABC transporters, indicating their role in detoxification process (Manabe et al., 2007). Similarly, 15 putative ABC transporter genes in S. pennellii were upregulated in response to imazapyr either at higher concentration or at both concentrations (Supplemental Table 5). However, three putative ABC transporter DEGs were repressed by imazapyr treatment, consistent with expression profiles of other acylsugar metabolic genes.

Unlike known and candidate acylsugar metabolic DEGs, amine oxidase-related defense DEGs did not show concentration-dependent downregulation in response to ALS inhibition (Supplemental Data Set 4: Sheet 2). The overall pattern of gene expression in response to imazapyr suggests that additional genes involved in acylsugar metabolism may be among these 171 DEGs. Of the 13 transcription factor DEGs in the initial list of 1,087, nine were present in the imazapyr treatment list (Table 3; Figure 2C). Putative S. lycopersicum orthologs of six of these nine DEGs have a trichome-enriched expression profile (Supplemental Table 4; Ning et al., 2015).

Identification of Coexpressed Gene Network Associated with Acylsugar Metabolism

Using expression profiles of S. pennellii genes in our 38 samples (29 from different accessions and nine from the imazapyr treatments), weighted gene correlation network analysis (WGCNA; Langfelder and Horvath, 2008) identified 41 coexpressed gene modules. Known and candidate acylsugar metabolic genes were clustered in the “darkorange” module (182 genes; Figure 5A; Supplemental Figure 8). One-third of the genes that responded to imazapyr (57/171) were also placed in this module (Supplemental Data Set 4: Sheet 3). We selected the top-50 most strongly connected genes (based on intramodular connectivity as determined by WGCNA) for each of the three ASAT genes and the invertase gene, and merged them to identify the most strongly interconnected gene network. This network included genes involved in BCAA/BCFA metabolism, SCFA synthesis (FAS components) and their activation (AAEs and the related peroxisomal membrane protein PMP22/Mpv17), acylsucrose metabolism (three ASATs, ASH2, ASH3 and the related carboxylesterase Sopen04g001210), as well as genes putatively encoding three ABC transporters, CCM proteins, six transcription factors (including three AP2-family transcription factors), and other proteins (Figure 5B).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

WGCNA (Langfelder and Horvath, 2008).

(A) Association of module eigengenes (MEs) with expression profiles of nine selected genes. WGCNA identified 41 modules; for each module, an ME was selected as representative of the expression profile for that module (first principle component). Numbers indicate SRCC values between MEs of two selected modules and nine genes’ expression profiles. P values associated with SRCC are given in parentheses (9E-23 = 9 × 10−23). Blue, white, and red colors indicate negative, zero, and positive correlations, respectively. Acylsugar metabolic genes (ASATs and invertase) and LOX-correlated genes were clustered in “darkorange” and “green” modules, respectively. NB-ARC-LRR indicates a recently characterized sequence (Sopen09g036080) that confers resistance against tospoviruses in S. lycopersicum and wild tomato species (Zhu et al., 2017). A complete list of WGCNA modules is given in Supplemental Figure 8.

(B) Simplified acylsugar metabolic gene network. Nodes and edges represent genes and intramodular connectivities, respectively. We used WGCNA to identify the top-50 genes most strongly connected to each of the three ASAT genes and the invertase gene. These four top-50 lists contained only 58 different genes.

However, not all genes in this network had an obvious connection to acylsugar metabolism. One gene on chromosome 11 (Sopen11g003320; UDP-glucose:catechin glucosyltransferase; Noguchi et al., 2008) and three sequential genes on chromosome 6 (Sopen06g034810, Sopen06g034820, and Sopen06g034830; myricetin methyltransferase (Schmidt et al., 2012; Kim et al., 2014) showed strong intramodular connectivity with acylsugar metabolic genes. These flavonoid metabolic genes were also strongly downregulated in response to imazapyr treatment (Supplemental Figure 9), and their putative orthologs in S. lycopersicum are preferentially expressed in trichomes (Ning et al., 2015). In addition to abundant acylsugar compounds, flavonoid compounds, such as methylated myricetin, have been reported in S. pennellii trichomes (McDowell et al., 2011). Together, these results suggest that metabolisms of acylsugar and flavonoid compounds are strongly connected.

Plant Defense Gene Network

We selected Sopen01g002520 (LOX) as a representative of the oxylipin-mediated defense system, due to its robust differential expression (Figure 4; Table 4), and WGCNA identified coexpressed genes in the “green” module. Enriched GO terms in this module include “defense response” (GO:0006952), “oxylipin biosynthetic process” (GO:0031408), and “glutathione metabolic process” (GO:0006749; Figure 6; Supplemental Figure 10). In the plant innate immune system, recognition and successful prevention of pathogen invasion rely on a series of molecular events, including pathogen detection at the cell surface by receptor kinases, a series of phosphorylation events with cytoplasmic kinases, nucleotide binding, release of intracellular calcium, protein ubiquitination, and transcriptional reprogramming (Couto and Zipfel, 2016). GO terms associated with these processes were significantly enriched in the “green” module (FDR = 8.83E−17 to 4.6E−02). Twenty of 104 expressed S. pennellii NB-ARC sequences and 23 of 150 leucine-rich repeat (LRR) sequences were clustered in the “green” module (Supplemental Data Set 5), including the recently characterized NB-ARC-LRR sequence Sw-5b (Sopen09g036080) that confers broad-spectrum resistance against American-type tospoviruses in S. lycopersicum and wild tomato species (Zhu et al., 2017).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Enrichment of Selected GO Terms Associated with WGCNA “Green” Module Genes.

Test set and reference set indicate “green” module genes and S. pennellii annotated sequences, respectively. Enriched sequences include those putatively involved in oxylipin metabolism and plant defense signaling (Couto and Zipfel, 2016). Enrichment analysis was performed using Blast2GO software (Fisher’s exact test). Only significant GO terms are shown (FDR < 0.05). A complete list of all significant GO terms, GO IDs, and associated FDR and P values is given in Supplemental Figure 10.

Functional Validation of Two Candidate Genes Involved in SCFA Biosynthesis

To our knowledge, there are no reports of genes involved in acylsugar SCFA synthesis. We selected two candidate genes for functional validation (Sopen08g002520 and Sopen12g004240; predicted to encode KAS III and KAS IV/KAS II-like enzymes, respectively). These genes have 235- and 60-fold, respectively, higher expression levels in isolated stem trichomes compared with underlying tissue (Figure 7A). In fact, Sopen08g002520 showed a SpASAT1-like trichome-enriched expression profile. We identified Solyc08g006560 and Solyc12g009260 as putative S. lycopersicum orthologs of Sopen08g002520 and Sopen12g004240, respectively, and these two tomato orthologs also show trichome-enriched expression profiles (90- and 103-fold, respectively; Supplemental Table 4; Ning et al. (2015)). Phylogenetic analyses supported these orthologous relationships and revealed distinct clades of KAS III and KAS II sequences in solanaceous species (Supplemental Figure 11; Supplemental Files 1 to 4).

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Functional Validation of Two Candidate KAS Genes.

(A) Transcript levels of Sopen08g002520 (KAS III) and Sopen12g004240 (KAS IV/II-like) in isolated stem trichomes and underlying tissues (shaved stems) were measured by RT-qPCR. Error bars indicate se (n = 5 individual plants).

(B) Acyl chain composition of acylsugars extracted from leaves of control plants (empty TRV vector) and two other groups of plants in which Sopen08g002520 and Sopen12g004240, respectively, were targeted for VIGS. VIGS of both Sopen08g002520 and Sopen12g004240 resulted in significant decrease in SCFAs (n-decanoic acid and n-dodecanoic acid). Relative abundance of acyl chains was measured by GC–MS. Asterisks denote statistical significance compared with control plants as determined by one-way ANOVA (Dunnett’s test; *P < 0.05, **P < 0.01, ***P < 0.001). Error bars indicate se (n = 10 individual plants).

We targeted these two S. pennellii loci in LA0716 for VIGS using tobacco rattle virus (TRV)-based silencing vectors, and acylsugar SCFA production was reduced by up to 40% in both n-C10 and n-C12 amount (P < 0.001; Figure 7B). VIGS resulted in significant downregulation of target genes (60% and 67% reduction in transcript levels for Sopen08g002520 and Sopen12g004240, respectively; Supplemental Figures 12A and 12B). Residual expression presumably is due to the fact that VIGS in S. pennellii is incomplete and inconsistent (Supplemental Figure 12C). Additionally, no morphologically distinct areas would allow for better area-selective metabolite profiling, because no significant morphological differences were observed between silenced and control plants (Supplemental Figure 12D). Nevertheless, our results indicate that the two candidate genes, which are preferentially expressed in trichomes, play important roles in acylsugar SCFA biosynthesis. This is consistent with the previous prediction that KAS enzyme(s) other than KAS I might be involved in SCFA production in S. pennellii acylsugars (Slocombe et al., 2008).

DISCUSSION

Acylsugars are powerful natural pesticides (Puterka et al., 2003), and increasing acylsugar-mediated insect pest resistance has been a target of various breeding programs in solanaceous crops (Bonierbale et al., 1994; Lawson et al., 1997; Alba et al., 2009; Leckie et al., 2012). The success of these breeding programs and bioengineering of acylsugar production will largely depend on unraveling the network of structural and regulatory genes. Here, we used transcriptomics approaches to identify candidate acylsugar and defense gene networks in S. pennellii.

Defense Pathway Genes

Enriched GO terms associated with pathogen recognition and downstream signaling cascades in the “green” module indicate an active defense network in low-acylsugar–producing accessions (Figure 6). Oxylipins (Halitschke and Baldwin, 2003; Prost et al., 2005) and SLPs (Figueiredo et al., 2014) have important roles in plant defense, and higher expression levels of these defense genes predominantly in low-acylsugar–producing accessions (Table 4) may indicate a compensatory mechanism for diminished defense activities of acylsugars in these accessions. Alternatively, expression of these defense genes may be suppressed in high-acylsugar–producing accessions due to the protection provided by acylsugars. Our data do not allow us to distinguish between these two possibilities, but it is important to note that all plants used in this study were grown in the same growth chamber, and no insects or symptoms of disease were apparent on any plant.

Most defense DEGs (both oxylipin-related and amine oxidase-related) had similar expression levels in the “MEDIUM” accession LA1302 and in low-acylsugar–producing accessions (Figure 4). Multidimensional scaling plots also showed that the overall gene expression profile of LA1302 is similar to that of low-acylsugar–producing accessions, but very different from that of high-acylsugar–producing accessions (Supplemental Figures 3A and 3B). These results are consistent with the prediction that LA1302 is more closely related to low-acylsugar–producing accessions (Shapiro et al., 1994). Despite the similarity of defense DEGs expression and the overall transcriptome profile between LA1302 and low-acylsugar–producing accessions, many DEGs with known and putative functions in acylsugar metabolism showed intermediate expression levels in LA1302 (Figure 2). Together, these results suggest that there is no strict linear relationship between acylsugar accumulation and expression of defense pathways.

Acylsugar Metabolic Genes

Our initial hypothesis for this work was that expression levels of previously reported acylsugar biosynthetic genes would be consistent with acylsugar levels in different accessions of S. pennellii, and that additional, unknown genes involved in acylsugar metabolism would be coordinately expressed with them. Our approach was validated by the identification of several DEGs with known functions in acylsugar biosynthesis, such as BCAA/BCFA metabolic genes and ASATs. Next, we used a chemical genetics approach to refine the list of acylsugar candidate genes by treating S. pennellii leaves with imazapyr to inhibit BCAA biosynthesis. This approach allowed us to separate the acylsugar metabolic DEGs from defense DEGs, which did not respond to imazapyr treatment (Table 5; Supplemental Data Set 4). Although we expected that genes related to synthesis of acylsugar BCFAs (derivatives of BCAAs) would be downregulated as a result of lack of substrate availability in response to imazapyr, our results showed that genes related to SCFA metabolism and other aspects of acylsugar metabolism (including known acylsucrose metabolic genes and candidate genes) also responded to imazapyr in a dose-dependent manner. It is important to note that the repression of acylsugar metabolic genes we observed after imazapyr treatment is distinct from the induction of general stress response genes that can result from amino acid starvation (Zhao et al., 1998).

Strong correlations among expression profiles of selected acylsugar metabolic genes (Figure 3) encouraged us to perform gene coexpression network analysis, which identified candidate gene networks involved in acylsugar metabolism and also plant defense (Figures 5 and 6; Supplemental Data Set 5). Additionally, data of Ning et al. (2015) indicate that putative S. lycopersicum orthologs of most of our acylsugar candidate genes, but not defense DEGs, are preferentially expressed in trichomes (Supplemental Table 4). Together, these results strongly support the candidacy of novel acylsugar metabolic genes that we have identified; candidate genes include those encoding KAS enzymes and other FAS components (for SCFA synthesis), AAE1 members (for activating SCFAs for ASATs), a PMP22/Mpv17 family peroxisomal membrane protein (for controlling permeability of peroxisomal membrane), ABC transporters (for acylsugar export), and CCM proteins (for providing carbon to support acylsugar production), including a glutathione s-transferase (for alleviating oxidative stress).

In plant primary metabolism, the condensation reactions of fatty acid biosynthesis are catalyzed by KAS enzymes of the FAS complex, and at least three different KAS isoenzymes (KAS I, KAS II, and KAS III) are involved in de novo fatty acid biosynthesis. We selected two candidate KAS genes that are preferentially, if not exclusively, expressed in S. pennellii trichomes, and knockdown of these two genes led to significant decrease in SCFA synthesis (Figure 7), demonstrating their role in acylsugar metabolism. Duplication of a primary metabolic gene followed by neofunctionalization (Ning et al., 2015) and spatial gene regulation have important functions in regulating trichome metabolism of acylsugar molecules, and KAS sequences also appear to follow that trend.

Regulation of Acylsugar Production

Both comparative transcriptomics and imazapyr treatment studies revealed that transcript levels for the two previously reported phase-2 acylglucose biosynthetic enzymes (UDP-Glc:FA GT and SCPL GAT) were not positively correlated with acylsugar levels (Tables 2 and 5). Recently, Leong et al. (2019) clearly demonstrated the role of the trichome-expressed invertase (SpASFF1) in acylglucose biosynthesis in S. pennellii trichomes. Their results contradict the involvement of UDP-Glc:FA GT and SCPL GAT in the biosynthesis of acylglucose molecules, and our results are consistent with this finding. Our results also suggest that a majority of the acylsugar metabolic pathway, from acyl chain synthesis to acylsugar secretion, represents a remarkable coexpressed genetic network. WGCNA revealed that six transcription factor genes form a strong gene coexpression network with many acylsugar and flavonoid metabolic genes. Downregulation of these transcription factor genes in response to imazapyr treatment parallels significant downregulation of most acylsugar and flavonoid metabolic genes (Table 5; Supplemental Data Set 4). Additionally, putative orthologs of these transcription factor genes in S. lycopersicum exhibit trichome-enriched expression, similar to many genes involved in acylsugar and flavonoid metabolism (Ning et al., 2015). Together, these results suggest important roles for these transcription factors in specialized metabolic pathways of S. pennellii trichomes.

METHODS

Plant Materials

Seeds of all Solanum pennellii accessions were obtained from the C.M. Rick Tomato Genetics Resource Center (University of California, Davis). Seeds were treated with a solution of 20% (v/v) commercial bleach (1.2% [w/v] sodium hypochlorite) for 20 min, and then germinated on moistened filter paper. Seedlings were transferred to soil at the cotyledon stage, and plants were grown in a growth chamber (24°C day/20°C night temperature, 150 mMol m−2 s−1 photosynthetically active radiation, 16-h photoperiod, and 75% humidity).

RNA Sequencing Library Preparation and Sequencing

Three biological replicates were used for each accession in the “LOW” versus “HIGH” comparison, and for the “MEDIUM” accession LA1302. Four biological replicates were used for each accession in the LA1920-versus-LA0716 comparison. Before RNA extraction, secreted acylsugars were removed from the leaf surface of 10-week–old plants by dipping them in ethanol for 2 to 3 s. Leaves were then immediately frozen in liquid nitrogen and stored at −80°C until further use. Total RNA was extracted from leaves using the RNAqueous Total RNA Isolation Kit (Thermo Fisher Scientific). Genomic DNA was removed using the TURBO DNA-Free Kit (Thermo Fisher Scientific). Quality of each RNA sample was analyzed with the Agilent 2200 TapeStation software A01.04 (Agilent Technologies). RNA sequencing (RNA-Seq) libraries were prepared from polyA+-selected RNA samples using the TruSeq RNA Library Preparation Kit v2 (Illumina), and then sequenced on the HiSeq 2500 v4 125-bp × 125-bp paired-end sequencing platform (High Output Mode; Illumina) following the manufacturer’s specifications. Sequencing was performed at the Texas A&M Genomics and Bioinformatics Service Center. Sequence cluster identification, quality prefiltering, base calling, and uncertainty assessment were done in real time using Illumina’s HCS v2.2.68 and RTA v1.18.66.3 software with default parameter setting. Sequencer.bcl basecall files were demultiplexed and formatted into .fastq files using bcl2fastq v2.17.1.14 script conFigureBclToFastq.pl. Sequencing reads were submitted to the National Center for Biotechnology Information Sequence Read Archive under accession number SRP136022.

Quality Control and Mapping of Sequencing Reads

Phred quality score distributions of sequencing reads were analyzed with the software FastQC v0.11.4 (Andrews, 2010). All RNA-Seq libraries had average Phred quality scores of >35, indicating >99.97% base call accuracy. To remove all the low-quality, ambiguous “N” nucleotides (Phred quality score = 2, ASCII code = #), the following settings were applied for quality control of RNA-Seq reads using the program Trimmomatic v0.35 (Bolger et al., 2014b): ILLUMINACLIP:TruSeq3-PE-2.fa:2:30:10, SLIDINGWINDOW = 1:14, MINLEN = 80. Using these settings, >98% of the sequencing reads could be retained, and it ensured that none of the nucleotides in the reads were ambiguous.

Processed reads were mapped to the S. pennellii genome v2.0 (LA0716; Bolger et al., 2014a) using TopHat2 v2.1.0 (Kim et al., 2013). The following parameters were used for the mapping process: -mate-inner-dist = 0, -mate-std-dev = 50, -read-realign-edit-dist = 0, -read-edit-dist = 4, -library-type = fr-unstranded, -read-mismatches = 4, -bowtie-n, -min-anchor-len = 8, -splice-mismatches = 0, -min-intron-length = 50, -max-intron-length = 50,000, -max-insertion-length = 3, -max-deletion-length = 3, -max-multihits = 20, -min-segment-intron = 50, -max-segment-intron = 50,000, -segment-mismatches = 2, -segment-length = 25. A summary of TopHat2 alignment results are given in Supplemental Table 2. Reads mapped to selected loci of interest were visualized with IGV (Thorvaldsdóttir et al., 2013).

DEG Analyses

Uniquely-mapped reads from TopHat2 were counted for all S. pennellii annotated genes using HTSeq-Count (Anders et al., 2015). The following parameters were used for the counting process: -f bam, -r name, -s no, -m union, -a 20. These count files were used for identification of DEGs using the program edgeR (Robinson et al., 2010). Genes with very low expression levels were filtered out before differential testing. For comparison between the “LOW” and “HIGH” groups, genes with >1 count per million (CPM) in at least six samples were retained for analysis after read counts were adjusted for transcript lengths. Library sizes were recalculated after the filtering process, and all the samples had post-filter library sizes of >99.70%. Tagwise dispersions were calculated (prior.df = 3), and genes were identified as differentially expressed when P values, corrected for multiple testing using Benjamini–Hochberg multiple testing correction (Benjamini and Hochberg, 1995), were <0.05 (FDR < 0.05), and had at least 2-fold expression differences. For comparison between LA1920 and LA0716, we used the same pipeline, except genes with >1 CPM in at least two samples were retained for analysis, and tagwise dispersions were calculated with prior.df = 9. In this analysis, all the samples showed >99.85% post-filter library sizes. Lists of DEGs as determined by edgeR (Robinson et al., 2010) are given in Supplemental Data Set 1.

Functional Annotation of DEGs

Transcript sequences of DEGs were extracted from the annotated gene models, and Blast2GO (Conesa et al., 2005) was used for the functional annotation of these transcript sequences. Sequences were searched (BLASTx; word size = 3, HSP length cutoff = 33) against the National Center for Biotechnology Information nonredundant database (subset Viridiplantae, taxa: 33,090) with e-value ≤ 1.0E−10, and the top-20 BLAST hits were reported. BLAST hits for each sequence were then mapped with the GO terms. Next, annotation was performed with the default parameter settings except the e-value threshold (e-value ≤ 1.0E−10). This process also generated the Enzyme Code numbers from the Kyoto Encyclopedia of Genes and Genomes pathway. After executing the whole functional annotation process (BLASTx, mapping, and annotation), transcript sequences were exported with sequence description and GO terms. Descriptions of transcript sequences were generated based on similarity levels (e-value and percentage of identity) to the subject genes, as determined by Blast2GO (Conesa et al., 2005).

Identification of Putative Orthologs and Their Trichome-Enriched Expression

Reciprocal BLAST was used to identify putative orthologs of S. pennellii genes in S. lycopersicum. An all-versus-all BLAST was performed between annotated proteins of S. pennellii v2.0 (Bolger et al., 2014a) and S. lycopersicum ITAG2.3 (Tomato Genome Consortium, 2012; Fernandez-Pozo et al., 2015) with the following parameters: minimum percentage identity = 70, minimum percentage query coverage = 50. Pairs with the lowest e-value for each BLASTp comparison were selected as putative orthologs. Gene expression profiles of S. pennellii putative orthologs in S. lycopersicum were obtained from Ning et al. (2015).

Gene Coexpression Network Analysis

Simple correlation analysis among expression profiles (fragments per kilobase of transcript per million mapped reads [FPKM] values in our 29 samples) of selected genes was performed using the “cor” function in the R programming language (R Core Team, 2014). Gene coexpression network analysis and identification of modules (clusters of strongly correlated genes) were performed using the R package WGCNA (Langfelder and Horvath, 2008). Normalized gene expression data (FPKM values) in our 38 samples were used as the input, and the 19,378 genes which passed our filtering criteria for minimum expression levels during “LOW” versus “HIGH” differential testing (Sopen12g004230 and Sopen12g004240 were merged) were selected for WGCNA analysis. A matrix of pairwise SRCCs between these 19,378 genes was created and then transformed into an adjacency (connection strength) matrix using “signed” network type and a soft threshold (β) value of 12. The value of β was determined based on the scale-free topology fit index (R2 > 0.85) and low mean connectivity. The weighted adjacency matrix was then converted to a topological overlap matrix (TOM), and genes were hierarchically clustered based on dissimilarity (dissTOM = 1 - TOM) of the topological overlap. Coexpression modules were defined as the branches of the clustering tree (dendrogram) using the following parameters in “cutreeDynamic” function: distM = disTOM, deepSplit = 2, minClusterSize = 50. A module eigengene (ME; representative of gene expression profile in each module) was defined as the first principle component in each module. A set of 41 modules was obtained. For each gene, a “module membership” was calculated as the correlation (SRCC) between expression profile of that gene and ME. Associations between external traits and modules were measured as determined by SRCC between trait data and MEs; “gene significance” values for each gene in a module of interest were calculated as the correlation (SRCC) between expression profile of a gene and trait data. The top-50 intramodular connections for selected genes were visualized with the program VisANT (Hu et al., 2005).

Inhibitor Study (S. pennellii Accession LA0716)

Compound leaves of S. pennellii LA0716 bearing five leaflets were used for the inhibitor study. Imazapyr (Sigma-Aldrich) was dissolved in distilled water, and administered as 0.1 mM and 1 mM solutions for 36 h. Cuttings were transferred to small trays containing imazapyr or control solution (deionized water) and placed under growth chamber conditions (24°C day/20°C night temperature, 150 mMol m−2 s−1 photosynthetically active radiation, 16-h photoperiod, and 75% humidity). Total RNA was extracted as described before, and RNA-Seq libraries were prepared using the TruSeq Stranded RNA LT Kit (Illumina). Libraries were sequenced on the HiSeq 4000 (Illumina) 150-bp × 150-bp paired-end sequencing platform at the Texas A & M Genomics and Bioinformatics Service Center, College Station. Raw reads were processed with the software Trimmomatic (Bolger et al., 2014b), and read-mapping to the S. pennellii genome was performed with the software Tophat2 (Kim et al., 2013) using the same parameters as described in the Quality Control and Mapping of Sequencing of Reads section, except with the following change: -library-type = fr-firststrand. Uniquely-mapped reads from TopHat2 were counted with HTSeq-Count (Anders et al., 2015) using the following options: -f bam, -r name, -s reverse, -m union, -a 20. DEGs (FDR < 0.05; log2FC > 1) were identified using the software edgeR (Robinson et al., 2010) with tagwise dispersions (prior.df = 12), and genes with >1 CPM in at least three samples were retained for the analysis. Summary results of the inhibitor study are given in Supplemental Data Set 6.

RT Quantitative PCR (S. pennellii Accession LA0716)

Stems of 10-week–old S. pennellii LA0716 plants were cut into small pieces and immediately frozen with liquid nitrogen. Stem trichomes were isolated by gentle scraping with a scalpel. Total RNA was extracted using the RNAqueous Total RNA Isolation kit (Thermo Fisher Scientific), followed by genomic DNA digestion using the TURBO DNA-free kit (Thermo Fisher Scientific), and 1 μg of total RNA was reverse-transcribed using the SuperScript IV VILO Master Mix (Thermo Fisher Scientific), following the manufacturer’s guidelines. Quantitative PCR (qPCR) was performed using SYBR Green Master Mix (Bio-Rad) with QuantStudio 6 Flex System, and the following PCR program was used: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Differential gene expression was measured using the comparative CT (ΔΔCT) method. Three housekeeping genes (actin, ubiquitin, and TATA-box binding protein) were assessed for their use as endogenous control(s), and final normalization was performed with the ubiquitin gene Sopen02g027600. For VIGS RT-qPCR analysis, similar-sized young leaves were used for RNA extraction, and RT-qPCR primers were designed outside the VIGS-targeted region. A list of primers used in this study is given in Supplemental Table 6.

VIGS and Metabolite Profiling (S. pennellii Accession LA0716)

VIGS constructs were designed using the SGN online tool (http://vigs.solgenomics.net/). Approximately 495-bp to 510-bp target fragments were cloned into the pTRV2-LIC vector (Dong et al., 2007). Agrobacterium tumefaciens strain GV3101 carrying pTRV1 and pTRV2 vectors were grown overnight at 28°C. Overnight cultures were centrifuged at 8,000g for 5 min at 4°C, and cells were resuspended in infiltration medium (10 mM of MES at pH 5.5, 10 mM of MgCl2, and 200 μM of acetosyringone). TRV1 and TRV2 suspensions were incubated at room temperature for 3 h, and then mixed in a 1:1 ratio with a final OD600 of 1. Infiltration of S. pennellii LA0716 was performed at the first true leaf stage, and all plants were grown for about six weeks in the same growth chamber mentioned before. A phytoene desaturase gene (Sopen03g041530) was used as a positive control to assess the efficiency of VIGS (Supplemental Figure 12C).

For the analysis of acylsugar acyl chain composition, secreted acylsugars were extracted from similar-sized young leaves by dipping them in ethanol. Ethanol was evaporated, and acylsugars were dissolved in n-heptane. Transesterification reaction was performed by adding 500 μL of 20% sodium ethoxide to 1 mL of heptane mixture at room temperature for 10 min with gentle vortexing. The reaction mixture was washed twice with 500 μL of saturated sodium chloride solution, and the heptane layer was used for gas chromatography-mass spectrometry (GC–MS) analysis using Trace GC Ultra with DSQII system (Thermo Fisher Scientific). A 30-m × 0.25-mm column with 0.25-μm thickness (DB-5MS; Agilent) was used as the GC column. A quantity of 1 μL of heptane extract was used for split injection (1:50) with injector temperature of 225°C, and the following GC temperature program was used: starting temperature of 30°C held for 2 min, increased at 15°C/min up to 300°C and held for 5 min. The carrier gas flow was 1.5 mL/min. Thermo Xcalibur (v3.0.63) and Qual Browser applications were used for data acquisition and analysis. Summary of statistical tests is given in Supplemental Data Set 7.

Phylogenetic Analyses

Predicted transit peptide sequences were removed from protein sequences before multiple sequence alignment using MUSCLE (Edgar, 2004). The evolutionary history was inferred by using the maximum likelihood method based on the JTT+G model (Jones et al., 1992; lowest Bayesian Information Criterion value). Bootstrap values were obtained from 1,000 replicates. All positions containing gaps and missing data were eliminated. Evolutionary analyses were conducted in the software MEGA7 (Kumar et al., 2016). Trees were drawn to scale, with branch lengths measured in the number of substitutions per site.

Accession Numbers

RNA-Seq reads used in this study were submitted to the National Center for Biotechnology Information Sequence Read Archive under accession numbers SRP136022 and SRP136034 (imazapyr study).

Supplemental Data

  • Supplemental Figure 1. Previous model of acylglucose biosynthesis.

  • Supplemental Figure 2. Enriched GO terms associated with the 1,087 common DEGs.

  • Supplemental Figure 3. Differential gene expression analyses between low- and high-acylsugar–producing accessions.

  • Supplemental Figure 4. Expression levels of 10 selected “housekeeping genes” in different biological groups.

  • Supplemental Figure 5. Biosynthesis of branched-chain and straight-chain acyl molecules.

  • Supplemental Figure 6. Single transcript produced from the locus containing Sopen12g004230-Sopen12g004240.

  • Supplemental Figure 7. Correlation among expression profiles of putative defense genes.

  • Supplemental Figure 8. WGCNA (Langfelder and Horvath, 2008).

  • Supplemental Figure 9. Selected o-methyltransferase genes on chromosome 6 involved in flavonoid metabolism in S. pennellii and S. lycopersicum.

  • Supplemental Figure 10. Enrichment of GO terms associated with WGCNA “green” module genes.

  • Supplemental Figure 11. Molecular phylogenetic analysis of KAS sequences from different solanaceous species.

  • Supplemental Figure 12. VIGS of two candidate genes.

  • Supplemental Table 1. Amount of acylsugars produced by different accessions of S. pennellii.

  • Supplemental Table 2. Read mapping (TopHat2 alignment) results with different accessions of S. pennellii.

  • Supplemental Table 3. Sequence similarity between Cuphea KAS IV/KAS II-like enzymes and Sopen12g004230-Sopen12g004240.

  • Supplemental Table 4. Trichome-enriched expression of candidate genes’ putative ortholog in S. lycopersicum.

  • Supplemental Table 5. Effect of imazapyr treatment on expression levels of ABC transporter genes.

  • Supplemental Table 6. List of primers used in this study.

  • Supplemental Data Set 1. Expression levels of all S. pennellii-annotated genes and results of differential gene expression analysis between low- and high-acylsugar–producing accessions.

  • Supplemental Data Set 2. SRCC among selected genes’ expression profiles in 29 biological samples (different accessions of S. pennellii).

  • Supplemental Data Set 3. Reciprocal best hits (RBH) between S. pennellii and S. lycopersicum annotated protein sequences.

  • Supplemental Data Set 4. Expression of selected genes in response to imazapyr treatment.

  • Supplemental Data Set 5. Results of WGCNA.

  • Supplemental Data Set 6. Summary of the imazapyr treatment study.

  • Supplemental Data Set 7. Summary of statistical tests in the VIGS study.

  • Supplemental File 1. Alignment of KAS III sequences from different solanaceous plants.

  • Supplemental File 2. Alignment of KAS IV/II-like sequences from different solanaceous plants.

  • Supplemental File 3. Phylogenetic relationships among KAS III sequences from different solanaceous plants (mtsx file).

  • Supplemental File 4. Phylogenetic relationships among KAS IV/II-like sequences from different solanaceous plants (mtsx file).

DIVE Curated Terms

The following phenotypic, genotypic, and functional terms are of significance to the work described in this paper:

  • AAE1 Gramene: AT1G20560

  • AAE1 Araport: AT1G20560

  • PMP22 Gramene: AT4G04470

  • PMP22 Araport: AT4G04470

Acknowledgments

We thank Charlie Johnson and the staff of the Texas A&M Genomics and Bioinformatics center for performing Illumina sequencing, Rodolfo Aramayo and Ricardo Perez for technical support on the BioGalaxy server, and Yohannes Rezenom from the Texas A&M Chemistry Mass Spectroscopy Center for his assistance with GC–MS analysis. Most of the sequence analysis was performed at the Texas A&M University High Performance Research Computer, where Michael Dickens provided invaluable assistance. We sincerely appreciate Alan Pepper’s (Department of Biology, Texas A&M University) assistance with multiple aspects of this project, from initial design to critical review of the article. This work was partially supported by the US Department of Agriculture (grant 2011-38821-30891).

AUTHOR CONTRIBUTIONS

S.M. designed experiments, prepared samples, analyzed data, and wrote the article; W.J. designed experiments, prepared samples, and reviewed the article; T.D.M. designed experiments, analyzed data, and wrote the article; all authors read and approved the final article.

Footnotes

  • www.plantcell.org/cgi/doi/10.1105/tpc.19.00552

  • The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors is: Thomas D. McKnight (mcknight{at}bio.tamu.edu).

  • ↵[OPEN] Articles can be viewed without a subscription.

  • Received July 19, 2019.
  • Revised September 30, 2019.
  • Accepted October 16, 2019.
  • Published October 18, 2019.

References

  1. ↵
    1. Alba, J.M.,
    2. Montserrat, M.,
    3. Fernández-Muñoz, R.
    (2009). Resistance to the two-spotted spider mite (Tetranychus urticae) by acylsucroses of wild tomato (Solanum pimpinellifolium) trichomes studied in a recombinant inbred line population. Exp. Appl. Acarol. 47: 35–47.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Anders, S.,
    2. Pyl, P.T.,
    3. Huber, W.
    (2015). HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 31: 166–169.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Andrews, S.
    (2010). FastQC: A quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk_?_/projects/fastqc/. Accessed 2011 October 6.
  4. ↵
    1. Balcke, G.U.,
    2. Bennewitz, S.,
    3. Bergau, N.,
    4. Athmer, B.,
    5. Henning, A.,
    6. Majovsky, P.,
    7. Jiménez-Gómez, J.M.,
    8. Hoehenwarter, W.,
    9. Tissier, A.
    (2017). Multi-omics of tomato glandular trichomes reveals distinct features of central carbon metabolism supporting high productivity of specialized metabolites. Plant Cell 29: 960–983.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Benjamini, Y.,
    2. Hochberg, Y.
    (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57: 289–300.
    OpenUrlCrossRef
  6. ↵
    1. Binder, S.
    (2010). Branched-chain amino acid metabolism in Arabidopsis thaliana. Arabidopsis Book 8: e0137.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Blauth, S.L.,
    2. Steffens, J.C.,
    3. Churchill, G.A.,
    4. Mutschler, M.A.
    (1999). Identification of QTLs controlling acylsugar fatty acid composition in an intraspecific population of Lycopersicon pennellii (Corr.) D’Arcy. Theor. Appl. Genet. 99: 373–381.
    OpenUrlCrossRef
  8. ↵
    1. Bolger, A., et al.
    (2014a). The genome of the stress-tolerant wild tomato species Solanum pennellii. Nat. Genet. 46: 1034–1038.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bolger, A.M.,
    2. Lohse, M.,
    3. Usadel, B.
    (2014b). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30: 2114–2120.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Bonierbale, M.W.,
    2. Plaisted, R.L.,
    3. Pineda, O.,
    4. Tanksley, S.D.
    (1994). QTL analysis of trichome-mediated insect resistance in potato. Theor. Appl. Genet. 87: 973–987.
    OpenUrlPubMed
  11. ↵
    1. Brosius, U.,
    2. Dehmel, T.,
    3. Gärtner, J.
    (2002). Two different targeting signals direct human peroxisomal membrane protein 22 to peroxisomes. J. Biol. Chem. 277: 774–784.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Chortyk, O.T.,
    2. Kays, S.J.,
    3. Teng, Q.
    (1997). Characterization of insecticidal sugar esters of Petunia. J. Agric. Food Chem. 45: 270–275.
    OpenUrlCrossRef
  13. ↵
    1. Chortyk, O.T.,
    2. Severson, R.F.,
    3. Cutler, H.C.,
    4. Sisson, V.A.
    (1993). Antibiotic activities of sugar esters isolated from selected Nicotiana species. Biosci. Biotechnol. Biochem. 57: 1355–1356.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Cona, A.,
    2. Rea, G.,
    3. Angelini, R.,
    4. Federico, R.,
    5. Tavladoraki, P.
    (2006). Functions of amine oxidases in plant development and defence. Trends Plant Sci. 11: 80–88.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Conesa, A.,
    2. Götz, S.,
    3. García-Gómez, J.M.,
    4. Terol, J.,
    5. Talón, M.,
    6. Robles, M.
    (2005). Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Couto, D.,
    2. Zipfel, C.
    (2016). Regulation of pattern recognition receptor signalling in plants. Nat. Rev. Immunol. 16: 537–552.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Dehesh, K.,
    2. Edwards, P.,
    3. Fillatti, J.,
    4. Slabaugh, M.,
    5. Byrne, J.
    (1998). KAS IV: A 3-ketoacyl-ACP synthase from Cuphea sp. is a medium chain specific condensing enzyme. Plant J. 15: 383–390.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Dong, Y.,
    2. Burch-Smith, T.M.,
    3. Liu, Y.,
    4. Mamillapalli, P.,
    5. Dinesh-Kumar, S.P.
    (2007). A ligation-independent cloning tobacco rattle virus vector for high-throughput virus-induced gene silencing identifies roles for NbMADS4-1 and -2 in floral development. Plant Physiol. 145: 1161–1170.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Edgar, R.C.
    (2004). MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32: 1792–1797.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Fan, P.,
    2. Miller, A.M.,
    3. Liu, X.,
    4. Jones, A.D.,
    5. Last, R.L.
    (2017). Evolution of a flipped pathway creates metabolic innovation in tomato trichomes through BAHD enzyme promiscuity. Nat. Commun. 8: 2080.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Fan, P.,
    2. Miller, A.M.,
    3. Schilmiller, A.L.,
    4. Liu, X.,
    5. Ofner, I.,
    6. Jones, A.D.,
    7. Zamir, D.,
    8. Last, R.L.
    (2016). In vitro reconstruction and analysis of evolutionary variation of the tomato acylsucrose metabolic network. Proc. Natl. Acad. Sci. USA 113: E239–E248.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Fernandez-Pozo, N., et al.
    (2015). The Sol Genomics Network (SGN)—from genotype to phenotype to breeding. Nucleic Acids Res. 43: D1036–D1041.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Figueiredo, A.,
    2. Monteiro, F.,
    3. Sebastiana, M.
    (2014). Subtilisin-like proteases in plant-pathogen recognition and immune priming: A perspective. Front. Plant Sci. 5: 739.
    OpenUrl
  24. ↵
    1. Fobes, J.F.,
    2. Mudd, J.B.,
    3. Marsden, M.P.
    (1985). Epicuticular lipid accumulation on the leaves of Lycopersicon pennellii (Corr.) D’Arcy and Lycopersicon esculentum Mill. Plant Physiol. 77: 567–570.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Ghangas, G.S.,
    2. Steffens, J.C.
    (1993). UDPglucose: Fatty acid transglucosylation and transacylation in triacylglucose biosynthesis. Proc. Natl. Acad. Sci. USA 90: 9911–9915.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Halitschke, R.,
    2. Baldwin, I.T.
    (2003). Antisense LOX expression increases herbivore performance by decreasing defense responses and inhibiting growth-related transcriptional reorganization in Nicotiana attenuata. Plant J. 36: 794–807.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hare, J.D.
    (2005). Biological activity of acyl glucose esters from Datura wrightii glandular trichomes against three native insect herbivores. J. Chem. Ecol. 31: 1475–1491.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Hawthorne, D.J.,
    2. Shapiro, J.A.,
    3. Tingey, W.M.,
    4. Mutschler, M.A.
    (1992). Trichome-borne and artificially applied acylsugars of wild tomato deter feeding and oviposition of the leafminer Liriomyza trifolii. Entomol. Exp. Appl. 65: 65–73.
    OpenUrlCrossRef
  29. ↵
    1. Hill, K.,
    2. Rhode, O.
    (1999). Sugar-based surfactants for consumer products and technical applications. Fett-Lipid 101: 25–33.
    OpenUrl
  30. ↵
    1. Hu, Z.,
    2. Mellor, J.,
    3. Wu, J.,
    4. Yamada, T.,
    5. Holloway, D.,
    6. Delisi, C.
    (2005). VisANT: Data-integrating visual framework for biological networks and modules. Nucleic Acids Res. 33: W352–W357.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Jackson, P.A.,
    2. Galinha, C.I.,
    3. Pereira, C.S.,
    4. Fortunato, A.,
    5. Soares, N.C.,
    6. Amâncio, S.B.,
    7. Pinto Ricardo, C.P.
    (2001). Rapid deposition of extensin during the elicitation of grapevine callus cultures is specifically catalyzed by a 40-kilodalton peroxidase. Plant Physiol. 127: 1065–1076.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Jones, D.T.,
    2. Taylor, W.R.,
    3. Thornton, J.M.
    (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8: 275–282.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Juvik, J.A.,
    2. Shapiro, J.A.,
    3. Young, T.E.,
    4. Mutschler, M.A.
    (1994). Acylglucoses from wild tomatoes alter behavior and reduce growth and survival of Helicoverpa zea and Spodoptera exigua (Lepidoptera, Noctuidae). J. Econ. Entomol. 87: 482–492.
    OpenUrlCrossRef
  34. ↵
    1. Kim, D.,
    2. Pertea, G.,
    3. Trapnell, C.,
    4. Pimentel, H.,
    5. Kelley, R.,
    6. Salzberg, S.L.
    (2013). TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14: R36.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Kim, J.,
    2. Matsuba, Y.,
    3. Ning, J.,
    4. Schilmiller, A.L.,
    5. Hammar, D.,
    6. Jones, A.D.,
    7. Pichersky, E.,
    8. Last, R.L.
    (2014). Analysis of natural and induced variation in tomato glandular trichome flavonoids identifies a gene not present in the reference genome. Plant Cell 26: 3272–3285.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Kroumova, A.B.,
    2. Zaitlin, D.,
    3. Wagner, G.J.
    (2016). Natural variability in acyl moieties of sugar esters produced by certain tobacco and other Solanaceae species. Phytochemistry 130: 218–227.
    OpenUrl
  37. ↵
    1. Kuai, J.P.,
    2. Ghangas, G.S.,
    3. Steffens, J.C.
    (1997). Regulation of triacylglucose fatty acid composition (uridine diphosphate glucose:fatty acid glucosyltransferases with overlapping chain-length specificity). Plant Physiol. 115: 1581–1587.
    OpenUrlAbstract
  38. ↵
    1. Kumar, S.,
    2. Stecher, G.,
    3. Tamura, K.
    (2016). MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for bigger datasets. Mol. Biol. Evol. 33: 1870–1874.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Langfelder, P.,
    2. Horvath, S.
    (2008). WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics 9: 559.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Lawson, D.M.,
    2. Lunde, C.F.,
    3. Mutschler, M.A.
    (1997). Marker-assisted transfer of acylsugar-mediated pest resistance from the wild tomato, Lycopersicon pennellii, to the cultivated tomato, Lycopersicon esculentum. Mol. Breed. 3: 307–317.
    OpenUrlCrossRef
  41. ↵
    1. Leckie, B.M.,
    2. De Jong, D.M.,
    3. Mutschler, M.A.
    (2012). Quantitative trait loci increasing acylsugars in tomato breeding lines and their impacts on silverleaf whiteflies. Mol. Breed. 30: 1621–1634.
    OpenUrlCrossRef
  42. ↵
    1. Leong, B.J.,
    2. Lybrand, D.B.,
    3. Lou, Y.-R.,
    4. Fan, P.,
    5. Schilmiller, A.L.,
    6. Last, R.L.
    (2019). Evolution of metabolic novelty: A trichome-expressed invertase creates specialized metabolic diversity in wild tomato. Sci. Adv. 5: eaaw3754.
    OpenUrlFREE Full Text
  43. ↵
    1. Li, A.X.,
    2. Eannetta, N.,
    3. Ghangas, G.S.,
    4. Steffens, J.C.
    (1999). Glucose polyester biosynthesis. Purification and characterization of a glucose acyltransferase. Plant Physiol. 121: 453–460.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Li, A.X.,
    2. Steffens, J.C.
    (2000). An acyltransferase catalyzing the formation of diacylglucose is a serine carboxypeptidase-like protein. Proc. Natl. Acad. Sci. USA 97: 6902–6907.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Liedl, B.E.,
    2. Lawson, D.M.,
    3. White, K.K.,
    4. Shapiro, J.A.,
    5. Cohen, D.E.,
    6. Carson, W.G.,
    7. Trumble, J.T.,
    8. Mutschler, M.A.
    (1995). Acylsugars of wild tomato Lycopersicon pennellii alters settling and reduces oviposition of Bemisia argentifolii (Homoptera, Aleyrodidae). J. Econ. Entomol. 88: 742–748.
    OpenUrlCrossRef
  46. ↵
    1. Luu, V.T.,
    2. Weinhold, A.,
    3. Ullah, C.,
    4. Dressel, S.,
    5. Schoettner, M.,
    6. Gase, K.,
    7. Gaquerel, E.,
    8. Xu, S.,
    9. Baldwin, I.T.
    (2017). O-acyl sugars protect a wild tobacco from both native fungal pathogens and a specialist herbivore. Plant Physiol. 174: 370–386.
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Maloney, G.S.,
    2. Kochevenko, A.,
    3. Tieman, D.M.,
    4. Tohge, T.,
    5. Krieger, U.,
    6. Zamir, D.,
    7. Taylor, M.G.,
    8. Fernie, A.R.,
    9. Klee, H.J.
    (2010). Characterization of the branched-chain amino acid aminotransferase enzyme family in tomato. Plant Physiol. 153: 925–936.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Manabe, Y.,
    2. Tinker, N.,
    3. Colville, A.,
    4. Miki, B.
    (2007). CSR1, the sole target of imidazolinone herbicide in Arabidopsis thaliana. Plant Cell Physiol. 48: 1340–1358.
    OpenUrlCrossRefPubMed
  49. ↵
    1. McDowell, E.T., et al.
    (2011). Comparative functional genomic analysis of Solanum glandular trichome types. Plant Physiol. 155: 524–539.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Moghe, G.D.,
    2. Leong, B.J.,
    3. Hurney, S.M.,
    4. Daniel Jones, A.,
    5. Last, R.L.
    (2017). Evolutionary routes to biochemical innovation revealed by integrative analysis of a plant-defense related specialized metabolic pathway. eLife 6.
  51. ↵
    1. Mueller, S.,
    2. Hilbert, B.,
    3. Dueckershoff, K.,
    4. Roitsch, T.,
    5. Krischke, M.,
    6. Mueller, M.J.,
    7. Berger, S.
    (2008). General detoxification and stress responses are mediated by oxidized lipids through TGA transcription factors in Arabidopsis. Plant Cell 20: 768–785.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Niebel, A.,
    2. De Almeida Engler, J.,
    3. Tire, C.,
    4. Engler, G.,
    5. Van Montagu, M.,
    6. Gheysen, G.
    (1993). Induction patterns of an extensin gene in tobacco upon nematode infection. Plant Cell 5: 1697–1710.
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Ning, J.,
    2. Moghe, G.D.,
    3. Leong, B.,
    4. Kim, J.,
    5. Ofner, I.,
    6. Wang, Z.,
    7. Adams, C.,
    8. Jones, A.D.,
    9. Zamir, D.,
    10. Last, R.L.
    (2015). A feedback-insensitive isopropylmalate synthase affects acylsugar composition in cultivated and wild tomato. Plant Physiol. 169: 1821–1835.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Noguchi, A.,
    2. Sasaki, N.,
    3. Nakao, M.,
    4. Fukami, H.,
    5. Takahashi, S.,
    6. Nishino, T.,
    7. Nakayama, T.
    (2008). cDNA cloning of glycosyltransferases from Chinese wolfberry (Lycium barbarum L.) fruits and enzymatic synthesis of a catechin glucoside using a recombinant enzyme (UGT73A10). J. Mol. Catal. B Enzym. 55: 84–92.
    OpenUrlCrossRef
  55. ↵
    1. Pérez-Castorena, A.L.,
    2. Martínez, M.,
    3. Maldonado, E.
    (2010). Labdanes and sucrose esters from Physalis sordida. J. Nat. Prod. 73: 1271–1276.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Prost, I., et al.
    (2005). Evaluation of the antimicrobial activities of plant oxylipins supports their involvement in defense against pathogens. Plant Physiol. 139: 1902–1913.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Puterka, G.J.,
    2. Farone, W.,
    3. Palmer, T.,
    4. Barrington, A.
    (2003). Structure-function relationships affecting the insecticidal and miticidal activity of sugar esters. J. Econ. Entomol. 96: 636–644.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Reumann, S.,
    2. Quan, S.,
    3. Aung, K.,
    4. Yang, P.,
    5. Manandhar-Shrestha, K.,
    6. Holbrook, D.,
    7. Linka, N.,
    8. Switzenberg, R.,
    9. Wilkerson, C.G.,
    10. Weber, A.P.,
    11. Olsen, L.J.,
    12. Hu, J.
    (2009). In-depth proteome analysis of Arabidopsis leaf peroxisomes combined with in vivo subcellular targeting verification indicates novel metabolic and regulatory functions of peroxisomes. Plant Physiol. 150: 125–143.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Robinson, M.D.,
    2. McCarthy, D.J.,
    3. Smyth, G.K.
    (2010). edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Schilmiller, A.L.,
    2. Charbonneau, A.L.,
    3. Last, R.L.
    (2012). Identification of a BAHD acetyltransferase that produces protective acyl sugars in tomato trichomes. Proc. Natl. Acad. Sci. USA 109: 16377–16382.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Schilmiller, A.L.,
    2. Gilgallon, K.,
    3. Ghosh, B.,
    4. Jones, A.D.,
    5. Last, R.L.
    (2016). Acylsugar acylhydrolases: Carboxylesterase-catalyzed hydrolysis of acylsugars in tomato trichomes. Plant Physiol. 170: 1331–1344.
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Schilmiller, A.L.,
    2. Moghe, G.D.,
    3. Fan, P.,
    4. Ghosh, B.,
    5. Ning, J.,
    6. Jones, A.D.,
    7. Last, R.L.
    (2015). Functionally divergent alleles and duplicated Loci encoding an acyltransferase contribute to acylsugar metabolite diversity in Solanum trichomes. Plant Cell 27: 1002–1017.
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Schmidt, A.,
    2. Li, C.,
    3. Jones, A.D.,
    4. Pichersky, E.
    (2012). Characterization of a flavonol 3-o-methyltransferase in the trichomes of the wild tomato species Solanum habrochaites. Planta 236: 839–849.
    OpenUrlCrossRefPubMed
  64. ↵
    1. Schütt, B.S.,
    2. Abbadi, A.,
    3. Loddenkötter, B.,
    4. Brummel, M.,
    5. Spener, F.
    (2002). Beta-ketoacyl-acyl carrier protein synthase IV: A key enzyme for regulation of medium-chain fatty acid synthesis in Cuphea lanceolata seeds. Planta 215: 847–854.
    OpenUrlPubMed
  65. ↵
    1. Shapiro, J.A.,
    2. Steffens, J.C.,
    3. Mutschler, M.A.
    (1994). Acylsugars of the wild tomato Lycopersicon pennellii in relation to geographic distribution of the species. Biochem. Syst. Ecol. 22: 545–561.
    OpenUrlCrossRef
  66. ↵
    1. Shockey, J.M.,
    2. Fulda, M.S.,
    3. Browse, J.
    (2003). Arabidopsis contains a large superfamily of acyl-activating enzymes. Phylogenetic and biochemical analysis reveals a new class of acyl-coenzyme a synthetases. Plant Physiol. 132: 1065–1076.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Singh, V.,
    2. Roy, S.,
    3. Giri, M.K.,
    4. Chaturvedi, R.,
    5. Chowdhury, Z.,
    6. Shah, J.,
    7. Nandi, A.K.
    (2013). Arabidopsis thaliana FLOWERING LOCUS D is required for systemic acquired resistance. Mol. Plant Microbe Interact. 26: 1079–1088.
    OpenUrlCrossRefPubMed
  68. ↵
    1. Slabaugh, M.B.,
    2. Leonard, J.M.,
    3. Knapp, S.J.
    (1998). Condensing enzymes from Cuphea wrightii associated with medium chain fatty acid biosynthesis. Plant J. 13: 611–620.
    OpenUrlCrossRefPubMed
  69. ↵
    1. Slocombe, S.P.,
    2. Schauvinhold, I.,
    3. McQuinn, R.P.,
    4. Besser, K.,
    5. Welsby, N.A.,
    6. Harper, A.,
    7. Aziz, N.,
    8. Li, Y.,
    9. Larson, T.R.,
    10. Giovannoni, J.,
    11. Dixon, R.A.,
    12. Broun, P.
    (2008). Transcriptomic and reverse genetic analyses of branched-chain fatty acid and acyl sugar production in Solanum pennellii and Nicotiana benthamiana. Plant Physiol. 148: 1830–1846.
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. R Core Team.
    (2014). R: A Language and Environment for Statistical Computing. (Vienna, Austria: R Foundation for Statistical Computing).
  71. ↵
    1. Thorvaldsdóttir, H.,
    2. Robinson, J.T.,
    3. Mesirov, J.P.
    (2013). Integrative Genomics Viewer (IGV): High-performance genomics data visualization and exploration. Brief. Bioinform. 14: 178–192.
    OpenUrlCrossRefPubMed
  72. ↵
    1. Tomato Genome Consortium.
    (2012). The tomato genome sequence provides insights into fleshy fruit evolution. Nature 485: 635–641.
    OpenUrlCrossRefPubMed
  73. ↵
    1. Walters, D.S.,
    2. Steffens, J.C.
    (1990). Branched chain amino acid metabolism in the biosynthesis of Lycopersicon pennellii glucose esters. Plant Physiol. 93: 1544–1551.
    OpenUrlAbstract/FREE Full Text
  74. ↵
    1. Zhao, J.,
    2. Williams, C.C.,
    3. Last, R.L.
    (1998). Induction of Arabidopsis tryptophan pathway enzymes and camalexin by amino acid starvation, oxidative stress, and an abiotic elicitor. Plant Cell 10: 359–370.
    OpenUrlAbstract/FREE Full Text
  75. ↵
    1. Zhu, M., et al.
    (2017). The intracellular immune receptor Sw-5b confers broad-spectrum resistance to tospoviruses through recognition of a conserved 21-amino acid viral effector epitope. Plant Cell 29: 2214–2232.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top

Table of Contents

Print
Download PDF
Email Article

Thank you for your interest in spreading the word on Plant Cell.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Candidate Gene Networks for Acylsugar Metabolism and Plant Defense in Wild Tomato Solanum pennellii
(Your Name) has sent you a message from Plant Cell
(Your Name) thought you would like to see the Plant Cell web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Candidate Gene Networks for Acylsugar Metabolism and Plant Defense in Wild Tomato Solanum pennellii
Sabyasachi Mandal, Wangming Ji, Thomas D. McKnight
The Plant Cell Jan 2020, 32 (1) 81-99; DOI: 10.1105/tpc.19.00552

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Candidate Gene Networks for Acylsugar Metabolism and Plant Defense in Wild Tomato Solanum pennellii
Sabyasachi Mandal, Wangming Ji, Thomas D. McKnight
The Plant Cell Jan 2020, 32 (1) 81-99; DOI: 10.1105/tpc.19.00552
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Extras

  • Peer review report

Jump to section

  • Article
    • Abstract
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • METHODS
    • DIVE Curated Terms
    • Acknowledgments
    • AUTHOR CONTRIBUTIONS
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

In this issue

The Plant Cell: 32 (1)
The Plant Cell
Vol. 32, Issue 1
Jan 2020
  • Table of Contents
  • Table of Contents (PDF)
  • Cover (PDF)
  • About the Cover
  • Index by author
View this article with LENS

More in this TOC Section

  • ARADEEPOPSIS, an Automated Workflow for Top-View Plant Phenomics using Semantic Segmentation of Leaf States
  • Alternative Crassulacean Acid Metabolism Modes Provide Environment-Specific Water-Saving Benefits in a Leaf Metabolic Model
  • Nonsense-Mediated RNA Decay Factor UPF1 Is Critical for Posttranscriptional and Translational Gene Regulation in Arabidopsis
Show more LARGE-SCALE BIOLOGY ARTICLES

Similar Articles

Our Content

  • Home
  • Current Issue
  • Plant Cell Preview
  • Archive
  • Teaching Tools in Plant Biology
  • Plant Physiology
  • Plant Direct
  • Plantae
  • ASPB

For Authors

  • Instructions
  • Submit a Manuscript
  • Editorial Board and Staff
  • Policies
  • Recognizing our Authors

For Reviewers

  • Instructions
  • Peer Review Reports
  • Journal Miles
  • Transfer of reviews to Plant Direct
  • Policies

Other Services

  • Permissions
  • Librarian resources
  • Advertise in our journals
  • Alerts
  • RSS Feeds
  • Contact Us

Copyright © 2021 by The American Society of Plant Biologists

Powered by HighWire