Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping

A matter of smell: The complex regulation of aroma production in melon

Melon fruit flavor is one of the most valuable traits for consumers. Aroma, formed by volatile organic compounds (VOCs), is a major component of flavor but has been neglected in breeding programs because of its complex regulation. Although the genetic regulation of VOCs biosynthesis is not fully understood, several advances have been recently achieved. VOCs originate from the degradation of fatty acids, aminoacids and terpenes, and the role of newly described enzymes, transcription factors and putative regulators is here discussed. Furthermore, ethylene plays a key role in fruit aroma production in melon, triggering the conversion of green-flavored aldehydes into fruity-flavored esters. A current challenge is to understand the ethylene-independent regulation of VOCs formation. Environmental conditions and human processing can also shape the melon volatile profile, and future research should focus on studying the effect of climate change in aroma formation.

A wild melon reference genome provides novel insights into the domestication of a key gene responsible for melon fruit acidity

KEY MESSAGE

A wild melon reference genome elucidates the genomic basis of fruit acidity domestication. Structural variants (SVs) have been reported to impose major effects on agronomic traits, representing a significant contributor to crop domestication. However, the landscape of SVs between wild and cultivated melons is elusive and how SVs have contributed to melon domestication remains largely unexplored. Here, we report a 379-Mb chromosome-scale genome of a wild progenitor melon accession "P84", with a contig N50 of 14.9 Mb. Genome comparison identifies 10,589 SVs between P84 and four cultivated melons with 6937 not characterized in previously analysis of 25 melon genome sequences. Furthermore, the population-scale genotyping of these SVs was determined in 1175 accessions, and 18 GWAS signals including fruit acidity, fruit length, fruit weight, fruit color and sex determination were detected. Based on these genotyped SVs, we identified 3317 highly diverged SVs between wild and cultivated melons, which could be the potential SVs associated with domestication-related traits. Furthermore, we identify novel SVs affecting fruit acidity and proposed the diverged evolutionary trajectories of CmPH, a key regulator of melon fruit acidity, during domestication and selection of different populations. These results will offer valuable resources for genomic studies and genetic improvement in melon.

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species

KEY MESSAGE

Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Clpf encodes pentatricopeptide repeat protein (PPR5) and regulates pink flesh color in watermelon (Citrullus lanatus L.)

KEY MESSAGE

The gene controlling pink flesh in watermelon was finely mapped to a 55.26-kb region on chromosome 6. The prime candidate gene, Cla97C06G122120 (ClPPR5), was identified through forward genetics. Carotenoids offer numerous health benefits; while, they cannot be synthesized by the human body. Watermelon stands out as one of the richest sources of carotenoids. In this study, genetic generations derived from parental lines W15-059 (red flesh) and JQ13-3 (pink flesh) revealed the presence of the recessive gene Clpf responsible for the pink flesh (pf) trait in watermelon. Comparative analysis of pigment components and microstructure indicated that the disparity in flesh color between the parental lines primarily stemmed from variations in lycopene content, as well as differences in chromoplast number and size. Subsequent bulk segregant analysis (BSA-seq) and genetic mapping successfully narrowed down the Clpf locus to a 55.26-kb region on chromosome 6, harboring two candidate genes. Through sequence comparison and gene expression analysis, Cla97C06G122120 (annotated as a pentatricopeptide repeat, PPR) was predicted as the prime candidate gene related to pink flesh trait. To further investigate the role of the PPR gene, its homologous gene in tomato was silenced using a virus-induced system. The resulting silenced fruit lines displayed diminished carotenoid accumulation compared with the wild-type, indicating the potential regulatory function of the PPR gene in pigment accumulation. This study significantly contributes to our understanding of the forward genetics underlying watermelon flesh traits, particularly in relation to carotenoid accumulation. The findings lay essential groundwork for elucidating mechanisms governing pigment synthesis and deposition in watermelon flesh, thereby providing valuable insights for future breeding strategies aimed at enhancing fruit quality and nutritional value.

T2T reference genome assembly and genome-wide association study reveal the genetic basis of Chinese bayberry fruit quality

Chinese bayberry (Myrica rubra or Morella rubra; 2n = 16) produces fruit with a distinctive flavor, high nutritional, and economic value. However, previous versions of the bayberry genome lack sequence continuity. Moreover, to date, no large-scale germplasm resource association analysis has examined the allelic and genetic variations determining fruit quality traits. Therefore, in this study, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar 'Zaojia' using PacBio HiFi long reads. The resulting 292.60 Mb T2T genome, revealed 8 centromeric regions, 15 telomeres, and 28 345 genes. This represents a substantial improvement in the genome continuity and integrity of Chinese bayberry. Subsequently, we re-sequenced 173 accessions, identifying 6 649 674 single nucleotide polymorphisms (SNPs). Further, the phenotypic analyses of 29 fruit quality-related traits enabled a genome-wide association study (GWAS), which identified 1937 SNPs and 1039 genes significantly associated with 28 traits. An SNP cluster pertinent to fruit color was identified on Chr6: 3407532 to 5 153 151 bp region, harboring two MYB genes (MrChr6G07650 and MrChr6G07660), exhibiting differential expression in extreme phenotype transcriptomes, linked to anthocyanin synthesis. An adjacent, closely linked gene, MrChr6G07670 (MLP-like protein), harbored an exonic missense variant and was shown to increase anthocyanin production in tobacco leaves tenfold. This SNP cluster, potentially a quantitative trait locus (QTL), collectively regulates bayberry fruit color. In conclusion, our study presented a complete reference genome, uncovered a suite of allelic variations related to fruit-quality traits, and identified functional genes that could be harnessed to enhance fruit quality and breeding efficiency of bayberries.

Genetic regulation of volatile production in two melon introgression line collections with contrasting ripening behavior

The importance of melon aroma in determining fruit quality has been highlighted in recent years. The fruit volatile profile is influenced by the type of fruit ripening. Non-climacteric fruits contain predominantly aldehydes, while climacteric fruits mainly produce esters. Several genes have been described to participate in volatile organic compounds (VOCs) biosynthesis pathways, but knowledge in this area is still incomplete. In this work we analysed the volatile profile of two reciprocal Introgression Line (IL) collections generated from a cross between 'Piel de Sapo' (PS) and 'Védrantais' (VED) melons, differing in their aroma profile and ripening behaviour. SPME GC-MS was performed to identify genes responsible for VOCs formation. More than 1000 QTLs for many volatiles were detected taken together both populations. Introgressions on chromosomes 3, 5, 6, 7 and 8 modified ester-aldehyde balance and were correlated to ripening changes in both genetic backgrounds. Some previously identified QTLs for fruit ripening might be involved in these phenotypes, such as ETHQV8.1 on chromosome 8 and ETHQV6.3 on chromosome 6. PS alleles on chromosomes 2, 6, 10 and 11 were found to increase ester content when introgressed in VED melons. Terpenes showed to be affected by several genomic regions not related to ripening. In addition, several candidate genes have been hypothesized to be responsible for some of the QTLs detected. The analysis of volatile compounds in two reciprocal IL collections has increased our understanding of the relationship between ripening and aroma and offers valuable plant material to improve food quality in melon breeding programs.

A high-density linkage map construction in guava (Psidium guajava L.) using genotyping by sequencing and identification of QTLs for leaf, peel, and pulp color in an intervarietal mapping population

Psidium guajava L. is an important fruit crop in the tropical and subtropical regions of the world. The advanced breeding methods are not employed for important commercial traits like peel and pulp color, seed hardiness, fruit size, etc., due to the scarcity of genome-wide molecular markers and high-density linkage maps. In this study, we employed single-nucleotide polymorphism (SNP) markers and identified quantitative trait loci (QTL) regions that are associated with color traits of leaf, peel, and pulp in the guava intervarietal mapping population. The mapping population was developed from the contrasting genotypes of fruit and leaf color. Variations in color among the segregating hybrids were recorded both visually and using a Color reader. A high-density linkage map of guava was constructed using the SNP markers from genotyping by sequencing (GBS) of 150 hybrid individuals of the cross 'Arka Poorna' (green) x 'Purple Local' (purple). The integrated linkage map consisted of 1426 SNPs mapped on 11 linkage groups (LG), spanning a total distance of around 730 cM with an average of 129.6 markers per LG. Through QTL analysis for color traits, a minor QTL region was identified for visually scored leaf color and peel color on LG1, whereas a major QTL was detected for pulp color in LG4. The Hunter color values (L* and, a*) also had major QTLs with overlapping marker intervals for leaf and peel colors, establishing the association of SNP markers to the trait. The QTLs harbored genes and transcription factors involved in lycopene and anthocyanin pigment biosynthesis. This is the first report of a high-density linkage map based on SNP markers in guava and QTL mapping for color characters in leaf, fruit peel and pulp. The genotyping information generated in this study can aid in genetic engineering and marker-assisted breeding in guava.

The effect of environmental factors on the genetic differentiation of Cucurbita ficifolia populations based on whole-genome resequencing

BACKGROUND

Cucurbita ficifolia is one of the squash species most resistant to fungal pathogens, and has especially high resistance to melon Fusarium wilt. This species is therefore an important germplasm resource for the breeding of squash and melon cultivars.

RESULTS

Whole-genome resequencing of 223 individuals from 32 populations in Yunnan Province, the main cucurbit production area in China, was performed and 3,855,120 single-nucleotide polymorphisms (SNPs) and 1,361,000 InDels were obtained. SNP analysis suggested that levels of genetic diversity in C. ficifolia were high, but that different populations showed no significant genetic differentiation or geographical structure, and that individual C. ficifolia plants with fruit rinds of a similar color did not form independent clusters. A Mantel test conducted in combination with geographical distance and environmental factors suggested that genetic distance was not correlated with geographical distance, but had a significant correlation with environmental distance. Further associations between the genetic data and five environmental factors were analyzed using whole-genome association analysis. SNPs associated with each environmental factor were investigated and genes 250 kb upstream and downstream from associated SNPs were annotated. Overall, 15 marker-trait-associated SNPs (MTAs) and 293 genes under environmental selection were identified. The identified genes were involved in cell membrane lipid metabolism, macromolecular complexes, catalytic activity and other related aspects. Ecological niche modeling was used to simulate the distribution of C. ficifolia across time, from the present and into the future. We found that the area suitable for C. ficifolia changed with the changing climate in different periods.

CONCLUSIONS

Resequencing of the C. ficifolia accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs identified in this study suggest that environmental factors mediated the formation of the population structure of C. ficifolia in China. These SNPs and Indels might also contribute to the variation in important pathways of genes for important agronomic traits such as yield, disease resistance and stress tolerance. Moreover, the genome resequencing data and the genetic markers identified from 223 accessions provide insight into the genetic variation of the C. ficifolia germplasm and will facilitate a broad range of genetic studies.

Regulation of climacteric fruit ripening in melon: recent advances and future challenges

Fruit ripening is a complex and highly regulated process where tomato and strawberry have been the model species classically used for studying climacteric and non-climacteric fleshy fruit ripening types, respectively. Melon has emerged as an alternative ripening model because climacteric and non-climacteric cultivars exist, which makes it possible to dissect the regulation of ripening using a genetic approach. Several quantitative trait loci that regulate climacteric fruit ripening have been identified to date, and their combination in both climacteric and non-climacteric genetic backgrounds resulted in lines with different ripening behaviors, demonstrating that the climacteric intensity can be genetically modulated. This review discusses our current knowledge of the physiological changes observed during melon climacteric fruit ripening such as ethylene production, fruit abscission, chlorophyll degradation, firmness, and aroma, as well as their complex genetic control. From pioneer experiments in which ethylene biosynthesis was silenced, to the recent genetic edition of ripening regulators, current data suggest that the climacteric response is determined by the interaction of several loci under quantitative inheritance. The exploitation of the rich genetic diversity of melon will enable the discovery of additional genes involved in the regulation of the climacteric response, ultimately leading to breeding aromatic melon fruits with extended shelf life.

Molecular Marker-Assisted Mapping, Candidate Gene Identification, and Breeding in Melon (Cucumis melo L.): A Review

Melon (Cucumis melo L.) is an important crop that is cultivated worldwide for its fleshy fruit. Understanding the genetic basis of a plant's qualitative and quantitative traits is essential for developing consumer-favored varieties. This review presents genetic and molecular advances related to qualitative and quantitative phenotypic traits and biochemical compounds in melons. This information guides trait incorporation and the production of novel varieties with desirable horticultural and economic characteristics and yield performance. This review summarizes the quantitative trait loci, candidate genes, and development of molecular markers related to plant architecture, branching patterns, floral attributes (sex expression and male sterility), fruit attributes (shape, rind and flesh color, yield, biochemical compounds, sugar content, and netting), and seed attributes (seed coat color and size). The findings discussed in this review will enhance demand-driven breeding to produce cultivars that benefit consumers and melon breeders.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

Identification of a novel locus C2 controlling canary yellow flesh color in watermelons

The flesh color of watermelon is an important trait that is determined by carotenoid composition and affects consumers' fruit desirability. Although a complete dominant control by C locus (Cllcyb) for canary yellow flesh (CY) over red flesh has been reported, red and CY colors frequently appear as a mixed pattern in the same flesh (incomplete canary yellow, ICY) in F1 and inbred lines carrying dominant C alleles. Therefore, we examined the genetic control of the mixed color pattern in ICY using whole-genome resequencing of three ICY (ICY group) and three CY inbred lines (CY group), as well as genetic linkage mapping of an F2 population. The segregation pattern in 135 F2 plants indicated that CY is controlled by a single locus (named C 2) dominant over ICY. The whole-genome resequencing of ICY and CY inbred lines revealed an ICY/CY-specific region of approximately 27.60-27.88 Mb on Chr. 2 that was polymorphic between the ICY and CY groups. Our genetic map, using nine cleaved amplified polymorphic sequence markers developed based on the single-nucleotide polymorphisms from the ICY/CY-specific region, confirmed that C 2 is located on Chr. 2 and cosegregated with the marker (M7) derived from a non-synonymous single-nucleotide polymorphism of the pentatricopeptide repeat (PPR) gene (ClPPR, Cla97C02G039880). Additionally, 27 watermelon inbred lines of ICY, CY, and red flesh were evaluated using previously reported Cllcyb (C locus)-based markers and our C 2 locus-linked ClPPR-based marker (M7). As a result, dominant alleles at the C 2 locus were required to produce CY, in addition to dominant alleles at the C locus, while a recessive homozygous genotype at the C locus gave the red flesh irrespective of the genotype at the C 2 locus. Using a ClPPR-based cleaved amplified polymorphic sequence developed in this study and Cllcyb-based markers, watermelon cultivars with CY, ICY, and red flesh could be successfully discerned, implying that the combined use of these markers will be efficient for marker-assisted selection of flesh color in watermelon breeding.

QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon

Melon (Cucumis melo L.) is an economically important Cucurbitaceae crop grown around the globe. The sweetness of melon is a significant factor in fruit quality and consumer appeal, and the soluble solids content (SSC) is a key index of melon sweetness. In this study, 146 recombinant inbred lines (RILs) derived from two oriental melon materials with different levels of sweetness containing 1427 bin markers, and 213 melon accessions containing 1,681,775 single nucleotide polymorphism (SNP) markers were used to identify genomic regions influencing SSC. Linkage mapping detected 10 quantitative trait loci (QTLs) distributed on six chromosomes, seven of which were overlapped with the reported QTLs. A total of 211 significant SNPs were identified by genome-wide association study (GWAS), 138 of which overlapped with the reported QTLs. Two new stable, co-localized regions on chromosome 3 were identified by QTL mapping and GWAS across multiple environments, which explained large phenotypic variance. Five candidate genes related to SSC were identified by QTL mapping, GWAS, and qRT-PCR, two of which were involved in hydrolysis of raffinose and sucrose located in the new stable loci. The other three candidate genes were involved in raffinose synthesis, sugar transport, and production of substrate for sugar synthesis. The genomic regions and candidate genes will be helpful for molecular breeding programs and elucidating the mechanisms of sugar accumulation.

Identification of chromosome region and candidate genes for canary-yellow flesh (Cyf) locus in watermelon (Citrullus lanatus)

Genetic control of fruit flesh color in watermelon is complex, and significant knowledge gaps still exist. In the present study, we investigated the genetic basis of canary-yellow flesh color in watermelon inbred line PI 635597 using a segregating population derived from a cross between PI 635597 and another inbred line, Cream of Saskatchewan (pale yellow flesh color). We showed that a single dominant gene controls the canary-yellow flesh color for the Cyf (canary-yellow flesh) trait. Bulk segregant analysis (BSA) and fine genetic mapping narrowed down the Cyf locus to a 79.62-kb region on chromosome 6, which harbors 10 predicted genes. Sequence variation analysis in the promoter and coding regions and gene expression analysis in both parental lines and selected watermelon accessions with diverse fruit flesh colors support Cla97C06G122050 (unknown protein) and Cla97C06G122120 (pentatricopeptide repeat) as predicted candidate genes for the Cyf locus. Marker-assisted selection and sequence alignment showed that the Cyf locus could differentiate canary-yellow flesh and pale-yellow flesh. Our results indicate that the Cyf locus might be responsible for canary-yellow flesh color and carotenoid accumulation levels.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Thorough Characterization of ETHQB3.5, a QTL Involved in Melon Fruit Climacteric Behavior and Aroma Volatile Composition

The effect of the QTL involved in climacteric ripening ETHQB3.5 on the fruit VOC composition was studied using a set of Near-Isogenic Lines (NILs) containing overlapping introgressions from the Korean accession PI 16375 on the chromosome 3 in the climacteric 'Piel de Sapo' (PS) genetic background. ETHQB3.5 was mapped in an interval of 1.24 Mb that contained a NAC transcription factor. NIL fruits also showed differences in VOC composition belonging to acetate esters, non-acetate esters, and sulfur-derived families. Cosegregation of VOC composition (23 out of 48 total QTLs were mapped) and climacteric ripening was observed, suggesting a pleiotropic effect of ETHQB3.5. On the other hand, other VOCs (mainly alkanes, aldehydes, and ketones) showed a pattern of variation independent of ETHQB3.5 effects, indicating the presence of other genes controlling non-climacteric ripening VOCs. Network correlation analysis and hierarchical clustering found groups of highly correlated compounds and confirmed the involvement of the climacteric differences in compound classes and VOC differences. The modification of melon VOCs may be achieved with or without interfering with its physiological behavior, but it is likely that high relative concentrations of some type of ethylene-dependent esters could be achieved in climacteric cultivars.

Two independent loss-of-function mutations in anthocyanidin synthase homeologous genes are responsible for the all-green phenotype of sweet basil

Sweet basil, Ocimum basilicum L., is an important culinary herb grown worldwide. Although basil is green, many landraces, breeding lines, and exotic cultivars have purple stems and flowers. This anthocyanin pigmentation is unacceptable in traditional Italian basil used for Pesto sauce production. In the current study, we aimed to resolve the genetics that underlines the different colors. We used the recently published sweet basil genome to map quantitative trait loci (QTL) for flower and stem color in a bi-parental F₂ population. It was found that the pigmentation is governed by a single QTL, harboring an anthocyanidin synthase (ANS) gene (EC 1.14.20.4). Further analysis revealed that the basil genome harbors two homeologous ANS genes, each carrying a loss-of-function mutation. ObANS1 carries a single base pair insertion resulting in a frameshift, and ObANS2 carries a missense mutation within the active site. In the purple-flower parent, ANS1 is functional, and ANS2 carries a nonsense mutation. The functionality of the ObANS1 active allele was validated by complementation assay in an Arabidopsis ANS mutant. Moreover, we have restored the functionality of the missense-mutated ObANS2 using site-directed activation. We found that the non-functional alleles were expressed to similar levels as the functional allele, suggesting polyploids invest futile effort in expressing non-functional genes, offsetting their advantageous redundancy. This work demonstrated the usefulness of the genomics and genetics of basil to understand the basic mechanism of metabolic traits and raise fundamental questions in polyploid plant biology.

Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo)

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F₂ populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

Fruit Morphology and Ripening-Related QTLs in a Newly Developed Introgression Line Collection of the Elite Varieties 'Védrantais' and 'Piel de Sapo'

Melon is an economically important crop with widely diverse fruit morphology and ripening characteristics. Its diploid sequenced genome and multiple genomic tools make this species suitable to study the genetic architecture of fruit traits. With the development of this introgression line population of the elite varieties 'Piel de Sapo' and 'Védrantais', we present a powerful tool to study fruit morphology and ripening traits that can also facilitate characterization or pyramidation of QTLs in inodorous melon types. The population consists of 36 lines covering almost 98% of the melon genome, with an average of three introgressions per chromosome and segregating for multiple fruit traits: morphology, ripening and quality. High variability in fruit morphology was found within the population, with 24 QTLs affecting six different traits, confirming previously reported QTLs and two newly detected QTLs, FLQW5.1 and FWQW7.1. We detected 20 QTLs affecting fruit ripening traits, six of them reported for the first time, two affecting the timing of yellowing of the rind (EYELLQW1.1 and EYELLQW8.1) and four at the end of chromosome 8 affecting aroma, abscission and harvest date (EAROQW8.3, EALFQW8.3, ABSQW8.3 and HARQW8.3). We also confirmed the location of several QTLs, such as fruit-quality-related QTLs affecting rind and flesh appearance and flesh firmness.

A multi-omics framework reveals strawberry flavor genes and their regulatory elements

Flavor is essential to consumer preference of foods and is an increasing focus of plant breeding programs. In fruit crops, identifying genes underlying volatile organic compounds has great promise to accelerate flavor improvement, but polyploidy and heterozygosity in many species have slowed progress. Here we use octoploid cultivated strawberry to demonstrate how genomic heterozygosity, transcriptomic intricacy and fruit metabolomic diversity can be treated as strengths and leveraged to uncover fruit flavor genes and their regulatory elements. Multi-omics datasets were generated including an expression quantitative trait loci map with 196 diverse breeding lines, haplotype-phased genomes of a highly-flavored breeding selection, a genome-wide structural variant map using five haplotypes, and volatile genome-wide association study (GWAS) with > 300 individuals. Overlaying regulatory elements, structural variants and GWAS-linked allele-specific expression of numerous genes to variation in volatile compounds important to flavor. In one example, the functional role of anthranilate synthase alpha subunit 1 in methyl anthranilate biosynthesis was supported via fruit transient gene expression assays. These results demonstrate a framework for flavor gene discovery in fruit crops and a pathway to molecular breeding of cultivars with complex and desirable flavor.

High-resolution genome-wide association study of a large Czech collection of sweet cherry (Prunus avium L.) on fruit maturity and quality traits

In sweet cherry (Prunus avium L.), quantitative trait loci have been identified for fruit maturity, colour, firmness, and size to develop markers for marker-assisted selection. However, resolution is usually too low in those analyses to directly target candidate genes, and some associations are missed. In contrast, genome-wide association studies are performed on broad collections of accessions, and assemblies of reference sequences from Tieton and Satonishiki cultivars enable identification of single nucleotide polymorphisms after whole-genome sequencing, providing high marker density. Two hundred and thirty-five sweet cherry accessions were sequenced and phenotyped for harvest time and fruit colour, firmness, and size. Genome-wide association studies were used to identify single nucleotide polymorphisms associated with each trait, which were verified in breeding material consisting of 64 additional accessions. A total of 1 767 106 single nucleotide polymorphisms were identified. At that density, significant single nucleotide polymorphisms could be linked to co-inherited haplotype blocks (median size ~10 kb). Thus, markers were tightly associated with respective phenotypes, and individual allelic combinations of particular single nucleotide polymorphisms provided links to distinct phenotypes. In addition, yellow-fruit accessions were sequenced, and a ~ 90-kb-deletion on chromosome 3 that included five MYB10 transcription factors was associated with the phenotype. Overall, the study confirmed numerous quantitative trait loci from bi-parental populations using high-diversity accession populations, identified novel associations, and genome-wide association studies reduced the size of trait-associated loci from megabases to kilobases and to a few candidate genes per locus. Thus, a framework is provided to develop molecular markers and evaluate and characterize genes underlying important agronomic traits.

An allelic variant in the ACS7 gene promotes primary root growth in watermelon

Gene mining in a C. lanatus × C. amarus population revealed one gene, ACS7, linked to primary root elongation in watermelon. Watermelon is a xerophytic crop characterized by a long primary root and robust lateral roots. Therefore, watermelon serves as an excellent model for studying root elongation and development. However, the genetic mechanism underlying the primary root elongation in watermelon remains unknown. Herein, through bulk segregant analysis we identified a genetic locus, qPRL.Chr03, controlling primary root length (PRL) using two different watermelon species (Citrullus lanatus and Citrullus amarus) that differ in their root architecture. Fine mapping revealed that xaa-Pro dipeptidase and 1-aminocyclopropane-1-carboxylate synthase 7 (ACS7) are candidate regulators of the primary root growth. Allelic variation in the delimited region among 193 watermelon accessions indicated that the long-root alleles might only exist in C. amarus. Interestingly, the discrepancy in PRL among the C. amarus accessions was clearly associated with a nonsynonymous single nucleotide polymorphism variant within the ACS7 gene. The ACS7 expression and ethylene levels in the primary root tips suggested that ethylene is a negative regulator of root elongation in watermelon, as supported by the application of 1-aminocyclopropane-1-carboxylate (ACC, the ethylene precursor) or 2-aminoethoxyvinyl glycine (AVG, an ACS inhibitor). To the best of our knowledge, these findings provide the first description of the genetic basis of root elongation in watermelon. The detected markers of the ACS7 gene will facilitate marker-assisted selection for the PRL trait to improve water and nutrient use efficacy in watermelon and beyond.

Genome-based high-resolution mapping of fusarium wilt resistance in sweet basil

Fusarium wilt of basil is a disease of sweet basil (Ocimum basilicum L.) plants caused by the fungus Fusarium oxysporum f. sp. basilici (FOB). Although resistant cultivars were released > 20 years ago, the underlying mechanism and the genes controlling the resistance remain unknown. We used genetic mapping to elucidate FOB resistance in an F₂ population derived from a cross between resistant and susceptible cultivars. We performed genotyping by sequencing of 173 offspring and aligning the data to the sweet basil reference genome. In total, 23,411 polymorphic sites were detected, and a single quantitative trait locus (QTL) for FOB resistance was found. The confidence interval was < 600 kbp, harboring only 60 genes, including a cluster of putative disease-resistance genes. Based on homology to a fusarium resistance protein from wild tomato, we also investigated a candidate resistance gene that encodes a transmembrane leucine-rich repeat - receptor-like kinase - ubiquitin-like protease (LRR-RLK-ULP). Sequence analysis of that gene in the susceptible parent vs. the resistant parent revealed multiple indels, including an insertion of 20 amino acids next to the transmembrane domain, which might alter its functionality. Our findings suggest that this LRR-RLK-ULP might be responsible for FOB resistance in sweet basil and demonstrate the usefulness of the recently sequenced basil genome for QTL mapping and gene mining.

Amino acids metabolism as a source for aroma volatiles biosynthesis

Aroma volatiles are essential for plant ecological fitness and reproduction. Plants produce and use volatiles to attract pollinators and seed dispersers, repel herbivores and recruit their natural enemies, and communicate with other plants. Amino acids and their biosynthetic intermediates play key roles as precursors for the biosynthesis of plant volatiles. Different plants utilize different strategies and biosynthetic pathways to meet their specific biological needs. This review focuses on the different biosynthetic pathways that plants utilize to form amino acid-derived aroma volatiles, emphasizing their common and unique aspects and stressing the importance of the limiting enzymes residing in the primary-specialized metabolism interface. We also briefly review how biotechnology has used this interface and point to promising future directions for improving the quality of agricultural produce and the production of key volatiles.

Identification of Chilling Accumulation-Associated Genes for Litchi Flowering by Transcriptome-Based Genome-Wide Association Studies

Litchi is an important Sapindaceae fruit tree. Flowering in litchi is triggered by low temperatures in autumn and winter. It can be divided into early-, medium-, and late-flowering phenotypes according to the time for floral induction. Early-flowering varieties need low chilling accumulation level for floral induction, whereas the late-flowering varieties require high chilling accumulation level. In the present study, RNA-Seq of 87 accessions was performed and transcriptome-based genome-wide association studies (GWAS) was used to identify candidate genes involved in chilling accumulation underlying the time for floral induction. A total of 98,155 high-quality single-nucleotide polymorphism (SNP) sites were obtained. A total of 1,411 significantly associated SNPs and 1,115 associated genes (AGs) were identified, of which 31 were flowering-related, 23 were hormone synthesis-related, and 27 were hormone signal transduction-related. Association analysis between the gene expression of the AGs and the flowering phenotypic data was carried out, and differentially expressed genes (DEGs) in a temperature-controlled experiment were obtained. As a result, 15 flowering-related candidate AGs (CAGs), 13 hormone synthesis-related CAGs, and 11 hormone signal transduction-related CAGs were further screened. The expression levels of the CAGs in the early-flowering accessions were different from those in the late-flowering ones, and also between the flowering trees and non-flowering trees. In a gradient chilling treatment, flowering rates of the trees and the CAGs expression were affected by the treatment. Our present work for the first time provided candidate genes for genetic regulation of flowering in litchi using transcriptome-based GWAS.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

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

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.

High-density genetic linkage mapping reveals low stability of QTLs across environments for economic traits in Eucalyptus

Introduction

Eucalyptus urophylla, E. tereticornis and their hybrids are the most important commercial forest tree species in South China where they are grown for pulpwood and solid wood production. Construction of a fine-scale genetic linkage map and detecting quantitative trait loci (QTL) for economically important traits linked to these end-uses will facilitate identification of the main candidate genes and elucidate the regulatory mechanisms.

Method

A high-density consensus map (a total of 2754 SNPs with 1359.18 cM) was constructed using genotyping by sequencing (GBS) on clonal progenies of E. urophylla × tereticornis hybrids. QTL mapping of growth and wood property traits were conducted in three common garden experiments, resulting in a total of 108 QTLs. A total of 1052 candidate genes were screened by the efficient combination of QTL mapping and transcriptome analysis.

Results

Only ten QTLs were found to be stable across two environments, and only one (qSG10Stable mapped on chromosome 10, and associated with lignin syringyl-to-guaiacyl ratio) was stable across all three environments. Compared to other QTLs, qSG10Stable explained a very high level of phenotypic variation (18.4-23.6%), perhaps suggesting that QTLs with strong effects may be more stably inherited across multiple environments. Screened candidate genes were associated with some transcription factor families, such as TALE, which play an important role in the secondary growth of plant cell walls and the regulation of wood formation.

Discussion

While QTLs such as qSG10Stable, found to be stable across three sites, appear to be comparatively uncommon, their identification is likely to be a key to practical QTL-based breeding. Further research involving clonally-replicated populations, deployed across multiple target planting sites, will be required to further elucidate QTL-by-environment interactions.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Research Progress of PPR Proteins in RNA Editing, Stress Response, Plant Growth and Development

RNA editing is a posttranscriptional phenomenon that includes gene processing and modification at specific nucleotide sites. RNA editing mainly occurs in the genomes of mitochondria and chloroplasts in higher plants. In recent years, pentatricopeptide repeat (PPR) proteins, which may act as trans-acting factors of RNA editing have been identified, and the study of PPR proteins has become a research focus in molecular biology. The molecular functions of these proteins and their physiological roles throughout plant growth and development are widely studied. In this minireview, we summarize the current knowledge of the PPR family, hoping to provide some theoretical reference for future research and applications.

A novel introgression line collection to unravel the genetics of climacteric ripening and fruit quality in melon

Introgression lines are valuable germplasm for scientists and breeders, since they ease genetic studies such as QTL interactions and positional cloning as well as the introduction of favorable alleles into elite varieties. We developed a novel introgression line collection in melon using two commercial European varieties with different ripening behavior, the climacteric cantalupensis 'Védrantais' as recurrent parent and the non-climacteric inodorus 'Piel de Sapo' as donor parent. The collection contains 34 introgression lines, covering 99% of the donor genome. The mean introgression size is 18.16 Mb and ~ 3 lines were obtained per chromosome, on average. The high segregation of these lines for multiple fruit quality traits allowed us to identify 27 QTLs that modified sugar content, altered fruit morphology or were involved in climacteric ripening. In addition, we confirmed the genomic location of five major genes previously described, which control mainly fruit appearance, such as mottled rind and external color. Most of the QTLs had been reported before in other populations sharing parental lines, while three QTLs (EAROQP11.3, ECDQP11.2 and FIRQP4.1) were newly detected in our work. These introgression lines would be useful to perform additional genetic studies, as fine mapping and gene pyramiding, especially for important complex traits such as fruit weight and climacteric ripening.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

Genetic dissection of aroma biosynthesis in melon and its relationship with climacteric ripening

Aroma is an essential trait in melon fruit quality, but its complexity and genetic basis are still poorly understood. The aim of this study was the identification of quantitative trait loci (QTLs) underlying volatile organic compounds (VOCs) biosynthesis in melon rind and flesh, using a Recombinant Inbred Line (RIL) population from the cross 'Piel de Sapo' (PS) × 'Védrantais' (VED), two commercial varieties segregating for ripening behavior. A total of 82 VOCs were detected by gas chromatography-mass spectrometry (GC-MS), and 166 QTLs were identified. The main QTL cluster was on chromosome 8, collocating with the previously described ripening-related QTL ETHQV8.1, with an important role in VOCs biosynthesis. QTL clusters involved in esters, lipid-derived volatiles and apocarotenoids were also identified, and candidate genes have been proposed for ethyl 3-(methylthio)propanoate and benzaldehyde biosynthesis. Our results provide genetic insights for deciphering fruit aroma in melon and offer new tools for flavor breeding.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Transition of aromatic volatile and transcriptome profiles during melon fruit ripening

Melon (Cucumis melo L.) is an important diploid crop with a wide variety of flavors due to its distinct aromatic volatile organic compounds (VOC). To understand the development of VOC profiles during fruit development, we performed metabolomic and transcriptomic analysis of two cantaloupe varieties over the course of fruit development. A total of 130 metabolites were detected in fruit samples, and 449014207 reads were mapped to the melon genome. A total of 4469 differentially expressed genes in fruits were identified and used to visualize the transition of VOC and transcriptomic profiles during the fruit development. A shift of VOC profiles in both varieties was observed from early-fruit profiles enriched in C5-C8 lipid-derived VOCs to late-fruit profiles abundant in C9 lipid-derived VOCs, apocarotenoids, and esters. The shift coincided with the expression of specific isoforms of lipid and carotenoid metabolizing enzymes as well as transcription factors involved in fruit ripening, metabolite regulation, and hormone signaling.

Transcriptomic analysis of a near-isogenic line of melon with high fruit flesh firmness during ripening

BACKGROUND

A near-isogenic line (NIL) of melon (SC10-2) with introgression in linkage group X was studied from harvest (at firm-ripe stage of maturity) until day 18 of postharvest storage at 20.5 °C together with its parental control ('Piel de Sapo', PS).

RESULTS

SC10-2 showed higher flesh firmness and whole fruit hardness but lower juiciness than its parental. SC10-2 showed a decrease in respiration rate accompanied by a decrease in ethylene production during ripening, both of which fell to a greater extent than in PS. The introgression affected 11 volatile organic compounds (VOCs), the levels of which during ripening were generally higher in SC10-2 than in PS. Transcriptomic analysis from RNA-Seq revealed differentially expressed genes (DEGs) associated with the effects studied. For example, 909 DEGs were exclusive to the introgression, and only 23 DEGs were exclusive to postharvest ripening time. Major functions of the DEGs associated with introgression or ripening time were identified by cluster analysis. About 37 genes directly and/or indirectly affected the delay in ripening of SC10-2 compared with PS in general and, more particularly, the physiological and quality traits measured and, probably, the differential non-climacteric response. Of the former genes, we studied in more detail at least five that mapped in the introgression in linkage group (LG) X, and 32 outside it.

CONCLUSION

There is an apparent control of textural changes, VOCs and fruit ripening by an expression quantitative trait locus located in LG X together with a direct control on them due to genes presented in the introgression (CmTrpD, CmNADH1, CmTCP15, CmGDSL esterase/lipase, and CmHK4-like) and CmNAC18. © 2020 Society of Chemical Industry.

High-Resolution Transcriptome Atlas and Improved Genome Assembly of Common Buckwheat, Fagopyrum esculentum

Common buckwheat (Fagopyrum esculentum) is an important non-cereal grain crop and a prospective component of functional food. Despite this, the genomic resources for this species and for the whole family Polygonaceae, to which it belongs, are scarce. Here, we report the assembly of the buckwheat genome using long-read technology and a high-resolution expression atlas including 46 organs and developmental stages. We found that the buckwheat genome has an extremely high content of transposable elements, including several classes of recently (0.5-1 Mya) multiplied TEs ("transposon burst") and gradually accumulated TEs. The difference in TE content is a major factor contributing to the three-fold increase in the genome size of F. esculentum compared with its sister species F. tataricum. Moreover, we detected the differences in TE content between the wild ancestral subspecies F. esculentum ssp. ancestrale and buckwheat cultivars, suggesting that TE activity accompanied buckwheat domestication. Expression profiling allowed us to test a hypothesis about the genetic control of petaloidy of tepals in buckwheat. We showed that it is not mediated by B-class gene activity, in contrast to the prediction from the ABC model. Based on a survey of expression profiles and phylogenetic analysis, we identified the MYB family transcription factor gene tr_18111 as a potential candidate for the determination of conical cells in buckwheat petaloid tepals. The information on expression patterns has been integrated into the publicly available database TraVA: http://travadb.org/browse/Species=Fesc/. The improved genome assembly and transcriptomic resources will enable research on buckwheat, including practical applications.

Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon

Domestication and improvement are two important stages in crop evolution. Melon (Cucumis melo L.) is an important vegetable crop with wide phenotypic diversity in many horticultural traits, especially fruit size, flesh thickness and aroma, which are likely the results of long-term extensive selection during its evolution. However, selective signals in domestication and improvement stages for these remarkable variations remain unclear. We resequenced 297 wild, landrace and improved melon accessions and obtained 2 045 412 high-quality SNPs. Population structure and genetic diversity analyses revealed independent and two-step selections in two subspecies of melon: ssp. melo and ssp. agrestis during melon breeding. We detected 233 (~18.35 Mbp) and 159 (~17.71 Mbp) novel potential selective signals during the improvement stage in ssp. agrestis and spp. melo, respectively. Two alcohol acyltransferase genes (CmAATs) unique to the melon genome compared with other cucurbit crops may have undergone stronger selection in ssp. agrestis for the characteristic aroma as compared with other cucurbits. Genome-wide association analysis identified eight fruit size and seven flesh thickness signals overlapping with selective sweeps. Compared with thin-skinned ssp. agrestis, thick-skinned ssp. melo has undergone a stronger selection for thicker flesh. In most melon accessions, CmCLV3 has pleiotropic effects on carpel number and fruit shape. Findings from this study provide novel insights into melon crop evolution, and new tools to advance melon breeding.

Melon Genetic Resources Characterization for Rind Volatile Profile

A melon core collection was analyzed for rind volatile compounds as, despite the fact that they are scarcely studied, these compounds play an important role in consumer preferences. Gas chromatography coupled to mass spectrometry allowed the detection of 171 volatiles. The high volatile diversity found was analyzed by Hierarchical Cluster Analysis (HCA), giving rise to two major clusters of accessions. The first cluster included climacteric and aromatic types such as Cantalupensis, Ameri, Dudaim and Momordica, rich in esters; the second one mainly included non-climacteric non-aromatic types such as Inodorus, Flexuosus, Acidulus, Conomon and wild Agrestis, with low volatiles content, specifically affecting esters. Many interesting accessions were identified, with different combinations of aroma profiles for rind and flesh, such as Spanish Inodorus landraces with low aroma flesh but rind levels of esters similar to those in climacteric Cantalupensis, exotic accessions sharing high contents of specific compounds responsible for the unique aroma of Dudaim melons or wild Agrestis with unexpected high content of some esters. Sesquiterpenes were present in rinds of some Asian Ameri and Momordica landraces, and discriminate groups of cultivars (sesquiterpene-rich/-poor) within each of the two most commercial melon horticultural groups (Cantalupensis and Inodorus), suggesting that the Asian germplasm is in the origin of specific current varieties or that this feature has been introgressed more recently from Asian sources. This rind characterization will encourage future efforts for breeding melon quality as many of the characterized landraces and wild accessions have been underexploited.

A Genome-Wide Analysis of the Pentatricopeptide Repeat (PPR) Gene Family and PPR-Derived Markers for Flesh Color in Watermelon (Citrullus lanatus)

Watermelon (Citrullus lanatus) is an economically important fruit crop grown for consumption of its large edible fruit flesh. Pentatricopeptide-repeat (PPR) encoding genes, one of the large gene families in plants, are important RNA-binding proteins involved in the regulation of plant growth and development by influencing the expression of organellar mRNA transcripts. However, systematic information regarding the PPR gene family in watermelon remains largely unknown. In this comprehensive study, we identified and characterized a total of 422 C. lanatus PPR (ClaPPR) genes in the watermelon genome. Most ClaPPRs were intronless and were mapped across 12 chromosomes. Phylogenetic analysis showed that ClaPPR proteins could be divided into P and PLS subfamilies. Gene duplication analysis suggested that 11 pairs of segmentally duplicated genes existed. In-silico expression pattern analysis demonstrated that ClaPPRs may participate in the regulation of fruit development and ripening processes. Genotyping of 70 lines using 4 single nucleotide polymorphisms (SNPs) from 4 ClaPPRs resulted in match rates of over 0.87 for each validated SNPs in correlation with the unique phenotypes of flesh color, and could be used in differentiating red, yellow, or orange watermelons in breeding programs. Our results provide significant insights for a comprehensive understanding of PPR genes and recommend further studies on their roles in watermelon fruit growth and ripening, which could be utilized for cultivar development of watermelon.

iTRAQ and RNA-Seq analyses revealed the effects of grafting on fruit development and ripening of oriental melon (Cucumis melo L. var. makuwa)

Rootstocks are among the primary factors that influence fruit yield and quality as well as melon development. To understand the differences in the molecular mechanisms and gene expression networks of fruit development between grafted and nongrafted plants in oriental melon, we performed a comprehensive analysis of the transcriptome and proteome dynamic gene/protein expression profiles during fruit development in oriental melon (Cucumis melo L. var. makuwa). Using pairwise comparisons between grafted and nongrafted samples by transcriptome analysis, we identified a large number of candidate genes involved in hormonal signaling pathways, transcription factors, resistance-related biosynthetic pathways and photosynthesis-related metabolic pathways. Many transcription factor-encoded genes were significantly more strongly expressed in the grafted samples, for example, AP2/ERF, C2H2, MYB, bHLH, and AUX/IAA, which are well-known participants in the regulation of developmental processes and hormonal signaling metabolism. Some differentially expressed genes (DEGs) were enriched in flavonoid biosynthesis and phenylpropanoid biosynthesis and determined plant resistance. In addition, some differentially expressed proteins (DEPs) were enriched in photosynthesis-related pathways, which could improve fruit quality and yield. Moreover, through weighted gene coexpression network analyses, we identified modules of coexpressed genes and hub genes specifically related to grafting for different fruit developmental stages. The results suggested that graft-related modules and hub genes were primarily associated with photosynthate metabolism and hormonal signaling pathways. The results obtained in this study provide a valuable resource for dissecting the role of candidate genes governing graft-related metabolism in oriental melon fruit, suggesting an interesting correlation with the effects of rootstock on fruit development.

The Chromosome-Scale Genome of Melon Dissects Genetic Architecture of Important Agronomic Traits

Comparative and evolutionary genomics analyses are the powerful tools to provide mechanistic insights into important agronomic traits. Here, we completed a chromosome-scale assembly of the "neglected" but vital melon subspecies Cucumis melo ssp. agrestis using single-molecule real-time sequencing, Hi-C, and an ultra-dense genetic map. Comparative genomics analyses identified two targeted genes, UDP-sugar pyrophosphorylase and α-galactosidase, that were selected during evolution for specific phloem transport of oligosaccharides in Cucurbitaceae. Association analysis of transcriptome and the DNA methylation patterns revealed the epigenetic regulation of sucrose accumulation in developing fruits. We constructed the melon recombinant inbred lines to uncover Alkaline/Neutral Invertase (CINV), Sucrose-Phosphatase 2 (SPP2), α-galactosidase, and β-galactosidase loci related to sucrose accumulation and an LRR receptor-like serine/threonine-protein kinase associated with gummy stem blight resistance. This study provides essential genomic resources enabling functional genomics studies and the genomics-informed breeding pipelines for improving the fruit quality and disease resistance traits.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

Using networks to identify and interpret natural variation

Studies on natural variation and network biology inherently work to summarize vast amounts of information and data. The combination of these two areas of study while creating datasets of immense complexity is critical to their mutual progress. Networks are necessary as a way to work to reduce the dimensionality inherent in natural variation with 100 s to 1000 s of genotypes. Correspondingly natural variation is essential for testing how networks may or may not be shared across individuals or species. Advances in this area of cross-fertilization including using networks directly as phenotypes and the use of networks to help in prioritizing candidate gene validation efforts. Interesting new observations on frequent presence-absence variation in gene content and adaptation is beginning to highlight the potential for natural variation in network presence-absence. This review attempts to delve into these new insights.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

Comparative Metabolomics and Molecular Phylogenetics of Melon (Cucumis melo, Cucurbitaceae) Biodiversity

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo, we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

A roadmap for research in octoploid strawberry

The cultivated strawberry (Fragaria × ananassa) is an allo-octoploid species, originating nearly 300 years ago from wild progenitors from the Americas. Since that time the strawberry has become the most widely cultivated fruit crop in the world, universally appealing due to its sensory qualities and health benefits. The recent publication of the first high-quality chromosome-scale octoploid strawberry genome (cv. Camarosa) is enabling rapid advances in genetics, stimulating scientific debate and provoking new research questions. In this forward-looking review we propose avenues of research toward new biological insights and applications to agriculture. Among these are the origins of the genome, characterization of genetic variants, and big data approaches to breeding. Key areas of research in molecular biology will include the control of flowering, fruit development, fruit quality, and plant-pathogen interactions. In order to realize this potential as a global community, investments in genome resources must be continually augmented.

The Genetics of Differential Gene Expression Related to Fruit Traits in Strawberry (Fragaria ×ananassa)

Octoploid strawberry (Fragaria ×ananassa) is a major specialty crop under intense annual selection for traits relating to plant vigor and fruit quality. Most functional validation experiments rely on transgenic or transient gene expression assays in the mature receptacle. These findings are not typically translatable to breeding without identifying a natural genetic source of transcript level variation, and developing reliable markers for selection in octoploids. Expression QTL (eQTL) analysis is a genetic/transcriptomic association approach for identifying sequence variants predicting differential expression. This eQTL study analyzed a wide array of mature receptacle-expressed genes, encompassing the majority of total mature receptacle transcript accumulation and almost all strawberry genes described in the literature. These results identified segregating genetic variants associated with the differential expression of hundreds of strawberry genes, many with known interest to breeders. Several of these eQTL pertain to published genes whose expression levels have been demonstrated to influence mature receptacle phenotypes. Many include key genes of the phenylpropanoid pathway, vitamin C, carotenoid, pectin, and receptacle carbohydrate/sugar metabolism. These subgenome-specific genetic markers may allow breeders to select for desired ranges of target gene expression. These results may also guide basic research efforts and facilitate the identification of causal genes underlying trait QTL.