Transcriptomic and metabolomic characterization of antibacterial activity of Melastoma dodecandrum

Antibacterial resistance poses a significant global threat, necessitating the discovery of new therapeutic agents. Plants are a valuable source of secondary metabolites with demonstrated anticancer and antibacterial properties. In this study, we reveal that Melastoma dodecandrum exhibits both bacteriostatic and bactericidal effects against Pseudomonas aeruginosa and Staphylococcus aureus. Treatment with plant extracts results in membrane damage and a reduction in P.aeruginosa swimming and swarming motility. A comparative analysis of bacterial transcriptomes exposed to M.dodecandrum extracts and four distinct antibiotics indicates that the extracts may trigger similar transcriptomic responses as triclosan, a fatty acid synthesis inhibitor. Activity-guided fractionation suggests that the antibacterial activity is not attributable to hydrolyzable tannins, but to unidentified minor compounds. Additionally, we identified 104 specialized metabolic pathways and demonstrated a high level of transcriptional coordination between these biosynthetic pathways and phytohormones, highlighting potential regulatory mechanisms of antibacterial metabolites in M.dodecandrum.


Introduction
Antimicrobial resistance (AMR) has been identified by the World Health Organization (WHO) as one of the top ten global health threats affecting humanity (Antimicrobial resistance).This problem is exacerbated by the lack of new effective antibacterials within clinical pipelines (Antimicrobial resistance).Breakthroughs in genomics and bioinformatics in the 90s enabled the identification of essential bacterial genes and targets.Consequently, pharmaceutical companies have attempted to generate new leads via target-based rational drug design.However, this strategy had limited success, partly due to an incomplete understanding of resistance mechanisms and factors influencing drug membrane permeability (Payne et al., 2007;Tommasi et al., 2015;Lewis, 2020).Alternatively, bioprospecting, the examination of natural sources for active bioproducts, may be considered in the search for new antibacterial compounds.To this end, the plant kingdom represents a rich source of plant secondary metabolites, many of which are known to have antimicrobial properties and have been used in traditional medicine to treat bacterial infections (Rıós and Recio, 2005;Frey and Meyers, 2010;Kessler and Kalske, 2018;Singh et al., 2018;Lal et al., 2020).
While bioprospecting may enable the discovery of more targetspecific and structurally complex active biological compounds, it comes with associated challenges (Cushnie et al., 2020).For example, there may be difficulties in isolating and identifying active compounds from biological extracts which may be further compounded if there are multiple active compounds that contribute synergistically to the therapeutic activity of the extract.Similarly, the subsequent elucidation of the mechanism of action of an active compound may also be a tedious process (Atanasov et al., 2021).Additionally, the further use and industrial production of plant metabolites may be hindered by a limited understanding of their biosynthetic pathways and that the active metabolites may only be produced under specific conditions in specific organs (Kumar et al., 2019;Atanasov et al., 2021).Multi-omics methods, consisting of a combination of transcriptomics, metabolomics, genomics and/or proteomics, have been utilized in the drug discovery process and have been recognized as a useful tool to identify new targets and the mechanisms of lead compounds (Goff et al., 2020).The multi-omics strategy may also be used to discover plant metabolic genes and pathways using co-expression analysis and metabolite-based genome-wide association studies (Zhan et al., 2022;Zhao and Rhee, 2022;Julca et al;Usadel et al., 2009;Lim et al., 2022).Taken together, such multi-omics approaches may represent the means to address challenges associated with bioprospecting and bridge the knowledge gap.
Members of the Melastoma genus, such as Melastoma malabathricum, Melastoma candidum D. Don, and Melastoma dodecandrum, have long been used in areas such as Malaysia, Taiwan, and China as traditional medicine for the treatment of dysentery, wounds, high blood pressure, diabetes, and skin diseases, amongst others (Wang et al., 2008;Joffry et al., 2012).More recently, studies have been carried out to scientifically evaluate M. malabathricum and M. candidum D. Don for acute toxicity, and antibacterial, antioxidant, and immunomodulatory activities (Alnajar et al., 2012;Wong et al., 2012;Che Omar et al., 2013;Mamat et al., 2013;Zheng et al., 2021).Similarly, efforts have been made to identify the chemical constituents of M. dodecandrum and its associated activities (Yang et al., 2014;Wang et al., 2017;Tong et al., 2019;Huang et al., 2021;Xu et al., 2023).However, limited studies have been carried out on evaluating the antimicrobial activity of the plant and identifying associated active compounds or mechanisms, with only a recent paper reporting moderate antibacterial activity of its ethanolic extract against diarrheagenic bacterial pathogens following an activity screen of 32 plants (Kudera et al., 2021).
Through the screenings conducted in this study, we have identified M. dodecandrum Lour.as a plant with antimicrobial and bactericidal effects on Pseudomonas aeruginosa.Phenotypic studies on P. aeruginosa indicate that the plant's active compounds likely result in membrane damage and affect bacteria motility.Additionally, a comparative transcriptomic analysis of P. aeruginosa suggests that out of four antibiotics that target various vulnerabilities of bacteria, the plant extract's mode of action overlaps with triclosan.Subsequent evaluation of fractions of the plant extract and activity testing of purified compounds suggests that an unknown compound is responsible for the antibacterial activity.By collecting samples of different organs, and carrying out transcriptomics studies and coexpression analysis, we investigated the biosynthetic potential of M.dodecandrum and showed that the plant likely synthesized many types of phenylpropanoids, with high levels of expression of sulfur-containing compounds, phenylpropanoids, and terpenoids.Finally, we showed co-expression patterns between the specialized metabolic pathways and hormones, suggesting that the pathways in the plant are transcriptionally linked.

Plant sampling
The flowers, fruits, leaves, stems, and roots of M. dodecandrum were sampled from the NTU Community Herb Garden.Leaves used solely for bioactivity screening were transported in dry ice and lyophilized using a freeze dryer (Labconco).Plant organs used for both bioactivity and transcriptomic studies were transported in liquid nitrogen.

Metabolite extraction
Plant material was ground in liquid nitrogen to a fine powder.One hundred mg of the powder was aliquoted into a screw cap tube.Three 3 mm steel beads and 500 μL of 80% HPLC-MS grade methanol were then added.The sample was subsequently homogenized for two cycles at 2000 rpm (PowerLyzer 24, Qiagen).The steel beads were removed, and the samples were placed on a shaker for 10 mins at 1500 rpm at room temperature in the dark.The samples were centrifuged at 10 min, 12000 × g.The supernatant was transferred into a new 2 mL Eppendorf tube and kept on ice.Five hundred μL of 80% HPLC-MS grade methanol was added to the remaining pellet, and the shaking and centrifuge step was repeated twice until 1.5 mL of supernatant was collected.The supernatant was dried completely using a centrifugal vacuum concentrator (Centrivap, Labconco) and dissolved in 100 -500 μL of DMSO before use.

Antimicrobial screening and activity testing
Dried methanol extracts derived from 100 mg of plant material or different plant organs were dissolved in 500 μL DMSO to give a stock solution of 200 mg/mL plant material.Overnight bacterial cultures were subcultured 1:10 in 10 mL of LB10 broth until the bacteria reached the exponential phase.Exponential phase bacteria were then diluted in Mueller-Hinton broth (MHB) to a final concentration of 5 × 10 5 CFU/mL, and 2 μL of the stock DMSO solution was added to 100 μL of MHB to a final concentration of 4 mg/mL starting plant material or 2% v/v DMSO.Evaluation of HPLC fractions was carried out at 2% v/v DMSO, with each fraction corresponding to ~10 mg/mL mass equivalent of the crude methanol extract.

Extraction, bioassay guide fractionation, and chemical characterization
Freeze-dried leaves of M. dodecandrum (11 g) were milled and extracted twice with 400 mL of methanol (MeOH).The combined MeOH extract was dried under reduced pressure to afford 3 g of dried MeOH extract.Crude MeOH extract (2 g) dissolved in 5 mL of MeOH was fractionated using a Sephadex LH20 open column (5 x 35 cm).Compounds were eluted sequentially with 1 L of MeOH, 50% aqueous acetone, and acetone to generate 13 fractions.Biological testing results showed activity in only one fraction (58K, 120 mg) which was eluted in the acetone eluent.50 mg of active fraction 58K was reconstituted in 15% aqueous MeOH, sonicated, and centrifuged.The insoluble pellet was collected (labeled as PPT, 5 mg) while the supernatant was injected into a Waters XTerra Prep MS C18 Column (10 x 300 mm) on an Agilent 1260 Infinity II Preparative HPLC system.The mobile phase consisted of water (A) and acetonitrile (B), both with 0.1% formic acid.The gradient elution started at 15% B for 5 mins, 15-30% B in 20 mins, 30-60% B in 25 mins, held at 60% for 10 mins, then increased from 60-100% B in 2 mins and washed at 100% B for 10 mins.The flow rate was 20 mL/min and the detection UV wavelength was at 254 nm. 8 fractions were generated from the preparative HPLC.All HPLC fractions, PPT and parent fraction 58K were submitted for biological evaluation.
Purified fractions were analyzed on an Agilent UPLC1290 coupled with a Quadrupole Time-of-Flight (Q-TOF) system.Separation was carried out with a reversed-phase C18 column (2.1 x 50 mm) at 0.5 mL/min, using a 10 mins linear gradient with 0.1% formic acid in both solvent A (water) and solvent B (acetonitrile).The typical QTOF operating parameters were as follows: negative ionization mode; sheath gas nitrogen flow, 12 L/ min at 275°C; drying gas nitrogen flow, 8 L/min at 250°C; nebulizer pressure, 30 psi; nozzle voltage, 1.5 kV; capillary voltage, 1.5 kV.Lock masses in negative ion mode: TFA anion at m/z 112.9856 and HP-0921 TFA adduct at m/z 1033.9881.1H NMR spectra were acquired on a Bruker DRX-400 NMR spectrometer with a 5-mm BBI Cryoprobe.

Determination of minimum inhibitory concentrations
Minimum inhibitory concentrations (MIC) were determined using broth microdilution methods as described in (Wiegand et al., 2008).Briefly, plant extracts or equivalent concentrations of DMSO were diluted in MHB to four times the desired final concentration.50 μL of the samples were added to 50 μL of MHB media on the first column of a 96-well plate, then serially diluted two-fold down the row.Untreated media-only growth control wells were also included.Overnight cultures of PAO1 were subcultured in LB10 medium and further incubated at 37°C with 200 rpm shaking.Exponential phase cells were then diluted in MHB to 1 × 10 6 CFU/mL.Fifty μL of the bacterial culture was added to each well to a final concentration of 5 × 10 5 CFU/mL, with the final concentration of plant extract falling between 0 -50 mg/mL (0 -5% DMSO).Sterility control wells, consisting of only MHB media, and blank control wells, containing 0 -50 mg/mL of plant extract and MHB media, were also included.The 96-well plate was incubated for 24 h at 37°C with OD readings taken at 600 nm every 10 mins using a microplate reader (Tecan M200).Readings from blank control wells were subtracted from wells with bacteria added to account for high background readings at high plant extract concentrations.MIC was similarly determined for 0 -512 μg/mL of the antibiotics and antimicrobial compounds triclosan, rifampicin, ceftazidime, and gentamicin.MIC was defined as the first concentration at which no growth is observed after 24 h incubation.
Antibacterial activities across experiments were compared using the area under the curve (AUC), which represents bacterial growth over a period of 18 h.The AUC is calculated using Graphpad prism (version 9.3.0).The same parameters (t = 0 -1080 mins) were used for all data.Eighteen hours were used for the AUC calculations even though OD data were collected over 24 h, as cell growth has already plateaued and increased noise was observed at later time points, possibly due to evaporation of the small volume of medium used.
The AUC values were plotted and compared using one-way ANOVA with multiple comparisons in Graphpad prism.

Time-kill kinetics assay
Overnight cultures of PAO1 were subcultured in LB10 medium and further incubated at 37°C with 200 rpm shaking.Exponential phase cells were diluted to 1 × 10 6 CFU/mL in MHB media and treated with 0 -2 × MIC of plant extract.Vehicle controls were treated with equivalent volumes of 0 -5% v/v DMSO.At each timepoint, samples were collected for serial dilution in PBS and CFU enumeration.Five μL of each dilution was spot plated onto LB10 agar plates and incubated at room temperature overnight.The first spots containing between 3 -30 colonies were counted.The assay was carried out independently twice with two technical replicates each time.

Evaluation of membrane integrity
The effect of plant extract on membrane integrity was evaluated using the LIVE/DEAD stain (L7012, Thermo Fischer) and microtiter plate reader (Tecan) as described in the commercial protocol.Briefly, sub-cultured exponential phase bacterial cells were washed in 1 × PBS (PBS).5 × 10 7 CFU/mL of bacterial cells diluted in PBS were subsequently treated with a non-inhibitory concentration of 0 -4 mg/mL of filtered plant extract or equivolume of DMSO for 10 mins.Samples were then treated for a further 30 mins with LIVE/DEAD stain readings taken with excitation at 485 nm and emission at 535 (red) and 635 nm (green).Dead cells were prepared via incubation of bacteria at 65°C for 30 mins.The percentage of live cells is determined using a standard curve generated through a mixture of live and dead cells and the calculation of the green/red fluorescence ratio as stated in the commercial protocol.When LIVE/DEAD stain is used in microscopy, 2 × 10 7 cells were seeded into a well of a poly-L lysine pre-treated ibidi 8-well chamber slide (80826, ibidi).The cells were allowed to attach for 10 mins before each well was washed and treated for 10 mins with 0 -16 mg/mL of plant extract or equivolume of 0 -1.6% v/v DMSO.Following treatment, each well was washed again and stained with LIVE/DEAD stain as described in the commercial protocol.At least three images per sample were taken using an epifluorescence microscope (Axio Observer Z1, Carl Zeiss) using the FITC filter for imaging SYTO 9 fluorescence and AF 568 filter for imaging propidium iodide (PI) staining.Images of fluorescence staining were analyzed using Image J.

Motility assays
Swimming and swarming assays with P. aeruginosa were carried out similarly to (Ha et al., 2014a;Ha et al., 2014b).Briefly, 0.3% w/v or 0.5% w/v of 1X M8 agar (1X M8 solution: 42.3 mM Na 2 HPO 4, 22 mM KH 2 PO 4, 8.56 mM NaCl) were prepared for use in swimming and swarming assay respectively.Plant extracts or equivalent volumes of DMSO were added to the well with the highest concentration of the extracts, then serially diluted 4-fold down the wells before the agar solidified to a final concentration of 0 -4 mg/mL plant extract.The agar was dried under laminar flow for 45 mins and used immediately after drying.Overnight cultures of PAO1 were used for inoculation.The agar was inoculated either by dipping an Eppendorf pipette tip into the culture and then stabbing it midway of the depth of the agar (swimming assay) or by spotting 1 μL of the culture on the surface of the agar at the centre of the well (swarming assay).Following inoculation, the plates were incubated at 30°C overnight and the images were taken 16 h after inoculation.Swimming and swarming areas were assessed using Image J by setting an appropriate scale and selecting the region of interest.The resulting values are plotted in Graphpad (9.3.0) and analyzed using one-way ANOVA with multiple comparisons against the control group.

RNA extraction and sequencing
Following M. dodecandrum sample collection and grinding into a fine powder (as described in 'Metabolite extraction'), RNA was extracted from 100mg of a sample using Spectrum ™ Total Plant RNA Kit (Sigma) Protocol A following manufacturer's instructions.Three biological replicates were collected for each adult organ, while single replicates were collected for the young flowers.To extract RNA from bacteria (PAO1), overnight cultures of PAO1 were first subcultured and grown to mid-log phase.Cells were adjusted to OD=0.4 and then treated with equivolume of 2X antibiotic dissolved in MHB medium to a final concentration of OD=0.2, MIC=1X (or 25 mg/mL plant extract and 128 μg/mL triclosan), and a final volume of 500 μL.The samples were incubated at 37°C statically for 30 mins.Following treatment, 1 mL of RNA protect reagent (76506, Qiagen) was added to 500 μL of the sample.The sample was treated as described in the commercial kit and the dried pellet was stored at -80°C.RNA was extracted using the RNeasy plus mini kit (Qiagen) as per the kit's instructions.For both plant and bacteria, quality control of all extracted RNA was carried out by Novogene (Singapore).Each sample was evaluated for its quantity, integrity, and purity using agarose gel electrophoresis, Nanodrop, and Agilent 2100 Bioanalyzer.Library construction was performed by Novogene where mRNA of plants was enriched from total RNA with oligo(-dT) magnetic beads.No ribosomal RNA depletion was performed for bacteria library construction due to high sequencing depth.The library was quantified with Qubit and real-time PCR, and then sequenced using Illumina NovaSeq 6000, with paired-end sequencing of 150 base pairs (bp) per read and a sequencing depth of approximately 60-70 million reads.

RNA sequencing analysis (P. aeruginosa)
Low-quality RNA-seq reads were removed, and the remaining reads were trimmed with fastp (v0.23.2) (Chen et al., 2018).Reads were then aligned to the PAO1 genome (NC_002516.2) and read counts were obtained using EDGE-pro (Magoc et al., 2013).Due to poor overall read alignment, one DMSO-treated sample (replicate 1) was discarded (data not shown).One-third of the reads from the remaining two DMSO-treated samples were removed and collated to form a new third replicate.To evaluate sample-level variability, read counts were transformed using variance-stabilizing transformation, followed by hierarchical clustering of samples using Euclidean distance as the distance metric (Love et al., 2014).Differentially expressed genes (DEGs) between each antibiotic-treated sample compared to the DMSO control were identified with read counts using DESeq2.Genes with jlog 2 fold change| ≥ 1 and Benjamini-Hochberg (BH) adjusted p-value< 0.01 were considered significant (Benjamini and Hochberg, 1995) for multiple-testing correction.Gene ontology (GO) enrichment analysis was performed via the GO webtool on each list of significant DEGs against the PAO1 GO biological process data set by applying Fisher's Exact test with BH correction (adjusted pv a l u e < 0 .0 1 ) ( M i e t a l ., 2 0 1 9 ; The Gene Ontology Consortium, 2021).

M.dodecandrum gene expression data preprocessing
Quality control and read trimming of raw RNA-seq data were performed with fastp (v0.23.2).To obtain normalized transcript abundances (transcripts per million; TPM) of the RNA-seq samples, the reads were pseudoaligned against the coding sequences of the M.dodecandrum genome (PRJCA005299 from the National Genomics Data Center) (Hao et al., 2022) and quantified using Kallisto (v0.46.1) (Bray et al., 2016).To evaluate sample-level variability, samples were clustered using Spearman's correlation as the distance metricd x,y = 1 − jr x,y j, where x and y are vectors comprising TPM expression profiles for all transcripts of a pair of RNA-seq samples respectively and r is the Spearman correlation coefficient.The final gene expression matrix was derived after filtering out genes with an average TPM of 0 across all organs.Genes were translated into protein sequences using the transeq module (EMBOSS v6.6.0.0) (Rice et al., 2000).
The following analyses were performed in M.dodecandrum

Enzyme annotation and pathway prediction
Enzyme annotation and pathway prediction were performed following the pipeline described by (Hawkins et al., 2021).Briefly, protein sequences were annotated with E2P2 (v4), pathway membership was predicted using PathoLogic in PathwayTools (v26) with default settings, and semi-automated validation of the final pathways was performed by SAVI (v3.1).All pathways assigned to the manual-to-validate list by SAVI were retained except for PWY0-501 and PWY-5723, following recommended validation steps by SAVI.From SAVI's list of pathways to add, PWY-5173 was discarded as it is a deprecated ID in the current version of MetaCyc used by PathwayTools.

Gene co-expression and pathway subcluster identification
To determine if genes within an SM pathway are more coexpressed than by chance, Pearson's correlation coefficient (PCC) for each gene-pair within a pathway was first calculated.To simulate random pathways, a list of pathway-labeled gene-pairs was generated.The pathway labels were then randomly shuffled amongst the gene-pairs.The PCC for each gene-pair was calculated, and the median PCC of each random pathway was extracted.This was repeated 100 times.
To identify subclusters of a pathway, genes within a pathway were hierarchically clustered based on their pattern of expression across samples using the dendrogram function from scipy with default parameters (where the distance threshold is 70% of the maximum dendrogram distance).A representative subcluster was only identified for a pathway if the subcluster consisted of genes that represented at least 75% of the annotated reactions in the pathway in the plant (degree of completeness).As some pathways have multiple genes catalyzing the same reaction, subclusters of differing sizes can have the same degree of completeness.In these cases, the largest subclusters will be chosen as representative subclusters.For pathways with no representative sub-cluster, the original pathway was the representative 'sub-cluster' in downstream analyses.

Co-expression analysis of connected pathways
To investigate patterns of inter-pathway coexpression within and between SM classes, we generated a pathway coexpression network.First, gene-pairs of all SM genes were calculated using the 'pearsonr' function from the scipy package.Then, p-values were corrected with the BH method (Benjamini and Hochberg, 1995).Gene-pairs with PCC ≥ 0.6 and BH-adjusted p-value (two-tailed)< 0.05 were considered significant.Subsequently, the gene nodes within the gene co-expression network were labeled with their respective SM pathways (as determined by pathway prediction).Multiple pathway labels will be assigned to genes that are predicted to be involved in multiple SM pathways.Labels of "Superpathways", defined as pathways that combine individual pathways, will only be assigned to a gene if the gene is not already labeled with their respective subpathways.
Next, to determine if the intra-pathway coexpression was better than random, the network was then subjected to the following permutation analysis: for each pathway, the number of intrapathway edges was counted.Pathway labels were then randomly shuffled and the number of intra-pathway edges was counted.The following empirical p-value was calculated for each pathway: where N is the number of permutations, I is an indicator function that takes a value of 0 or 1 if the statement is false or true, respectively.This was repeated 10 000 times.Pathways with BH-adjusted p-values of< 0.05 were considered to be coexpressed (gene members of the particular pathway are coexpressed with each other).Inter-pathway coexpression between pathway-pairs were also determined similarly by counting inter-pathway edges that connect genes across any given two pathways.
To generate a pathway coexpression network, edges were used to connect SM pathways (represented as nodes) that were coexpressed with each other.Nodes were named according to the SM classes each pathway belongs to.Each node has a suffix of "_X" where X represents a single pathway from an SM class.e.g., if a PWY-A belongs to Y number of SM classes, then there will be Y number of PWY-A nodes in the network ("PWY-A_1, "PWY-A_2", …, "PWY-A_Y") while each edge joins 2 significantly coexpressed pathways.

Melastoma dodecandrum Lour. shows antimicrobial activities
To assess the antimicrobial activity of M. dodecandrum, we sampled roots, petioles, leaves, vegetative branches, stems, fruits, and flowers from the Nanyang Technological University herb garden in Singapore (Figure 1A), and tested the antibacterial activity of methanol extracts derived from these organs.Bacterial growth inhibition studies indicated that the plant has antibacterial activity against gram-negative P. aeruginosa (Figure 1B; Table S2) and gram-positive S. aureus (Figure 1C; Table S2).When tested at 2% v/v (~4 mg/mL extract) against P. aeruginosa, the stems showed the highest activity with a 33.13% (ANOVA, p< 0.0001) inhibition in growth over 18 h as compared to the DMSO controls, followed by flowers (32.56%), fruits (27.9%), leaves (15.54%), and roots (11.26%) (Figure 1B).Similarly, the stems showed the highest activity against S. aureus with 50.21% growth inhibition (p< 0.05) (Figure 1C).The vegetative branch, fruit, roots, and leaves showed moderate activity in a descending order ranging from 30.02 -24.27% growth inhibition.In comparison, the flowers and petiole showed the lowest activity against S. aureus at 10.67 and 14.15% growth inhibition, respectively (Figure 1C).
Next, we evaluated the minimum inhibitory concentration (MIC), defined as the concentration at which no growth is observed after overnight incubation.While leaves did not have the strongest activity amongst all organs (Figure 1B), they were the most abundant and easiest to collect.Hence, leaf methanol extract was used to evaluate the MIC of the extract on P. aeruginosa via a broth microdilution method.Complete growth inhibition was achieved when the extract was used at 50 mg/mL (Figures 1E, F, Table S4).The inhibitory effect was concentration dependent, with > 50% growth inhibition by 1.563 mg/mL of the extract compared to treatment with equivalent volumes of DMSO (Figures 1E, F, Tables S4, S5).While 5% v/v DMSO treatment resulted in longer lag times, growth resumed by 720 min (Figure 1D).
The leaf methanol extract was bactericidal at 50 mg/mL, with complete cell death within four hours after treatment (Figure 1G, Table S6).In contrast, with the addition of equivalent 5% v/v DMSO, P. aeruginosa CFU recovered by 24 h (7.00 × 10 9 CFU/mL) to levels comparable to that of untreated samples (5.65 × 10 9 CFU/ mL) (Figure 1G).A lower cell count was measured from 0 -8 h with treatment with 25 mg/mL extract (Figure S1), but the bacterial density reached untreated control levels by 24 h.Finally, P. aeruginosa treated with 4 mg/mL of the plant extract or equivalent concentrations of DMSO showed similar growth across all time points (Figure S1), suggesting that a higher concentration of the plant extract is needed to achieve growth inhibition of bacteria.

M.dodecandrum extracts cause membrane damage and reduced swimming and swarming motility
To elucidate how the M. dodecandrum extracts affect the bacterial viability, we used LIVE/DEAD staining to assess potential membrane damage.The LIVE/DEAD stain consists of two nucleic acid staining dyes, propidium iodide (PI) and SYTO9.The green fluorescing SYTO9 can enter all cells, while the red fluorescing PI can only enter cells with compromised or damaged cytoplasmic membranes.PI, however, has a stronger affinity to DNA and hence can displace SYTO9 if present (Stiefel et al., 2015).While the kit is more accurately said to distinguish between cells with intact or compromised cell membranes, it is more often used to infer populations of living or dead cells under the assumption that all membrane-compromised cells are likely dead.
DMSO, when used at equivalent volumes of 0 -0.8% v/v and a treatment time of 10 minutes, did not result in increased uptake of PI (Figure 2A, nearly all cells are green), and the proportion of live cells was calculated to be at ~100% across all concentrations of DMSO (Figure 2B).In contrast, 10-minute treatment with 0 -8 mg/ mL of the plant extract resulted in a concentration-dependent increase in PI uptake and a corresponding reduction in the proportion of live cells from ~100% in untreated samples to 13.6% in 8 mg/mL treated samples (Figures 2A, B, nearly all cells are red, Table S7).
Antibacterial compounds can also negatively affect the motility of bacteria, effectively disrupting their pathogenesis (Erhardt, 2016).Microscopy revealed reduced motility of cells treated with M. dodecandrum plant extract compared to the DMSO control groups, even at low concentrations of 1 mg/mL, where most P. aeruginosa cells appeared alive following 10 mins of treatment (Video 1 and 2 for DMSO control and extract treatment, respectively).To quantify and evaluate the effect of plant extract treatment on P. aeruginosa motility, swarming and swimming motility assays were carried out (Figures 2C-F).Swarming motility was significantly impaired at 0.0625 mg/mL (p< 0.05) (Figures 2C, D; Table S7), while swimming motility was impaired at 0.25 mg/mL (p< 0.0001) (Figures 2E, F; Table S7).These results indicated that M. dodecandrum leaf extracts exhibit a multipronged effect on P. aeruginosa, inhibiting viability, growth, and motility.

Comparative transcriptomic analysis of antibiotic-and extract-induced gene expression changes in P. aeruginosa
Antibiotics can stimulate or depress gene expression in bacteria (Shitikov et al., 2022).Studying gene expression changes in responses to antibiotics can help the understanding of the mode of action of the various drug classes on bacterial adaptive ability, physiology, and metabolism.Since different antibiotics can elicit specific transcriptional responses (Hesketh et al., 2011), we hypothesized that gene expression analysis of P. aeruginosatreated M. dodecandrum extracts and several antibiotics would allow us to propose the mode of action of the plant extracts.
To see if transcriptome changes caused by our plant extracts are similar to the changes caused by common antibiotics, we performed a comparative RNA sequencing analysis.We treated P. aeruginosa with rifampicin (RNA synthesis inhibitor), triclosan (fatty acid synthesis inhibitor), ceftazidime (cell wall synthesis inhibitor), gentamicin (protein synthesis inhibitor), M. dodecandrum extracts and DMSO control (Table S8).Notably, rRNA depletion was not needed to obtain a high number (on average ~20 million) of reads mapping to protein-coding genes when increasing the sequencing depth to ~40 million reads (Table S9), indicating that this expensive step can be skipped if sequencing depth can be adjusted.

Concentrations of extracts and antimicrobial compounds used
for RNA-sequencing were selected based on results from microtiter plate growth curves (Figure S2).While 25 mg/mL of the extract did not fully inhibit P. aeruginosa growth by 24 h (Figures 1E, H), it was selected over 50 mg/mL as it inhibited growth to a similar extent within 18 h of growth without strong inhibition of growth associated with higher DMSO concentrations at 50 mg/mL (Figure 1F).Rifampicin, ciprofloxacin, and gentamicin were used at 1 × MIC of 32, 0.5, and 2 μg/mL respectively (Figure S2).Lastly, no growth inhibition was observed for triclosan over the range of concentration tested.Higher concentrations of triclosan treatment were not evaluated due to the precipitation of the compound in the media.As such, the highest concentration of 128 μg/mL, where some growth inhibition was observed, was used.DMSO was adjusted to the same concentration of 2.5% v/v for all samples to enable comparison with the control group.The resulting RNA-sequencing data was used to generate the gene expression matrix (Table S10).
the hierarchical clustering analysis of the 18 RNAseq samples revealed that M. dodecandrum extracts show the highest Euclidean distance to all other samples (Figure 3A, M. dodecandrum RNA-seq samples form an outgroup).Similarly, rifampicin and triclosan also showed unique transcriptomic responses, while DMSO control, ceftazidime, and gentamicin displayed a similar response (Figure 3A).These responses could be explained by the number of significantly differentially up-(red) and down-regulated genes (Figure 3B, BH-adj.p-value< 0.01).M.dodecandrum elicited the highest number of up-( 712) and down-regulated (757) genes, followed by rifampicin (331 up, 283 down) and triclosan (250, 135)(Table S11).Conversely, ceftazidime (40, 13) and gentamicin (99, 13) caused the mildest transcriptomic response.Overall, we observed that the M.dodecandrum extracts showed some similarities of the differentially expressed genes with rifampicin and triclosan (Figure S3).
To see which biological processes are differentially expressed, we performed a Gene Ontology enrichment analysis (see methods).We observed that the DEGs caused by triclosan resulted in the highest number of expressed GO terms (Figure S4).These terms comprised many metabolic, biosynthetic, respiration, and macromolecule localization processes.All 8 GO terms enriched in the plant extract-treated samples overlap with triclosan and are mostly metabolic and biosynthetic processes (translation, peptide, amide, protein, and nitrogen compounds).Conversely, rifampicin only induced changes in the type VI secretion system (Figure S4).Interestingly, gentamicin caused changes in terms involved in the high-temperature response (folding, response to heat, and temperature stimulus), while ceftazidime responses were external stimulus and DNA damage (Figure S4).Taken together, these results indicate the M.dodecandrum extracts may have an overlapping mode of action with triclosan (fatty acid synthesis inhibitor).

The antibacterial activity is conferred by unidentified metabolites
To identify the antibacterial compounds in M.dodecandrum, we performed a bioassay-guided purification on the MeOH extract from the leaves.The extract was further fractionated by preparative reversed-phase HPLC and yielded 8 subfractions (a-h).Antimicrobial activity was observed in the Parent fraction and PPT, however, no activity was found in the a-h fractions (Figure 4A), indicating that the activity is caused by more than one synergistically-acting metabolites, or that the active metabolite was diluted during the fractionation.
To study which metabolites might be present in the active fractions, we compared the active material (Parent and PPT) to the three inactive fractions f-h, using 1 H NMR and LC-HRMS.HR-ESI-MS (High-resolution electrospray ionization mass spectrometry) of fractions f-h showed a major m/z of 1402.1149[z=2]; 784.0776 [z=2] and 935.0793 [z=2], respectively (Figure 4B).A dominant fragment ion observed at m/z 301 indicated the existence of a hexahydroxydiphenoyl group (HHDP), which was characteristic of ellagitannins (Era et al., 2020).The MS/MS profiles were also in good agreement with that of hydrolyzable tannin casuarinin (m/z 935.0805 [z=1])(Figure 4C) (Chang et al., 2019).Coincidentally, fraction h presented the same m/z 935.0805, but a doubly charged [M-2H] 2-.Therefore, the mass of h was calculated to be 1872.1738and the molecular formula was determined as C 82 H 56 O 52 .From the predicted formula, fractions h, f (C 123 H 82 O 78 ), and g (C 68 H 50 O 44 ) were putatively identified as nobotanin analogs.The 1 H NMR spectra (Figure 4D) of the major compound in fraction H also resembles spectra containing nobotanins (Chang et al., 2019).However, as none of the eight fractions a-h is active, it is unlikely that the tannins are significantly conferring antibacterial activity.
The overlay of 1 H NMR spectra (Figure 4D) showed that the precipitation of the Parent fraction using 15% MeOH/H 2 O removed most of the hydrolysable tannins.indicating a complex mixture consisting of phenolic compounds (9 -11 ppm) and other unidentified compounds that could be responsible for the antibacterial activity.
Next, we investigated the types of metabolites produced by the detected SM pathways.The majority of the compounds belonged to phenylpropanoids, of which flavonoids are most abundant (Figure 5A).Terpenoids, and more specifically diterpenoids biosynthesis comes in a close second.These observations are concurrent with previous studies which report these compounds in the larger M.dodecandrumcae family (Serna and Martıńez, 2015).
To propose which specialized metabolites are most abundant in M.dodecandrum, we analyzed the gene expression of the SM pathways.Since SM pathways are thought to be under strong transcriptional control (Mutwil, 2020), we set out to investigate the gene expression patterns in M.dodecandrum.To this end, we isolated RNA from the same samples used to measure antibacterial activity (Figure 1A) and constructed a gene expression atlas for M.dodecandrum (Figure S5, Table S14).We then ranked the SM pathways within each organ based on its average TPM expression.To identify the 20 highest and lowest-ranked pathways, we ordered them according to the sum of their ranks in each organ (Figure 5B).The top 5 pathways included a sulfur-containing compound, two terpenoids, a phenylpropanoid, and a hormone synthesis pathway.There are some pathways that show organ-specific expression (e.g.young-flower-specific nitrogen-containing secondary compound, PWY-6442).
Since SM pathways expand and diversify by gene duplications (Hofberger et al., 2013), the predicted SM pathways might likely represent several separate, related pathways biosynthesizing related metabolites.For example, several enzymes implicated in flavonoid biosynthesis (PWY-7897, Figure 5C) form an aerial root-specific cluster (green clade).This suggests that these genes biosynthesize a flavonoid in aerial roots, which might differ from flavonoids produced in flowers (first five enzymes, pink and brown clade, Figure 5C).To investigate whether these functional sub-clusters are frequently found within the predicted pathways, we calculated the gene expression correlation (captured by Pearson Correlation Coefficient, PCC) within a pathway.We assumed that the median PCC should be high if the enzymes within a pathway are involved in the biosynthesis of the same metabolite.Surprisingly, while some of the pathways showed a high average PCC (blue line, Figure S6), many pathways showed average PCC values like randomly chosen gene pairs (orange line, average PCC ~0.25).Conversely, using tissuespecific gene expression data to identify the clusters containing >75% of annotated enzymatic reactions within a pathway resulted in a much higher average PCC value for nearly all pathways (green line, Figure S6D).This result shows that sequence similarity-based approaches such as PathwayTools can be supported by gene expression analyses to identify genes from tissue-specific pathways.
Transcriptional wiring of specialized metabolism in M.dodecandrum M.dodecandrum might express over 100 SM pathways (Table S12) and metabolites (Figure 4A).Since these pathways should have specific functions, we expect the functionally-related pathways to be expressed at the same time and place.To identify such patterns, we investigate which pathways are more connected in the coexpression network than expected by chance.We first set a PCC threshold of ≥0.6, as the PCC values obtained from the randomized pathway assignments were all<0.6 (Figure S5, randomized pathway assignment confidence interval is shown in orange).We then counted the number of co-expression edges connecting any two pathways and calculated which pathways were more connected than expected by chance (see methods).Interestingly, we observed high connectivity between the different pathways, where, e.g., most phenylpropanoid, terpenoid, and hormone biosynthesis pathways were significantly connected, indicating transcriptional coordination between the pathways (Figure 6A).Between 93 pathways, there were 251 significant connections.
To take a closer look at the connections that are found within and across the pathways, we investigated the gene co-expression networks of two connected pathways.The first pathway involves the biosynthesis of the plant hormone ethylene, which has numerous roles in plant development and stress responses (Chang, 2016).Ethylene is synthesized from L-methionine by three enzymatic steps comprising methionine adenosyltransferase, 1-aminocyclopropane-1carboxylate synthase (synthase), and 1-aminocyclopropane-1carboxylate oxidase (oxidase) (Chang, 2016) Interestingly, the coexpression network of the ethylene biosynthesis pathway shows several groups of co-expressed synthases and oxidases (Figure 6B, left pathway, blue edges connect co-expressed genes), suggesting that M.dodecandrum contains several ethylene biosynthesis pathways.The other pathway represents the biosynthesis of phytosterol terpenoids, which comprise multiple steps that convert cycloartenol to stigmasterol, crinosterol, and brassicasterol in Arabidopsis thaliana.This pathway in M.dodecandrum contains several dehydrogenases, methyltransferases, isomerases, demethylases, and desaturases that modify the cycloartenol (or a related compound) to other sterols (Figure 6B, right pathway).While PWY-2541 has been characterized in Arabidopsis, it is likely that this pathway produces other, related compounds in M.dodecandrum.Within each pathway, we also observed subclusters of genes that are co-expressed with the same gene(s) in the other pathway.While the analysis does not reveal any clear organ-specific clusters like in Figure 4C, this could suggest separate clusters of SM pathway expansion where subclustering occurs not only within a pathway but with other pathways too.The co-expression edges found between the ethylene and phytosterol pathways indicate that the two pathways are transcriptionally positively coordinated.This suggests that an increase in ethylene might result in increased biosynthesis of phytosterols.This is in line with a recent study that showed that exogenous ethylene led to a profound change in the ratio of stigmasterol to sitosterol (Markowski et al., 2022).

Discussion
In this study, we showed that M. dodecandrum shows antibacterial activities against gram-positive and -negative bacteria (Figure 1).The plant is known to produce a diverse range of chemical compounds and metabolites, including flavonoids, tannins and ellagitannins, phenylpropanoids, longchain fatty acids, aromatic acids, terpenoids, steroids, alkaloids, glycosides and monosaccharides (Zheng et al., 2021).Plant specialized metabolites can exert their antibacterial effects through various mechanisms.This includes but is not limited to, inhibition of bacterial nucleic acid and protein synthesis, inhibition of enzymes and respiration or metabolism, disruption of cell wall synthesis or cell membrane permeability or function, influencing the expression of virulence genes or chelation of essential metal ions such as iron (Rıós and Recio, 2005;Payne et al., 2007;Frey and Meyers, 2010;Tommasi et al., 2015;Kessler and Kalske, 2018;Singh et al., 2018;Lewis, 2020).
Results from LIVE/DEAD staining, in conjunction with viable CFU counts, suggest that the M. dodecandrum extract results in membrane damage (Figure 2), butviable cell counting suggests that there is no significant difference or drop in CFU compared to the control groups at 25 mg/mL (Figure 1H) or in viable cell counts in time-kill studies involving samples treated with 4 mg/mL w/v of the plant extract within an hour of treatment (Figure 1I).Recent studies suggest intermediate living/dead states in populations of bacteria stained with the LIVE/DEAD kit (Cushnie et al., 2020).Furthermore, there is evidence that, in some cases, PI staining may underestimate the viability of bacterial cells (Rosenberg et al., 2019).Taken together, it is likely that the increased PI staining following extract treatment is partly due to membrane damage and not actual cell death.
In addition to membrane damage, treatment with M. dodecandrum extract leads to reduced swimming and swarming motility in P. aeruginosa (Figures 2C-F, Video 1).Both forms of motility are flagella-mediated (Kazmierczak et al., 2015).The effects on motility may be related to the extract's effects on the bacterial membranes, as flagella rotation is driven by ion flux and membrane-active antibiotics that result in the alteration of membrane and changes in cation permeability have been suggested to reduce flagella activity and in turn, motility in P. fluorescens (Faust and Doetsch, 1971;Manson et al., 1980).Similarly, cationic peptides that target cell membranes may also affect flagella integrity, and some cationic peptides can reduce swimming and swarming motility in E. coli (Faust and Doetsch, 1971).Unlike swimming motility, swarming motility is a much more complex and tightly regulated adaptation process involving the upregulation of virulence gene expression and antibiotic resistance (Overhage et al., 2008).In addition to functional flagella, swarming is also dependent on quorum sensing and type IV villi (Köhler et al., 2000;Deźiel et al., 2003).Several plant metabolites have been reported to inhibit swarming via various mechanisms.For example, ginseng extract reduces swarming by reducing the production of QS signals in P. aeruginosa, while some flavonoids act as allosteric inhibitors of QS receptors (Wu et al., 2011;Paczkowski et al., 2017;Rütschlin and Böttcher, 2020).Some phenolic plant compounds or tannins can affect swarming but not swimming motility by affecting surfactant production (O' May and Tufenkji, 2011;Rütschlin and Böttcher, 2020).Beyond plant metabolites, antibiotics such as gentamicin and ciprofloxacin have also been reported to reduce swarming either via disruption of QS or through SOS responses that prevented polar chemosensory array formation (Irazoki et al., 2017;Rütschlin and Böttcher, 2020).
Comparison of gene expression responses of P.aeruginosatreated with M.dodecandrum extracts revealed that the plant extracts elicited Co-expression network analysis reveals a transcriptional association between specialized metabolic pathways and hormones.(A) Nodes represent pathways, while edges connect significantly connected pathways.Edges within the same pathway type are colored blue, while edges connecting different pathway types are colored gray.Node colors represent different pathway types.(B) Ethylene (left) and phytosterol (right) biosynthesis pathways.Nodes represent enzymes, while edges connect enzymes correlated with r>0.6.Edges within the same pathway are blue, while edges connecting two pathways are gray.
the most unique and highest number of differentially expressed compared to four antibiotics (Figures 3A, B).This is not surprising, as the plant extract contains hundreds of different compounds that might synergistically affect different aspects of P.aeruginosabiology.Surprisingly, the observed transcriptomics responses were not aligned with the primary mode of action of the antibiotics (Figure S4).For example, while rifampicin inhibits RNA biosynthesis, the observed transcriptomic changes affect the unrelated type VI secretion system (Figure S4).While these results indicate that transcriptomics might not be suitable to reveal the primary protein target of antibiotics, the downstream secondary effects can still be inferred (O'Rourke et al., 2020).Antibiotics with classical mechanism-of-actions can be clustered using their gene expression profiles (O'Rourke et al., 2020).Thus the overlapping differentially expressed biological processes between M.dodecandrum and triclosan could represent some of the similar downstream transcriptomic effects caused by triclosan (fatty acid biosynthesis inhibitor), which is in line with the observed membrane damage phenotypes (Figure S4).However, more time points and concentrations of the extracts and antibiotics should be compared to differentiate between primary and downstream effects.
To identify the metabolite(s) responsible for the antibacterial activity, we performed activity-guided fractionation.However, none of the purified fractions showed activity (Figure 5A), suggesting that either the activity is caused by multiple metabolites that together cause the activity or that the active metabolite is poorly soluble.The latter is more likely, as the insoluble material showed antibacterial activity (Figure 5A).The 1 H NMR analysis of the partially purified Parent fraction indicated the presence of hydrolysable tannins (Figure 5B), which are known to show broad antibacterial activities (Tomiyama et al., 2016), and other health-promoting activities due to their anti-oxidant properties (Abu Zarin et al., 2016).However, a deeper analysis using different fractionation methods showed that the activity was not due to the major hydrolysable tannins, but some other unidentified minor compounds.A larger scale of purification is needed to identify the active compounds.
The genomic and transcriptomic analysis revealed that M. dodecandrum contains many phenylpropanoid and terpenoid biosynthetic pathways (Figure 4A), and we were able to identify the most expressed pathways producing a sulfur-containing compound, two terpenoids, a phenylpropanoid (Figure 4B).Interestingly, within the predicted pathways, we observed several groups of highly co-expressed enzymes (Figures 4C, D), suggesting that these groups might represent several related pathways active in the different parts of the plant.Since enzymes biosynthesizing a given metabolite are expected to be expressed in the same cell, tissue and organ (Delli-Ponti et al., 2020), pathway inference approaches should be augmented with gene expression data.With >300,000 RNA-seq samples publicly available for hundreds of species (Julca et al., 2023), and several approaches that can be used to infer the biosynthetic pathways (Delli-Ponti et al., 2020), our and other groups can expect to provide more accurate inventories of specialized metabolism pathways.
Finally, we observed many co-expressed metabolic and hormonal pathways.This indicates that the specialized metabolic pathways are transcriptionally and functionally connected and that they are under hormonal control.This is in line with the existing knowledge linking many hormones to the biosynthesis of specialized metabolites (Chen et al., 2019;Markowski et al., 2022).Since co-expression analysis can identify hormonepathway pairs (Figure 6A), our analyses can generate testable hypotheses of how plants' hormones and metabolism are wired.
In conclusion, the high connectivity observed between various pathways, such as phenylpropanoid, terpenoid, and hormone biosynthesis pathways in M.dodecandrum, highlights the intricate transcriptional coordination between these pathways, suggesting that they work together to regulate key aspects of plant development, stress responses, and secondary metabolite production, ultimately contributing to the plant's adaptive capabilities and overall survival.Future studies could investigate how modulating plant hormones could affect the production of anti-bacterial metabolites.

FIGURE 1
FIGURE 1Antibacterial activities of M. dodecandrum.(A) Various organs of M. dodecandrum Lour were collected for evaluating antimicrobial activity.The antibacterial activity of each plant organ against (B) P. aeruginosa PAO1 and (C) S. aureus 25923 was assessed using 2% v/v of plant methanol extract (~4 mg/mL).The area under the growth curve (AUC) over 18 hours was measured and plotted.GC = Growth control (untreated bacteria in media).Boxplots represent data from at least two biological replicates, with differences in activity evaluated using one-way ANOVA with multiple comparisons against the 2% v/v DMSO control.*, P,< 0.05; ****, P,<0.0001.(D-G) Antibacterial activity of methanol plant extract from M. dodecandrum leaves was further evaluated in P.aeruginosa at higher concentrations of 0-50 mg/mL (D) and equivalent DMSO concentrations of 0-5% v/v (E) over 24 hours.The AUC over 18 hours was plotted in (F), the open and closed circles represent data from two biological replicates while error bars represent the standard deviation of the mean.(G) The time-kill activity of 50 mg/mL of methanol plant extract from M. dodecandrum leaves or equivalent volumes of DMSO was assessed over 24 hours, with data points from each biological replicate plotted.

FIGURE 2
FIGURE 2 Viability, swimming and swarming analysis if P.aeruginosa.(A) Confocal microscopy images of P. aeruginosa PAO1 cells treated with equivolume of 0.1% (top left), 0.4% (top middle), or 1.6% v/v (top right) DMSO and 1 mg/mL (bottom left), 4 mg/mL (bottom middle), or 16 mg/mL (bottom right) M. dodecandrum extract and stained with LIVE/DEAD stain.Green and red cells indicate live and dead cells, respectively.(B) Percentage of live cells following treatment with 0-8 mg/mL of M. dodecandrum leaf extract or an equal volume of DMSO, as determined using LIVE/DEAD stain and quantified using a microtiter plate reader.Two replicates were performed.(C) Swarming area experiment over 0-4 mg/mL M. dodecandrum extract ranges.(D) Normalized area values from the swarming area experiments.Two replicates were performed.(E) Swimming area experiment over 0-4 mg/mL M. dodecandrum extract ranges.F) Normalized area values from the swimming area experiments.Each data point (B, D, F) represents data from at least two independent experiments, with error bars representing the standard deviation of the mean.* = p ≤ 0.05, ** = p ≤ 0.01, **** = p ≤ 0.0001.
FIGURE 3 Gene expression responses of P.aeruginosa to M. dodecandrum extracts and four antibiotics.(A) Hierarchical clustering of RNA sequencing samples for DMSO, M.dodecandrum extract, and four antibiotics.Each sample type is represented by three replicates, indicated by its suffix.(B) The number of significantly differentially expressed genes (BH-adjusted p-value< 0.01), where up-and down-regulated genes are indicated in red and blue bars, respectively.
FIGURE 4 Bioassay results for HPLC fractions A-H, Parent fraction, and insoluble pellet (PPT).(A) Antibacterial activity of eight HPLC fractions(A-H), Parent, insoluble pellet (PPT) and DMSO control.(B) MS/MS profiles of the three fractions (F-H).(C) Fragmentation patterns of casuarinin.(D) Overlay 1 H NMR spectra of fractions H, PPT, and Parent, acquired in acetone-d 6 .
FIGURE 5Identification of specialized metabolic pathways in M. dodecandrum.(A) The number and types of putative specialized metabolic pathways identified by PathwayTools.(B) Identification of the top 5 most and least expressed pathways from each organ.The scale represents the rank of a pathway (white = lowest rank, dark blue = highest rank).The rows represent pathways, while the columns indicate the sampled organs (as in Figure1A).The lower rank (lighter color) indicates a higher expression.(C) Clustering analysis of flavonoid biosynthesis pathway (PWY-7897).