Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2021.657156Plant ScienceOriginal ResearchAn R2R3-MYB Transcription Factor Positively Regulates the Glandular Secretory Trichome Initiation in Artemisia annua L.QinWeiXieLihuiLiYongpengLiuHangFuXueqingChenTiantianHassaniDanialLiLingSunXiaofenTangKexuan*Joint International Research Laboratory of Metabolic and Developmental Sciences, Key Laboratory of Urban Agriculture (South) Ministry of Agriculture, Plant Biotechnology Research Center, Fudan-SJTU-Nottingham Plant Biotechnology R&D Center, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
Edited by: Kazufumi Yazaki, Kyoto University, Japan
Reviewed by: Lei Zhang, Second Military Medical University, China; Takashi Aoyama, Kyoto University, Japan
*Correspondence: Kexuan Tang, kxtang@sjtu.edu.cn
This article was submitted to Plant Metabolism and Chemodiversity, a section of the journal Frontiers in Plant Science
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Artemisia annua L. is known for its specific product “artemisinin” which is an active ingredient for curing malaria. Artemisinin is secreted and accumulated in the glandular secretory trichomes (GSTs) on A. annua leaves. Earlier studies have shown that increasing GST density is effective in increasing artemisinin content. However, the mechanism of GST initiation is not fully understood. To this end, we isolated and characterized an R2R3-MYB gene, AaMYB17, which is expressed specifically in the GSTs of shoot tips. Overexpression of AaMYB17 in A. annua increased GST density and enhanced the artemisinin content, whereas RNA interference of AaMYB17 resulted in the reduction of GST density and artemisinin content. Additionally, neither overexpression lines nor RNAi lines showed an abnormal phenotype in plant growth and the morphology of GSTs. Our study demonstrates that AaMYB17 is a positive regulator of GSTs’ initiation, without influencing the trichome morphology.
Malaria is a mosquito-borne infectious disease caused by the Plasmodium species which targets human red blood cells (Wilson et al., 2011). According to the World Health Organization (WHO), malaria is responsible for 228 million cases, including 405,000 deaths worldwide in 2018 (World Health Organization [WHO], 2019). Artemisinin-based combination therapy (ACT) is considered the most efficient treatment to cure malaria (World Health Organization [WHO], 2019). Artemisinin is a bioactive compound, synthesized and accumulated in the glandular trichome of a traditional Chinese herb, Artemisia annua L. (Graham et al., 2010; Tan et al., 2015).
Trichomes are unicellular or multicellular structures derived from aerial epidermal cells (Duke and Paul, 2013). Their morphological diversity, different number of cells, and their potential in secondary metabolism distinguish them from each other (Werker, 2000; Serna and Martin, 2006). Depending on secondary metabolism capability, they are divided into glandular trichomes and non-glandular trichomes (Huchelmann et al., 2017). The powerful biosynthetic secreting and accumulating capabilities of glandular trichomes for the production of plants’ secondary metabolites including terpenoids, flavonoid, alkaloids, polysaccharides, polyphenols, and phenylpropanoids have made them an exciting research topic in recent years for plant metabolic engineering strategies (Wagner et al., 2004; Fridman et al., 2005; Gershenzon and Dudareva, 2007; Xie et al., 2008; Tattini et al., 2010; Weinhold and Baldwin, 2011). Many efforts have been dedicated to promoting the accumulation of secondary metabolites in plants, and it is feasible to improve secondary metabolic productivity by increasing glandular trichome density (Tissier, 2012). Unlike glandular trichomes, non-glandular trichomes are not efficient in secondary metabolism (Soetaert et al., 2013; Huchelmann et al., 2017). Besides, both non-glandular and glandular trichomes play an essential role in the defense against abiotic and biotic stress, such as UV light, water absorption, salt stress, and herbivores (Mauricio and Rausher, 1997; Serna and Martin, 2006). Therefore, systematic investigations on the molecular basis of trichome initiation are needed.
There are two types of trichomes on the surface of A. annua leaves, glandular secretory trichomes (GSTs) and T-shaped trichomes (TSTs), both have multicellular structures (Duke and Paul, 2013; Xiao et al., 2016). Artemisinin is specifically synthesized and accumulated in GSTs. Thus, promoting GST initiation is a potential strategy for increasing the artemisinin content in A. annua. Even though the mechanism of unicellular trichome initiation in Arabidopsis has been studied thoroughly, the mechanism of multicellular trichome initiation might vary and needs further investigation (Payne et al., 1999; Martin, 2006; Maes et al., 2011; Yang et al., 2011).
In plants, transcription factors regulate a variety of biological processes (Mitsuda and Ohme-Takagi, 2009). Many TF families including the well-known MYB family were reported to regulate trichome initiation. MYBs are a family of proteins that contain the conserved MYB DNA-binding domain (Romero et al., 1998). MYB proteins are classified into four classes according to the number of adjacent repeats (one, two, three, or four): 4R-MYB, R1R2R3-type MYB (3R-MYB), 1R-MYB/MYB-related, and R2R3-MYB (Rosinski and Atchley, 1998; Jin and Martin, 1999; Stracke et al., 2001; Dubos et al., 2010). The R2R3-MYBs are the largest subfamily of MYB TFs containing an N terminal DNA-binding domain and a C terminus activation or repression domain (Martin and Paz-Ares, 1997; Dubos et al., 2010). In recent studies, the R2R3-MYBs have been proven to regulate primary and secondary metabolism, plant development, cell fate and identity, response to biotic and abiotic stresses, and light and hormone signaling pathways (Martin and Paz-Ares, 1997; Dubos et al., 2010; Brockington et al., 2013). Based on the conserved amino-acid sequence motifs present at the carboxy terminal to the MYB domain, R2R3-type MYB factors have been categorized into 22 subgroups (Martin and Paz-Ares, 1997; Kranz et al., 2010). Former studies have reported that subgroup 9 of MYB TFs, including MIXTA or MIXTA-like, are essential in plant cellular differentiation, especially in epidermal cells differentiation (Brockington et al., 2013). For instance, the first MIXTA characterized is the snapdragon Antirrhinum majus MIXTA (AmMIXTA) which controls formation of the conical shape of petal epidermal cells (Noda et al., 1994). Moreover, some other MIXTA genes in A. majus, such as AmMYBML1, are also described as having a similar function (Oshima et al., 2013). AtMYB106 and AtMYB16 which also belong to subgroup 9 regulate trichome branch formation and epidermal cell morphology in Arabidopsis (Jakoby et al., 2008; Oshima et al., 2013). Furthermore, AaMIXTA1, belonging to subgroup 9, positively regulates glandular trichome initiation and cuticle biosynthesis in A. annua (Shi et al., 2018). In Populus euphratica, overexpression of an AtMYB106 homologous gene, PtaMYB186, can increase trichome density (Plett et al., 2010). Nevertheless, the molecular mechanism of MYB-involved trichome initiation needs further study.
In this study, an R2R3-MYB gene, AaMYB17, was isolated from 13 MYBs which are highly and specifically expressed in the meristem. Overexpression of AaMYB17 increased the GST density, whereas RNA interference of AaMYB17 resulted in the decrease of GST density in A. annua. Taken together, we identified a new TF which could promote artemisinin accumulation by positively regulating the GST density in A. annua leaves.
Materials and MethodsPlant Materials and Growth Conditions
The A. annua we used in this study was “Huhao 1,” which originated from Chongqing, China and has been planted and selected in Shanghai, China for several years (Shen et al., 2016). The plants were grown under a 16 h/8 h light/dark photoperiod and 65% relative humidity at 25 ± 2°C. Nicotiana benthamiana, used for transient transformation, was grown at 25 ± 2°C under a 16-h light photoperiod.
Phylogenetic Tree and Amino Acid Sequence Alignment
Myeloblastosis (MYB) TFs highly and specifically expressed in the meristem were identified from our unpublished transcriptome data (Shi et al., 2018). The phylogenetic tree was analyzed by Mega5 software (Tamura et al., 2011). Protein sequences of MIXTA or MIXTA-like TFs from other species were downloaded from the National Center for Biotechnology Information (NCBI). Amino acid sequence alignment of AaMYB17 and MIXTA/MIXTA-like proteins from other species was performed with Genedoc (Nicholas, 1997).
RNA Isolation and Reverse Transcription
RNA of different tissues and leaves from different phyllotaxis of A. annua was extracted using the RNA prep Pure Plant Kit following the manufacturer’s instructions (Tiangen, Beijing, China). Different tissues (flower bud, leaves, flower, stem, and root) of the wild-type plants were collected from 6-month-old A. annua grown in a glasshouse. In addition, the leaf samples were gathered from leaf 0 (meristem), leaf 1 (first leaf below meristem), leaf 2, leaf 3, leaf 4, leaf 5, leaf 9, and leaf 16, counting from the apical top of the main stem (Yan et al., 2016). RNA samples were reverse transcribed into cDNA using the PrimeScript II RT Master Mix (Takara, Dalian, China).
Quantitative Real-Time PCR (qRT-PCR)
Quantitative real-time PCR was performed on a Roche LightCycler 96 real-time PCR machine (Roche, Basel, Switzerland) using SuperReal PreMix Plus SYBR-Green (Tiangen Biotech, China). qRT-PCR was conducted as described previously (Shen et al., 2016), and the relative expression levels were calculated as described previously (Livak and Schmittgen, 2002). The experiments were performed using three biological replicates. All the primers used in the qRT-PCR are listed in Supplementary Table 1.
Transformation of <italic>A. annua</italic>
The full-length cDNA sequence of AaMYB17 was amplified using the cDNA of the A. annua meristem through PCR using KOD plus DNA polymerase (Toyobo, Osaka, Japan). It was further cloned into a pHB vector under a double CaMV35S promoter to generate pHB-CaMV35S:AaMYB17-YFP:NOS with the YFP fused to the C-terminal of AaMYB17 (Mao et al., 2005). A 269 bp AaMYB17 fragment was recombined into the phellsgate12 vector via a gateway LR recombination reaction (Invitrogen) to construct the AaMYB17-RNAi vector. A 2185 bp AaMYB17 promoter fragment was cloned for the construction of pCAMBIA 1391Z-PMYB17. These constructs were introduced into Agrobacterium tumefaciens strain EHA105, following Agrobacterium-mediated transformation of A. annua as described previously (Shen et al., 2016).
Subcellular Localization of AaMYB17
The full-length ORF of the AaMYB17 gene without the terminator codon was inserted into the pHB-YFP expression vector under the CaMV35S promoter to form a pHB-AaMYB17–YFP fusion protein. Then the plasmid and p19 protein were introduced into the A. tumefaciens strain GV3101 for N. benthamiana leaf transient expression (Sparkes et al., 2006). The fluorescent signals were observed 60–72 h after infiltration using a TCS SP5-II confocal laser microscopy (Leica Microsystems, Wetzlar, Germany). Three biological repeats were performed to verify these results.
GUS Staining Assay
A 2185 bp AaMYB17 promoter fragment was cloned using KODFX (Toyobo, Japan) and inserted into the pCAMBIA1391Z vector which carries the GUS gene. This construction was transformed into A. annua. The GUS assay of the transgenic plants (T1) was performed as previously described (Jefferson et al., 1987). Transgenic plants were stained in a GUS staining solution and incubated at 37°C in the dark overnight. After GUS staining, ethyl alcohol was used to remove chlorophyll.
Glandular Trichome Density Counting
The mature leaves (leaf 9, the ninth leaf below the meristem) of A. annua plants (T0) grown in the glasshouse were selected to count the density of glandular trichomes. Each leaf was imaged by a ×5 objective using fluorescence microscopy (Olympus, Tokyo, Japan). The ImageJ program1 was used to measure the leaf area and the number of glandular trichomes, as previously described (Cheng et al., 2014). Three different leaves of each independent plant at the same position were selected to count the trichome numbers.
Scanning Electron Microscopy (SEM)
Mature leaves (leaf 9, the ninth leaf below the meristem) of A. annua plants (T0) grown in the glasshouse were selected and treated following the previously reported method (Singh et al., 2016). Leaves were imaged with a Hitachi (Hitachi Ltd., Tokyo, Japan) S-3400N scanning electron microscope.
Artemisinin Content Measurement
Leaves of 5-month-old A. annua plants (T0) grown in the glasshouse were gathered and dried in 50°C for 24 h to measure the artemisinin content. The dried leaves were ground into powder samples, and 1 g of each sample was extracted by methanol and treated ultrasonically twice (55 HZ, 30 min). The artemisinin content was measured by the Waters Alliance 2695 HPLC system (Milford, MA, United States) using high-performance liquid chromatography (HPLC) as described previously (Shen et al., 2016). Three biological repeats were measured for each sample.
Dual-Luciferase (Dual-LUC) Assay
PHB-AaMYB17 was transformed into A. tumefaciens strain GV3101 to act as an effector. The empty pHB vector was used as a control. The promoters of ADS, CYP71AV1, DBR2, ALDH1, AaKCS5, AaCER1, AaCYP77A1, AaCYP86A1, and AaABCG12 were cloned into the vector pGREEN II0800 to act as the reporters. GV3101 strains harboring the indicated combinations of effectors and reporters were co-infiltrated into 6-week-old N. benthamiana leaves. The leaves of N. benthamiana were gathered after 24 h of cultivation in dark conditions and 24 h of cultivation in light conditions. After being quick-frozen in liquid nitrogen, the samples were subjected to firefly LUC and REN activities analysis using the Dual-Luciferase® Reporter Assay System (Promega, United States). Three biological repeats were measured for each sample.
ResultsIdentification and Characterization of <italic>AaMYB17</italic>
To further study the function of MYB proteins in trichome initiation, we selected 13 MYBs which are highly and specifically expressed in the meristem where trichome initiation occurs using the GST transcriptome database (Shi et al., 2018) and other published databases (Graham et al., 2010). Phylogenetic analysis of these 13 candidate MYBs revealed that contig133232 clustered with the MIXTA and MIXTA-like TFs from the other species (Figure 1A). Moreover, contig133232 was found to share a highly conserved R2R3-MYB domain by means of an amino acid alignment (Figure 1B). The characteristic of “HMAQWESARxEAEAxLxMDS” demonstrated that contig133232 belonged to subgroup 9 of R2R3-MYBs (Stracke et al., 2001; Brockington et al., 2013). We further named contig133232 as AaMYB17 which is homologous with AtMYB17 in Arabidopsis thaliana and selected it as a candidate for further study.
Identification and characterization of AaMYB17. (A) Phylogenetic tree of MYB proteins expressed in young leaf trichomes, flower bud trichomes, and shoots of A. annua, and some MYBs involved in organ development and cell differentiation from other species. The subgroup 9 MYBs described previously are highlighted in red, and the candidate contig is starred. Sequences were downloaded from the TAIR and NCBI databases. (B) Amino acid sequence alignment of AaMYB17 and MIXTA/MIXTA-like proteins from other species were performed with Genedoc (Nicholas, 1997). The conserved MYB domain and subgroup 9 motif are represented by a line above the sequence. The candidate is starred.
Expression Profile of <italic>AaMYB17</italic>
qRT-PCR was performed using cDNA from different tissues to analyze the transcript level of AaMYB17 and investigate its expression pattern. The results indicated that AaMYB17 was highly expressed in the young leaf, bud, flower, and especially in the shoot, where GSTs were abundant (Figure 2A). The expression of AaMYB17 in different leaves followed a rapid descending pattern associated with leaf aging (Figure 2B). To further explore the tissue-specific expression pattern of AaMYB17, we cloned a 2186-bp sequence of the AaMYB17 promoter to generate a pAaMYB17-GUS plasmid, and transformed it into A. annua. GUS staining of the transgenic plants showed that AaMYB17 was specifically expressed in shoot tip GSTs (Figures 2C–E and Supplementary Figure 1).
Expression profile of AaMYB17. (A) Relative expression of AaMYB17 in young leaf, old leaf, shoot, bud, flower, root, and stem of A. annua by quantitative real-time polymerase chain reaction. (B) Relative expression of AaMYB17 in A. annua leaves of different developmental ages. β-actin was used as an internal control in A and B. Error bars represent the standard deviation (n = 3). (C–E) GUS staining of transgenic A. annua plants was observed in the GSTs of shoot tip. Bars: (C) 100 μm; (D) 50 μm, and (E) 20 μm.
Subcellular Localization of AaMYB17
To investigate the subcellular localization of AaMYB17, a yellow fluorescent protein (YFP) was fused to the N-terminus of AaMYB17. The YFP fluorescence of 35s: AaMYB17-YFP was observed in the nucleus of N. benthamiana epidermal cells (Figure 3A). However, the fluorescence of control was observed from the whole cell of N. benthamiana (Figure 3B). These results demonstrated that AaMYB17 was localized in the nucleus, which is also consistent with its role as a TF.
The subcellular localization of AaMYB17 in leaves of N. benthamiana. (A) Localization of AaMYB17 protein fused with YFP. (B) Localization of empty vector. Bars = 40 μm.
AaMYB17 Positively Regulates GST Initiation in <italic>A. annua</italic>
To further clarify the molecular function of AaMYB17 in A. annua, we generated AaMYB17-overexpressed transgenic plants driven by the cauliflower mosaic virus (CaMV) 35S promoter. Quantitative RT-PCR was performed and revealed that AaMYB17 was significantly overexpressed in transgenic plants (Figure 4D). In three independent overexpression lines, the number of GSTs on the adaxial leaf side was increased 1.3–1.6-fold compared with the check control (Figures 4A,B,F). As expected, the artemisinin content of the OE-AaMYB17 lines was increased from 8 to 15 mg g–1 DW compared to the wild-type plants (Figure 4G). On the other hand, AaMYB17-suppressed RNAi lines were generated under the CaMV 35S promoter and AaMYB17 expression decreased significantly (Figure 4E). In these four independent RNAi lines, the number of GSTs on the adaxial leaf side was decreased 1.5- to 2.3-fold compared with check control (Figures 4A,C,F). The artemisinin content of the AaMYB17-RNAi lines was decreased from 8 to 6 mg g–1 DW compared to the control (Figure 4G). Meanwhile, the shape of GSTs and TSTs was observed using SEM, and there was no difference between transgenic plants and wild-type plants (Figures 5A–I). Furthermore, neither the overexpression nor the RNAi of AaMYB17 affected the growth of transgenic plants (Supplementary Figure 2).
The expression of AaMYB17 affects A. annua glandular trichome initiation. The glandular trichomes on the adaxial side of mature leaves derived from (A) wild-type (CK) plants, (B)OE-AaMYB17 transgenic A. annua plants, and (C)AaMYB17-RNAi transgenic A. annua plants (bars, 200 μm). Expression level of AaMYB17 in the A. annua leaves in (D)OE-AaMYB17 transgenic A. annua plants and (E)AaMYB17-RNAi transgenic A. annua plants. (F) Glandular trichomes density of mature leaves derived from CK and AaMYB17 transgenic plants. (G) High-performance liquid chromatography (HPLC) analysis of the artemisinin content (mg g–1 leaf DW) in transgenic A. annua plants. Data are given as means ± SD (n = 3) (**P < 0.01; Dunnett’s test).
Scanning electron microscope (SEM) analysis of the surface of leaves. (A) GSTs and T-shape trichomes (TSTs) on the adaxial sides of mature leaves derived from wild-type plants, (B)OE-AaMYB17 transgenic A. annua plants, and (C)AaMYB17-RNAi transgenic A. annua plants. (D) The morphologies of GST on the leaf of wild-type plants, (E)OE-AaMYB17 transgenic A. annua plants, and (F)AaMYB17-RNAi transgenic A. annua plants. (G) The morphologies of TST on the leaf of wild-type plants, (H)OE-AaMYB17 transgenic A. annua plants, and (I)AaMYB17-RNAi transgenic A. annua plants. Bars: (A–C) 1.00 mm; (D–F) 50 μm; and (G–I) 100 μm.
The dual-LUC assay was performed to detect whether AaMYB17 activated the expression of the key enzyme genes in the artemisinin biosynthesis pathway, such as ADS, CYP71AV1, DBR2, and ALDH1. It showed that AaMYB17 had no significant influence on them (Supplementary Figure 3). These results indicated that AaMYB17 regulated trichome initiation positively.
Discussion
Former studies revealed that MIXTA/MIXTA-like TFs play an important role in the regulation of trichome initiation or cell development. A. majus MIXTA (AmMIXTA) as the first MIXTA gene to be reported regulates the development and initiation of the conical cell shape of the petal epidermis (Noda et al., 1994; Glover et al., 1998). AtMYB16 and AtMYB106 of A. thaliana are also reported to regulate trichome development (Baumann et al., 2007; Jakoby et al., 2008). In Gossypium hirsuta, GhMYB25 and GhMYB25-like genes regulate early cotton fiber and trichome development (Machado et al., 2009; Walford et al., 2011). These studies demonstrate that MIXTA/MIXTA-like genes play very crucial roles in trichome development and initiation. In this study, we characterized MIXTA TF AaMYB17. By generating overexpression transgenic A. annua plants, we found that AaMYB17 could increase GST density and artemisinin content compared to the control (Figures 4A,B,F,G). On the contrary, GST density and artemisinin content decreased in RNAi transgenic plants (Figures 4A,C,F,G).
Glandular trichomes are special structures derived from the epidermal cells of many plants (Werker, 2000). They have the potential to be a powerful secondary metabolic factory (Tissier, 2012). Furthermore, the various compounds in glandular trichomes are potential active ingredients for many applications, including curing diseases, in fragrance, in killing pests, etc. (Dixon, 2001; Schilmiller et al., 2008; Tissier, 2012; Huchelmann et al., 2017). In A. annua, one of the most well-known compounds is artemisinin which is synthesized in GSTs (Olofsson et al., 2011). It is feasible that artemisinin content could be increased by increasing GST density.
From the previous studies, R2R3-MYB and HD-ZIP IV TFs play very important roles in glandular trichome initiation in A. annua. AaMYB1 is the first R2R3-MYB found to regulate GST initiation positively in A. annua (Matías-Hernández et al., 2017). Moreover, AaMIXTA1 is found as a positive regulator in GST initiation (Shi et al., 2018). As AaMYB17 is homologous with AaMIXTA1, we tried to determine whether AaMYB17 and AaMIXTA1 have the same functions. qRT-PCR was performed to analyze the expression level of cutin- and wax-related synthase genes AaCYP77A1, AaCYP86A1, AaABCG12, AaKCS5, and AaCER1, which are activated by AaMIXTA1 significantly (Shi et al., 2018) in AaMYB17 overexpression and RNAi transgenic A. annua plants. The results indicated that there were no apparent differences (Supplementary Figure 4). The dual-LUC assay results in N. benthamiana revealed that AaMYB17 had no significant influence on AaCYP77A1, AaCYP86A1, AaABCG12, AaKCS5, and AaCER1 (Supplementary Figure 5). These results indicate that AaMYB17 differs from AaMIXTA1 in cuticle biosynthesis. Furthermore, two HD-ZIP IV TFs, AaHD1 and AaHD8, were found to positively regulate GST initiation (Yan et al., 2016, 2018). These studies indicate that R2R3-MYB and HD-ZIP IV TFs have a powerful function in GST initiation, and there may be more unknown R2R3-MYB and HD-ZIP IV TFs involved in GST initiation.
In A. thaliana, there is only one type of trichome which is non-glandular and unicellular. It is thought that unicellular trichomes and multicellular trichomes form a distinct pathway (Yang and Ye, 2013). In A. annua, there are two types of multicellular trichomes, GSTs and TSTs. This might provide a reference for multicellular trichome development. In this study, an R2R3-MYB transcription factor, AaMYB17, which positively regulates GST initiation, was identified. Our study expands the knowledge of the molecular mechanism of multicellular trichome initiation.
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
WQ and KT designed the research. WQ, LX, YL, HL, and TC carried out the expression analysis, vector construction, transgenic plant generation, subcellular localization, SEM, and dual-luciferase. WQ drafted the manuscript. DH, YL, LX, LL, XS, and KT revised the manuscript. All authors approved the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Funding. This research was supported by the National Key R&D Program of China (2018YFA0900600), grants from the Bill & Melinda Gates Foundation (OPP1199872), the National Natural Science Foundation of China (18Z103150043), and the China National Key Research and Development Program (2017ZX09101002-003-002).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.657156/full#supplementary-material
ReferencesBaumannK.Perez-RodriguezM.BradleyD.VenailJ.MartinC. (2007). Control of cell and petal morphogenesis by R2R3 MYB transcription factors.Development1341691–1701. 10.1242/dev.0283617376813BrockingtonS. F.Alvarez-FernandezR.LandisJ. B.AlcornK.WalkerR. H.ThomasM. M. (2013). Evolutionary analysis of the MIXTA gene family highlights potential targets for the study of cellular differentiation.Mol. Biol. Evol.30526–540. 10.1093/molbev/mss26023188591ChengY.CaoL.WangS.LiY.ZhouY. (2014). Analyses of plant leaf cell size, density and number, as well as trichome number using cell counter plugin.Bio Protoc.4:e1165.DixonR. A. (2001). Natural products and plant disease resistance.Nature411843–847. 10.1038/3508117811459067DubosC.StrackeR.GrotewoldE.WeisshaarB.MartinC.LepiniecL. C. (2010). MYB transcription factors in Arabidopsis.Trends Plant Sci.15573–581.DukeS. O.PaulR. N. (2013). Development and fine structure of the glandular trichomes of Artemisia annua L.Int. J. Plant Sci.154107–118. 10.1086/297096FridmanE.WangJ.IijimaY.FroehlichJ. E.GangD. R.OhlroggeJ. (2005). Metabolic, genomic, and biochemical analyses of glandular trichomes from the wild tomato species lycopersicon hirsutum identify a key enzyme in the biosynthesis of methylketones.Plant Cell171252–1267. 10.1105/tpc.104.02973615772286GershenzonJ.DudarevaN. (2007). The function of terpene natural products in the natural world.Nat. Chem. Biol.3:408. 10.1038/nchembio.2007.517576428GloverB. J.Perez-RodriguezM.MartinC. (1998). Development of several epidermal cell types can be specified by the same MYB-related plant transcription factor.Development125:3497.GrahamI. A.BesserK.BlumerS.BraniganC. A.CzechowskiT.EliasL. (2010). The genetic map of Artemisia annua L. identifies loci affecting yield of the antimalarial drug artemisinin.Science327328–331. 10.1126/science.118261220075252HuchelmannA.BoutryM.HachezC. (2017). Plant glandular trichomes: natural cell factories of high biotechnological interest.Plant Physiol.1756–22. 10.1104/pp.17.0072728724619JakobyM. J.FalkenhanD.MaderM. T.BrininstoolG.WischnitzkiE.PlatzN. (2008). Transcriptional profiling of mature Arabidopsis trichomes reveals that NOECK encodes the MIXTA-like transcriptional regulator MYB106.Plant Physiol.1481583–1602. 10.1104/pp.108.12697918805951JeffersonR.KavanaghT.BevanM. (1987). GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants.EMBO J.63901–3907. 10.1002/j.1460-2075.1987.tb02730.xJinH.MartinC. (1999). Multifunctionality and diversity within the plant MYB-gene family.Plant Mol. Biol.41577–585.KranzH. D.DenekampM.GrecoR.JinH.LeyvaA.MeissnerR. C. (2010). Towards functional characterisation of the members of theR2R3-MYBgene family from Arabidopsis thaliana.Plant J.16263–276. 10.1046/j.1365-313x.1998.00278.x9839469LivakK. J.SchmittgenT. D. (2002). Analysis of Relative gene expression data using real-time quantitative PCR.Methods25402–408. 10.1006/meth.2001.126211846609MachadoA.WuY.YangY.LlewellynD. J.DennisE. S. (2009). The MYB transcription factor GhMYB25 regulates early fibre and trichome development.Plant J.5952–62. 10.1111/j.1365-313x.2009.03847.x19309462MaesL.Van NieuwerburghF. C.ZhangY.ReedD. W.PollierJ.Vande CasteeleS. R. (2011). Dissection of the phytohormonal regulation of trichome formation and biosynthesis of the antimalarial compound artemisinin in Artemisia annua plants.New Phytol.189176–189. 10.1111/j.1469-8137.2010.03466.x20874804MaoJ.ZhangY.SangY.LiQ.YangH. (2005). A role for Arabidopsis cryptochromes and COP1 in the regulation of stomatal opening.Proc. Natl. Acad. Sci.U.S.A.10212270–12275. 10.1073/pnas.050101110216093319MartinC.Paz-AresJ. (1997). MYB transcription factors in plants.Trends Genet.1367–73. 10.1016/s0168-9525(96)10049-4MartinS. C. (2006). Trichomes: different regulatory networks lead to convergent structures.Trends Plant Sci.11274–280. 10.1016/j.tplants.2006.04.00816697247Matías-HernándezL.JiangW.YangK.TangK.PelazS. (2017). AaMYB1 and its orthologue AtMYB61 affect terpene metabolism and trichome development in Artemisia annua and Arabidopsis thaliana.Plant J.90520–534. 10.1111/tpj.1350928207974MauricioR.RausherM. D. (1997). Experimental Manipulation of putative selective agents provides evidence for the role of natural enemies in the evolution of plant defense.Evolution511435–1444. 10.1111/j.1558-5646.1997.tb01467.x28568624MitsudaN.Ohme-TakagiM. (2009). Functional analysis of transcription factors in Arabidopsis.Plant Cell Physiol.501232–1248. 10.1093/pcp/pcp07519478073NicholasK. B. (1997). GeneDoc: analysis and visualization of genetic variation, EMBNEW.Embnew News4:14.NodaK. I.GloverB. J.LinsteadP.MartinC. (1994). Flower colour intensity depends on specialized cell shapecontrolled by a Myb-related transcription factor.Nature369:661. 10.1038/369661a08208293OlofssonL.EngströmA.LundgrenA.BrodeliusP. E. (2011). Relative expression of genes of terpene metabolism in different tissues of Artemisia annua L.BMC Plant Biol.11:45. 10.1186/1471-2229-11-4521388533OshimaY.ShikataM.KoyamaT.OhtsuboN.MitsudaN.Ohme-TakagiM. (2013). MIXTA-like transcription factors and WAX INDUCER1/SHINE1 coordinately regulate cuticle development in Arabidopsis and Torenia fournieri.Plant Cell251609–1624. 10.1105/tpc.113.11078323709630PayneT.ClementJ.ArnoldD.LloydA. (1999). Heterologous myb genes distinct from GL1 enhance trichome production when overexpressed in Nicotiana tabacum.Development126671–682.PlettJ. M.WilkinsO.CampbellM. M.RalphS. G.ReganS. (2010). Endogenous overexpression of Populus MYB186 increases trichome density, improves insect pest resistance, and impacts plant growth.Plant J.64419–432. 10.1111/j.1365-313x.2010.04343.x20807210RomeroI.FuertesA.BenitoM. J.MalpicaJ. M.LeyvaA.Paz-AresJ. (1998). More than 80R2R3-MYB regulatory genes in the genome of Arabidopsis thaliana.Plant J.14273–284. 10.1046/j.1365-313x.1998.00113.x9628022RosinskiJ. A.AtchleyW. R. (1998). Molecular evolution of the Myb family of transcription factors: evidence for polyphyletic origin.J. Mol. Evol.4674–83. 10.1007/pl000062859419227SchilmillerA. L.LastR. L.PicherskyE. (2008). Harnessing plant trichome biochemistry for the production of useful compounds.Plant J.54702–711. 10.1111/j.1365-313x.2008.03432.x18476873SernaL.MartinC. (2006). Trichomes: different regulatory networks lead to convergent structures.Trends Plant Sci.11274–280.ShenQ.LuX.YanT.FuX.LvZ. (2016). The jasmonate-responsive AaMYC2 transcription factor positively regulates artemisinin biosynthesis in Artemisia annua.New Phytol.2101269–1281. 10.1111/nph.1387426864531ShiP.FuX.ShenQ.LiuM.PanQ.TangY. (2018). The roles of AaMIXTA1 in regulating the initiation of glandular trichomes and cuticle biosynthesis in Artemisia annua.New Phytol.217261–276. 10.1111/nph.1478928940606SinghN. D.KumarS.DaniellH. (2016). Expression of β-glucosidase increases trichome density and artemisinin content in transgenic Artemisia annua plants.Plant Biotechnol. J.141034–1045. 10.1111/pbi.1247626360801SoetaertS. S.Van NesteC. M.VandewoestyneM. L.HeadS. R.GoossensA.Van NieuwerburghF. C. (2013). Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua.BMC Plant Biol.13:220. 10.1186/1471-2229-13-22024359620SparkesI. A.RunionsJ.KearnsA.HawesC. (2006). Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants.Nat. Protoc.12019–2025. 10.1038/nprot.2006.28617487191StrackeR.WerberM.WeisshaarB. (2001). The R2R3-MYB gene family in Arabidopsis thaliana.Curr. Opin. Plant Biol.4447–456. 10.1016/s1369-5266(00)00199-0TamuraK.PetersonD.PetersonN.StecherG.NeiM.KumarS. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.Mol. Biol. Evol.282731–2739. 10.1093/molbev/msr12121546353TanH.XiaoL.GaoS.LiQ.ChenJ.XiaoY. (2015). TRICHOME AND ARTEMISININ REGULATOR 1 is required for trichome development and artemisinin biosynthesis in Artemisia annua.Mol. Plant81396–1411. 10.1016/j.molp.2015.04.00225870149TattiniM.GravanoE.PinelliP.MulinacciN.RomaniA. (2010). Flavonoids accumulate in leaves and glandular trichomes of Phillyrea latifolia exposed to excess solar radiation.New Phytol.14869–77. 10.1046/j.1469-8137.2000.00743.xTissierA. (2012). Glandular trichomes: what comes after expressed sequence tags?Plant J. Cell Mol. Biol.7051–68. 10.1111/j.1365-313x.2012.04913.x22449043WagnerG. J.WangE.ShepherdR. W. (2004). New approaches for studying and exploiting an old protuberance, the plant trichome.Ann. Bot.933–11. 10.1093/aob/mch01114678941WalfordS. A.WuY.LlewellynD. J.DennisE. S. (2011). GhMYB25-like: a key factor in early cotton fibre development.Plant J.65785–797. 10.1111/j.1365-313x.2010.04464.x21235650WeinholdA.BaldwinI. T. (2011). Trichome-derived O-acyl sugars are a first meal for caterpillars that tags them for predation.Proc. Natl. Acad. Sci. U.S.A.1087855–7859. 10.1073/pnas.110130610821518882WerkerE. (2000). Trichome diversity and development.Adv. Bot. Res.311–35. 10.1016/s0065-2296(00)31005-9WilsonM. E.AnuK.SakariJ. T. (2011). Review of cases with the emerging fifth human malaria parasite, Plasmodium knowlesi.Clin. Infect. Dis.521356–1362. 10.1093/cid/cir18021596677World Health Organization [WHO] (2019). World Malaria Report 2019.Geneva: World Health Organization.XiaoL.TanH.ZhangL. (2016). Artemisia annua glandular secretory trichomes: the biofactory of antimalarial agent artemisinin.Sci. Bull.6126–36. 10.1007/s11434-015-0980-zXieZ.KapteynJ.GangD. R. (2008). A systems biology investigation of the MEP/terpenoid and shikimate/phenylpropanoid pathways points to multiple levels of metabolic control in sweet basil glandular trichomes.Plant J.54349–361. 10.1111/j.1365-313x.2008.03429.x18248593YanT.ChenM.ShenQ.LiL.FuX.PanQ. (2016). HOMEODOMAIN PROTEIN 1 is required for jasmonate-mediated glandular trichome initiation in Artemisia annua.New Phytol.2131145–1155. 10.1111/nph.1420527659595YanT.LiL.XieL.ChenM.ShenQ.PanQ. (2018). A novel HD-ZIP IV/MIXTA complex promotes glandular trichome initiation and cuticle development in Artemisia annua.New Phytol.218567–578. 10.1111/nph.1500529377155YangC.YeZ. (2013). Trichomes as models for studying plant cell differentiation.Cell Mol. Life Sci.701937–1948. 10.1007/s00018-012-1147-622996257YangC.LiH.ZhangJ.LuoZ.GongP.ZhangC. (2011). A regulatory gene induces trichome formation and embryo lethality in tomato.Proc. Natl. Acad. Sci.10811836–11841. 10.1073/pnas.110053210821730153