Identification of key gene networks controlling organic acid and sugar metabolism during watermelon fruit development by integrating metabolic phenotypes and gene expression profiles

Abstract

The organoleptic qualities of watermelon fruit are defined by the sugar and organic acid contents, which undergo considerable variations during development and maturation. The molecular mechanisms underlying these variations remain unclear. In this study, we used transcriptome profiles to investigate the coexpression patterns of gene networks associated with sugar and organic acid metabolism. We identified 3 gene networks/modules containing 2443 genes highly correlated with sugars and organic acids. Within these modules, based on intramodular significance and Reverse Transcription Quantitative polymerase chain reaction (RT-qPCR), we identified 7 genes involved in the metabolism of sugars and organic acids. Among these genes, Cla97C01G000640, Cla97C05G087120 and Cla97C01G018840 (r² = 0.83 with glucose content) were identified as sugar transporters (SWEET, EDR6 and STP) and Cla97C03G064990 (r² = 0.92 with sucrose content) was identified as a sucrose synthase from information available for other crops. Similarly, Cla97C07G128420, Cla97C03G068240 and Cla97C01G008870, having strong correlations with malic (r² = 0.75) and citric acid (r² = 0.85), were annotated as malate and citrate transporters (ALMT7, CS, and ICDH). The expression profiles of these 7 genes in diverse watermelon genotypes revealed consistent patterns of expression variation in various types of watermelon. These findings add significantly to our existing knowledge of sugar and organic acid metabolism in watermelon.

References

[1] Faostat, F. Available online: http://www.fao.org/faostat/en/#data.QC (accessed on January 2018) (2017).

[2] Guo, S. et al. Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits. Nat. Genet. 51, 1-8 (2019).

[3] Akashi, K. et al. Spatial accumulation pattern of citrulline and other nutrients in immature and mature watermelon fruits. J. Sci. Food Agric. 97, 479-487 (2017).

[4] Hayashi, T. et al. l-Citrulline and l-arginine supplementation retards the progression of high-cholesterol-diet-induced atherosclerosis in rabbits. Proc. Natl Acad. Sci. USA 102, 13681-13686 (2005).

[5] Collins, J. K. et al. Watermelon consumption increases plasma arginine concentrations in adults. Nutrition 23, 261-266 (2007).

[6] Saminathan, T. et al. Metagenomic and metatranscriptomic analyses of diverse watermelon cultivars reveal the role of fruit associated microbiome in carbohydrate metabolism and ripening of mature fruits. Front. Plant Sci. 9, 4 (2018).

[7] Wechter, W. P. et al. Gene expression in developing watermelon fruit. BMC Genomics 9, 275 (2008).

[8] Borsani, J. et al. Carbon metabolism of peach fruit after harvest: changes in enzymes involved in organic acid and sugar level modifications. J. Exp. Bot. 60, 1823-1837 (2009).

[9] Zhu, Q. et al. Comparative transcriptome analysis of two contrasting watermelon genotypes during fruit development and ripening. BMC Genomics 18, 3 (2017).

[10] Montero, T. M., Mollá, E. M., Esteban, R. M. & López-Andréu, F. J. Quality attributes of strawberry during ripening. Sci. Horticulturae 65, 239-250 (1996).

[11] Jiang, C. C., Fang, Z. Z., Zhou, D. R., Pan, S. L. & Ye, X. F. Changes in secondary metabolites, organic acids and soluble sugars during the development of plum fruit cv. ‘Furongli’ (Prunus salicina Lindl). J. Sci. Food Agric. 99, 1010-1019 (2019).

[12] Jawad, U. M. et al. Expression pattern of sugars and organic acids regulatory genes during watermelon fruit development. Sci. Horticulturae 265, 109102 (2020).

[13] Miron, D. & Schaffer, A. A. Sucrose phosphate synthase, sucrose synthase, and invertase activities in developing fruit of Lycopersicon esculentum Mill. and the sucrose accumulating Lycopersicon hirsutum Humb. and Bonpl. Plant Physiol. 95, 623-627 (1991).

[14] Granot, D., David-Schwartz, R. & Kelly, G. Hexose kinases and their role in sugar-sensing and plant development. Front. Plant Sci. 4, 44 (2013).

[15] D'Ambrosio, C. et al. Proteomic analysis of apricot fruit during ripening. J. Proteom. 78, 39-57 (2013).

[16] dos Santos, R. S. et al. Genetic regulation and the impact of omics in fruit ripening. Plant Omics 8, 78 (2015).

[17] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 9, 559 (2008).

[18] Bai, Y., Dougherty, L., Cheng, L., Zhong, G.-Y. & Xu, K. Uncovering co-expression gene network modules regulating fruit acidity in diverse apples. BMC genomics 16, 612 (2015).

[19] Zhang, Q. et al. Transcriptional regulatory networks controlling taste and aroma quality of apricot (Prunus armeniaca L.) fruit during ripening. BMC genomics 20, 1-15 (2019).

[20] Li, H. et al. MdMYB8 is associated with flavonol biosynthesis via the activation of the MdFLS promoter in the fruits of Malus crabapple. Horticulture Res. 7, 1-13 (2020).

[21] Gao, L. et al. Comparative transcriptome analysis reveals key genes potentially related to soluble sugar and organic acid accumulation in watermelon. PLoS ONE 13, e0190096 (2018).

[22] Bartolozzi, F., Bertazza, G., Bassi, D. & Cristoferi, G. Simultaneous determination of soluble sugars and organic acids as their trimethylsilyl derivatives in apricot fruits by gas-liquid chromatography. J. Chromatogr. A 758, 99-107 (1997).

[23] Shi, C.-Y. et al. Citrus PH5-like H+-ATPase genes: identification and transcript analysis to investigate their possible relationship with citrate accumulation in fruits. Front. Plant Sci. 6, 135 (2015).

[24] Wang, N. et al. Transcriptomic analysis of red-fleshed apples reveals the novel role of MdWRKY11 in flavonoid and anthocyanin biosynthesis. J. Agric. Food Chem. 66, 7076-7086 (2018).

[25] Anders, S., Pyl, P. T. & Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166-169 (2015).

[26] Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L. & Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621 (2008).

[27] Andino, G. K., Gribskov, M., Anderson, D. L., Evans, J. D. & Hunt, G. J. Differential gene expression in Varroa jacobsoni mites following a host shift to European honey bees (Apis mellifera). BMC Genomics 17, 926 (2016).

[28] Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc.: Ser. B (Methodol.) 57, 289-300 (1995).

[29] Anders, S. & Huber, W. Differential expression analysis for sequence count data. Nature Preceedings 11, R106 (2010).

[30] Mao, X., Cai, T., Olyarchuk, J. G. & Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787-3793 (2005).

[31] Zhang, B. & Horvath, S. A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 4, 17 (2005).

[32] Kariyanna, B. et al. Identification of suitable reference genes for normalization of RT-qPCR data in eggplant fruit and shoot borer (Leucinodes orbonalis Guenée). Biologia 75, 289-297 (2020).

[33] Kong, Q. et al. Identification of suitable reference genes for gene expression normalization in qRT-PCR analysis in watermelon. PloS ONE 9, e90612 (2014).

[34] Gao, L. et al. Comparative transcriptome analysis reveals key genes potentially related to soluble sugar and organic acid accumulation in watermelon. PLoS ONE 13, e0190096 (2018).

[35] Li, M., Feng, F. & Cheng, L. Expression patterns of genes involved in sugar metabolism and accumulation during apple fruit development. PLoS ONE 7, e33055 (2012).

[36] Smeekens, S. Sugar-induced signal transduction in plants. Annu. Rev. Plant Biol. 51, 49-81 (2000).

[37] Lombardo, V. A. et al. Metabolic profiling during peach fruit development and ripening reveals the metabolic networks that underpin each developmental stage. Plant Physiol. 157, 1696-1710 (2011).

[38] Famiani, F., Battistelli, A., Moscatello, S., Cruz-Castillo, J. G. & Walker, R. P. The organic acids that are accumulated in the flesh of fruits: occurrence, metabolism and factors affecting their contents-a review. Rev. Chapingo Ser. Horticultura 21, 97-128 (2015).

[39] Yao, Y.-X. et al. Molecular cloning of three malic acid related genes MdPEPC, MdVHA-A, MdcyME and their expression analysis in apple fruits. Sci. Horticulturae 122, 404-408 (2009).

[40] Ren, Y. et al. A tonoplast sugar transporter underlies a sugar accumulation QTL in watermelon. Plant Physiol. 176, 836-850 (2018).

[41] Peng, Q. et al. Functional analysis reveals the regulatory role of PpTST1 encoding tonoplast sugar transporter in sugar accumulation of peach fruit. Int. J. Mol. Sci. 21, 1112 (2020).

[42] Li, S.-j et al. The Citrus transcription factor, CitERF13, regulates citric acid accumulation via a protein-protein interaction with the vacuolar proton pump, CitVHA-c4. Sci. Rep. 6, 20151 (2016).

[43] De Angeli, A. et al. The vacuolar channel VvALMT9 mediates malate and tartrate accumulation in berries of Vitis vinifera. Planta 238, 283-291 (2013).

[44] Etienne, A., Génard, M., Lobit, P., Mbeguié-A-Mbeguié, D. & Bugaud, C. What controls fleshy fruit acidity? A review of malate and citrate accumulation in fruit cells. J. Exp. Bot. 64, 1451-1469 (2013).

[45] Etienne, C. et al. Candidate genes and QTLs for sugar and organic acid content in peach [Prunus persica (L.) Batsch]. Theor. Appl. Genet. 105, 145-159 (2002).

[46] Weber, A. et al. Identification, purification, and molecular cloning of a putative plastidic glucose translocator. Plant Cell 12, 787-801 (2000).

[47] Martinoia, E., Massonneau, A. & Frangne, N. Transport processes of solutes across the vacuolar membrane of higher plants. Plant Cell Physiol. 41, 1175-1186 (2000).

[48] Wang, H. X. et al. A Golgi-localized hexose transporter is involved in heterotrimeric G protein-mediated early development in Arabidopsis. Mol. Biol. Cell 17, 4257-4269 (2006).

[49] Chen, L.-Q. et al. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science 335, 207-211 (2012).

[50] Chardon, F. et al. Leaf fructose content is controlled by the vacuolar transporter SWEET17 in Arabidopsis. Curr. Biol. 23, 697-702 (2013).

[51] Klemens, P. A. et al. Overexpression of a proton-coupled vacuolar glucose exporter impairs freezing tolerance and seed germination. N. phytologist 202, 188-197 (2014).

[52] Poschet, G. et al. A novel Arabidopsis vacuolar glucose exporter is involved in cellular sugar homeostasis and affects the composition of seed storage compounds. Plant Physiol. 157, 1664-1676 (2011).

[53] Terrier, N., Deguilloux, C., Sauvage, F.-X., Martinoia, E. & Romieu, C. Proton pumps and anion transport in Vitis vinifera: the inorganic pyrophosphatase plays a predominant role in the energization of the tonoplast. Plant Physiol. Biochem. 36, 367-377 (1998).

[54] Ma, B. et al. Genes encoding aluminum-activated malate transporter II and their association with fruit acidity in apple. The Plant Genome 8, 1-14 (2015).

[55] Bai, Y. et al. A natural mutation-led truncation in one of the two aluminum-activated malate transporter-like genes at the Ma locus is associated with low fruit acidity in apple. Mol. Genet. genomics 287, 663-678 (2012).

[56] Berüter, J., Feusi, M. E. S. & Rüedi, P. Sorbitol and sucrose partitioning in the growing apple fruit. J. Plant Physiol. 151, 269-276 (1997).

[57] Ludewig, F. & Flügge, U.-I. Role of metabolite transporters in source-sink carbon allocation. Front. Plant Sci. 4, 231 (2013).