Integrated analysis of the transcriptome and metabolome of purple and green leaves of Tetrastigma hemsleyanum reveals gene expression patterns involved in anthocyanin biosynthesis.

To gain better insight into the regulatory networks of anthocyanin biosynthesis, an integrated analysis of the metabolome and transcriptome in purple and green leaves of Tetrastigma hemsleyanum was conducted. Transcript and metabolite profiles were archived by RNA-sequencing data analysis and LC-ESI-MS/MS, respectively. There were 209 metabolites and 4211 transcripts that were differentially expressed between purple and green leaves. Correlation tests of anthocyanin contents and transcriptional changes showed 141 significant correlations (Pearson correlation coefficient >0.8) between 16 compounds and 14 transcripts involved in the anthocyanin biosynthesis pathway. Some novel genes and metabolites were discovered as potential candidate targets for the improvement of anthocyanin content and superior cultivars.

Introduction

Tetrastigma hemsleyanum Diels et Gilg (T. hemsleyanum, belonging to the family Vitaceae), also known as “Sanyeqing” (SYQ), is distributed in tropical to subtropical areas in Asia, mainly in southern China, such as Zhejiang, Guizhou and Guangxi provinces. The entire herb and its root tubers have been used as a broad-spectrum antibiotic material for the treatment of fever and sore throat in China for a long time. Previous findings have demonstrated that the principal and functional components of SYQ, such as polysaccharides and polyphenols are beneficial for health [1,2].

Flavonoids such as anthocyanins and flavanols are the most dominant and important components in SYQ [3,4]. These compounds are beneficial for human health and are an important group of pigments that colour the leaves and flowers of many plants [5]. Anthocyanins are also involved in biotic and abiotic stress defence responses [6,7]. It is vital to elucidate the flavonoid biosynthesis and regulatory pathways in SYQ. Although several genes in the anthocyanin regulation pathway have been identified [8], the unique regulation mechanism in SYQ remains unclear.

The colour of plant leaves are mainly determined by the composition and concentration of anthocyanins [9], products of the branched flavonoid biosynthesis pathway. The accumulation of anthocyanins is regulated by the expression of multiple genes, including chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3'H), flavonoid 3′,5′-hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS) and flavonoid 3-O-glucosyltransferase (UFGT)[10].

In the past decade, RNA sequencing has rapidly become an efficient approach to analyse the function of genes in a high-throughput way [11]. Transcripts with low abundance and unknown transcripts cannot be identified [12]. Liquid chromatography / mass spectrometry (LC/MS) has advanced metabolomics by enabling the discovery ofa large number of compounds compared with traditional chemical analysis [13,14]. Recently, the integrated analysis of metabolic and gene expression has been widely used in the exploitation network and correlation between metabolites and genes [15-17].

In this study, we explored the regulatory networks of flavonoid and anthocyanin biosynthesis in green (RG) and purple leaves (PL) of SYQ at the transcriptome and metabolome levels. The accumulation of different types of anthocyanin and the expression of related genes were investigated. A connection network was constructed to highlight the regulatory genes associated with specific metabolites. Our findings provide insights into the accumulation mechanism of SYQ leaf colour pigments and the regulation of anthocyanin biosynthesis.

Materials and methods

Plant materials and growth conditions

The green (RG) and purple leaf (PL) genotypes (S1 Fig) of Tetrastigma hemsleyanum (Sanyeqing) were collected from Guangxi and Hunan provinces, respectively. They were grown in the plant garden of Hangzhou Academy of Agricultural Sciences (Hangzhou, Zhejiang Province, China). Leaves at the third node away from the top were collected, frozen in liquid nitrogen and stored at -80°C for RNA and metabolite isolation.

RNA isolation and transcriptome sequencing

Total RNA was isolated with TRIzol reagent (Invitrogen, USA) according to the manufacturer's protocol. The transcriptome sequence library was constructed using NEBNext Ultra RNA Library Prep Kits for Illumina (NEB, USA). Sequencing was performed on an Illumina HiSeq 2500 platform (Novogene, China). The reads were aligned to the genome using TopHat (2.0.9) software after removing the reads containing adapter or poly-N and low-quality reads from the raw data. Total number of reads per kilobase per million reads (RPKM) of each gene was calculated based on the length of the gene and the counts of reads mapped to this gene. GO annotation was implemented using Blast2GO software. KOBAS (2.0) software was used for KEGG enrichment analysis of differentially expressed genes.

Metabolic profiling

The samples were freeze-dried and crushed into power for metabolite isolation. One hundred milligrams of each sample were extracted with 1.0 ml 70% aqueous methanol at 4°C overnight. The extracts were absorbed (CNWBOND Carbon-GCB SPE Cartridge, 250 mg, 3 ml; ANPEL, Shanghai, China, www.anpel.com.cn/cnw) and filtered (SCAA-104, 0.22μm pore size; ANPEL, Shanghai, China, http://www.anpel.com.cn/) before LC-MS analysis.

The sample extracts were analysed using an LC-ESI-MS/MS system (HPLC, Shim-pack UFLC SHIMADZU CBM30A system, www.shimadzu.com.cn; MS, Applied Biosystems 6500 Q TRAP, www.appliedbiosystems.com.cn). A Waters ACQUITY UPLC HSS T3 C18 (1.8 μm, 2.1 mm*100 mm) was used for compound separation. The analytical conditions were set as follows: solvent system, water (0.04% acetic acid): acetonitrile (0.04% acetic acid); gradient program,95:5V/V at 0min, 5:95V/V at 11.0min, 5:95V/V at 12.0min, 95:5V/V at 12.1min, 95:5V/V at 15.0 min; flow rate, 0.40 ml/min; temperature, 40°C; injection volume, 2 μl. Mass data acquisition was performed in both positive and negative modes using the following parameters: ion source, turbo spray; source temperature, 500°C; ion spray voltage (IS), 5500 V; ion source gas I (GSI), gas II(GSII), and curtain gas (CUR) of 55, 60, and 25.0 psi, respectively; collision gas(CAD) was high. The mass fragmentations were compared to the HMDB (http://www.hmdb.ca), METLIN (http://metlin.scrippps.edu) and KEGG (http://kegg.jp) databases. The obtained data were used by SIMCA-P V12.0.0 Demo (Umetric, Umea, Sweden) for principal components analysis (PCA) and partial least-squares discriminant analysis (PLS-DA).

Verification of candidate genes by quantitative real-time PCR (qRT-PCR)

The expression levels of transcripts involved in anthocyanin biosynthesis pathways were validated by qPCR using the same RNA samples used for sequencing. The primers were designed using beacon designer 7.8 software, and cDNAs were synthesized using SuperScript™ III First-Strand Synthesis SuperMix for qRT-PCR (Invitrogen). Power SYBR® Green PCR Master Mix (Applied Biosystems) was selected for the identification of the PCR products on a CFX384 Real-time PCR system (Bio-Rad). Three replicates were performed for each sample, and the primers of each sample are listed in S1 Table.

Results

Metabolic differences in green leaves and purple leaves

To compare the metabolite composition in the two accessions, LC-MS was used for metabolite identification. In total, 597 metabolites with known structures were detected in RG and PL. Differentially expressed proteins were defined as those with a variable importance for projection (VIP) value >1.5 compared to RG (Fig 1). Among the detected metabolites, 111 metabolites were upregulated and 98 metabolites were downregulated. The top 10 most differentially accumulated metabolites are listed in Table 1. The top 10 upregulated metabolites were mainly anthocyanins and flavone C-glycosides. Compared to RG, cyanidin 3-O-rutinoside showed the highest fold change in PL, with a log2 fold change (FC) (PL/RG) value of 20.59. Of the 10 most downregulated metabolites, quinic acid O-glucuronic acid had the lowest log2 fold change (FC) (PL/RG) value of -21.73.

Fig 1

VIP (variable importance for projection) plot of differentially accumulated metabolites in PL vs RG.

Red dots are upregulated metabolites, and green dots are downregulated metabolites.

Table 1

List of top 10 up- and downregulation metabolites in PL vs RG.

Differentially expressed proteins were defined as those with a variable importance for projection (VIP) value >1.5.

Index	Compounds	Class	RG1	RG2	RG3	PL1	PL2	PL3	VIP	Fold_Change	Log2FC
pme1773	Cyanidin 3-O-rutinoside (Keracyanin)	Anthocyanins	9	9	9	19500000	11500000	11700000	1.363135	1581481.481	20.59285
pmb0675	C-pentosyl-apigenin O-p-coumaroylhexoside	Flavone C-glycosides	9	9	9	4880000	18500000	819000	1.351613	896259.2593	19.77356
pme3038	5-oxoproline	Amino acid derivatives	9	9	9	2880000	7300000	2470000	1.362247	468518.5185	18.83775
pme1793	Pelargonin	Anthocyanins	9	9	9	2760000	2500000	3770000	1.363342	334444.4444	18.35141
pmb0837	Procyanidin A3	Proanthocyanidins	9	9	9	3160000	1610000	2300000	1.362756	261851.8519	17.99839
pmb3074	3-O-p-Coumaroyl quinic acid	Quinate and its derivatives	9	9	9	1450000	1170000	1030000	1.363438	135185.1852	17.04458
pme3300	Tricetin	Flavone	9	9	9	2510000	311000	819000	1.355341	134814.8148	17.04062
pmb0808	6,7-dihydroxycoumarin 6-O-quinic acid	Coumarins	9	9	9	1450000	9	1460000	0.958286	107778.1111	16.7177
pme1816	Neochlorogenic acid (5-O-Caffeoylquinic acid)	Quinate and its derivatives	9	9	9	773000	868000	1050000	1.363533	99666.66667	16.60482
pmb0749	O-Feruloyl quinic acid	Quinate and its derivatives	9	9	9	173000	546000	1510000	1.355143	82555.55556	16.33308
pmb3058	Quinic acid O-glucuronic acid	Quinate and its derivatives	33800000	30800000	29200000	9	9	9	1.363641	2.88E-07	-21.7282
pme0398	Chlorogenic acid (3-O-Caffeoylquinic acid)	Quinate and its derivatives	22800000	17200000	16900000	9	9	9	1.363516	4.75E-07	-21.007
pmb0639	8-C-hexosyl-apigenin O-hexosyl-O-hexoside	Flavone C-glycosides	14500000	8060000	24600000	9	9	9	1.362386	5.73E-07	-20.7362
pma6647	C-hexosyl-chrysoeriol O-hexoside	Flavone C-glycosides	6340000	5840000	6310000	9	9	9	1.363675	1.46E-06	-19.3854
pmb0652	C-hexosyl-apigenin O-pentoside	Flavone C-glycosides	6300000	5470000	5240000	9	9	9	1.363606	1.59E-06	-19.265
pmb1108	Luteolin 6-C-hexoside 8-C-hexosyl-O-hexoside	Flavone C-glycosides	1530000	382000	11500000	9	9	9	1.346659	2.01E-06	-18.9221
pmb0623	6-C-hexosyl chrysoeriol O-hexoside	Flavone C-glycosides	2360000	1200000	5440000	9	9	9	1.360537	0.000003	-18.3466
pma6516	C-hexosyl-apigenin O-hexosyl-O-hexoside	Flavone C-glycosides	2710000	1850000	2170000	9	9	9	1.363427	4.01E-06	-17.9273
pmb0629	Chrysoeriol 6-C-hexoside	Flavone C-glycosides	2590000	1940000	1210000	9	9	9	1.362641	4.70E-06	-17.6977
pmb4777	4-Hydroxy-7-methoxycoumarin-beta-rhamnoside	Coumarins	584000	856000	851000	9	9	9	1.363425	1.18E-05	-16.3727

All differentially accumulated metabolites were subjected to KEGG analysis (Fig 2).Metabolites participating in flavonoid biosynthesis, flavone and flavanol biosynthesis, tryptophan metabolism and biosynthesis of alkaloids derived from the shikimate pathway were predominantly enriched.

Fig 2

KEGG pathway enrichment analysis based on the differentially accumulated metabolites in SYQ purple and green leaves.

Principal component analysis (PCA) was employed to identify the differences in metabolite profiles among samples (Fig 3). The results showed that the first principal component (PC1, 52.98% of the total variables) was clearly separated in the RG and PL samples, indicating that the accumulation patterns of metabolites were different in green and purple leaves.

Fig 3

PCA score plot of SYQ purple and green leaves.

Transcriptome analysis in green leaves and purple leaves

After removing the low-quality reads, approximately 30,000,000 reads of each sample were generated. De novo transcriptome assembly using Trinity built 205,387 transcripts and 100,540 unigenes for the samples, with an average transcript length of 1151 bp. The unigenes were mapped to the NR database and then applied to the Pfam database for annotation.

Using a false discovery rate (FDR) <0.01, log2FC>1 as threshold values, 4211 genes were found to be differentially expressed in PL vs. RG. Among them, 2035 genes were upregulated and 2176 genes were downregulated (Fig 4).

Fig 4

Volcano plot of differentially regulated genes in PL vs RG.

The red dots are upregulated genes and green dots are downregulated genes.

Differentially expressed genes potentially involved in anthocyanin biosynthesis were identified. These included CHS, CHI, F3H, F3'H, F3'5'H, DFR, ANS and UFGT.

The enriched GO (Gene Ontology) terms for DEGs were analysed (Fig 5). In biology process (BP) categories, the most significant terms were metabolic process, cellular process, and single-organism process. In molecular function (MF) and cellular component (CC) categories, binding and catalytic activity, and cell and cell part were the most enriched, respectively.

Fig 5

GO classification of all DEGs in purple and green SYQ leaves.

Upregulated genes (A) and downregulated genes (B).

The identified genes were mapped to 105 KEGG reference pathways (Fig 6). These pathways included flavonoid biosynthesis (ko00941), purine metabolism (ko00230) and so on, and plant pathogen interaction (ko04626) and carbon metabolism (ko01200) had the highest enrichment factor values. In the flavonoid biosynthesis pathway, most genes were upregulated in purple leaves.

Fig 6

KEGG pathway enrichment analysis based on all DEGs of SYQ purple and green leaves.

Upregulated genes (A) and downregulated genes (B).

Transcription factors are essential regulators in anthocyanin biosynthesis. A total of 176 TFs were identified as either up- or downregulated between purples and green leaves. Of the most extensively studied TFs, one MYB and six bHLH were significantly different between purple and green leaves. MYB was highly expressed in purple leaves, four bHLH were upregulated, and two were downregulated.

Network analysis of metabolites and transcripts in SYQ

The identified anthocyanins, their relevant compounds and the related genes were mapped onto their corresponding position in anthocyanin pathway. The results indicated that the compounds in this pathway were different in purple and green leaves. The abundances of the transcripts and composition of the compounds are shown via heatmap (Fig 7).

Fig 7

Biosynthetic pathway of anthocyanins.

This pathway is constructed based on the KEGG pathway. Each coloured cell represents the normalized intensity of each compound ion according to the colour scale (three biological replicates×two cultivars, n = 6). The expression levels of chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3'H), flavonoid 3′,5′-hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS) and flavonoid 3-O-glucosyltransferase (UFGT) are shown.

As shown in Fig 7, pelargonidin, dihydrokaempferol, and pelargonidin 3-O-glucoside were significantly accumulated in purple leaves, and dihydromyricetin, cyanidin, and cyanidin 3-glucoside had higher levels in green leaves. The majority of the transcripts were promoted in purple leaves, although the expression levels of CHS and CHI were inhibited. In our results, the expression levels of the transcripts and metabolites in the pelargonidin biosynthesis pathway were increased.

To gain a better understanding of the regulatory network of the leaf colour in the two accessions, the Pearson correlation test was performed for the metabolites and the transcripts. In total, 14 genes and 16 metabolites involved in flavonoid and anthocyanin pathways were subjected to Pearson correlation analysis (Fig 8). There were 141 significant correlation combinations between the genes and metabolites that had a Pearson correlation coefficient >0.8 and p value<0.05.

Fig 8

Connection network between regulatory genes and anthocyanin-related metabolites.

The networks between metabolites and transcripts were visualized with Cytoscape software (version 2.8.2).

The expression levels of the key genes in the anthocyanin biosynthesis pathway were analysed by real-time quantitative PCR (qPCR). As shown in Fig 9, UGFT, ANS, DFR, F3’5’H, and F3H accumulated at a much higher level in the purple leaves, whereas the expression of CHI was reduced in purple leaves. The expression pattern of the genes involved in the anthocyanin biosynthesis pathway revealed by qPCR was consistent with the RNA-seq results.

Fig 9

qPCR analysis of anthocyanin biosynthesis-related genes in green- and purple-leaf SYQ cultivars.

Discussion

In this study, the differences in transcripts and metabolites in green and purple leaves were compared. Using ultra-performance liquid chromatography and tandem mass spectrometry, 597 metabolites were detected (Fig 1). As shown by hierarchical sample clustering, there was a clear separation between purple and green leaves (S2 Fig). In the differential accumulated metabolites, six metabolites annotated to anthocyanins were up-accumulated in purple leaves, while only two anthocyanin compounds were inhibited.

To investigate the genes regulating the differential pigmentation of SYQ leaves, transcriptome profiling was engaged. Anthocyanin accumulation corresponds to the expression of the genes involved in the biosynthesis pathway [5,18]. Our study showed that most of the key genes in the anthocyanin biosynthesis pathway were upregulated in purple leaves (Fig 7), indicating that the anthocyanin pathway was promoted in purple leaves. This finding was consistent with the gene expression results of the key genes in anthocyanin biosynthesis pathway (Fig 9). The largest log2FC (purple/green) value was found in F3H, which had a value of 5.869. F3H catalyses the conversion of flavanone into dihydroflavanol[19]. The upregulation of F3H facilitates the accumulation of anthocyanins.

There are many types of anthocyanins; cyanidin, delphinidin and pelargonidin have been recognized as the most common anthocyanins in plants. Different distribution of anthocyanins results in different leaf colours in fruits and leaves. For instance, the major type of anthocyanins in berries is cyanidin [20]; red pigment is generated from the accumulation of pelargonidin [21], and delphinidin makes plants look purple. In our results, the biosynthesis pathway of pelargonidin-3-glucoside was promoted at both the transcriptional and metabolic levels. The genes and metabolites work together to regulate the production of anthocyanins.

Transcription factors including MYB, bHLH and WD40 play vital roles in the regulation of flavonoid and anthocyanin pathway [22-24]. Only MYB had a higher expression level in purple leaves, and different bHLHs transcripts were differentially expressed in different leaves. There are 8 WRKYs that showed higher expression levels in purple leaves. As suggested previously, WRKY may be involved in the regulation of anthocyanin biosynthesis. More anthocyanin was accumulated in the WRKY overexpressed lines compared to wild type in Arabidopsis [25].

Polysaccharides, phenolic acids and flavanols are the main compounds in SYQ with biological activities that can be used to cure some diseases. Although there are abundant bioactive compounds in SYQ root tubers [26], the shootsmay be more practical for use due to their easy access. In our results, the composition of the bioactive compounds in purple leaves was higher than in green leaves. Purple leaves may contain higher levels of antioxidant activity compared to green leaves. Therefore, purple leaves with higher anthocyanin contents are a superior resource for the improvement of SYQ quality.

Network analysis of the metabolites and transcripts was conducted. There were 141 significant correlation combinations between 14 genes and 16 metabolites. Transcript c90950.graph_c0 (CHS) was correlated with 14 metabolites, the most of those investigated (Fig 8). CHS is the initial enzyme of flavonoid biosynthesis; it catalyses the synthesis of naringenin chalcone from 4-coumaroyl-CoA and malonyl-CoA. In different plant species, the expression pattern of CHS of red tissues and green tissues are different [27]. The expression level of CHS was found to be significantly upregulated in red tomato fruits, while CHS expression was similar between pigmented and non-pigmented tomato mutants [28]. Data obtained indicated that the metabolites and genes work co-ordinately to regulate anthocyanin biosynthesis. Regulatory genes that were highly correlated with the accumulation of metabolites were identified; these could provide new insights into the regulatory mechanism of anthocyanin biosynthesis.

Conclusions

We analysed the regulatory network of anthocyanin biosynthesis integrating the metabolome and transcriptome in the purple and green leaves of SYQ. Correlation analysis of the metabolites and transcripts involved in anthocyanin biosynthesis pathway was conducted. A regulatory network of metabolites and transcripts related to leaf colour could be a resource for the exploitation of the mechanism of anthocyanin regulation in SYQ.

Phenotypes of purple-leaf SYQ (A-PL) and green-leaf SYQ (B-RG).

Red and yellow arrows indicate the leaves shown in S1 Fig. Yellow arrows indicate the position where the leaves were collected for further experiments.

(TIFF)

Click here for additional data file.

Hierarchical clustering of differentially accumulated metabolites in PL and RG sample.

Normalized contents of different compounds are represented as colours ranging from green (-2) to red (2).

(TIFF)

Click here for additional data file.

Sequence information for primers.

(XLSX)

Click here for additional data file.

Gene accession numbers and sequences of DEGs.

(XLSX)

Click here for additional data file.

2 Jan 2020

PONE-D-19-31044

Integrated analysis of transcriptome and metabolome reveals alterations in expression of genes involved in anthocyanin biosynthesis in purple and green leaves of Tetrastigma hemsleyanum

PLOS ONE

Dear Mr. Ruan,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

We would appreciate receiving your revised manuscript by Feb 16 2020 11:59PM. When you are ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter.

To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

Please include the following items when submitting your revised manuscript:

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). This letter should be uploaded as separate file and labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. This file should be uploaded as separate file and labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. This file should be uploaded as separate file and labeled 'Manuscript'.

Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out.

We look forward to receiving your revised manuscript.

Kind regards,

Zhong-Hua Chen, Ph.D.

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

http://www.journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and http://www.journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. We suggest you thoroughly copyedit your manuscript for language usage, spelling, and grammar. If you do not know anyone who can help you do this, you may wish to consider employing a professional scientific editing service.

Whilst you may use any professional scientific editing service of your choice, PLOS has partnered with both American Journal Experts (AJE) and Editage to provide discounted services to PLOS authors. Both organizations have experience helping authors meet PLOS guidelines and can provide language editing, translation, manuscript formatting, and figure formatting to ensure your manuscript meets our submission guidelines. To take advantage of our partnership with AJE, visit the AJE website (http://learn.aje.com/plos/) for a 15% discount off AJE services. To take advantage of our partnership with Editage, visit the Editage website (www.editage.com) and enter referral code PLOSEDIT for a 15% discount off Editage services.  If the PLOS editorial team finds any language issues in text that either AJE or Editage has edited, the service provider will re-edit the text for free.

Upon resubmission, please provide the following:

The name of the colleague or the details of the professional service that edited your manuscript

A copy of your manuscript showing your changes by either highlighting them or using track changes (uploaded as a *supporting information* file)

A clean copy of the edited manuscript (uploaded as the new *manuscript* file)

3. In the Introduction to your manuscript, please remove any claims regarding the medicinal benefits of Tetrastigma hemsleyanum extracts which are not substantiated by robust scientific and clinical evidence.

4. In your Data Availability statement, you have not specified where the minimal data set underlying the results described in your manuscript can be found. PLOS defines a study's minimal data set as the underlying data used to reach the conclusions drawn in the manuscript and any additional data required to replicate the reported study findings in their entirety. All PLOS journals require that the minimal data set be made fully available. For more information about our data policy, please see http://journals.plos.org/plosone/s/data-availability.

Upon re-submitting your revised manuscript, please upload your study’s minimal underlying data set as either Supporting Information files or to a stable, public repository and include the relevant URLs, DOIs, or accession numbers within your revised cover letter. For a list of acceptable repositories, please see http://journals.plos.org/plosone/s/data-availability#loc-recommended-repositories. Any potentially identifying patient information must be fully anonymized.

Important: If there are ethical or legal restrictions to sharing your data publicly, please explain these restrictions in detail. Please see our guidelines for more information on what we consider unacceptable restrictions to publicly sharing data: http://journals.plos.org/plosone/s/data-availability#loc-unacceptable-data-access-restrictions. Note that it is not acceptable for the authors to be the sole named individuals responsible for ensuring data access.

We will update your Data Availability statement to reflect the information you provide in your cover letter.

5. PLOS requires an ORCID iD for the corresponding author in Editorial Manager on papers submitted after December 6th, 2016. Please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field. This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager. Please see the following video for instructions on linking an ORCID iD to your Editorial Manager account: https://www.youtube.com/watch?v=_xcclfuvtxQ

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: Please modify the manuscript according to the revised MS.

The study described in their manuscript, authors explored the regulatory networks of anthocyanin biosynthesis from the metabolome and transcriptome, the plant materials were very interesting with the green and purple leaves. Transcript and metabolite profiles were integrated to find the differential metabolites and genes, and their correlation for anthocyanin biosysthesis.

The author should discribe the background of the two accessions of SYQ, and their ralation.

I comment on related issues and hope the authors find these helpful to polish their work before final publication.

Reviewer #2: This manuscript applied metabolome and transcriptome to reveal the alteration of regulatory network in anthocyanin biosynthesis in SYQ. It is meaningful to discover some novel genes and metabolites for the improvement of anthocyanin content and superior cultivars. However, there are some problems on writing of this manuscript.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files to be viewed.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org. Please note that Supporting Information files do not need this step.

Submitted filename: MS revised.docx

Click here for additional data file.

Submitted filename: reviews-Zhigang Han.docx

Click here for additional data file.

22 Jan 2020

Responses to reviewers’ comments and suggestion

First of all, we are deeply indebted to your careful examination and valuable comments for further improvement of this manuscript. Accordingly, we revised the MS according to your comments and suggestion. The detail revisions are as follows.

From Vicky Stabler:

1.Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming.

Au:Thank you for your comment.We have re-checked the style andfile naming of our manuscript.

2.We suggest you thoroughly copyedit your manuscript for language usage, spelling, and grammar. If you do not know anyone who can help you do this, you may wish to consider employing a professional scientific editing service.

Au: We employed American Journal Experts (AJE) to polish our manuscript. The new version was submitted.

3.In the Introduction to your manuscript, please remove any claims regarding the medicinal benefits of Tetrastigmahemsleyanum extracts which are not substantiated by robust scientific and clinical evidence.

Au: Relevant information was deleted in the revised version.

4.In your Data Availability statement, you have not specified where the minimal data set underlying the results described in your manuscript can be found.

Au: A new supplementary table (Table S2) of the sequences of all DEGs was added.

5.PLOS requires an ORCID iD for the corresponding author in Editorial Manager on papers submitted after December 6th, 2016.

Au: The ORCID iD was validated in the submission system.

Reviewer #1: Please modify the manuscript according to the revised MS.

1.The author should describe the background of the two accessions of SYQ, and their relation.

Au: Relevant information was added in the revised version.

2. I comment on related issues and hope the authors find these helpful to polish their work before final publication.

Au: Thank you for your careful examination and valuable comments. In the revised version, the spelling and grammatical errors were corrected.

Reviewer #2: This manuscript applied metabolome and transcriptome to reveal the alteration of regulatory network in anthocyanin biosynthesis in SYQ. It is meaningful to discover some novel genes and metabolites for the improvement of anthocyanin content and superior cultivars. However, there are some problems on writing of this manuscript.

-Minor issues

1. Line 35-36, “Previous findings have demonstrated that the principal and functional components of SYQ, such as flavonoids, polysaccharide and polyphenol, have various effect on anti-cancer, liver protection, and anti-inflammatory”, I think polyphenol include flavonoids.

Au:Thank you for your comment. It was corrected in revised version.

Line 40, “Anthocyanins also involved in biotic and abiotic stress defending responses”, rephrase this sentence.

Au: We rewrote it in the revised version.

Line 43-44, “The color of the plants leaves is mainly determined by the composition and amount of the anthocyanins, the anthocyanins are the products of the branched flavonoid biosynthesis pathway”, “amount” is not proper in this sentence.

Au:Sorry we made an error. “amount” was changed to “concentration” in the revised version.

Line 50, “The transcripts in low abundance or some unrevealed ever can be identified”, rephrase this sentence.

Au: This sentence was rewritten in revised version.

Line 51-52, “Advanced in metabolomics with the help of liquid chromatography / mass spectrometry (LC/MS) has been a powerful tool to discover a large number of compounds compared with traditional chemical analysis”, “Advanced in metabolomics”? Rephrase this sentence.

Au:Thank you for your comment. “Advanced” was changed to “Advances” in the revised version.

Line 203, check and rephrase this sentence.

Au: Relevant correction was made in revised version.

Line 220-221, “Although the bioactive compounds in root tubers of SYQ is abundant, the aerial parts of SYQ may be more practical due to its easily access”, what are the aerial parts? And why are these parts easily access?

Au: “aerial parts” was changed to “shoots” in the revised version.

Line 223-224, check and rephrase this sentence.

Au: The correction was made in the revised version.

Submitted filename: Responses to reviewers.docx

Click here for additional data file.

24 Feb 2020

Integrated analysis of the transcriptome and metabolome of purple and green leaves of Tetrastigma hemsleyanum reveals gene expression patterns involved in anthocyanin biosynthesis

PONE-D-19-31044R1

Dear Dr. Ruan,

We are pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it complies with all outstanding technical requirements.

Within one week, you will receive an e-mail containing information on the amendments required prior to publication. When all required modifications have been addressed, you will receive a formal acceptance letter and your manuscript will proceed to our production department and be scheduled for publication.

Shortly after the formal acceptance letter is sent, an invoice for payment will follow. To ensure an efficient production and billing process, please log into Editorial Manager at https://www.editorialmanager.com/pone/, click the "Update My Information" link at the top of the page, and update your user information. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, you must inform our press team as soon as possible and no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

With kind regards,

Zhong-Hua Chen, Ph.D.

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #2: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: (No Response)

Reviewer #2: (No Response)

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

26 Feb 2020

PONE-D-19-31044R1

Integrated analysis of the transcriptome and metabolome of purple and green leaves of Tetrastigma hemsleyanum reveals gene expression patterns involved in anthocyanin biosynthesis

Dear Dr. Ruan:

I am pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please notify them about your upcoming paper at this point, to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

For any other questions or concerns, please email plosone@plos.org.

Thank you for submitting your work to PLOS ONE.

With kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Zhong-Hua Chen

Academic Editor

PLOS ONE