Genome-Wide Single Nucleotide Polymorphism Discovery and the Construction of a High-Density Genetic Map for Melon (Cucumis melo L.) Using Genotyping-by-Sequencing

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

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.

Genetic analysis and QTL mapping of yield and fruit traits in bitter gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an economically important vegetable crop grown in tropical parts of the world. In this study, a high-density linkage map of M. charantia was constructed through genotyping-by-sequencing (GBS) technology using F2:3 mapping population generated from the cross DBGy-201 × Pusa Do Mausami. About 2013 high-quality SNPs were assigned on a total of 20 linkage groups (LGs) spanning over 2329.2 CM with an average genetic distance of 1.16 CM. QTL analysis was performed for six major yield-contributing traits such as fruit length, fruit diameter, fruit weight, fruit flesh thickness, number of fruits per plant and yield per plant. These six quantitative traits were mapped with 19 QTLs (9 QTLs with LOD > 3) using composite interval mapping (CIM). Among 19 QTLs, 12 QTLs derived from 'Pusa Do Mausami' revealed a negative additive effect when its allele increased trait score whereas 7 QTLs derived from 'DBGy-201' revealed a positive additive effect when its allele trait score increased. The phenotypic variation (R2%) elucidated by these QTLs ranged from 0.09% (fruit flesh thickness) on LG 14 to 32.65% (fruit diameter) on LG 16 and a total of six major QTLs detected. Most QTLs detected in the present study were located relatively very close, maybe due to the high correlation among the traits. This information will serve as a significant basis for marker-assisted selection and molecular breeding in bitter gourd crop improvement.

Quantitative Trait Loci Mapping for Earliness, Fruit, and Seed Related Traits Using High Density Genotyping-by-Sequencing-Based Genetic Map in Bitter Gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an important vegetable crop having numerous medicinal properties. Earliness and yield related traits are main aims of bitter gourd breeding program. High resolution quantitative trait loci (QTLs) mapping can help in understanding the molecular basis of phenotypic variation of these traits and thus facilitate marker-assisted breeding. The aim of present study was to identify genetic loci controlling earliness, fruit, and seed related traits. To achieve this, genotyping-by-sequencing (GBS) approach was used to genotype 101 individuals of F₄ population derived from a cross between an elite cultivar Punjab-14 and PAUBG-6. This population was phenotyped under net-house conditions for three years 2018, 2019, and 2021. The linkage map consisting of 15 linkage groups comprising 3,144 single nucleotide polymorphism (SNP) markers was used to detect the QTLs for nine traits. A total of 50 QTLs for these traits were detected which were distributed on 11 chromosomes. The QTLs explained 5.09-29.82% of the phenotypic variance. The highest logarithm of the odds (LOD) score for a single QTL was 8.68 and the lowest was 2.50. For the earliness related traits, a total of 22 QTLs were detected. For the fruit related traits, a total of 16 QTLs and for seed related traits, a total of 12 QTLs were detected. Out of 50 QTLs, 20 QTLs were considered as frequent QTLs (FQ-QTLs). The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for these traits in bitter gourd improvement program.

High-Density GBS-Based Genetic Linkage Map Construction and QTL Identification Associated With Yellow Mosaic Disease Resistance in Bitter Gourd (Momordica charantia L.)

Yellow mosaic disease (YMD) in bitter gourd (Momordica charantia) is a devastating disease that seriously affects its yield. Although there is currently no effective method to control the disease, breeding of resistant varieties is the most effective and economic option. Moreover, quantitative trait locus (QTL) associated with resistance to YMD has not yet been reported. With the objective of mapping YMD resistance in bitter gourd, the susceptible parent "Punjab-14" and the resistant parent "PAUBG-6" were crossed to obtain F₄ mapping population comprising 101 individuals. In the present study, the genotyping by sequencing (GBS) approach was used to develop the genetic linkage map. The map contained 3,144 single nucleotide polymorphism (SNP) markers, consisted of 15 linkage groups, and it spanned 2415.2 cM with an average marker distance of 0.7 cM. By adopting the artificial and field inoculation techniques, F_4:5 individuals were phenotyped for disease resistance in Nethouse (2019), Rainy (2019), and Spring season (2020). The QTL analysis using the genetic map and phenotyping data identified three QTLs qYMD.pau_3.1, qYMD.pau_4.1, and qYMD.pau_5.1 on chromosome 3, 4, and 5 respectively. Among these, qYMD.pau_3.1, qYMD.pau_4.1 QTLs were identified during the rainy season, explaining the 13.5 and 21.6% phenotypic variance respectively, whereas, during the spring season, qYMD.pau_4.1 and qYMD.pau_5.1 QTLs were observed with 17.5 and 22.1% phenotypic variance respectively. Only one QTL qYMD.pau_5.1 was identified for disease resistance under nethouse conditions with 15.6% phenotypic variance. To our knowledge, this is the first report on the identification of QTLs associated with YMD resistance in bitter gourd using SNP markers. The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for disease resistance.

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.

Genome-wide characterization of two-component system (TCS) genes in melon (Cucumis melo L.)

To better understand cytokinin signaling in melon (Cucumis melo L.), one of the most important fruit crops in the Cucurbitaceae family, we identified and characterized melon two-component system (TCS) genes in this study. The results showed that there were 51 genes encoding putative TCS proteins in melon, and these TCS genes were classified into 3 subgroups, with 17 HK(L)s (histidine kinase/histidine-kinase like; 9 HKs and 8 HKLs), 9 HPs (histidine phosphotransfer proteins; 6 authentic and 3 pseudo), and 25 RRs (response regulators; 8 Type-A, 11 Type-B and 6 pseudo). The identity values of these cytokinin signaling proteins were revealed by analyzing their conserved motifs, domains and amino acid sequences. By analyzing TCS genes in different plant species, we found that melon HK(L)s, HPs and RRs had closer phylogenetic relationships with cucumber genes than with the genes of other plants, and the expansion of melon cytokinin signaling genes might be attributed to segmental duplication events. Analysis of the putative promoter regions (2-kb upstream regions of the start codon) revealed the enrichment of stress- and hormone-response cis-elements. The involvement of these putative TCS genes in melon cytokinin signaling was further supported by qRT-PCR data.

Genetic Analysis and Fine Mapping of a Spontaneously Mutated Male Sterility Gene in Brassica rapa ssp. chinensis

Male sterility has been widely used in hybrid seed production in Brassica, but not in B. rapa ssp. chinensis, and genetic models of male sterility for this subspecies are unclear. We discovered a spontaneous mutant in B. rapa ssp. chinensis A series of progeny tests indicated that male sterility in B. rapa ssp. chinensis follows a three-allele model with BrMsa , BrMsb , and BrMsc The male sterility locus has been mapped to chromosome A07 in BC1 and F2 populations through genotyping by sequencing. Fine mapping in a total of 1,590 F2 plants narrowed the male sterility gene BrMs to a 400 kb region, with two SNP markers only 0.3 cM from the gene. Comparative gene mapping shows that the Ms gene in B. rapa ssp. pekinensis is different from the BrMs gene of B. rapa ssp. chinensis, despite that both genes are located on chromosome A07. Interestingly, the DNA sequence orthologous to a male sterile gene in Brassica napus, BnRf, is within 400 kb of the BrMs locus. The BnRf orthologs of B. rapa ssp. chinensis were sequenced, and one KASP marker (BrMs_indel) was developed for genotyping based on a 14 bp indel at intron 4. Cosegregation of male sterility and BrMs_indel genotypes in the F2 population indicated that BnRf from B. napus and BrMs from B. rapa are likely to be orthologs. The BrMs_indel marker developed in this study will be useful in marker-assisted selection for the male sterility trait.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

Genomic and phenotypic evaluation of rice susceptible check TN1 collected in Taiwan

BACKGROUND

Taichung Native 1 (TN1), a variety of rice (Oryza sativa L.) developed in Taiwan, has played a key role in the green revolution of this major staple crop because of its semi-dwarf characteristics. Due to its susceptibility, it has been used as a susceptibility indicator in rice insect and pathogen resistance studies worldwide. While within-variety differences have been reported for agronomic traits in other rice varieties, no study has addressed the within-variety consistency of pathogen and insect susceptibility of TN1, which would influence the result interpretation of plant-pest interaction studies. Therefore, the objective of this study was to evaluate the genomic consistency and to assess a range of agronomic and insect susceptibility traits in three representative accessions of TN1 in Taiwan.

RESULTS

Among these three accessions, two were identical across 43,325 genome-wide single nucleotide polymorphisms (SNPs) while the third one differed at four SNPs. Of the three accessions of TN1, there were minor differences in seed length, seed breadth, length/width ratio, number of leaves and tillers, and number of unfilled seeds. Besides, there was no effect on relative growth rate of Cnaphalocrocis medinalis larvae fed on the three accession sources. Furthermore, there is no different on plant susceptibility among these three accessions against C. medinalis and Nilaparvata lugens.

CONCLUSION

Our study indicates that it is appropriate to use TN1 in Taiwan to test for rice insect susceptibility as it yields consistent results.

QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map

BACKGROUND

Melon shows a broad diversity in fruit morphology and quality, which is still underexploited in breeding programs. The knowledge of the genetic basis of fruit quality traits is important for identifying new alleles that may be introduced in elite material by highly efficient molecular breeding tools.

RESULTS

In order to identify QTLs controlling fruit quality, a recombinant inbred line population was developed using two commercial cultivars as parental lines: "Védrantais", from the cantalupensis group, and "Piel de Sapo", from the inodorus group. Both have desirable quality traits for the market, but their fruits differ in traits such as rind and flesh color, sugar content, ripening behavior, size and shape. We used a genotyping-by-sequencing strategy to construct a dense genetic map, which included around five thousand variants distributed in 824 bins. The RIL population was phenotyped for quality and morphology traits, and we mapped 33 stable QTLs involved in sugar and carotenoid content, fruit and seed morphology and major loci controlling external color of immature fruit and mottled rind. The median confidence interval of the QTLs was 942 kb, suggesting that the high density of the genetic map helped in increasing the mapping resolution. Some of these intervals contained less than a hundred annotated genes, and an integrative strategy combining gene expression and resequencing data enabled identification of candidate genes for some of these traits.

CONCLUSION

Several QTLs controlling fruit quality traits in melon were identified and delimited to narrow genomic intervals, using a RIL population and a GBS-based genetic map.

Construction of a genome-anchored, high-density genetic map for melon (Cucumis melo L.) and identification of Fusarium oxysporum f. sp. melonis race 1 resistance QTL

Four QTLs and an epistatic interaction were associated with disease severity in response to inoculation with Fusarium oxysporum f. sp. melonis race 1 in a recombinant inbred line population of melon. The USDA Cucumis melo inbred line, MR-1, harbors a wealth of alleles associated with resistance to several major diseases of melon, including powdery mildew, downy mildew, Alternaria leaf blight, and Fusarium wilt. MR-1 was crossed to an Israeli cultivar, Ananas Yok'neam, which is susceptible to all of these diseases, to generate a recombinant inbred line (RIL) population of 172 lines. In this study, the RIL population was genotyped to construct an ultra-dense genetic linkage map with 5663 binned SNPs anchored to the C. melo genome and exhibits the overall high quality of the assembly. The utility of the densely genotyped population was demonstrated through QTL mapping of a well-studied trait, resistance to Fusarium wilt caused by Fusarium oxysporum f. sp. melonis (Fom) race 1. A major QTL co-located with the previously validated resistance gene Fom-2. In addition, three minor QTLs and an epistatic interaction contributing to Fom race 1 resistance were identified. The MR-1 × AY RIL population provides a valuable resource for future QTL mapping studies and marker-assisted selection of disease resistance in melon.

A high-density linkage map and QTL mapping of fruit-related traits in pumpkin (Cucurbita moschata Duch.)

Pumpkin (Cucurbita moschata) is an economically worldwide crop. Few quantitative trait loci (QTLs) were reported previously due to the lack of genomic and genetic resources. In this study, a high-density linkage map of C. moschata was structured by double-digest restriction site-associated DNA sequencing, using 200 F2 individuals of CMO-1 × CMO-97. By filtering 74,899 SNPs, a total of 3,470 high quality SNP markers were assigned to the map spanning a total genetic distance of 3087.03 cM on 20 linkage groups (LGs) with an average genetic distance of 0.89 cM. Based on this map, both pericarp color and strip were fined mapped to a novel single locus on LG8 in the same region of 0.31 cM with phenotypic variance explained (PVE) of 93.6% and 90.2%, respectively. QTL analysis was also performed on carotenoids, sugars, tuberculate fruit, fruit diameter, thickness and chamber width with a total of 12 traits. 29 QTLs distributed in 9 LGs were detected with PVE from 9.6% to 28.6%. It was the first high-density linkage SNP map for C. moschata which was proved to be a valuable tool for gene or QTL mapping. This information will serve as significant basis for map-based gene cloning, draft genome assembling and molecular breeding.

A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo

The melon (Cucumis melo) genome and genetic maps with hundreds to thousands of single nucleotide polymorphism markers were recently released. However, a high-resolution genetic map was lacking. Gummy stem blight (Gsb) is a destructive disease responsible for considerable economic losses during melon production. We herein describe the development of an ultra-dense genetic map consisting of 12,932 recombination bin markers covering 1,818 cM, with an average distance of 0.17 cM between adjacent tags. A comparison of the genetic maps for melon, watermelon, and cucumber revealed chromosome-level syntenic relationships and recombination events among the three Cucurbitaceae species. Our genetic map was useful for re-anchoring the genome scaffolds of melon. More than 92% assembly was anchored to 12 pseudo-chromosomes and 90% of them were oriented. Furthermore, 1,135 recombination hotspots revealed an unbalanced recombination rate across the melon genome. Genetic analyses of the Gsb-resistant and -susceptible lines indicated the resistance phenotype is mediated by a single dominant gene. We identified Gsb-resistance gene candidates in a 108-kb region on pseudo-chromosome 4. Our findings verify the utility of an ultra-dense genetic map for mapping a gene of interest, and for identifying new disease resistant genes.