- 1College of Forestry and Landscape Architecture, South China Agricultural University, Guangzhou, China
- 2Guangdong Provincial Key Laboratory of Plant Molecular Breeding, South China Agricultural University, Guangzhou, China
- 3Guangzhou Flower Research Center, Guangzhou, China
- 4Mid-Florida Research and Education Center, Environmental Horticulture Department, Institute of Food and Agricultural Sciences, University of Florida, Apopka, FL, United States
Modern anthuriums, Anthurium andraeanum (Hort.) are among the most popular flowering plants and widely used for interior decoration. Their popularity is largely attributed to the exotic spathes with different colors. Previous studies have reported color development in red spathe cultivars, but limited information is available on key genes regulating white and green colored spathes. This study analyzed anthocyanin, chlorophyll, and carotenoid contents as well as transcript differences in spathes of eight cultivars that differed in spathe colors ranging from red to white and green. Results showed that increased expression of a transcription factor AaMYB2 was associated with elevated levels of anthocyanin in spathes, but decreased expression of AaMYB2 and increased expression of AaLAR (leucoanthocyanidin reductase) and AaANR (anthocyanidin reductase) were accompanied with the accumulation of colorless proanthocyanidin, thus the white spathe. As to the green colored spathe, chlorophyll content in the green spathe cultivar was substantially higher than the other cultivars. Correspondingly, transcripts of chlorophyll biosynthesis-related genes AaHemB (porphobilinogen synthase) and AaPor (protochlorophyllide oxidoreductase) were highly upregulated but almost undetectable in white and red spathes. The increased expression of AaHemB and AaPor was correlated with the expression of transcription factor AaMYB124. Subsequently, qRT-PCR analysis confirmed their expression levels in nine additional cultivars with red, white, and green spathes. A working model for the formation of white and green spathes was proposed. White colored spathes are likely due to the decreased expression of AaMYB2 which results in increased expression of AaLAR and AaANR, and the green spathes are attributed to AaMYB124 enhanced expression of AaHemB and AaPor. Further research is warranted to test this working model.
Introduction
Anthurium Schott is the largest genus in the family Araceae and comprises more than 1,000 species (Croat, 1992). Among them, A. andraeanum Linden ex Andre and A. scherzerianum Schott were produced as cut flower and potted flowering plants, respectively (Kamemoto and Kuehnle, 1996; Anthura, 2002). With the efforts in anthurium breeding, increasing numbers of anthurium cultivars have been developed through the hybridization of A. andraeanum with other species or cultivars (Marco van Herk et al., 1998; Anthura, 2002). These commercial cultivars are complex interspecific hybrids (Elibox and Umaharan, 2008) and are collectively referred to as A. andraeanum (Hort.) or modern anthuriums to distinguish them from A. andraeanum Linden ex Andre (Kamemoto and Kuehnle, 1996).
Anthurium andraeanum (Hort.) cultivars have some distinct characteristics: (1) Spathes vary greatly in color, size, shape, and texture (Elibox and Umaharan, 2008); (2) they have more compact growth forms and can be used as either cut flower or potted flowering plants (Chen et al., 1999); and (3) spathes have extended longevity (Kamemoto and Kuehnle, 1996; Chen et al., 2002). As a result, these cultivars have quickly gained popularity in the world floral market. The turnover of potted anthurium was €9.94 million in 1997 and increased to €23.22 million in 2001 in Dutch auctions (Anthura, 2002), and subsequently the turnover increased to €50 million (Rikken, 2010). In China, the wholesale value of anthuriums ranked second after orchids, and more than 20 million potted anthuriums were sold in 2013, of which 6.5 million were sold during the Chinese New Year (Hua, 2014).
The increasing popularity is mainly attributed to their exotic shape and colorful spathes (Henny and Chen, 2003). The spathe is actually a modified leaf that initially covers a cylindrical shaped spadix on which both male and female flowers are found. The primary colors of spathes are red including bright red, dark red, pink, orange, brown, and coral as well as white. To increase colorations, cultivars with green, brown, and bi-color spathes were developed (Kamemoto and Kuehnle, 1996). Since color is by far the most important characteristic of ornamental plants (Halloran and Kuehnle, 1998; Chen, 2021), there is a continuous effort to develop more colorful spathes in anthurium breeding using genetic manipulation (Gopaulchan et al., 2013; Teixeira da Silva et al., 2015; Li et al., 2016; Li et al., 2019).
Flower color is primarily determined by pigments (Clegg and Durbin, 2000). There are four major pigments in plants: anthocyanins, betalains, carotenoids, and chlorophyll. Among them, anthocyanins, carotenoids, and chlorophyll are implicated in anthurium spathe color, and their concentrations and ratios define the spathe color and its intensity (Lin et al., 2022). Betalains are restricted to the core Caryophyllales and do not occur in other plants including anthurium (Timoneda et al., 2019).
Anthocyanins are the largest group of pigments ranging from pink, red, purple to blue, which is synthesized through the flavonoid pathway (Pucker and Selmar, 2022). Cinnamate-4-hydroxylase (C4H) plays a key role in flavanol accumulation during the early stage (Mizutani et al., 1993; Ren et al., 2019). Chalcone synthase (CHS) converts 4-coumaroyl-CoA and malonyl-CoA to naringenin chalcone (Feinbaum and Ausubel, 1988; Franken et al., 1991). Then naringenin chalcone is isomerized to naringenin through chalcone isomerase (CHI) (Grotewold and Peterson, 1994). The formation of dihydrokaempferol and dihydroquercetin from naringenin is catalyzed by flavanone 3-hydroxylase (F3H) and flavonoid 3′-hydroxylase (F3’H), respectively (Sparvoli et al., 1994; Brugliera et al., 1999). The conversion of dihydroflavonols into leucoanthocyanidins is catalyzed by dihydroflavonol 4-reductase (DFR) (Helariutta et al., 1993). Leucoanthocyanidins are channeled into anthocyanidins by anthocyanidin synthase (ANS) (Saito et al., 1999); and meanwhile leucoanthocyanidin reductase (LAR) converts leucoanthocyanidins to catechin (Tanner et al., 2003; Akagi et al., 2009). Anthocyanidin then may be converted to epicatechin by anthocyanidin reductase (ANR) (Xie et al., 2004; Akagi et al., 2009). Both catechin and epicatechin are colorless. Colored anthocyanidins were produced by flavonoid 3-O-glucosyl transferase (UFGT) (Boss et al., 1996; Fan et al., 2022). Anthocyanin O-methyltransferase (AOMT) is responsible for the methylation of cyanidin glycosides and plays an important role in purple coloration (Du et al., 2015). Flavonoids are synthesized in the cytosol and then transported to the vacuole. Glutathione S-transferase (GST) participates in flavonoid transport (Mueller et al., 2000; Kitamura et al., 2012; Sun et al., 2012).
Chlorophylls and carotenoids are photosynthetic pigments, which are biosynthesized in plastids from metabolic precursors derived from the methylerythritol 4-phosphate (MEP) pathway (Rodriguez-Concepcion and Boronat, 2002). Chlorophylls are green pigments. A series of important enzymes are involved in chlorophyll metabolism, such as porphobilinogen synthase (HemB) and protochlorophyllide reductase (Por) (Zhou et al., 2022). Carotenoids range from colorless to yellow, orange, and red; and a large number of genes participates in carotenoid biosynthesis, including phytoene desaturase gene (PDS), β-ring carotene hydroxylase gene (CrtR-b), zeaxanthin epoxidase gene (ZEP), and lutein deficient 5 gene (LUT5) (Nisar et al., 2015).
Transcription factors (TF), like R2R3-MYB and members of the basic helix-loop-helix (bHLH) family, are the most important classes of transcriptional regulators controlling secondary metabolism, developmental processes, and signal transduction (Yan et al., 2021). Anthocyanin and proanthocyanidin are secondary metabolites. The MYB-TF family has been found to act favorably in regulation of the pigment pathway by binding directly to the promoter of pigment biosynthesis genes (Aharoni et al., 2001; Ithal and Reddy, 2004; Jiang et al., 2021; Yang et al., 2022). Flavonoid pathway is also regulated by a combination model of transcription factors, including members of MYB, bHLH, and WD-repeat (WDR) families (Gonzalez et al., 2008; Petroni and Tonelli, 2011; Li, 2014). In Arabidopsis, TT2 (MYB123), TT8 (bHLH042), and TTG1 (WDR) are shown to participate in the control of proanthocyanidin biosynthesis (Baudry et al., 2004; Baudry et al., 2006; Peng et al., 2020). Overexpression of MYB75/PAP1 and MYB90/PAP2, whose sequences are similar to MYB113 and MYB114, respectively, results in the accumulation of anthocyanin (Borevitz et al., 2000; Gonzalez et al., 2008). MYB11, MYB12, and MYB111 control flavonol accumulation in different parts of the Arabidopsis seedling (Mehrtens et al., 2005; Stracke et al., 2007). MYB TFs activate or suppress anthocyanin biosynthesis genes, which either improve or inhibit anthocyanin accumulation (Li et al., 2022). In Anthurium, it has been reported that AaMYB2 is highly expressed in the red, pink, and purple spathes of cultivars but barely identifiable in cultivars with white and green spathes (Osorio-Guarín et al., 2021; Liu et al., 2023). AaMYB3 interacts with AabHLH1 to regulate proanthocyanidin (PA) accumulation, and several structure genes were co-expressed with AaMYB3 in red spathes of anthuriums (Li et al., 2019).
Previous studies of anthuriums have been largely focused on specific gene regulations in red colored spathes (Li et al., 2016, Li et al., 2019; Lin et al., 2022), but there have been no reports on key genes in regulation of white and green colorations. Additionally, there has been no information about the interaction between TFs and structure genes in regulation of these two colors. In this study, we hypothesized that the white colored spathe could be due to the increased expression of LAR and ANR, leading to the production of colorless PA, while the cause of the green spathe could be attributed to the elevated expression of genes related to chlorophyll biosynthesis, and TFs, likely some members of the MYB family could be implicated in the regulation of those key structure genes.
The objectives of this study were to analyze transcript differences among cultivars with different colored spathes through RNA-Seq technology, identify key genes implicated in spathe coloration, test our hypotheses on the occurrence of green and white spathes, and elucidate the relationship between TFs and key structure genes in regulation of the two colors in modern anthurium cultivars.
Materials and methods
Plant materials
Eight cultivars of A. andraeanum (Hort.) were grown in 20-cm plastic containers filled with a soilless substrate composed of 60% peat mixed with coconut chaff and perlite (STANLEY, Linyi, China) in a shaded greenhouse under a maximum photosynthetic radiation of 200 μmol m-2 s-1 in Guangzhou Flower Research Center, Guangzhou, Guangdong Province, China. Daytime and nighttime temperatures in the shaded greenhouse were maintained from 25 to 28°C and 19 to 21°C, respectively with a relative humidity ranging from 70 to 90%. Plants were fertilized and watered as described by Chen et al. (2012). Pinching of flower buds was used to promote vegetative growth in the initial flowering stage. Subsequently, plants entered their bloom stage.
Spathes of these cultivars varied in color ranging from red to white. For transcriptomic analysis of spathe colorations, three plants were randomly selected per cultivar, one spathe was collected from each plant at stage 7 of spathe development, i.e., peduncle fully extended but spathe showing full coloration without being open (Supplementary Figure S1). Thus, there were three biological replicates per cultivar. To analyze anthocyanin, chlorophyll, and carotenoid levels in fully open spathes, spathes at their developmental stage of 10 (Supplementary Figure S1), i.e., spathes were fully expanded but flowers were not dehisced (Collette et al., 2004), were collected (about 200 mg) from three randomly selected plants per cultivar. Collected spathes were immediately extracted with appropriate solutions mentioned below.
RNA extraction and transcriptomic analysis of spathes
Total RNA was extracted from collected spathes. RNA concentrations were measured using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies,CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA) and also checked by electrophoresis using 1% agarose gel (Supplemental Table S1). The mRNA was enriched using magnetic bead with oligo (dT) and used for cDNA synthesis. The first strand cDNA was synthesized using random hexamer primer and M-MuLV reverse transcriptase (Dingguo Changsheng, Inc. Beijing, China). The second strand cDNA synthesis was subsequently performed using DNA polymerase I with RNase H. cDNA fragments were preferentially selected from 150 to 200 bp in length. PCR amplification was performed for 13 cycles as recommended. The library preparation was performed according to the method described by Liu et al. (2021). The cDNA libraries were sequenced by Illumina Novaseq™ system at NextOmics Bio-Tech. Co. (Wuhan, China). Paired-end sequencing of the polyA enriched library was performed on an Illumina Novaseq 6000 platform generating 150 nt reads.
RNA-Seq data were analyzed to identify key structural genes and transcription factors involved in flavonoid biosynthesis and chlorophyll metabolism. Trinity software (version 2.4.0, Supplementary Table S2) was used to assemble the short reads (Grabherr et al., 2011). Trinity commands were as follows: Trinity –seqType fq –max_memory 120G –CPU 30 –min_kmer_cov 2 –full_cleanup –left clean.R1.fq.gz –right clean.R2.fq.gz –output trinity.out. Other parameters were default. Based on the sequence similarity principle, TIGR Gene Indices clustering tools (TGICL, version 2.1, Supplementary File S1) was further clustered to eliminate redundancy (Pertea et al., 2003). CDS and protein sequences were predicted by Transdecoder software (version 5.7.0, https://github.com/TransDecoder/TransDecoder/releases, default parameters); then, a unigene sequence set of 8 cultivars (24 samples) was obtained (Supplementary Table S3). The analysis on the transcriptome assembly to assess completeness was conducted using BUSCO (version 5.4.2) (Manni et al., 2021). The lineage dataset is embryophyta_odb10.
Unigene functions were annotated based on the following database: the NCBI COG/KOG database (http://www.ncbi.nlm.nih.gov/COG/), Kyoto Encyclopedia of Gene and Genomes database (KEGG, https://www.genome.jp/kegg/), non-redundant protein database (NR, ftp://ftp.ncbi.nih.gov/blast/db/), SwissProt protein database (https://www.uniprot.org/), and Gene Ontology (GO, http://geneontology.org/). Identification of flavonoid biosynthesis genes was annotated using KIPEs (Knowledge-based Identification of Pathway Enzymes, https://pbb-tools.de/KIPEs/) based on transcript sequences (Pucker et al., 2020). MYB were also predicted with a dedicated tool MYB annotator (Automatic annotation of MYBs, https://pbb-tools.de/MYB_annotator/) (Pucker, 2022). bHLH were predicted with a dedicated tool bHLH annotator (Automatic annotation of the bHLH gene family in plants, https://pbb-tools.de/bHLH_annotator/) (Thoben and Pucker, 2023). The identifications of flavonoid biosynthesis genes, MYB and bHLH were listed in Supplementary File S2.
Gene expression levels were estimated by RSEM software (version 1.1.12) and presented as FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) (Li and Dewey, 2011). Differentially expressed genes (DEGs) were analyzed between different groups by DESeq2 software (version 1.38.3) (Love et al., 2014). DEGs were detected and screened between the eight comparisons, using fold change ≥ 2 and false discovery rate (FDR) ≤ 0.05. The “pheatmap” R package (version 4.2.2, https://www.r-project.org/) was used to determine DEGs abundance, and the scale of heatmap for DEGs were converted from the FPKM was Log2 of the fold change. The enrichment analysis was performed in GO term and KEGG pathway.
For validation of candidate genes associated with the spathe coloration, red spathe cultivars Te Lun Sa (TLS), 2016, and A302; green spathe cultivars Kai Xin Guo (KXG), A166, and A086; and white spathe cultivars Bai Ma (BM), A231BAI, and A168 at stage 7 (Supplementary Figure S1) were harvested, immediately immersed in liquid nitrogen, and stored at -80°C for qRT-PCR analysis. These nine cultivars were grown in the same conditions as the aforementioned eight cultivars.
Extraction and quantification of anthocyanin, chlorophyll, and carotenoid levels
Anthocyanin is red in acid solution. Its color is in proportion to the anthocyanin content. The absorption peak of anthocyanin acid solution is 530, with which molar extinction coefficient is 4.62 × 104. Thus, ultraviolet-visible (UV-Vis) spectrophotometer can be used to measure the anthocyanin content.
To analyze anthocyanin, collected 0.2 g spathes from eight cultivars were extracted in 10 mL of ethanol: hydrochloride (99:1, v/v) at 32°C for 4 h. The sample was centrifuged at 5,000 rpm for 5 min (Centrifuge 5804 R, Eppendorf, Shanghai, China), and the supernatant was stored at 4°C. The residue was extracted with the extraction solution 1 to 2 times until the supernatant turned colorless. Absorbencies were measured at 530 nm and 650 nm wavelengths with an UV-Vis spectrophotometer (UV4800, Unico, Shanghai, China). The relative content of anthocyanin was calculated according to the method (Yang et al., 2022):
For quantification of total chlorophyll and carotenoid, fresh spathe samples were ground to powder in a mortar with liquid nitrogen, extracted with 80% acetone, filtered into a test tube, and finally extracted with ethyl acetate. The suspension was centrifuged for 5 min at 2,000 rpm. The absorbance of the supernatant was recorded at 663 nm, 645 nm, and 440 nm, respectively. The total contents of chlorophyll (Chl) and carotenoid were calculated according to the following formula (Zhou et al., 2013; Vaculík et al., 2015):
Proanthocyanidins extraction and determination
The UV/Vis spectrophotometry was used to evaluate the proanthocyanidins (PA) content. About 200 mg spathe was used for extraction. Soluble PAs from fresh anthurium spathes were extracted and measured by Innovabio Testing Technology Co. (Nanjing, Jiangsu). After ethanol precipitation, the PA content was determined by normal butyl alcohol-hydrochloric acid method at 546 nm. Proanthocyanidins (Cas: 4852-22-6, Yuanye, Shanghai, China) was used as a standard for PA quantification. The content of PA was calculated based on the standard curve according to the method (Li X. et al., 2021).
Cyanidin and pelargonidin contents measurement
About 100 mg spathe was ground to powder and extracted by ethyl alcohol: distilled water: hydrochloric acid (2:1:1, v/v/v) for supernatant preparation. SCIEX Qtrap6500 mass spectrometer System (AB, Massachusetts, USA) coupled to a LC-30AD high performance liquid chromatograph (Shimadzu, Kyoto, Japan) was used for liquid chromatography-tandem mass spectrometry (LC-MS/MS). Cyanidin chloride (Cas: 528-58-5), and pelargonidin chloride (Cas: 134-04-3, Yuanye, Shanghai, China) was used as a standard for cyanidin and pelargonidin quantification, respectively. The standard substance was dissolved in hydrochloric acid-methanol solution as single standard stock solution. LC-MS/MS analysis was performed by Innovabio Co. (Nanjing, Jiangsu) according to the method of Tomas (2022).
RNA isolation for qRT-PCR analysis
Spathes at stage 7 (Supplementary Figure S1) harvested from the above-mentioned cultivars were ground to a powder with a mortar and a pestle in liquid nitrogen. RNA was extracted using the RaPure Plant RNA Kit (Magen, Guangzhou, China) following the instructions of the manufacturer. Total RNA was quantified using NanodropND1000 UV/Vis spectrophotometer (Thermo Scientific, Walthman, MA), and the purity was assessed by the absorbance ratios of 260/280 nm and 260/230 nm. The integrity of the purified RNA was confirmed with 2% agarose gel electrophoresis. The first-strand cDNA was synthesized from 1 μg of total RNA using a TIANGEN FastKing cDNA Kit (TIANGEN, Beijing, China) according to the manufacturer’s instructions. The 20 μL cDNA samples were diluted to a 1:5 ratio using 80 μL RNase-free water for all gene amplification reactions.
qRT-PCR reactions were performed using SYBR Green detection chemistry by running on 96 well-plates with a CFX connect Real-Time PCR System (BIO-RAD, Shanghai, China) with three replicates for each sample. Primers for each target were designed using Primer3 web service (https://primer3.ut.ee/), and the sequences of all primers are listed in Supplemental Table S4. qRT-PCR analysis data were averaged by three biological replicates and calculated using the 2-△△Ct method (Livak and Schmittgen, 2001). All expression levels were normalized with internal control, AaGAPDH (Li S. et al., 2021).
Prediction of MYB binding motifs annotated on promoters
The upstream sequences (5,000 bp) of AaLAR (GenBank accession no. OR241534), AaANR (GenBank accession no. OR241535), AaHemB (GenBank accession no. OR241536), and AaPor (GenBank accession no. OR241537) genes were obtained from the Anthurium Genome Database of A. andraenum cv. Xiaojiao-Texana. Potential binding motifs for TFs of MYB family on AaLAR, AaANR, AaHemB, and AaPor promoter regions were predicted by the online website JASPAR2022 (http://jaspar.genereg.net/) (Fornes et al., 2019). The TF binding sites and motifs were further screened by a relative profile score threshold of 90%. The binding sites on the positive strand were presented and the common MYB motifs in the anthocyanin or chlorophyll biosynthesis pathway were marked. As the motifs in the JASPAR database were only included in the model species, Arabidopsis thaliana, Oryza sativa, and Zea mays, the MYB proteins in A. andraeanum were aligned together by Clustalw for conducting phylogenetic analysis using the neighbor-joining method in MEGA version 7.0 to determine the genetic relationship (Supplementary Figures S2, S3).
Statistical analysis
Data for pigment contents and gene relative expression levels were statistically analyzed using SPSS (SPSS statistics 19, IBM Corp., Armonk, NY, USA). If significance occurred among cultivars, means were separated using Fisher’s LSD at P< 0.05 level. Data are presented as mean ± S.E. with either 3 or 5 replicates.
Results
Spathe pigment contents of eight cultivars
Spathes of eight anthurium cultivars had eight different colors, including ‘Rosa’ (Ros) with red spathe; ‘Wu Yang’ (WY), orange spathe; ‘Sante’ (San), coral spathe; ‘Guang Hua Zi Yun’ (GHZ), purple spathe; ‘Pink Champion’ (Pin), pink spathe; ‘Sonate’ (Son), light pink; ‘Midori’ (Mid), green spathe; and ‘Acropolis’ (Acr) (Figure 1A). Anthocyanin contents in spathes of Ros and GHZ were the highest, 14.0 and 15.0 ΔA/g·FW, respectively, which were substantially greater than the others (Figure 1B). Anthocyanin content in Acr was the lowest, only 0.5 ΔA/g·FW. The other cultivars had anthocyanin ranging from 1.0 to 3.8 ΔA/g·FW, suggesting that anthocyanin contents were associated with redness of spathes. The main component of anthocyanins in GHZ was cyanidin, up to 0.25 mg·g-1·FW, rather than pelargonidin, which was only 0.025 mg·g-1·FW (Supplementary Figure S4A). Proanthocyanidin (PA) contents of Ros and San were highest, 11.4 and 14.5 mg·g-1·FW, respectively, following by Son and Mid. PA content in Acr was about 5.27 mg·g-1·FW (Supplementary Figure S4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g001.jpg)
Figure 1 Spathe phenotype and pigment analysis of eight anthurium cultivars. Phenotypes of anthurium spathes (A), the total content of anthocyanin (B), chlorophyll (C), and carotenoid (D) in spathes of eight cultivars are shown as means with corresponding standard error (n = 5). Different letters above bars indicate significant differences analyzed by Fisher’s LSD test at P< 0.05 level.
The Chl and carotenoid levels were the highest in Mid at 0.16 ΔA/g·FW and 0.08 ΔA/g·FW, respectively (Figures 1C, D). Notably, the contents of both Chl and carotenoid in the other seven cultivars were not significantly different, ranging from 0.02 to 0.04 ΔA/g·FW for Chl, and 0.002 to 0.02 ΔA/g·FW for carotenoid, respectively. These results showed that the pigment composition and content in green and white spathes were different from those of the other colored spathes.
Transcriptome sequencing and annotation
Raw data or reads generated from the transcriptome sequencing of 24 samples (8 cultivars, each with three biological replicates) ranged from 40,964,478 to 63,961,814. After stringent quality control checks and filtering, a total of 1,947,453,054 paired-end clean reads were obtained with Q30 values higher than 92% in each sample library. The data volume per library was greater than 6 Gb. The clean reads were spliced and assembled using Trinity. N50 length of the Trinity assembly ranged from 1610 to 1799 bp (Supplementary Table S2). Trinity sequence of all samples were merged. According to the sequence similarity principle, the merged Trinity assembly sequences were clustered for a set of unigene sequences, which resulted in a total of 62,013 unigenes. The length of assembled unigenes varied from 220 to 21,222 bp. The average length of the unigenes was 2,038.45 bp, and the N50 was 2,519 bp (Supplementary Figure S5A, Supplementary Table S3). These data showed that the throughput and sequencing quality were high enough for further analysis. The complete BUSCOs of the transcriptome assembly was 1387, accounting for 85.9% (Supplementary Figure S6). In this study, the tissue used for RNA-Seq was spathes only. Some genes were hardly detectable in the spathe, suggesting that some genes were not expressed in this special organ. Thus, it was reasonable for the incompleteness of this transcriptome assembly.
A total of 48,582 coding sequences (CDS) were predicted, and their lengths varied from 261 to 15,300 bp with the majority ranging from 972 to 1,206 bp (Supplementary Figure S5B). These results implied that the reads derived from transcriptome sequencing had high quality and reliability. Further, 33,392 unigenes were assigned to 23 categories in the KOG database (Supplementary Figure S5C). Among them, 6,806 unigenes, “Posttranslational modification, protein turnover, chaperones” represented the largest group with 3,453 unigenes, followed by “Signal transduction mechanisms” (2,858 unigenes)”, and Transcription (1,965 unigenes). In the NR database, the identified unigenes were used to detect their homologous genes in other species by blast, and the results showed that the ratio of homologous genes in Elaeis guineensis, Phoenix dactylifera, Nelumbo nucifera, Musa acuminata subsp. Malaccensis, and Ananas comosus were 25.60%, 21.40%, 8.53%, 5.92%, and 4.59%, respectively (Supplementary Figure S5D).
Differentially expressed genes among different colored spathes
Transcriptome analysis identified differentially expressed genes (DEGs) in different colored spathes based on their FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) levels. To identify the key DEGs involved in different colorations, seven pair-wise comparisons were conducted which included Ros vs. Acr; WY vs. Acr; San vs. Acr; Mid vs. Acr; GHZ vs. Acr; Pin vs. Acr; and Son vs. Acr (Supplementary Figure S7A). Among these comparisons, the greatest abundance of DEGs (29,588) was found in Pin and Acr libraries, of which 13,891 and 15,769 genes were up-regulated and down-regulated, respectively. Conversely, the lowest abundance of DEGs (14,015) were recorded in GHZ and Acr libraries (Supplementary Figure S7A). Gene functional annotation was carried out to understand the function of DEGs. According to the GO annotation, a total of 51,469 DEGs were spread across 1,123 terms consisting of three domains: molecular functions accounting for 57.52%, cellular components 12.54%, and biological processes 29.94%. The most frequent terms under the molecular function were protein binding, ATP binding, protein kinase activity, and nucleic acid binding, DNA binding, and zinc ion binding (Supplementary Figure S7B). To further examine the metabolic process of the DEGs, the KEGG pathway for DEGs in seven comparison pairs were conducted. The top 20 enriched pathways were listed in all comparable groups (Supplementary Figure S7C). From the results of KEGG enrichment analysis, the pathways for flavonoid biosynthesis, and porphyrin and chlorophyll metabolism were enriched in the comparisons of Ros vs. Acr; GHZ vs. Acr; WY vs. Acr; and Mid vs. Acr (Supplementary Figures S7C, S8).
Differentially expressed genes in the flavonoid biosynthesis pathway
The DEGs associated with the flavonoid biosynthesis were analyzed through the examination of KEGG enrichment and gene annotations, from which 56 DEGs in the flavonoid pathway were identified, including nine MYB genes and seven bHLH genes regulating flavonoid biosynthesis, and one C4H, one 4-coumarate-CoA ligase-like (4CL), five CHS, four CHI, one F3H, two F3’H, one flavonol synthase (FLS), four DFR, six LAR, one ANS, three anthocyanidin 3-O-glucosyltransferase (3GT), two anthocyanidin 5,3-O-glucosyltransferase (5,3GT), four ANR, two caffeoyl-CoA O-methyltransferase (CCoAOMT), and three glutathione S-transferase F12 (also named TT19) in seven unrepetitive comparison pairs of Ros, WY, San, GHZ, Pin, Son, Mid, and Acr. According to the expression level of DEGs shown in heatmap (Figure 2), MYBs were highly expressed in Ros, WY, San, and GHZ compared to their expression in the other cultivars, of which MYB2 Unigene32682 exhibited notable changes with the cultivars or the content of anthocyanin in spathes (Figure 2, Supplementary Figure S9). The results implied that MYB genes could play an important role in regulation of anthocyanin biosynthesis in anthuriums. For CHI, its expression levels in Pin, Son, Mid, and Acr were relatively lower than in other cultivars (Figure 2), indicating that more naringenin chalcones was transformed into naringenin in Ros, WY, San, and GHZ. The expression level of F3’H in Son was significantly higher than other spathes, suggesting that more dihydroquercetin (DHQ) was synthesized in Son than the other spathes. F3H gene expression in San was the highest among eight anthurium cultivars, indicating that naringenin in San was catalyzed into more dihydrokaempferol (DHK). Higher levels of DFR expression occurred in Ros, WY, and GHZ, converting more dihydroquercetin/dihydrokaempferol to leuco-cyanidin/pelargonidin. The expression of ANS was increased in WY, San, and Son, with cyanidin and pelargonidin accumulation. The gene expression of 5,3GT was highly detected in WY and GHZ, accompanied with 3,5-O-glucoside product. Downregulation of a decorating enzyme CCoAOMT gene seemed to activate the anthocyanin biosynthesis in red- and purple-colored spathes. Anthocyanin transport carrier TT19 was induced in Ros, WY, San, GHZ, Pin, and Son, supporting anthocyanin accumulation and especially cyanidin accumulation. The increased expression of the LAR gene, such as Unigene59237 (Figure 2), could result in converting leucocyanidin and leucopelargonidin to catechin and afzelechin, respectively. The highly expressed ANR gene, especially Unigene7781, could make cyanine and pelargonidin flow into epicatechin and epiafzelechin, respectively. The results showed that the expression of LAR and ANR were significantly different in green and white spathe cultivars compared with the six other cultivars. The increased expression of LAR and ANR could reduce the catalytic substrate to synthesize anthocyanins, resulting in accumulation of proanthocyanidins in spathes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g002.jpg)
Figure 2 DEGs involved in the flavonoid metabolism and their expression levels. Heatmaps were constructed based on Fragments per kilobase per million mapped reads (FPKM) of eight anthurium cultivars with different spathe colors. Color bar: Log2 (fold changes).
Transcript factors and key structural genes related to flavonoid biosynthesis
Except for differential expression of flavonoid biosynthesis genes, we also analyzed if MYB-TFs were implicated in the regulation of some of the genes. A MYB-TF (Unigene32682, MYB2), a LAR gene (Unigene59237), and a ANR gene (Unigene7781) were chosen from DEGs to analyze their expressions in the eight cultivars (Figure 2). The transcript level of MYB2 was significantly lower in the green or white spathed cultivar compared with the color intensified ones. Conversely, in comparison to the color intensified spathes, green and white spathes had significantly high levels of LAR and ANR transcripts (Figure 3A, Supplementary Figure S10).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g003.jpg)
Figure 3 Molecular analysis of key genes in the flavonoid biosynthesis pathway. (A) Expression analysis of AaMYB2, AaLAR, and AaANR in the spathe tissue of six cultivars. Expression levels were normalized based on the expression of the GADPH gene. Data are presented as means of three biological replicates and error bar shows standard error (n = 5). Different letters above bars indicate significant differences determined by Fisher’s LSD test at P< 0.05 level. (B) Prediction of MYB2 transcription factor binding sites in the promoters of AaLAR and AaANR. The binding sites were shown as red ovals (positive strand).
Promoter analysis of LAR and ANR genes showed that the promoter of LAR contained two motifs bound by MYB2, and the promoter of ANR also had one motif bound by MYB2 (Figure 3B, Supplementary Figure S2), implying that the expression levels of LAR and ANR could be likely regulated by the TF-MYB2 (Figure 3B). In addition, expression analysis revealed that the expression pattern of MYB2 was opposite to that of LAR and ANR, suggesting MYB2 might negatively regulate the expression of LAR and ANR (Figure 3A).
bHLH family members were also analyzed. Expression levels of bHLH113 and bHLH49 increased in purple spathe GHZ, while bHLH62 and bHLH57 were highly detected in red spathe Ros (Figure 2), indicating bHLH might be involved in the positive regulation of anthocyanin biosynthesis.
Transcription factors and key structural genes related to chlorophyll biosynthesis
Comparative analysis of DEGs showed that the expression levels of a set of MYB-TFs, including MYB124 (Unigene12793) were significantly higher in Mid (green spathe) than in the others, which could suggest that those MYB-TFs may also participate in the regulation of Chl synthesis (Figure 4A). Thus, the DEGs in the chlorophyll biosynthesis pathway were further analyzed, and data indicated that the expressions levels of HemB (Unigene51685) and Por (Unigene38255) were significantly greater in Mid than in the other cultivars (Figure 4B, Supplementary Figure S10). To explore the expression pattern of the DEGs, qRT-PCR was conducted with three genes, MYB124, HemB, and Por, involved in the biosynthesis of Chl. The results showed that the relative expression levels of the three genes were similar to the FPKM values (Figure 4C), suggesting that MYB124 was positively correlated with the expression of HemB and Por in green spathe. These data implied that MYB124 could promote Chl synthesis by regulating HemB and Por genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g004.jpg)
Figure 4 Molecular analysis of key genes in the chlorophyll metabolism pathway. (A) Heatmap of MYBs related to chlorophyll metabolism. (B) Expression analysis of AaMYB124, AaHemB, and AaPor in the spathe tissue of six cultivars. Expression levels were normalized to the expression of the GADPH gene. Data are presented as means of three biological replicates and error bar shows standard error (n = 5). Different letters above bars indicate significant differences analyzed by Fisher’s LSD at P< 0.05 level. (C) DEGs involved in chlorophyll metabolism and their expression levels. Each colored cell represents the average FPKM value standardized by the Z-Score of each gene. (D) Prediction of MYB124 transcription factors binding sites in the promoters of AaHemB and AaPor. The binding sites were shown in red (MYB124) and purple ovals (positive strand).
To further investigate the regulatory relationships among MYB124, HemB, and Por, we analyzed the promoter sequences of HemB and Por, respectively. The analysis showed that upstream of HemB and Por contained MYB124-bounding motifs, respectively (Figure 4D, Supplementary Figure S3). The results supported the above analysis that TF-MYB124 enhanced the expression of the downstream genes, HemB and Por and promoted Chl accumulations, thus the green spathe in cultivar Mid.
Differentially expressed genes in the carotenoid biosynthesis pathway
Since Mid with green spathe was loaded with carotenoids (Figure 1D), we took a closer look at gene expression in carotenoid biosynthesis pathway. Five unigenes were identified, including one PDS (Unigene54090), two LUT5 (Unigene14810 and Unigene58790), one CrtR-b (Unigene54067), and one ZEP (Unigene8259). Upstream gene phytoene desaturase (PDS), which is associated with carotenoid content, was detected to be highly expressed in Mid. Increased expression of lutein deficient (LUT5) and β -carotene hydroxylase (CrtR-b) in Mid suggested more lutein accumulation (Supplementary Figure S11). Meanwhile, the level of zeaxanthin epoxidase (ZEP) was significantly greater in Mid than other cultivars, indicating that ZEP may be involved in carotenoid synthesis through the epoxidation of β-carotene and β-cryptoxanthin (Supplementary Figure S11). Thus, increased carotenoid content in Mid might be caused by these elevated expression genes in carotenoid biosynthesis pathway.
qRT-PCR analysis of candidate genes
To further confirm key genes involved in the change of spathe colors, we used qRT-PCR to analyze the expression pattern of seven candidate genes in nine additional cultivars with red, green, and white spathes at S7 stage (Figure 5, Supplementary Figure S1). As expected, the expression levels of LAR and ANR in green or white spathes were significantly higher than in red spathes. In contrast, MYB2 gene exhibited the opposite expression pattern of LAR and ANR (Figure 5, Supplementary Figure S9). Therefore, these results confirmed that the expression of MYB2 was negatively correlated with the expression of LAR and ANR.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g005.jpg)
Figure 5 Expression analysis of candidate genes in nine additional anthurium cultivars with different spathe colors. (A) Spathe phenotypes of red cultivars TLS, 2016, A302, green cultivars KXG, A166, A086, and white cultivars BM, A231BAI, and A168 at early bud (S7 stage). (B) The transcript level of AaMYB2, AaMYB124, AaLAR, AaANR, AaHemB, and AaPor in the spathe tissue of nine cultivars. Expression levels were normalized based on the expression of the GADPH gene. Data are presented as means of three biological replicates and the error bar shows standard error (n = 5). Different letters above bars indicate significant differences analyzed by Fisher’s LSD test at P< 0.05 level.
HemB and Por genes were highly expressed in green spathes but at very low levels in red and white spathes (Figure 5). Such expression levels were closely associated with Chl contents in the respective spathes (Figure 1C). Additionally, MYB124 was highly correlated with HemB or Por expression and the total Chl content, suggesting that HemB and Por may be regulated by MYB124 in regulation of Chl biosynthesis.
Taken together, these results showed that the white colored spathe was probably due to the increased expression of LAR and ANR, resulting in the production of colorless PA. The green-colored spathe was likely attributed to the elevated expression of HemB and Por, leading to chlorophyll biosynthesis. While MYB2 could negatively regulate the expression of LAR and ANR, and MYB124 enhanced the expression of HemB and Por.
Discussion
Modifying spathe colors has been a major emphasis in anthurium breeding since the 1950s (Kamemoto and Nakaone, 1955; Kamemoto and Kuehnle, 1996). Later, Elibox and Umaharan (2008) reported that three dominant genes, O, R, and M control spathe colors. Plants with O_R_ genotypes produce colorful spathes, while those with ooR_ or O_rr give rise to white spathes. Elibox and Umaharan (2008) proposed that R and O encoded transcription factors regulating CHS, F3H, ANS, and DFR expression, respectively, and M encoded F3’H (flavonoid 3’-hydroxylase) that catalyzes dihydrokaempferol to dihydroquercetin. Subsequent analyses of selected genes at mRNA and protein levels indicated that F3H and ANS were probably controlled by a transcription factor encoded by R locus, and CHS could be controlled by a different mechanism (Gopaulchan et al., 2014). These studies provide classical genetic and molecular information on anthocyanin biosynthesis in modern anthuriums. With the advent of RNA-Seq technology, Li et al. (2016) isolated an R2R3-MYB gene, MYB2 in A. andraeanum (Hort.) and found that MYB2 plays an important role in spathe accumulation of anthocyanin by positively regulating CHS, ANS, and F3’H expressions. Furthermore, Li et al. (2019) transformed AaMYB3 from A. andraeanum (Hort.) into tobacco and found that AaMYB3 upregulated NtDFR and NtANS. Although the present study was primarily focused on white and green spathes, our analysis of eight cultivars with different colors showed that cultivars with red colored spathes (Ros, WY, and GHZ) generally had higher contents of anthocyanin. The expression of MYB2 gene in Ros, WY, and GHZ was mostly higher than other cultivars (Figures 1, 2). Our results agree with the report of Li et al. (2016) that MYB genes are critical for anthocyanin accumulation in reddish spathes of anthuriums. Phylogenetic analysis suggests that AaMYB2 was highly similar to AtMYB113. AtMYB113 of Arabidopsis has been mentioned in activating anthocyanin synthesis (Gonzalez et al., 2008). Thus, MYB2 gene was selected as a target gene as the subsequent analysis. The present study explored the molecular basis of white colored spathes. In the flavonoid pathway, anthocyanidins can be catalyzed by ANR to produce epicatechin, and leucoanthocyanidins are catalyzed by LAR to synthesize catechin (Akagi et al., 2009), which are important steps in the biosynthesis of colorless proanthocyanidins (Bogs et al., 2007; Nie and Stürzenbaum, 2019). Our data showed that anthocyanin was barely detectable, and chlorophyll and carotenoid contents were lower in white spathe (Acr). Correspondingly, the expression of AaMYB2 in white spathes was extremely lower, but the expression levels of AaLAR and AaANR were higher in white spathes (Figures 2, 3). These data suggest that the decreased expression of AaMYB2 is correlated with the activation of LAR and ANR expressions. Increasing expression of LAR and ANR genes could enhance proanthocyanin production, and in turn alter the competition between anthocyanin and proanthocyanidin biosynthesis, resulting in white colored spathes. Similar results were reported in other crops as well. For example, transcript levels of LAR and ANR genes were higher in white than red fruits in strawberry (Salvatierra et al., 2013). The transcript accumulation of MmANR gene leads to the abundant catechin and epicatechin in white flowers in Michelia maudiae (Lang et al., 2019). Overexpression of ANR gene promotes the proanthocyanidins biosynthesis in tobacco flowers (Han et al., 2012). These reports concur with our results that increasing expression of LAR and ANR genes is associated with organ discoloration.
We also investigated molecular underpinnings of green colored spathes. Similar to the white spathe, anthocyanidins were hardly detected in green spathes (Figure 1). This could be attributed to the lower expression of MYB2 and increased expression of LAR and ANR (Figure 3), resulting in an extremely low anthocyanin content. In addition, our data showed that the expression level of AaMYB124 was higher, which was associated with the increased expression of HemB and Por and elevated contents of chlorophyll (Figures 4, 5). These results suggest that AaMYB124 may promote chlorophyll synthesis by positively regulating HemB and Por expression in green spathe. It was worth to pay more attention to the transcription factor AaMYB124 in enhancing the transcription of chlorophyll metabolism genes thus green spathe. Therefore, the formation of green spathes is not only dependent on the accumulation of a large amount of chlorophyll via upregulating HemB and Por by AaMYB124, but also related to the high expression level of AaLAR and AaANR, which limits anthocyanin biosynthesis.
Considering the above results, we proposed a working model to illustrate the formation of red, white, and green spathes in A. andraeanum (Hort.) (Figure 6). AaMYBs play important roles in regulation of pigment biosynthesis. Increased expression of AaMYB2 enhances the expression of AaCHS, AaF3’H, and AaANS and anthocyanin biosynthesis, resulting in the formation of reddish spathes. On the other hand, the decreased expression of AaMYB2 is associated with the increased expression of AaLAR and AaANR, which could block anthocyanin synthesis and promote biosynthesis of proanthocyanidins, leading to spathes with white color. Under the decreased expression of AaMYB2, the formation of green colored spathes is attributed to the increased expression of AaMYB124 as it promotes the expression of both HemB and Por, thus chlorophyll biosynthesis and the greenish spathes. This working model slightly differs from the report of Gopaulchan et al. (2014) that F3H and ANS were controlled by a TF encoded by R locus, but CHS could be controlled by a different mechanism. However, this model concurs with the report of Li et al. (2016) that AaMYB2 enhances the accumulation of anthocyanin in spathes by positively regulating CHS, ANS, and F3’H. Li et al. (2019) also mentioned that AaMYB3 was associated with NtLAR and NtANR in tobacco, resulting in a lighter color by reducing anthocyanin level and enhancing proanthocyanidin accumulation. As far as is known, this is the first attempt to elucidate the likely mechanisms underlying the formation of white and green colored spathes in modern anthuriums, i.e., their formations are regulated by AaMYBs: the decreased expression of AaMYB2 and the increased expressions of AaMYB124. Further research will be conducted to test this model.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1208226/fpls-14-1208226-HTML/image_m/fpls-14-1208226-g006.jpg)
Figure 6 A working model for the formation of red, white, and green spathes in A. andraeanum (Hort.). AaMYB3 interacts with AabHLH1 to regulate proanthocyanidin accumulation. AaMYB3 was coexpressed with AaCHS, AaF3H, AaDFR, AaANS, AaLAR, and AaANR in the developing red spathe. The different expression of AaUFGT might explain the formation of anthocyanin-loss mutant. Increased expression of AaMYB2 enhances the expression of AaCHS, AaF3’H, and AaANS as well as anthocyanin accumulation in reddish spathes. Meanwhile, decreased expression of AaMYB2 activates AaLAR and AaANR expression and promotes biosynthesis of proanthocyanidins, leading to spathes with white color. Under the suppressed expression of AaMYB2, two other TFs AaMYB124 upregulated HemB and Por expression, thus chlorophyll accumulation and spathe with green coloration.
Our effort on elucidating the molecular basis behind the spathe coloration is intended to improve modern anthuriums through genetic manipulation. There were about 14 reported anthurium transformations in the literature, but none of transformed genes are TFs (Teixeira da Silva et al., 2015). On the other hand, understanding the role of TFs and genetic transformation of the TFs could substantially alter plant phenotypes. For example, transformation of a member of R2R3-MYB gene VvMYBA1 from grape plants into grape, tobacco, and Ficus lyrata produced purple-leaved grape and purple flowered tobacco (Li et al., 2011) as well purple and variegated leaved fiddle fig plants (Zhao et al., 2013). With the availability of anthurium transformation systems and a better understanding of TFs involved in the pigment biosynthesis, genetic transformation or gene editing of TFs could be a viable approach for modification of spathe colors in modern anthuriums.
Conclusions
This study tested anthocyanin, chlorophyll, and carotenoid contents in different colored spathes of eight cultivars and analyzed their transcript differences using RNA-Seq. Anthocyanin was generally higher in red spathes depending on the cultivar, the elevated level of chlorophyll occurred only in green spathes, but both anthocyanin and chlorophyll were extremely lower in white spathes. Thus, RNA-Seq analysis was focused on flavonoid biosynthesis and chlorophyll metabolism pathways. Results showed that the decreased expression of AaMYB2 was associated with increased expression of AaLAR and AaANR, which led to the conversion of anthocyanidins to colorless proanthocyanidin, thus white spathes. Meanwhile, the decreased expression of AaMYB2 was accompanied with increased expression of AaMYB124, and their expression upregulated AaHemB and AaPor, key genes in chlorophyll metabolism, thus green spathes. A working model for the formation of red, white, and green spathes in modern anthuriums was proposed, which will be verified in further studies.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number of Ros (SAMN34164551), WY (SAMN34164548), San (SAMN34164547), GHZ (SAMN34164553), Pin (SAMN34164550), Son (SAMN34164549), Mid (SAMN34164552), and Acr (SAMN35719425) available at: https://submit.ncbi.nlm.nih.gov/about/sra/, and AaLAR (OR241534), AaANR (OR241535), AaHemB (OR241536), and AaPor (OR241537) available at: https://submit.ncbi.nlm.nih.gov/about/bankit/.
Author contributions
Writing-original draft, JL and QT; formal analysis, JL and LZ; funding acquisition, MY; methodology, ZY; data curation, ZY; resources, QX; investigation, LZ; project administration, XZ; interpretation, writing-review and editing, X-QZ and JC; supervision, HG. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by Key-Area Research and Development Program of Guangdong Province (2022B0202080004), Project of Guangzhou Municipal Science and Technology Bureau (202206010136), Project of Rural Revitalization Strategy of Guangdong Province (2022-NPY-01-006), and Project of Guangzhou Municipal Bureau of Agriculture and Rural Affairs (2112422, 22103647, 23105372).
Acknowledgments
The authors would like to thank Ms. Terri A. Mellich for critical review of this 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1208226/full#supplementary-material
References
Aharoni, A., De Vos, C. H., Wein, M., Sun, Z., Greco, R., Kroon, A., et al. (2001). The strawberry FaMYB1 transcription factor suppresses anthocyanin and flavonol accumulation in transgenic tobacco. Plant J. 28, 319–332. doi: 10.1046/j.1365-313x.2001.01154.x
Akagi, T., Ikegami, A., Suzuki, Y., Yoshida, J., Yamada, M., Sato, A., et al. (2009). Expression balances of structural genes in shikimate and flavonoid biosynthesis cause a difference in proanthocyanidin accumulation in persimmon (Diospyros kaki Thunb.) fruit. Planta 230, 899–915. doi: 10.1007/s00425-009-0991-6
Baudry, A., Caboche, M., Lepiniec, L. (2006). TT8 controls its own expression in a feedback regulation involving TTG1 and homologous MYB and BHLH factors, allowing a strong and cell-specific accumulation of flavonoids in Arabidopsis thaliana. Plant J. 46, 768–779. doi: 10.1111/j.1365-313X.2006.02733.x
Baudry, A., Heim, M. A., Dubreucq, B., Caboche, M., Weisshaar, B., Lepiniec, L. (2004). TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant J. 39, 366–380. doi: 10.1111/j.1365-313X.2004.02138.x
Bogs, J., Jaffé, F. W., Takos, A. M., Walker, A. R., Robinson, S. P. (2007). The grapevine transcription factor VvMYBPA1 regulates proanthocyanidin synthesis during fruit development. Plant Physiol. 143, 1347–1361. doi: 10.1104/pp.106.093203
Borevitz, J. O., Xia, Y., Blount, J., Dixon, R. A., Lamb, C. (2000). Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12, 2383–2394. doi: 10.1105/tpc.12.12.2383
Boss, P. K., Davies, C., Robinson, S. P. (1996). Expression of anthocyanin biosynthesis pathway genes in red and white grapes. Plant Mol. Biol. 32, 565–569. doi: 10.1007/BF00019111
Brugliera, F., Barri-Rewell, G., Holton, T. A., Mason, J. G. (1999). Isolation and characterization of a flavonoid 3’-hydroxylase cDNA clone corresponding to the Ht1 locus of Petunia hybrida. Plant J. 19, 441–451. doi: 10.1046/j.1365-313x.1999.00539.x
Chen, J. (2021). Ornamental plant research inaugural editorial. Ornam. Plant Res. 1, 1–2. doi: 10.48130/opr-2021-0001
Chen, J., Henny, R. J., McConnell, B. D. (2002). “Development of new foliage plant cultivars,” in Trends in New Crops and New Uses. Eds. Janick, J., Whipkey, A. (Portland, OR: Timber Press, Inc), 446–452.
Chen, J., Henny, R. J., Robinson, C. A., Mellich, T., Caldwell, R. D. (1999). Potted anthurium: An interior-flowering foliage plant. Proc. Fla. State Hortic. 112, 280–281. doi: 10.32473/edis-ep159-2003
Chen, J., McConnell, D. B., Henny, R. J., Everitt, K. C. (2012). Cultural Guidelines for Commercial Production of Interiorscape Anthurium. (Gainesville, FL: IFAS Extension Bulletin, University of Florida). ENH956, 3.
Clegg, M. T., Durbin, M. L. (2000). Flower color variation: A model for the experimental study of evolution. Proc. Natl. Acad. Sci. U.S.A. 97, 7016–7023. doi: 10.1073/pnas.97.13.7016
Collette, V. E., Jameson, P. E., Schwinn, K. E., Umaharan, P., Davies, K. M. (2004). Temporal and spatial expression of flavonoid biosynthetic genes in flowers of Anthurium andraeanum. Physiol. Plant 122, 297–304. doi: 10.1111/j.1399-3054.2004.00402.x
Croat, T. B. (1992). Species diversity of Araceae in Colombia: preliminary survey. Ann. Mo. Bot. Gard. 79, 17–28. doi: 10.2307/2399806
Du, H., Wu, J., Ji, K., Zeng, Q., Bhuiya, M., Su, S., et al. (2015). Methylation mediated by an anthocyanin, o-methyltransferase, is involved in purple flower coloration in Paeonia. J. Exp. Bot. 66, 6563–6577. doi: 10.1093/jxb/erv365
Elibox, W., Umaharan, P. (2008). Inheritance of major spathe colors in Anthurium andraeanum Hort. is determined by three major genes. HortScience 43, 787–791. doi: 10.21273/HORTSCI.43.3.787
Fan, W., Li, B., Tian, H., Li, X., Ren, H., Zhou, Q. (2022). Metabolome and transcriptome analysis predicts metabolism of violet-red color change in Lilium bulbs. J. Sci. Food Agric. 102, 2903–2915. doi: 10.1002/jsfa.11631
Feinbaum, R. L., Ausubel, F. M. (1988). Transcriptional regulation of the Arabidopsis thaliana chalcone synthase gene. Mol. Cell. Biol. 8, 1985–1992. doi: 10.1128/mcb.8.5.1985-1992.1988
Fornes, O., Castro-Mondragon, J. A., Khan, A., van der Lee, R., Zhang, X., Richmond, P. A., et al. (2019). Jaspar 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 48, D87–D92. doi: 10.1093/nar/gkz1001
Franken, P., Niesbach-Klosgen, U., Weydemann, U., Marechal-Drouard, L., Saedler, H., Wienand, U. (1991). The duplicated chalcone synthase genes C2 and Whp (white pollen) of Zea mays are independently regulated; Evidence for translational control of whp expression by the anthocyanin intensifying gene in. EMBO J. 10, 2605–2612. doi: 10.1002/j.1460-2075.1991.tb07802.x
Gonzalez, A., Zhao, M., Leavitt, J. M., Lloyd, A. M. (2008). Regulation of the anthocyanin biosynthetic pathway by the TTG1/BHLH/MYB transcriptional complex in Arabidopsis seedlings. Plant J. 53, 814–827. doi: 10.1111/j.1365-313X.2007.03373.x
Gopaulchan, D., Lennon, A. M., Umaharan, P. (2013). Identification of reference genes for expression studies using quantitative RT-PCR in spathe tissue of Anthurium andraeanum (Hort.). Sci. Hortic. 153, 1–7. doi: 10.1016/j.scienta.2013.01.024
Gopaulchan, D., Umaharan, P., Lennon, A. M. (2014). A molecular assessment of the genetic model of spathe color inheritance in Anthurium andraeanum (Hort.). Planta 239, 695–705. doi: 10.1007/s00425-013-2007-9
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Grotewold, E., Peterson, T. (1994). Isolation and characterization of a maize gene encoding chalcone flavonone isomerase. Mol. Gen. Genet. 242, 1–8. doi: 10.1007/BF00277341
Halloran, J. M., Kuehnle, A. R. (1998). What do anthurium buyers want in their flowers? Results of a market survey. University of Hawaii CTAHR cooperative extension services series: EFS 27 (Honolulu: University of Hawaii Press), 1–4.
Han, Y., Vimolmangkang, S., Soria-Guerra, R. E., Korban, S. S. (2012). Introduction of apple ANR genes into tobacco inhibits expression of both CHI and DFR genes in flowers, leading to loss of anthocyanin. J. Exp. Bot. 63, 2437–2447. doi: 10.1093/jxb/err415
Helariutta, Y., Elomaa, P., Kotilainen, M., Seppanen, P., Teeri, T. H. (1993). Cloning of cDNA coding for dihydroflavonol-4-reductase (DFR) and characterization of DFR expression in the corollas of Gerbera hybrida var. Regina (compositae). Plant Mol. Biol. 22, 183–193. doi: 10.1007/BF00014927
Henny, R. J., Chen, J. (2003). “Cultivar development of ornamental foliage plants,” in Plant Breeding Reviews, ed. Janick, J. (Hoboken, NJ: John Wiley and Sons, Inc.) 23, 245–290. doi: 10.1002/9780470650226.CH6
Hua, X. (2014). Flower market research report in Chinese New Year in 2013. China Flowers Hortic. 13, 24–28.
Ithal, N., Reddy, A. R. (2004). Rice flavonoid pathway genes, OsDFR and OsANS, are induced by dehydration, high salt and ABA, and contain stress responsive promoter elements that interact with the transcription activator, OsC1-MYB. Plant Sci. 166, 1505–1513. doi: 10.1016/j.plantsci.2004.02.002
Jiang, S., Sun, Q., Zhang, T., Liu, W., Wang, N., Chen, X. (2021). MdMYB114 regulates anthocyanin biosynthesis and functions downstream of MdbZIP4-like in apple fruit. J. Plant Physiol. 257, 153353. doi: 10.1016/j.jplph.2020.153353
Kamemoto, H. S., Kuehnle, A. R. (1996). Breeding Anthuriums in Hawaii (Honolulu: University of Hawaii Press).
Kamemoto, H. S., Nakaone, H. Y. (1955). Improving anthuriums through breeding. Hawaii Farm Sci. 3, 4–5.
Kitamura, S., Akita, Y., Ishizaka, H., Narumi, I., Tanaka, A. (2012). Molecular characterization of an anthocyanin-related glutathione s-transferase gene in cyclamen. J. Plant Physiol. 169, 636–642. doi: 10.1016/j.jplph.2011.12.011
Lang, X., Li, N., Li, L., Zhang, S. (2019). Integrated metabolome and transcriptome analysis uncovers the role of anthocyanin metabolism in Michelia maudiae. Int. J. Genomics 2019, 4393905. doi: 10.1155/2019/4393905
Li, S. (2014). Transcriptional control of flavonoid biosynthesis: fine-tuning of the MYB-BHLH-WD40 (MBW) complex. Plant Signal. Behav. 9, e27522. doi: 10.4161/psb.27522
Li, B., Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323. doi: 10.1007/s00299-016-2025-8
Li, Z. T., Dhekney, S. A., Gray, D. J. (2011). Use of the VvMYBA1 gene for non-destructive quantification of promoter activity via color histogram analysis in grapevine (Vitis vinifera) and tobacco. Transgenic Res. 20, 1087–1097. doi: 10.1007/s11248-010-9482-6
Li, S., Li, J., Li, X., Guan, Y., Chen, M., Zhu, J. (2021). Comparative transcriptome analysis reveals molecular regulators underlying pluripotent cell induction and callus formation in Anthurium andraeanum “Alabama”. In Vitro Cell. Dev. Biol. Plant 57, 235–247. doi: 10.1007/s11627-020-10138-0
Li, X., Liu, J., Chang, Q., Zhou, Z., Han, R., Liang, Z. (2021). Antioxidant and antidiabetic activity of proanthocyanidins from Fagopyrum dibotrys. Molecules 26, 2417. doi: 10.3390/molecules26092417
Li, C., Qiu, J., Yang, G., Huang, S., Yin, J. (2016). Isolation and characterization of a R2R3-MYB transcription factor gene related to anthocyanin biosynthesis in the spathes of Anthurium andraeanum (Hort.). Plant Cell Rep. 35, 2151–2165. doi: 10.1007/s00299-016-2025-8
Li, C., Qiu, J., Huang, S., Yin, J., Yang, G. (2019). AaMYB3 interacts with AabHLH1 to regulate proanthocyanidin accumulation in Anthurium andraeanum (Hort.)-another strategy to modulate pigmentation. Hortic. Res. 6, 14–14. doi: 10.1038/s41438-018-0102-6
Li, C., Yu, W., Xu, J., Lu, X., Liu, Y. (2022). Anthocyanin biosynthesis induced by MYB transcription factors in plants. Int. J. Mol. Sci. 23, 11701. doi: 10.3390/ijms231911701
Lin, F., Chen, S., Lin, K., Chen, C., Yao, F., Zhong, L. (2022). Integrated small RNA profiling and degradome analysis of Anthurium andraeanum cultivars with different-colored spathes. J. Plant Res. 135, 609–626. doi: 10.1007/s10265-022-01394-9
Liu, H., Liu, Q., Chen, Y., Zhu, Y., Zhou, X., Li, B. (2023). Full-length transcriptome sequencing provides insights into flavonoid biosynthesis in Camellia nitidissima petals. Gene 850, 146924. doi: 10.1016/j.gene.2022.146924
Liu, X., Zhao, F., Liu, H., Xie, Y., Zhao, D., Li, C. (2021). Transcriptomics and metabolomics reveal the adaption of Akkermansia muciniphila to high mucin by regulating energy homeostasis. Sci. Rep. 11, 9073. doi: 10.1038/s41598-021-88397-z
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Manni, M., Berkeley, M. R., Seppey, M., Zdobnov, E. M. (2021). BUSCO: assessing genomic data quality and beyond. Curr. Protoc. 1, e323. doi: 10.1002/cpz1.323
Marco van Herk, I. M., van Koppen, I. M., Smeding, I. S., van der Elzen, I. C., van Rosmalen, N., van Dijk, I. J., et al. (1998). Cultivation Guide Anthurium (Bleiswijk: Anthura).
Mehrtens, F., Kranz, H., Bednarek, P., Weisshaar, B. (2005). The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 138, 1083–1096. doi: 10.1104/pp.104.058032
Mizutani, M., Ward, E., DiMaio, J., Ohta, D., Ryals, J., Sato, R. (1993). Molecular cloning and sequencing of a cDNA encoding mung bean cytochrome P450 (P450C4H) possessing cinnamate 4-hydroxylase activity. Biochem. Biophys. Res. Commun. 190, 875–880. doi: 10.1006/bbrc.1993.1130
Mueller, L. A., Goodman, C. D., Silady, R. A., Walbot, V. (2000). AN9, a petunia glutathione s-transferase required for anthocyanin sequestration, is a flavonoid-binding protein. Plant Physiol. 123, 1561–1570. doi: 10.1104/pp.123.4.1561
Nie, Y., Stürzenbaum, S. (2019). Proanthocyanidins of natural origin: molecular mechanisms and implications for lipid disorder and aging-associated diseases. Adv. Nutr. 10, 464–478. doi: 10.1093/advances/nmy118
Nisar, N., Li, L., Lu, S., Khin, N., Pogson, B. (2015). Carotenoid metabolism in plants. Mol. Plant 8 (1), 68–82. doi: 10.1016/j.molp.2014.12.007
Osorio-Guarín, J. A., Gopaulchan, D., Quanckenbush, C., Lennon, A. M., Umaharan, P., Cornejo, O. E. (2021). Comparative transcriptomic analysis reveals key components controlling spathe color in Anthurium andraeanum (Hort.). PloS One 16, e261364. doi: 10.1371/journal.pone.0261364
Peng, Y., Thrimawithana, A. H., Cooney, J. M., Jensen, D. J., Espley, R. V., Allan, A. C. (2020). The proanthocyanin-related transcription factors MYBC1 and WRKY44 regulate branch points in the kiwifruit anthocyanin pathway. Sci. Rep. 10, 14161. doi: 10.1038/s41598-020-70977-0
Pertea, G., Huang, X., Liang, F., Antonescu, V., Sultana, R., Karamycheva, S., et al. (2003). TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large est datasets. Bioinformatics 19, 651–652. doi: 10.1093/bioinformatics/btg034
Petroni, K., Tonelli, C. (2011). Recent advances on the regulation of anthocyanin synthesis in reproductive organs. Plant Sci. 181, 219–229. doi: 10.1016/j.plantsci.2011.05.009
Pucker, B. (2022). Automatic identification and annotation of MYB gene family members in plants. BMC Genomics 23, 220. doi: 10.1186/s12864-022-08452-5
Pucker, B., Reiher, F., Schilbert, H. M. (2020). Automatic identification of players in the flavonoid biosynthesis with application on the biomedicinal plant croton tiglium. Plants 9, 1103. doi: 10.3390/plants9091103
Pucker, B., Selmar, D. (2022). Biochemistry and molecular basis of intracellular flavonoid transport in plants. Plants 11, 963. doi: 10.3390/plants11070963
Ren, J., Liu, Z., Chen, W., Xu, H., Feng, H. (2019). Anthocyanin degrading and chlorophyll accumulation lead to the formation of bicolor leaf in ornamental kale. Int. J. Mol. Sci. 20, 603. doi: 10.3390/ijms20030603
Rikken, M. (2010). The European Market for Fair and Sustainable Flowers and Plants (Belgium: Trade for Development Centre, Belgian Development Agency).
Rodriguez-Concepcion, M., Boronat, A. (2002). Elucidation of the methylerythritol phosphate pathway for isoprenoid biosynthesis in bacteria and plastids. A metabolic milestone achieved through genomics. Plant Physiol. 130, 1079–1089. doi: 10.1104/pp.007138
Saito, K., Kobayashi, M., Gong, Z., Tanaka, Y., Yamazaki, M. (1999). Direct evidence for anthocyanidin synthase as a 2-oxoglutarate-dependent oxygenase: molecular cloning and functional expression of cDNA from a red forma of Perilla frutescens. Plant J. 17, 181–189. doi: 10.1046/j.1365-313x.1999.00365.x
Salvatierra, A., Pimentel, P., Moya-Leon, M. A., Herrera, R. (2013). Increased accumulation of anthocyanins in Fragaria chiloensis fruits by transient suppression of FcMYB1 gene. Phytochemistry 90, 25–36. doi: 10.1016/j.phytochem.2013.02.016
Sparvoli, F., Martin, C., Scienza, A., Gavazzi, G., Tonelli, C. (1994). Cloning and molecular analysis of structural genes involved in flavonoid and stilbene biosynthesis in grape (Vitis vinifera L.). Plant Mol. Biol. 24, 743–755. doi: 10.1007/BF00029856
Stracke, R., Ishihara, H., Huep, G., Barsch, A., Mehrtens, F., Niehaus, K., et al. (2007). Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 50, 660–677. doi: 10.1111/j.1365-313X.2007.03078.x
Sun, Y., Li, H., Huang, J. (2012). Arabidopsis TT19 functions as a carrier to transport anthocyanin from the cytosol to tonoplasts. Mol. Plant 5, 387–400. doi: 10.1093/mp/ssr110
Tanner, G. J., Francki, K. T., Abrahams, S., Watson, J. M., Larkin, P. J., Ashton, A. R. (2003). Proanthocyanidin biosynthesis in plants. Purification of legume leucoanthocyanidin reductase and molecular cloning of its cDNA. J. Biol. Chem. 278, 31647–31656. doi: 10.1074/jbc.M302783200
Teixeira da Silva, J. A., Dobránszki, J., Zeng, S., Winarto, B., Lennon, A. M., Jaufeerally-Fakim, Y., et al. (2015). Genetic transformation and molecular research in Anthurium: progress and prospects. Plant Cell Tissue Organ Cult. 123, 205–219. doi: 10.1007/s11240-015-0832-1
Thoben, C., Pucker, B. (2023). Automatic annotation of the bHLH gene family in plants. BioRxiv 05, 539087. doi: 10.1101/2023.05.02.539087
Timoneda, A., Feng, T., Sheehan, H., Walker-Hale, N., Pucker, B., Lopez-Nieves, S., et al. (2019). The evolution of betalain biosynthesis in caryophyllales. New Phytol. 224, 71–85. doi: 10.1111/nph.15980
Tomas, M. (2022). Effect of dietary fiber addition on the content and in vitro bioaccessibility of antioxidants in red raspberry puree. Food Chem. 375, 131897. doi: 10.1016/j.foodchem.2021.131897
Vaculík, M., Pavlovič, A., Lux, A. (2015). Silicon alleviates cadmium toxicity by enhanced photosynthetic rate and modified bundle sheath’s cell chloroplasts ultrastructure in maize. Ecotoxicol. Environ. Saf. 120, 66–73. doi: 10.1016/j.ecoenv.2015.05.026
Xie, D. Y., Sharma, S. B., Dixon, R. A. (2004). Anthocyanidin reductases from Medicago truncatula and Arabidopsis thaliana. Arch. Biochem. Biophys. 422, 91–102. doi: 10.1016/j.abb.2003.12.011
Yan, H., Pei, X., Zhang, H., Li, X., Zhang, X., Zhao, M., et al. (2021). MYB-mediated regulation of anthocyanin biosynthesis. Int. J. Mol. Sci. 22, 3103. doi: 10.3390/ijms22063103
Yang, J., Chen, Y., Xiao, Z., Shen, H., Li, Y., Wang, Y. (2022). Multilevel regulation of anthocyanin-promoting R2R3-MYB transcription factors in plants. Front. Plant Sci. 13. doi: 10.3389/fgene.2022.936051
Zhao, J., Li, Z. T., Chen, J., Henny, R. J., Gray, D. J., Chen, J. (2013). Purple-leaved Ficus lyrata plants produced by overexpressing a grapevine VvMybA1 gene. Plant Cell Rep. 32, 1783–1793. doi: 10.1007/s00299-013-1491-5
Zhou, S., Hu, Z., Zhu, M., Zhang, B., Deng, L., Pan, Y., et al. (2013). Biochemical and molecular analysis of a temperature-sensitive albino mutant in kale named “White dove”. Plant Growth Regul. 71, 281–294. doi: 10.1007/s10725-013-9829-0
Keywords: anthocyanin, Anthurium andraeanum, spathe color, flavonoid, chlorophyll, MYB transcription factor
Citation: Li J, Tan Q, Yi M, Yu Z, Xia Q, Zheng L, Chen J, Zhou X, Zhang X-Q and Guo H-R (2023) Identification of key genes responsible for green and white colored spathes in Anthurium andraeanum (Hort.). Front. Plant Sci. 14:1208226. doi: 10.3389/fpls.2023.1208226
Received: 18 April 2023; Accepted: 17 August 2023;
Published: 06 September 2023.
Edited by:
He Huang, Beijing Forestry University, ChinaCopyright © 2023 Li, Tan, Yi, Yu, Xia, Zheng, Chen, Zhou, Zhang and Guo. 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.
*Correspondence: Xiaoyun Zhou, xyzh28@126.com; Xiang-Qian Zhang, aacrav@163.com; He-Rong Guo, guohe-rong@scau.edu.cn
†These authors have contributed equally to this work