Differential gene expression in leaves and roots of Hydrangea serrata treated with aluminium chloride

Hydrangea serrata, also knowen as the Japanese tea hortensia, is known for its sweet taste and health properties of bevarages produced from this plant. The H. serrata 3,4-dihydroisocoumarins, hydrangenol and phyllodulcin harbour a variety of biological activities and pharmacological properties. Therefore, a detailed understanding of dihydroisocoumarin biosynthesis in H. serrata is of major interest. Their biosynthesis is assumed to be enhanced by elicitors and mediated by polyketide synthases like in cases of phenylpropanoid derived phytoalexins. A de-novo transcriptome assembly of leaves and roots from the aluminium chloride treatment group versus the control group alongside with annotation was generated. Secondary plant metabolites were analysed by LC-MS. It revealed that a terpene synthase and a triterpenoid synthase gene as well as lignin biosynthesis encoding genes were upregulated in roots. Many genes for transporters, glycosyl, and other transferases as well as glycosylases were found to be differentially expressed in both organs. As no differentially expressed polyketide synthase gene homolog was found, the relative leaf and root 3,4-dihydroisocoumarin content was analysed by LC-MS measurement. Although Hydrangea species are known for their aluminium detoxification using phenylpropanoid-derived compounds, the levels of 3,4- dihydroisocoumarins were not enhanced. In this metabolite analysis, an organ- specific accumulation profile of hydrangenol, phyllodulcin, hydrangeic acid and their mono- and di-glycosides was figured out.


Introduction
The Japanese mountain Hydrangea (Hydrangea serrata) is an important cultural and medical plant due to the content of specialised metabolites.A traditional tea is prepared from the leaves.Health properties of H. serrata "Oamacha" and closely related cultivars are mainly due to 3,4-dihydroisocoumarins (DHICs) including hydrangenol (HYD), phyllodulcin (PHY) and stilbene carboxylic acids like hydrangeic acid (Hashimoto et al., 1979;Liu et al., 2013).They share similar positive health effects and anti-microbial properties with the group of isocoumarins (Saddiqa et al., 2017).In contrast to the mostly bitter-tasting isocoumarins, PHY tastes sweet and is traditionally consumed in tea and beverages (Yamato and Hashigaki, 1979).Within the genus Hydrangea, H. serrata the varieties "Oamacha" and "Amagi Amacha" (classified as a subspecies of H. macrophylla Thunbergii) is the only natural sweetener-producing member (Çicȩk et al., 2018).The bioconversion of PHY-glycoside to the aglycon as high-intensity sweetener from Hydrangea macrophylla Thunbergii leaves is established (Jung et al., 2016).
PHY and HYD occur as glycosides or aglycons and have a structure and synthesis similar to stilbenoids or flavonoids (Asahina and Asano, 1929;Kućand Rush, 1985;Liu et al., 2013).HYD and PHY are proposed to be derived from hydrangeic acid, a stilbene carboxylic acid (Akiyama et al., 1999;Eckermann et al., 2003;Schröder and Schröder, 1990).Like for other phenylpropanoids, the precursors phenylalanine and tyrosine are derived from the shikimate pathway (D'Auria and Gershenzon, 2005;Herrmann and Weaver, 1999).These aromatic amino acids are converted via phenylpropanoid pathway into p-coumarate-CoA (Hahlbrock and Scheel, 1989).An elongation of this intermediate with malonyl-CoA occurs 3× in addition to ring closure and is catalysed by a polyketide synthase (PKS) (Ibrahim and Towers, 1962).A type III PKS subgroup, the CHS family, includes chalcone synthases and stilbene synthases.Together they produce a large group of plant phenylpropanoids including many phytoalexins from specialised metabolism.Phytoalexins play important roles in plant defence, stress responses, or detoxification (Francini et al., 2019;Gundlach et al., 1992;Herrmann and Weaver, 1999).Their biosynthesis can be induced externally by abiotic stressors such as UV radiation, high light, drought, salinity, and toxic compounds such as metals and biotic stressors including certain signal molecules (elicitors) released upon attack by pathogenic organisms or herbivores (Adrian et al., 1996;Baier et al., 1999;Matern, 1991).The phytoalexin spectrum is huge.Terpenoids that are composed of isoprene subunits are non-phenylpropanoid derived phytoalexins (Chen et al., 2019;Ebel, 1986;Kućand Rush, 1985;Li et al., 2015).
Hydrangea species are commonly known for their robustness against aluminium and also for their decorative flowers.H. serrata and H. macrophylla display a variation of sepal colour in response to soil pH change.The flowers are pink in alkaline and turn blue in acidic soil.Aluminium ions are complexed with two phenylpropanoids.This hydrangea blue complex is deposited in the sepals as means of detoxification (Ito et al., 2019;Schreiber et al., 2011).AlCl 3 is toxic to plants according to other studies (Rahman and Upadhyaya, 2021).It changes gene expression upon stress response and especially effects the roots (Zhou et al., 2016).The response to treatment with AlCl 3 on global changes in transcript levels was analysed in H. macrophylla which produces no DHICs (Chen et al., 2015).Other phenylpropanoids are inducible by AlCl 3 or CuCl 2 treatment (Adrian et al., 1996;Hanawa et al., 1992;Hossain et al., 2023;Schmidlin et al., 2008;Wang et al., 2022).The focus of our study was to observe whether phytoalexins or phenylpropanoids like DHICs in H. serrata can be induced by AlCl 3 treatment.Phytoalexins in hydrangeas were not described.Some monoterpenoids, the iridoids-glycosides from hydrangeas were reviewed by Gousiadou et al. (2016).Secoiridoid glycosides and hydrangeosides A and C were isolated from H. serrata leaves (Inouye et al., 1980;Kikuchi et al., 2008).
In this study, we aim to describe changes in the H. serrata transcriptome under abiotic stress and compare the effects on roots and leaves.Aluminium chloride-induced alterations in leaf and root transcripts were analyzed to provide insights into specialized metabolism-related gene activity.Utilizing the Illumina NextSeq2000 technique, we employed a deep sequencing approach, which offers a solid foundation for further data analysis and comparisons using the established dataset.The H. serrata transcriptome data from this study includes differentially expressed sequences of terpenoid synthases and genes involved in monolignol biosynthesis.This dataset can be used to search for novel genes and gene functions.Additionally, our findings suggest that H. serrata is unlikely to produce phenylpropanoids in response to aluminium stress.

Plant material
H. serrata (H.macrophylla sups.serrata "Oamacha") plants of about 40 cm height (GartenBaumschule TIMM, Reppenstedt, Germany) were cultivated in the greenhouse at the Bielefeld University.The identity of the plants was confirmed by the production of characteristic secondary metabolites (see results) and by the morphology.According to the leaf and flower morphology the analysed species is H. serrata.We follow the taxonomical description of the Royal Botanical Garden Kew (https://powo.science.kew.org).Ten individual plants were purchased as one cultivar and divided into two groups.Five plants in the control group were treated daily with 100 mL of tap water for 7 days.The treatment group consisted of the other five plants, which were daily watered with 100 mL of 100 mM AlCl 3 (Carl Roth, Karlsruhe, Germany) solution for 7 days.The plants were kept in 12 l containers with common potting soil.From every group, a pool sample of vital leaves was collected with N = 8 per plant.In parallel, a pool sample of about 2 g vital and rinsed roots were collected per individual plant of every group.All pool samples were immediately frozen in liquid nitrogen and stored at −80°C.

Total RNA extraction
The total RNA was extracted from 95 to 99 mg of frozen sample using the Spectrum Plant Total RNA Kit (Sigma-Aldrich, Taufkirchen, Germany) according to the supplier's protocol.Sixteen samples were prepared, consisting of four biological replicates for each condition from leaves and roots.On-Column DNase I Digestion Set from Sigma-Aldrich (Taufkirchen, Germany) was used according to the supplier's instructions.For RNA elution nuclease-free water was used (Merck, Darmstadt, Germany).The RNA quality was validated using a 2,100 Bioanalyzer system (Agilent, Santa Clara, CA, USA).Samples with RNA integrity number above 8 went into library construction.

Library construction and cDNA sequencing
Sequencing libraries were constructed from 3 µg of total RNA, following the TrueSeq v2 protocol (Illumina, Cambridge, UK).Paired-end sequencing was performed on an Illumina NextSeq2000 instrument at the Sequencing Core Facility of the Center for Biotechnology (CeBiTec) at Bielefeld University.A total of 24 samples libraries were applied on two lanes of a P3 flow cell.The paired-end read length was 150 bp and the total amount was 1.3 Mrd reads per sample.

Data processing and de-novo transcriptome assembly using Trinity
Illumina paired-end reads were processed analogous to the best practises guideline for a de-novo assembly using Trinity (Freedman and Weeks 2020;Grabherr et al., 2011).All computational operations were calculated using the de.NBI cloud (https:// www.denbi.de/cloud).First, read quality was examined by FastQC v0.11.3 (Andrews, 2010).Removal of erroneous k-mers was made by applying rCorrector v1.0.4 (Song and Florea, 2015).Adapter and low-quality regions were trimmed with TrimGalore! v0.6.6 (Krueger, 2015).Trimmed reads were obtained and a total number of 424,302,023 bases was assembled using Trinity v2.12.0 in genome independent-and paired-end mode, including normalisation to targeted maximum achievable coverage.

Functional annotation and differential expression analysis
The transcriptome assembly quality was rated via (BUSCO v3.0.2) (Simão et al., 2015) using the reference dataset of embryophyte_odb9 containing 1,440 Benchmarking Universal Single-Copy Orthologs (BUSCOs) of 30 species.TransDecoder (Haas et al., 2013) was used to identify open reading frames (ORFs) and translate ORFs to peptide sequences with a minimum length of 100 amino acids.Genes were annotated by transferring the A. thaliana Araport11 (Cheng et al., 2016) functional annotation to the H. serrata gene models as described before (Schilbert et al., 2021).The quantitative transcriptome data were obtained from read counts per contiq.Therefore, Kallisto v.0.44 (Bray et al., 2016) was applied using default parameters to quantify transcript abundances and to generate count tables.The generated count tables were subjected to differential gene expression analysis via DESeq2 to obtain the log 2 fold change for each gene (Love et al., 2014).Differentially expressed genes (DEGs) with a log 2 fold change >2 and adjusted p-value <0.1 between leaves versus leaves_ AlCl 3 and roots versus roots_ AlCl 3 were identified.Identified DEGs were visualized by MapMan as follows: The functional annotations were sorted by biological categories of the proteins using Mercator4 v2.0 (Schwacke et al., 2019) and mapped together with the quantitative transcriptome data onto specialised metabolite pathways in MapMan 3.5.1R2(Usadel et al., 2009).Pathway map files were adjusted to the biological categories using GIMP 2.10.18.The DEGs were displayed colour scaled on pathway maps and summarized in tables (Supplementary Material).The biological context was obtained from protein category bins using Mercator4 and MapMan.DEGs from the selected biological categories of interest were visualised in Venn diagrams using http:// bioinformatics.psb.ugent.be/software.

Results
H. serrata plants of about 40 cm height were grown in potting soil under optimal greenhouse conditions.The plants were treated with AlCl 3 solution or with tap water as control.Leaves and roots were harvested, RNA was extracted and cDNA libraries were sequenced using the Illumina technology.The number of resulting 580,770 contigs includes isoforms, allelic variants and other sequences.The application of BUSCO revealed 96.2% complete BUSCO genes in the de-novo H. serrata transcriptome assembly.Among the complete BUSCOs, 10.4% were single-copy genes whereas 85.8% were duplicated.Moreover, 1.6% of all BUSCO genes were fragmented and only 2.2% were missing.
A differential expression analysis between AlCl 3 -treated organs and corresponding controls was made using DESeq.As displayed in volcano plots, a higher number of genes were differentially expressed in AlCl 3 -treated leaves compared to AlCl 3 -treated roots (Figure 1).
A total of 119,756 representative peptide sequences were annotated via reciprocal best BLAST hits (RBHs) and best BLAST hits (BBHs) using Arabidopsis thaliana data.A total of 75,100 annotated sequences were sorted in biological categories plus 42,432 enzyme classification categorized in more than 50 bins using Mercator 4.0.Each bin contained second and third order subclasses or more for detailed specification.

Biological categorisation of differential expressed genes in response to aluminium chloride treatment
Categories of functional annotations were summarized using Mercator4.Tables 1, 2 show the strongest regulated genes, filtered for a log 2 change <4>, and their annotations.Both H. serrata tissues showed transcriptional changes in the biological categories of solute transport, phytohormone action, external stimuli response, and specialised metabolism (Figures 2, 3).The most annotations were categorized into enzyme classification and solute transport followed by protein modification and vesicle trafficking (Figure 2).Many genes encoding for glycosyltransferases and oxidoreductases were regulated in both tissues.A higher number of glycosylases encoding genes were regulated in the roots (Figure 2; Tables 1, 2).
Among DEGs from AlCl 3 -treated roots (roots_AlCl 3 ) versus control with a log 2 fold change >4, several genes encode for enzymes such as oxidoreductases (six cases), glycosyltransferases (two cases), glycosylases (six cases), and one methylsterol monooxygenase from phytosterol metabolism.Additionally, differentially expressed sequences specialised metabolism include genes encoding terpene synthase and triterpenoid synthase.Other genes were annotated as transporter as well as ion channels.Downregulated genes with log 2 fold change <−4 were annotated as glycosylases (five cases), one glycosyltransferase, hexosyltransferases (two cases), and acyltransferases (two cases).Furthermore, downregulated genes encoded for iron cation sensor protein as well as many transporters.A polymeric acetyl-CoA carboxyltransferase subunit from the fatty acid metabolism encoding transcript was significantly downregulated (Table 1).
Genes for aluminium or general heavy metal ion response mechanisms like transporters were similarly differentially expressed, like Chen et al., 2015 reported for H. macrophylla.One interesting exception is the aluminium-activated malate transporter (QUAC/ALMT) which was slightly downregulated in roots (Supplementary Material).Also, in AlCl-treated roots of H. serrata the iron cation sensor protein HRZ/BRUTUS, ferric-citric complex transporter, and metal-citrate complex transporter FRD were downregulated.
DEGs in AlCl 3 -treated leaves (leaves_AlCl 3 ) versus control with log 2 fold change >4 was functionally annotated as transporters, transmembrane proteins (two cases), one acetyl-CoA binding protein, one salicylic acid binding protein, one glycosylase, oxidoreductases (three cases), one hexosyltransferase and one glycosyltransferase (Table 2).Downregulated genes with log 2 fold change <−4 encoded for other transporters, calcium cation channel, oxidoreductases (four cases), one glycosylase, and one ROS-induced xanthine dehydrogenase (Table 2).The biological categorisation of functional annotated DEGs were mapped onto pathways of the specialised metabolism.These MapMan images show which gene sets were regulated under the presence of AlCl 3 .The most significantly regulated gene sets covered the shikimate and MVA pathway, terpenoid, lignan, and lignin biosyntheses (Figures 4, 5).More of these genes were upregulated in roots_AlCl 3 than in leaves_AlCl 3 (Figure 4).One PAL isoform was downregulated in the roots (−2).Additionally, in the leaves an isoprenyl diphosphate synthase (−3.2) was downregulated, a triterpenoid synthase (6.1), a terpene synthase (4.1), and an MEP pathway reductase (2.3) were upregulated.Obviously, AlCl 3 induced stress activated genes of the terpenoid pathway and genes encoding for enzymes of the lignin biosynthesis.

LC-MS analysis of organ specific dihydroisocoumarin levels
The Illumina NextSeq2000 cDNA sequencing was followed by a general transcriptome analysis.In DHIC biosynthesis or among PKS homologues, no differential gene expression was observed.It was the aim of this study to decipher differences on DHIC level.No significant alteration as a response to aluminium chloride treatment, in terms of DHICs was detected, in the LC-MS analysis of the methanolic extracts.The metabolite analysis shows an organ specific accumulation profile of DHICs and mono-as well as diglycosidic versions (Table 3).Also, a small elicitor study using yeast extract solution on H. serrata leaves, as well as a drought experiment and jasmonic acid treatment neither enhanced HYD nor PHY in the LC-MS measurement of leaf or root extracts (data not shown).Interestingly, the plants showed signs of stress response in specialised metabolites.One sign was the accumulation of anthocyanins, which was visible to the eye in the experiment that combined high light and drought.Another observation was the formation of necrosis in leaves treated with yeast-elicitor solution (data not shown).The DHIC content was not affected in the described experiments (data not shown).

Discussion
Illumina NextSeq 2000 sequencing of cDNAs is a powerful technique to gain insight into differences in gene expression related to plant specialised metabolism under abiotic stress.This study is based on paired-end reads which were sequenced in two technical replicates from sixteen cDNA libraries, consisting of four biological replicates for each condition applied on two H. serrata organs.A total of 380,000 different transcripts was proposed as genes by Trinity out of 580,770 assembled contigs from H. serrata leaf and root sequencing data.For roots from H. macrophylla 256,127 unigenes of 401,215 contigs were reported (Chen et al., 2015).This transcriptome study examined two tissues under aluminium influence and we saw that a higher number of genes were expressed in leaves than in roots.AlCl 3 , in some studies, has been reported to be toxic to plants (Rahman and Upadhyaya, 2021).It changed gene expression as a stress response and especially effected the roots (Zhou et al., 2016).A response of H. serrata under AlCl 3 treatment detected as DEGs was found in roots and additionally in leaves.Similar genes encoding transporter proteins which were found in H. serrata were described for H. macrophylla roots and studies on other plants (Chen et al., 2015;Wu et al., 2023;Zhang et al., 2019;Zhou et al., 2016).
The transcriptional changes under aluminium chloride treatment were visualized using the MapMan software.Combined with the high-quality functional annotation via data from A. thaliana and a Mercator4 summarization, the DEGs were categorized according to the biological context of their functional annotation.Signs of DHIC biosynthesis were not visible in differential gene expression against our expectation.In contrast, genes encoding for enzymes of the terpenoid biosynthesis and cell wall lignification were upregulated.
From the context of specialised metabolism many genes were downregulated in roots and only a few were upregulated related to isoflavonol, anthocyanin, carotene, and betaine synthesis.Additionally, lignin-and lignan biosynthesis related genes were upregulated in both tissues.Enhanced lignin synthesis is known as a part of stress defence in plants (Moura et al., 2010;Sasaki et al., 1996;Vanholme et al., 2012).Metal ions cause oxidative stress as ROS in plants.The data show evidence that ROS detoxifying mechanisms were activated or enhanced.Genes from the sectors of glutathione synthesis, sulfur containing amino acids, and redox homeostasis were upregulated in roots, but not in leaves.Tissue-specific alterations in roots_ AlCl 3 versus control roots showed that genes encoding for genes related to betain, tocopherol, and phenylpropanoid metabolisms were downregulated.Genes related to MVA, MEP, non-MVA, carotenoid, and shikimate pathway were upregulated.Phenylpropanoid synthesis-involved genes were less expressed than expected.It is remarkable that no homologue of Type III PKS was among DEGs under AlCl 3 treatment.Moreover, genes from terpenoid and lignin metabolism were upregulated.Upregulated DEGs from leaves_ AlCl 3 versus control leaves belong to MVA, shikimate, tocopherol, and minor to phenylpropanoid pathways.In addition to various phenylpropanoids, the terpenoids are known to have phytoalexin functions (Chen et al., 2019;Li et al., 2015).
The external stimulation for H. serrata DHICs seems difficult because little is known about the biological purpose in the plant.Studies on highlight or shading experiments did not effectuate a change in the DHIC content in leaves (Moll et al., 2021).Also, herbivore attacks did not enhance the DHIC content (Ujihara et al., 1995).Aluminium chloride is another stressor that did not affect the DHIC biosynthesis.A recent study shows that neither ethylene nor salicylic acid had an influence on the HYD and PHY accumulation in H. serrata leaves.The HYD and PHY contents in leaves were enhanced due to a methyl jasmonic acid treatment (Preusche et al., 2022).So, methyl jasmonic acid seems to be the right elicitor for the induction of the DHIC biosynthesis and it needs further investigation in the discovery of genes involved in the DHIC biosynthesis and the biological assignment of DHICs.
Important bioactive terpenoids were isolated from hydrangeas, but they were not investigated as phytoalexins.The isolated terpenoids from H. macrophylla were identified as secoiridoid Counts of differentially expressed genes per biological category in leaves treated with aluminium chloride (AlCl 3 ) compared to control leaves and in roots treated with AlCl 3 compared to control roots.Functional annotated proteins were sorted in biologicals categories using Mercator4.The assigned functional protein categories were selected for specialised metabolism and the corresponding amount of DEGs having a log 2 fold change >2; <−2 and adjusted p-value <0.1 was counted.DHIC aglycons (H-aglyc and P-aglyc), mono-(H-glyc, P-glyc) and di-glycosides (H-2glyc, P-2glyc) occur in specific contents which differ between leaves and roots.Hydrangenol (H) and one hydrangenol isobar, putatively hydrangeic acid (H Acid), were more abundant in roots.The content of phyllodulcin (P-aglyc and P-glyc) was high in leaves but very low in roots.Both diglycosides (H-2glyc, P-2glyc) were detected in leaves as well as in roots.Four biological replicates were analysed for each condition and organ.The plants were treated with tap water in the control group and with 100 mM AlCl 3 -solution in the treatment group for 5 days.UHPLC-ESI-QToF-MS measurement was performed after methanolic DHIC extraction.
glycosides, loganin, secologanin, secologanic acid, and sweroside (Gousiadou et al., 2016;Yoshikawa et al., 1994).Secoiridoid glycosides and hydrangeosides A and C were isolated from H. serrata leaves (Inouye et al., 1980;Kikuchi et al., 2008).Other terpenoids then described in literature may exist as inducible phytoalexins.The use of terpenoids as phytoalexins in response to aluminium chloride in H. serrata would explain that less transcriptional changes related to phenylpropanoid metabolism were detected in this study.Even this extensive transcriptome dataset for this study on H. serrata was generated, its analysis in terms of DHIC was limited.The most limiting factor is the rareness of 3,4-DHICs in Hydrangea species and in plant kingdom at all and therefore a lack of suitable reference enzymes in forms of annotated PKS.

Conclusion
Many aspects about the main mechanisms of aluminium chloride detoxification in hydrangea which were described in other studies were observed in this study as well.It is important to mention that transporter and cell wall stability related genes were differentially expressed as a similar study on H. macrophylla showed.This work was more focussed on the specialised metabolites and especially in DHICs.Interestingly, H. serrata in contrast to H. macrophylla showed changes in the terpenoid biosynthesis as response to aluminium chloride.In terms of DHIC biosynthesis little references on genetic level are available to fully explore the differential gene expression.

FIGURE 1
FIGURE 1Volcano plots showing differentially expressed genes (DEGs) in leaves treated with aluminium chloride (AlCl 3 ) compared to control leaves as well as in AlCl 3 -treated roots compared to control roots.DEGs with a log 2 fold chance >2 and <−2 and an adjusted p-value of <0.05 are marked in orange.FDR is the false discovery rate of the adjusted p-value.The individual H. serrata plants were treated with tap water in the control group and with 100 mM AlCl 3 -solution in the treatment group for 5 days.Each condition group consists of four biological replicates.

FIGURE 2
FIGURE 2Differentially expressed genes (DEGs) are displayed for selected biologicals categories based on Mercator4 protein assignment to functional context.External stimuli response and solute transport in roots treated with aluminium chloride (AlCl 3 ) compared to control roots as well as in leaves treated with AlCl 3 compared to control leaves and specialised metabolism in roots treated with AlCl 3 compared to control roots and in leaves treated with AlCl 3 compared to control leaves.Both Venn diagrams show tissue specific DEG amounts, log 2 fold change >2; <−2 and an adjusted p-value <0.1 from aluminium chloride (AlCl 3 ) treated roots and leaves compared to the respective control.DEGs with same expression level are in the overlap.

FIGURE 4
FIGURE 4Pathway map for H. serrata specialised metabolism and differential expressed genes (DEGs) from roots treated with aluminium chloride (AlCl 3 ) compared to control roots and from leaves treated with AlCl 3 compared to control leaves.DEGs are shown as coloured squares in a MapMan-based image.Each square stands for one DEG or isoform.The genes were mapped on the pathway image based on biological categories of the encoded proteins which were summarized using Mercator4.The false colour scale represents log 2 fold change >0 in red and <−0 in blue for adjusted p-value <1.Each condition group consists of four biological replicates from H. serrata.The plants were treated with tap water in the control group and with 100 mM AlCl 3 -solution in the treatment group for 5 days.

FIGURE 5
FIGURE 5Pathway map of H. serrata specialised metabolism and differential expressed genes (DEGs) from and roots compared to leaves aluminium chloride (AlCl 3 ) treatment group.Coloured squares are shown in MapMan-based image.Each square stands for one differentially regulated gene or isoform.The genes were mapped on the pathway image based on biological categories of the encoded proteins which were summarized using Mercator4.The false colour scale represents log 2 fold change >2 in red and <−2 in blue for p-value adjusted <0.05.

FIGURE 6
FIGURE 6Overview of H. serrata metabolism and its involved differentially expressed genes (DEGs) from roots compared to leaves.DEGs are shown as coloured squares in a MapMan-based image.Each square stands for one DEG or isoform.The genes were mapped on the pathway image based on biological categories of the encoded proteins which were summarized using Mercator4.The.The false colour scale represents log 2 fold change >2 in red and <−2 in blue for p-value adjusted <0.05.

TABLE 1
Functional annotations of differential expressed genes from AlCl 3 -treated H. serrata roots versus control with a log 2 fold change.

TABLE 2
Functional annotations of differential expressed genes from AlCl 3 -treated H. serrata leaves versus control with a log 2 fold change.

TABLE 3
Relative content of dihydroisocoumarins (DHICs) in extracts from aluminium chloride (AlCl 3 ) treated and control Hydrangea serrata leaves and roots.
*Not confirmed by standard.