Integration of genomics, transcriptomics and metabolomics identifies candidate loci underlying fruit weight in loquat

Abstract

Fruit weight is an integral part of fruit quality and directly influences the commodity value and economic return of fruit crops. Despite the importance of fruit weight, its underlying molecular mechanisms remain understudied, especially for perennial fruit tree crops such as cultivated loquat (Eriobotrya japonica Lindl.). Auxin is known to regulate fruit development, but its role and metabolism during the development of loquat fruit remain obscure. In this study, we used a multi-omics approach, integrating whole-genome resequencing-based quantitative trait locus (QTL) mapping with an F₁ population, population genomics analysis using germplasm accessions, transcriptome analysis, and metabolic profiling to identify genomic regions potentially associated with fruit weight in loquat. We identified three major loci associated with fruit weight, supported by both QTL mapping and comparative genomic analysis between small- and large-fruited loquat cultivars. Comparison between two genotypes with contrasting fruit weight performance by transcriptomic and metabolic profiling revealed an important role for auxin in the regulation of fruit development, especially at the fruit enlargement stage. The multi-omics approach identified homologs of ETHYLENE INSENSITIVE 4 (EjEIN4) and TORNADO 1 (EjTRN1) as promising candidates for the control of fruit weight. Three single nucleotide polymorphism (SNP) markers were also closely associated with fruit weight. Results from this study provide insights into the genetic and metabolic controls of fruit weight in loquat from multiple perspectives. The candidate genomic regions, genes, and sequence variants will facilitate our understanding of the molecular basis of fruit weight and lay a foundation for future breeding and manipulation of fruit weight in loquat.

References

[1] Mauxion JP, Chevalier C, Gonzalez N . Complex cellular and molecular events determining fruit size. Trends Plant Sci. 2021; 26: 1023-38.

[2] Li YH, Zhang Z, Sun GM . Changes in cell number and cell size during pineapple (Ananas comosus L.) fruit development and their relationship with fruit size . Aust J Bot. 2010; 58: 673-8.

[3] Rawandoozi ZJ, Hartmann TP, Carpenedo S et al. Identification and characterization of QTLs for fruit quality traits in peach through a multi-family approach. BMC Genomics. 2020; 21: 1-18.

[4] Ozga JA, Reinecke DM . Hormonal interactions in fruit development. J Plant Growth Regul. 2003; 22: 73-81.

[5] Zheng WW, Kim YJ, Oh SM, Chun IJ . Anatomical analysis of fruit development of different-sized ‘Hongro’ and ‘Fuji’ apple fruits. Hortic Environ Biotechnol. 2009; 50: 160-5.

[6] Pattison RJ, Catalá C . Evaluating auxin distribution in tomato (Solanum lycopersicum) through an analysis of the PIN and AUX/LAX gene families . Plant J. 2012; 70: 585-98.

[7] Devoghalaere F, Doucen T, Guitton B et al. A genomics approach to understanding the role of auxin in apple (malus x domestica) fruit size control. BMC Plant Biol. 2012; 12: 1-15.

[8] Bu H, Yu W, Yuan H et al. Endogenous auxin content contributes to larger size of apple fruit. Front Plant Sci. 2020; 11: 592540. https://www.frontiersin.org/articles/10.3389/fpls.2020.592540/full.

[9] Agusti M, Almela V, Aznar M et al. Satsuma mandarin fruit size increased by 2, 4-DP. HortScience. 1994; 29: 279-81.

[10] Liao X, Li M, Liu B et al. Interlinked regulatory loops of ABA catabolism and biosynthesis coordinate fruit growth and ripening in woodland strawberry. Proc Natl Acad Sci U S A. 2018; 115: E11542-50.

[11] Wang H, Wu T, Liu J et al. PbGA20ox2 regulates fruit set and induces parthenocarpy by enhancing GA4 content . Front Plant Sci. 2020; 11: 113.

[12] Matsuo S, Kikuchi K, Fukuda M et al. Roles and regulation of cytokinins in tomato fruit development. J Exp Bot. 2012; 63: 5569-79.

[13] Xin T, Zhang Z, Li S et al. Genetic regulation of ethylene dosage for cucumber fruit elongation. Plant Cell. 2019; 31: 1063-76.

[14] Nitsch L, Kohlen W, Oplaat C et al. ABA-deficiency results in reduced plant and fruit size in tomato. J Plant Physiol. 2012; 169: 878-83.

[15] Fu FQ, Mao WH, Shi K et al. A role of brassinosteroids in early fruit development in cucumber. J Exp Bot. 2008; 59: 2299-308.

[16] Perrot-Rechenmann C, Napier RM . Auxins. Vitam Horm. 2005; 72: 203-33.

[17] Zhao Y, Christensen SK, Fankhauser C et al. A role for flavin monooxygenase-like enzymes in auxin biosynthesis. Science. 2001; 291: 306-9.

[18] Staswick PE, Tiryaki I, Rowe ML . Jasmonate response locus JAR1 and several related Arabidopsis genes encode enzymes of the firefly luciferase superfamily that show activity on jasmonic, salicylic, and indole-3-acetic acids in an assay for adenylation. Plant Cell. 2002; 14: 1405-15.

[19] Vanneste S, Friml J . Auxin: a trigger for change in plant development. Cell. 2009; 136: 1005-16.

[20] Tiwari SB, Hagen G, Guilfoyle T . The roles of auxin response factor domains in auxin-responsive transcription. Plant Cell. 2003; 15: 533-43.

[21] Lavy M, Estelle M . Mechanisms of auxin signaling. Development. 2016; 143: 3226-9.

[22] Tan X, Calderon-Villalobos LIA, Sharon M et al. Mechanism of auxin perception by the TIR1 ubiquitin ligase. Nature. 2007; 446: 640-5.

[23] Leyser O. Auxin signaling. Plant Physiol. 2018; 176: 465-79.

[24] Giovannoni JJ . Genetic regulation of fruit development and ripening. Plant Cell. 2004; 16: S170-80.

[25] Su L, Bassa C, Audran C et al. The auxin SlIAA17 transcriptional repressor controls fruit size via the regulation of endoreduplication-related cell expansion. Plant Cell Physiol. 2014; 55: 1969-76.

[26] Mounet F, Moing A, Kowalczyk M et al. Downregulation of a single auxin efflux transport protein in tomato induces precocious fruit development. J Exp Bot. 2012; 63: 4901-17.

[27] Guilfoyle TJ, Hagen G . Auxin response factors. Curr Opin Plant Biol. 2007; 10: 453-60.

[28] Ren Z, Li Z, Miao Q et al. The auxin receptor homologue in Solanum lycopersicum stimulates tomato fruit set and leaf morphogenesis. J Exp Bot. 2011; 62: 2815-26.

[29] De Jong M, Wolters-Arts M, Feron R et al. The Solanum lycopersicum auxin response factor 7 (SlARF7) regulates auxin signaling during tomato fruit set and development. Plant J. 2009; 57: 160-70.

[30] De Jong, Wolters-Arts M, Schimmel BCJ et al. Solanum lycopersicum AUXIN RESPONSE FACTOR 9 regulates cell division activity during early tomato fruit development . J Exp Bot. 2015; 66: 3405-16.

[31] Nardozza S, Cooney J, Boldingh HL et al. Phytohormone and transcriptomic analysis reveals endogenous cytokinins affect kiwifruit growth under restricted carbon supply. Meta. 2020; 10: 23.

[32] Frary A, Nesbitt TC, Grandillo S et al. fw2.2: a quantitative trait locus key to the evolution of tomato fruit size. Science. 2000; 289: 85-8.

[33] Cong B, Liu J, Tanksley SD . Natural alleles at a tomato fruit size quantitative trait locus differ by heterochronic regulatory mutations. Proc Natl Acad Sci U S A. 2002; 99: 13606-11.

[34] Liu J, Cong B, Tanksley SD . Generation and analysis of an artificial gene dosage series in tomato to study the mechanisms by which the cloned quantitative trait locus fw2.2 controls fruit size. Plant Physiol. 2003; 132: 292-9.

[35] Zygier S, Chaim AB, Efrati A et al. QTLs mapping for fruit size and shape in chromosomes 2 and 4 in pepper and a comparison of the pepper QTL map with that of tomato. Theor Appl Genet. 2005; 111: 437-45.

[36] Illa-Berenguer E, Van Houten J, Huang Z, van der Knaap E . Rapid and reliable identification of tomato fruit weight and locule number loci by QTL-seq. Theor Appl Genet. 2015; 128: 1329-42.

[37] Chakrabarti M, Zhang N, Sauvage C et al. A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc Natl Acad Sci U S A. 2013; 110: 17125-30.

[38] Li Q, Chakrabarti M, Taitano NK et al. Differential expression of SlKLUH controlling fruit and seed weight is associated with changes in lipid metabolism and photosynthesis-related genes. J Exp Bot. 2021; 72: 1225-44.

[39] Mu Q, Huang Z, Chakrabarti M et al. Fruit weight is controlled by cell size regulator encoding a novel protein that is expressed in maturing tomato fruits. PLoS Genet. 2017; 13: e1006930. https://journals.plos.org/plosgenetics/article?id=10.1371/journal.pgen.1006930.

[40] Potts SM, Khan MA, Han Y et al. Identification of quantitative trait loci (QTLs) for fruit quality traits in apple. Plant Mol Biol Rep. 2014; 32: 109-16.

[41] Duan N, Bai Y, Sun H et al. Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement. Nat Commun. 2017; 8: 1-11.

[42] Wu J, Li L-T, Li M et al. High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers. J Exp Bot. 2014; 65: 5771-81.

[43] Li X, Singh J, Qin M et al. Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus) . Plant Biotechnol J. 2019; 17: 1582-94.

[44] Li Y, Cao K, Zhu G et al. Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. Genome Biol. 2019; 20: 1-18.

[45] Shi P, Xu Z, Zhang S et al. Construction of a high-density SNP-based genetic map and identification of fruit-related QTLs and candidate genes in peach [ Prunus persica (L.) Batsch] . BMC Plant Biol. 2020; 20: 1-16.

[46] Doligez A, Bertrand Y, Farnos M et al. New stable QTLs for berry weight do not colocalize with QTLs for seed traits in cultivated grapevine (Vitis vinifera L.) . BMC Plant Biol. 2013; 13: 1-16.

[47] Mengist MF, Bostan H, Young E et al. High-density linkage map construction and identification of loci regulating fruit quality traits in blueberry. Hortic Res. 2021; 8: 1-13.

[48] Jiang S, An H, Luo J et al. Comparative analysis of transcriptomes to identify genes associated with fruit size in the early stage of fruit development in Pyrus pyrifolia . Int J Mol Sci. 2018; 19: 2342.

[49] Pei MS, Cao S-H, Wu L et al. Comparative transcriptome analyses of fruit development among pears, peaches, and strawberries provide new insights into single sigmoid patterns. BMC Plant Biol. 2020; 20: 1-15.

[50] Muñoz-Espinoza C, Di Genova, Correa J et al. Transcriptome profiling of grapevine seedless segregants during berry development reveals candidate genes associated with berry weight. BMC Plant Biol. 2016; 16: 1-17.

[51] Sharpe RH . Loquat: botany and horticulture. Hortic Rev. 2010; 23: 233.

[52] Su W, Jing Y, Lin S et al. Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe. Proc Natl Acad Sci U S A. 2021; 118: e2101767118.

[53] Liu Y, Zhang W, Xu C, Li X . Biological activities of extracts from loquat (Eriobotrya japonica Lindl.): a review . Int J Mol Sci. 2016; 17: 1983.

[54] Agusti M, Juan M, Almela V, Gariglio N . Loquat fruit size is increased through the thinning effect of naphthaleneacetic acid. Plant Growth Regul. 2000; 31: 167-71.

[55] Amorós A, Zapata P, Pretel MT et al. Role of naphthalene acetic acid and phenothiol treatments on increasing fruit size and advancing fruit maturity in loquat. Sci Hortic. 2004; 101: 387-98.

[56] Reig C, Mesejo C, Martínez-Fuentes A, Agustí M . Synthetic auxin 3, 5, 6-TPA increases fruit size of loquat (Eriobotrya japonica Lindl.) by reducing cell turgor pressure . Sci Hortic. 2016; 210: 213-9.

[57] Su W, Zhu Y, Zhang L et al. The cellular physiology of loquat (Eriobotrya japonica Lindl.) fruit with a focus on how cell division and cell expansion processes contribute to pome morphogenesis . Sci Hortic. 2017; 224: 142-9.

[58] Su W, Shao Z, Wang M et al. EjBZR1 represses fruit enlargement by binding to the EjCYP90 promoter in loquat . Hortic Res. 2021; 8: 1-14.

[59] Zhao C, Guo Y, Li S et al. Correlation and genetic analysis of fruit traits in Ninghaibai × Dafang loquat F1 hybrid population. J Fruit Sci. 2021; 38: 11.

[60] Jiang S, An H, Xu F, Zhang X . Chromosome-level genome assembly and annotation of the loquat (Eriobotrya japonica) genome . GigaScience. 2020; 9: giaa015.

[61] Liu D, Ma C, Hong W et al. Construction and analysis of high-density linkage map using high-throughput sequencing data. PLoS One. 2014; 9: e98855. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0098855.

[62] Cnops G, Neyt P, Raes J et al. The TORNADO1 and TORNADO2 genes function in several patterning processes during early leaf development in Arabidopsis thaliana . Plant Cell. 2006; 18: 852-66.

[63] Li Z, Chen H-W, Li Q-T et al. Three SAUR proteins SAUR76, SAUR77 and SAUR78 promote plant growth in Arabidopsis. Sci Rep. 2015; 5: 1-19.

[64] Liu S, Zhang Y, Feng Q et al. Tomato AUXIN RESPONSE FACTOR 5 regulates fruit set and development via the mediation of auxin and gibberellin signaling. Sci Rep. 2018; 8: 1-16.

[65] Hao YJ, Zhang Z, Kitashiba H et al. Molecular cloning and functional characterization of two apple S-adenosylmethionine decarboxylase genes and their different expression in fruit development, cell growth and stress responses. Gene. 2005; 350: 41-50.

[66] Cui X, Ge C, Wang R et al. The BUD2 mutation affects plant architecture through altering cytokinin and auxin responses in Arabidopsis. Cell Res. 2010; 20: 576-86.

[67] Kim JH, Choi D, Kende H . The AtGRF family of putative transcription factors is involved in leaf and cotyledon growth in Arabidopsis. Plant J. 2003; 36: 94-104.

[68] Beltramino M, Ercoli MF, Debernardi JM et al. Robust increase of leaf size by Arabidopsis thaliana GRF3-like transcription factors under different growth conditions. Sci Rep. 2018; 8: 1-13.

[69] Pattison RJ, Csukasi F, Catalá C . Mechanisms regulating auxin action during fruit development. Physiol Plant. 2014; 151: 62-72.

[70] Fenn MA, Giovannoni JJ . Phytohormones in fruit development and maturation. Plant J. 2021; 105: 446-58.

[71] Fukuda S, Ishimoto K, Terakami S et al. Genetic mapping of the loquat canker resistance gene pse-c in loquat (Eriobotrya japonica) . Sci Hortic. 2016; 200: 19-24.

[72] Fukuda S, Nagano Y, Matsuguma K et al. Construction of a high-density linkage map for bronze loquat using RAD-Seq. Sci Hortic. 2019; 251: 59-64.

[73] Hu Z, Deng G, Mou H et al. A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo . DNA Res. 2018; 25: 1-10.

[74] Shi TX, Li R-Y, Zheng R et al. Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum) . BMC Genomics. 2021; 22: 1-12.

[75] Gisbert AD, Martínez-Calvo J, Llácer G et al. Development of two loquat [ Eriobotrya japonica (Thunb.) Lindl.] linkage maps based on AFLPs and SSR markers from different Rosaceae species . Mol Breed. 2009; 23: 523-38.

[76] Fukuda S, Ishimoto K, Sato S et al. Genetic mapping of the loquat canker resistance locus in bronze loquat (Eriobotrya deflexa) . Tree Genet Genomes. 2014; 10: 875-83.

[77] Mejía N, Gebauer M, Muñoz L et al. Identification of QTLs for seedlessness, berry size, and ripening date in a seedless x seedless table grape progeny. Am J Enol Vitic. 2007; 58: 499-507.

[78] Costantini L, Battilana J, Lamaj F et al. Berry and phenology-related traits in grapevine (Vitis vinifera L.): from quantitative trait loci to underlying genes . BMC Plant Biol. 2008; 8: 1-17.

[79] Müller A, Düchting P, Weiler EW . A multiplex GC-MS/MS technique for the sensitive and quantitative single-run analysis of acidic phytohormones and related compounds, and its application to Arabidopsis thaliana . Planta. 2002; 216: 44-56.

[80] Agustí M, Gariglio N, Castillo A et al. Effect of the synthetic auxin 2, 4-DP on fruit development of loquat. Plant Growth Regul. 2003; 41: 129-32.

[81] Ulmasov T, Hagen G, Guilfoyle TJ . Dimerization and DNA binding of auxin response factors. Plant J. 1999; 19: 309-19.

[82] Audran-Delalande C, Bassa C, Mila I et al. Genome-wide identification, functional analysis and expression profiling of the aux/IAA gene family in tomato. Plant Cell Physiol. 2012; 53: 659-72.

[83] Guilfoyle TJ, Ulmasov T, Hagen G . The ARF family of transcription factors and their role in plant hormone-responsive transcription. Cell Mol Life Sci. 1998; 54: 619-27.

[84] Liscum E, Reed J . Genetics of aux/IAA and ARF action in plant growth and development. Plant Mol Biol. 2002; 49: 387-400.

[85] Ishida T, Thitamadee S, Hashimoto T . Twisted growth and organization of cortical microtubules. J Plant Res. 2007; 120: 61-70.

[86] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[87] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[88] Van Ooijen J. MapQTL^® 6, Software for the Mapping of Quantitative Trait Loci in Experimental Populations of Diploid Species . Wageningen: Kyazma B.V.; 2009.

[89] Baum BR . PHYLIP: phylogeny inference package. Version 3.2. Joel Felsenstein. Q Rev Biol. 1989; 64: 164-6.

[90] Purcell S, Neale B, Todd-Brown K et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81: 559-75.

[91] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[92] Kim D, Paggi JM, Park C et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15.

[93] Pertea M, Pertea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[94] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550.

[95] Liu Q, Chen T-T, Xiao D-W et al. OsIAGT1 is a glucosyltransferase gene involved in the glucose conjugation of auxins in rice . Rice (NY). 2019; 12: 1-13.

[96] Fischer S, Brunk BP, Chen F et al. Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups. Curr Protoc Bioinformatics. 2011; 35: 1-19.

[97] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[98] Jiang Y, Zhu Y, Zhang L et al. EjTFL1 genes promote growth but inhibit flower bud differentiation in loquat . Front Plant Sci. 2020; 11: 576.