QTL mapping and genomic analyses of earliness and fruit ripening traits in a melon recombinant inbred lines population supported by de novo assembly of their parental genomes

Abstract

Earliness and ripening behavior are important attributes of fruits on and off the vine, and affect quality and preference of both growers and consumers. Fruit ripening is a complex physiological process that involves metabolic shifts affecting fruit color, firmness, and aroma production. Melon is a promising model crop for the study of fruit ripening, as the full spectrum of climacteric behavior is represented across the natural variation. Using Recombinant Inbred Lines (RILs) population derived from the parental lines “Dulce” (reticulatus, climacteric) and “Tam Dew” (inodorus, non-climacteric) that vary in earliness and ripening traits, we mapped QTLs for ethylene emission, fruit firmness and days to flowering and maturity. To further annotate the main QTL intervals and identify candidate genes, we used Oxford Nanopore long-read sequencing in combination with Illumina short-read resequencing, to assemble the parental genomes de-novo. In addition to 2.5 million genome-wide SNPs and short InDels detected between the parents, we also highlight here the structural variation between these lines and the reference melon genome. Through systematic multi-layered prioritization process, we identified 18 potential polymorphisms in candidate genes within multi-trait QTLs. The associations of selected SNPs with earliness and ripening traits were further validated across a panel of 177 diverse melon accessions and across a diallel population of 190 F1 hybrids derived from a core subset of 20 diverse parents. The combination of advanced genomic tools with diverse germplasm and targeted mapping populations is demonstrated as a way to leverage forward genetics strategies to dissect complex horticulturally important traits.

References

[1] LeRoy P, Lyon CB . Inheritance studies on duration of developmental stages in crosses within the genus Lycopersicon. J Agric Res. 1941; 63: 129-48.

[2] Buckler ES, Holland JB, Bradbury PJ et al. The genetic architecture of maize flowering time. Science. 2009; 325: 714-8.

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

[4] Tucker G, Yin X, Zhang A et al. Ethylene† and fruit softening. Food Quality and Safety. 2017; 1: 253-67.

[5] Barry CS, Giovannoni J . Ethylene and fruit ripening. J Plant Growth Regul. 2007; 26: 143-59.

[6] Yang SF, Hoffman NE . Ethylene biosynthesis and its regulation in higher plants. Annu Rev Plant Physiol. 1984; 35: 155-89.

[7] Klee HJ, Giovannoni JJ . Genetics and control of tomato fruit ripening and quality attributes. Annu Rev Genet. 2011; 45: 41-59.

[8] Balague C, Watson CF, Turner AJ et al. Isolation of a ripening and wound-induced cDNA from Cucumis melo L. encoding a protein with homology to the ethylene-forming enzyme. Eur J Biochem. 1993; 212: 27-34.

[9] Lasserre E, Bouquin T, Hernandez JA et al. Structure and expression of three genes encoding ACC oxidase homologs from melon (Cucumis melo L.). Mol Gen Genet. 1996; 251: 81-90.

[10] Miki T, Yamamoto M, Nakagawa H et al. Nucleotide sequence of a cDNA for 1-aminocyclopropane-1-carboxylate synthase from melon fruits. Plant Physiol. 1995; 107: 297-8.

[11] Zheng X, Wolff D, Crosby K . Genetics of ethylene biosynthesis and restriction fragment length polymorphisms (RFLPs) of ACC oxidase and synthase genes in melon (Cucumis melo L.). Theor Appl Genet. 2002; 105: 397-403.

[12] Sato-Nara K, Yuhashi KI, Higashi K et al. Stage- and tissue-specific expression of ethylene receptor homolog genes during fruit development in muskmelon. Plant Physiol. 1999; 120: 321-30.

[13] Pech JC, Bouzayen M, Latché A . Climacteric fruit ripening: ethylene-dependent and independent regulation of ripening pathways in melon fruit. Plant Sci. 2008; 175: 114-20.

[14] Pereira L, Santo Domingo M, Ruggieri V et al. Genetic dissection of climacteric fruit ripening in a melon population segregating for ripening behavior. Hortic Res. 2020; 7: 187.

[15] Perpiñá G, Cebolla-Cornejo J, Esteras C et al. ‘MAK-10’: a long shelf-life Charentais breeding line developed by Introgression of a genomic region from Makuwa melon. HortScience. 2017; 52: 1633-8.

[16] Périn C, Gomez-Jimenez MC, Hagen L et al. Molecular and genetic characterization of a non-climacteric phenotype in melon reveals two loci conferring altered ethylene response in fruit. Plant Physiol. 2002; 129: 300-9.

[17] Vegas J, Garcia-Mas J, Monforte AJ . Interaction between QTLs induces an advance in ethylene biosynthesis during melon fruit ripening. Theor Appl Genet. 2013; 126: 1531-44.

[18] Moreno E, Obando JM, Dos-Santos N et al. Candidate genes and QTLs for fruit ripening and softening in melon. Theor Appl Genet. 2008; 116: 589-602.

[19] Perpiñá G, Esteras C, Gibon Y et al. A new genomic library of melon introgression lines in a cantaloupe genetic background for dissecting desirable agronomical traits. BMC Plant Biol. 2016; 16: 154.

[20] Ríos P, Argyris J, Vegas J et al. ETHQV6.3 is involved in melon climacteric fruit ripening and is encoded by a NAC domain transcription factor. Plant J. 2017; 91: 671-83.

[21] Saladié M, Cañizares J, Phillips MA et al. Comparative transcriptional profiling analysis of developing melon (Cucumis melo L.) fruit from climacteric and non-climacteric varieties. BMC Genomics. 2015; 16: 440.

[22] Leida C, Moser C, Esteras C et al. Variability of candidate genes, genetic structure and association with sugar accumulation and climacteric behavior in a broad germplasm collection of melon (Cucumis melo L.). BMC Genet. 2015; 16: 1-17.

[23] Nimmakayala P, Tomason YR, Abburi VL et al. Genome-wide differentiation of various melon horticultural groups for use in GWAS for fruit firmness and construction of a high resolution genetic map. Front Plant Sci. 2016; 7: 1437.

[24] Song YH, Ito S, Imaizumi T . Flowering time regulation: photoperiod- and temperature-sensing in leaves. Trends Plant Sci. 2013; 18: 575-83.

[25] Singh D, Debnath P, Sane AP et al. Expression of the tomato WRKY gene, SlWRKY23, alters root sensitivity to ethylene, auxin and JA and affects aerial architecture in transgenic Arabidopsis. Physiol Mol Biol Plants. 2020; 26: 1187-99.

[26] Zhang L, Chen L, Yu D . Transcription factor WRKY75 interacts with DELLA proteins to affect flowering. Plant Physiol. 2018; 176: 790-803.

[27] Lindhout P, Van Heusden S, Pet G et al. Perspectives of molecular marker assisted breeding for earliness in tomato. Euphytica. 1994; 79: 279-86.

[28] Kemble JM, Gardner RG . Inheritance of shortened fruit maturation in the cherry tomato Cornell 871213-1 and its relation to fruit size and other components of earliness. J Am Soc Hortic Sci. 2019; 117: 646-50.

[29] Nakano H, Kobayashi N, Takahata K et al. Quantitative trait loci analysis of the time of floral initiation in tomato. Sci Hortic. 2016; 201: 199-210.

[30] Monforte AJ, Oliver M, Gonzalo MJ et al. Identification of quantitative trait loci involved in fruit quality traits in melon (Cucumis melo L.). Theor Appl Genet. 2004; 108: 750-8.

[31] Garcia-Mas J, Benjak A, Sanseverino W et al. The genome of melon (Cucumis melo L.). Proc Natl Acad Sci. 2012; 109: 11872-7.

[32] Castanera R, Ruggieri V, Pujol M et al. An improved melon reference genome with single-molecule sequencing uncovers a recent burst of transposable elements with potential impact on genes. Front Plant Sci. 2020; 10: 1-10.

[33] Ruggieri V, Alexiou KG, Morata J et al. An improved assembly and annotation of the melon (Cucumis melo L.) reference genome. Sci Rep. 2018; 8: 1-9.

[34] Zhao G, Lian Q, Zhang Z et al. A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits. Nat Genet. 2019; 51: 1607-15.

[35] Liu S, Gao P, Zhu Q et al. Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon. Plant Biotechnol J. 2020; 18: 2545-58.

[36] Yano R, Nonaka S, Ezura H . Melonet-DB, a grand RNA-Seq gene expression atlas in melon (Cucumis melo L.). Plant Cell Physiol. 2018; 59: e4.

[37] Zheng Y, Wu S, Bai Y et al. Cucurbit genomics database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops. Nucleic Acids Res. 2019; 47: D1128-36.

[38] Pereira L, Ruggieri V, Pérez S et al. QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map. BMC Plant Biol. 2018; 18: 1-17.

[39] Oren E, Tzuri G, Dafna A et al. High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo. Theor Appl Genet. 2020; 133: 1927-45. https://doi.org/10.1007/s00122-020-03567-3.

[40] Branham SE, Daley J, Levi A, a. QTL mapping and marker development for tolerance to Sulfur Phytotoxicity in melon (Cucumis melo). Front Plant Sci. 2020; 11: 1-9.

[41] Badouin H, Gouzy J, Grassa CJ et al. The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature. 2017; 546: 148-52.

[42] Haberer G, Kamal N, Bauer E et al. European maize genomes highlight intraspecies variation in repeat and gene content. Nat Genet. 2020; 52: 950-7.

[43] Song JM, Guan Z, Hu J et al. Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus. Nat Plants. 2020; 6: 34-45. https://doi.org/10.1038/s41477-019-0577-7.

[44] Liu Y, Du H, Li P et al. Pan-genome of wild and cultivated soybeans. Cell. 2020; 182: 162-176.e13.

[45] Gao L, Gonda I, Sun H et al. The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat Genet. 2019; 51: 1044-51.

[46] Hufford MB, Seetharam AS, Woodhouse MR et al. De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes. Science. 2021; 373: 655-62.

[47] Bayer PE, Golicz AA, Scheben A et al. Plant pan-genomes are the new reference. Nat Plants. 2020; 6: 914-20. https://doi.org/10.1038/s41477-020-0733-0.

[48] Alonge M, Wang X, Benoit M et al. Major impacts of widespread structural variation on gene expression and crop improvement in tomato. Cell. 2020; 182: 145-161.e23.

[49] Sanseverino W, Hénaff E, Vives C et al. Transposon insertion, structural variations and SNPs contribute to the evolution of the melon genome. Mol Biol Evol. 2015; 32: 2760-74. https://doi.org/10.1093/molbev/msv152.

[50] Demirci S, Fuentes RR, van Dooijeweert W et al. Chasing breeding footprints through structural variations in Cucumis melo and wild relatives. G3 Genes|Genomes|Genetics. 2021; 11: 1689-99.

[51] Zhang H, Li X, Yu H et al. A high-quality melon genome assembly provides insights into genetic basis of fruit trait improvement. iScience. 2019; 22: 16-27.

[52] Yano R, Ariizumi T, Nonaka S et al. Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression. Commun Biol. 2020; 3: 1-13.

[53] Freilich S, Lev S, Gonda I et al. Systems approach for exploring the intricate associations between sweetness, color and aroma in melon fruits. BMC Plant Biol. 2015; 15: 71. https://doi.org/10.1186/s12870-015-0449-x.

[54] Cingolani P, Platts A, Wang LL et al. A program for annotating and predicting the effects of single nucleotide polymorphisms. Fly. 2012; 6: 80-92.

[55] Galpaz N, Gonda I, Shem-Tov D et al. Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping. Plant J. 2018; 94: 169-91.

[56] Dafna A, Halperin I, Oren E et al. Underground Heterosis for melons yield. J Exp Bot. 2021; 72: 6205-18.

[57] Oren E, Tzuri G, Vexler L et al. Multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon. J Exp Bot. 2019; 70: 3781-94.

[58] Flynn JM, Hubley R, Goubert C et al. RepeatModeler2 for automated genomic discovery of transposable element families. Proc Natl Acad Sci. 2020; 117: 9451-7.

[59] Shumate A, Salzberg SL . Liftoff: accurate mapping of gene annotations. Bioinformatics. 2021; 37: 1639-43.

[60] Song B, Sang Q, Wang H et al. Complement genome annotation lift over using a weighted sequence alignment strategy. Front Genet. 2019; 10: 1-10.

[61] Nattestad M, Schatz MC . Assemblytics: a web analytics tool for the detection of variants from an assembly. Bioinformatics. 2016; 32: 3021-3.

[62] Fernández-Trujillo JP, Obando-Ulloa JM, Martínez JA et al. Climacteric and non-climacteric behavior in melon fruit. Postharvest Biol Technol. 2008; 50: 125-34.

[63] Jaganathan D, Bohra A, Thudi M et al. Fine mapping and gene cloning in the post-NGS era: advances and prospects. Theor Appl Genet. 2020; 133: 1791-810.

[64] Stacey D, Fauman EB, Ziemek D et al. ProGeM: a framework for the prioritization of candidate causal genes at molecular quantitative trait loci. Nucleic Acids Res. 2019; 47: e3-13.

[65] Nguyen CV, Vrebalov JT, Gapper NE et al. Tomato GOLDEN2-LIKE transcription factors reveal molecular gradients that function during fruit development and ripening. Plant Cell. 2014; 26: 585-601.

[66] Pereira-Garcia L Genetic dissection of fruit quality and ripening traits in melon. 2018. doctoral dissertation, Universitat Autonoma de Barcelona (UAB). 2018.

[67] Gachomo EW, Jimenez-Lopez JC, Baptiste LJ et al. GIGANTUS1 (GTS1), a member of Transducin/WD40 protein superfamily, controls seed germination, growth and biomass accumulation through ribosome-biogenesis protein interactions in Arabidopsis thaliana. BMC Plant Biol. 2014; 14: 37.

[68] Wu HC, Bulgakov VP, Jinn TL . Pectin methylesterases: cell wall remodeling proteins are required for plant response to heat stress. Front Plant Sci. 2018; 871: 1-21.

[69] Schwihla M, Korbei B . The beginning of the end: initial steps in the degradation of plasma membrane proteins. Front Plant Sci. 2020; 11: 680.

[70] Jayakodi M, Padmarasu S, Haberer G et al. The barley pan-genome reveals the hidden legacy of mutation breeding. Nature. 2020; 588: 284-9.

[71] Bariah I, Keidar-Friedman D, Kashkush K . Identification and characterization of largescale genomic rearrangements during wheat evolution. PLoS One. 2020; 15: e0231323.

[72] Gur A, Tzuri G, Meir A et al. Genome-wide linkage-disequilibrium mapping to the candidate gene level in melon (Cucumis melo). Sci Rep. 2017; 7: 9770.

[73] Cahaner A, Hillel J . Estimating heritability and genetic correlation between traits from generations F2 and F 3 of self-fertilizing species: a comparison of three methods. Theor Appl Genet. 1980; 58: 33-8.

[74] Bradbury PJ, Zhang Z, Kroon DE et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633-5.

[75] Taylor J, Butler D . R package ASMap : efficient genetic linkage map construction and diagnosis. J Stat Softw. 2017; 79: 1-29. https://doi.org/10.18637/jss.v079.i06.

[76] Fulton TM, Chunwongse J, Tanksley SD . Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Mol Biol Report. 1995; 13: 207-9.

[77] Wick RR, Judd LM, Gorrie CL et al. Completing bacterial genome assemblies with multiplex MinION sequencing. Microb Genom. 2017; 3: 1-7.

[78] Kolmogorov M, Yuan J, Lin Y et al. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019; 37: 540-6.

[79] Vaser R, Sović I, Nagarajan N et al. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 2017; 27: 737-46.

[80] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

[81] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34: 3094-100.

[82] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[83] Alonge M, Soyk S, Ramakrishnan S et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 2019; 20: 224.

[84] Shen W, Le S, Li Y et al. SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One. 2016; 11: 1-11.

[85] Mikheenko A, Prjibelski A, Saveliev V et al. Versatile genome assembly evaluation with QUAST-LG. Bioinformatics. 2018; 34: i142-50.

[86] Waterhouse RM, Seppey M, Simão FA et al. BUSCO applications from quality assessments to gene prediction and Phylogenomics. Mol Biol Evol. 2018; 35: 543-8.

[87] Marçais G, Delcher AL, Phillippy AM et al. MUMmer4: A fast and versatile genome alignment system. PLoS Comput Biol. 2018; 14: 1-14.

[88] Lyons E, Freeling M . How to usefully compare homologous plant genes and chromosomes as DNA sequences. Plant J. 2008; 53: 661-73.

[89] Haug-Baltzell A, Stephens SA, Davey S et al. SynMap2 and SynMap3D: Web-based whole-genome synteny browsers. In: Bioinformatics. 2017; 33: 2197-8.

[90] Vaser R, Adusumalli S, Leng SN et al. SIFT missense predictions for genomes. Nat Protoc. 2016; 11: 1-9.

[91] Choi Y, Chan AP . PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015; 31: 2745-7.

[92] Broman KW, Wu H, Sen S et al. R/QTL: QTL mapping in experimental crosses. Bioinformatics. 2003; 19: 889-90.