Genetic mapping reveals a marker for yellow skin in watermelon (Citrullus lanatus L.)

Uncovering the genetic mechanism of rind color trait in watermelon using fine mapping and comparative transcriptomic analysis

The rind color of watermelon fruit is a significant trait that directly affects consumer acceptability. However, the genetic regulatory mechanisms underlying rind color remain poorly understood. In this study, we crossed two differentiated watermelon lines (K2Q "female parent line with a light green rind" and K2S "male parent line with a dark green rind") and developed segregated F2 mapping populations. The dynamic development of rind color was observed by identifying the critical period for color transformation as occurring between 7 and 14 days after pollination (DAP). Genetic segregation analysis indicated that a single dominant gene regulates the major genetic locus (ClRC) associated with the dark green rind trait. Whole-genome BSA-sequencing (BSA-seq) and fine mapping analysis exposed the delimited ClRC locus to a 37.52 kb region on chromosome 08 (Chr08), comprising five genes. The pairwise sequence comparisons analysis of the parental lines revealed the single major gene (Cla97C08G161570), which encodes a 2-phytyl-1,4-β-naphthoquinone methyltransferase protein, exhibiting one non-synonymous type single nucleotide polymorphism (nsSNP) at candidate site (Chr8:27994761, C-G). The real-time quantitative polymerase chain reaction (RT-qPCR) verified the higher expression level of the K2S line on the 14 DAP than that of the K2Q line. The analysis of comparative transcriptomes (RNA-sequencing) identified a total of 940 differentially expressed genes (DEGs) associated with rind coloration in the two parental lines at three dynamic stages of development (0, 7, and 14 DAP). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed key genes (C01G023430, C04G071470, C09G165830, C07G128820, C08G148460, and C08G155040) that share the same pathway as the Cla97C08G161570 gene and exhibited high levels of differential expression trend. Further, RT-qPCR verified that these genes display the same expression pattern as the Cla97C08G161570 gene, and expression levels in the dark green rind lines were significantly higher than those in the light green rind lines, suggesting the significant role in modulating the pigmentation activity.

Identification and characterization of ClAPRR2, a key candidate gene controlling watermelon stripe color

The stripe color of watermelon is a vital commercial trait and is the focus of attention of consumers and researchers. However, the genetic determinants of watermelon stripe color are incompletely understood. Based on the results of preliminary localization studies, we constructed a large-capacity F2 generation population (710 plants) using light-green striped ZXG1555 and green-striped Cream of Saskatchewan (COS) watermelon strains as parental lines for fine mapping. Genes controlling stripe color were located in an 85.284 kb region on chromosome 9, which contained five candidate genes. Combined with parental phenotypes, chlorophyll contents of rinds and stripes were assayed. Gene sequence alignment and transcriptional level analysis of parental lines predicted Cla97C09G175170 (encoding a two-component response regulator-like protein, APRR2) as the best candidate gene for stripe color trait. Two SNPs in the ClAPRR2 coding region caused amino acid substitutions, but were not located in the conserved domain, while a 12 bp insertion caused premature translation termination and a 35 amino acid deletion in the conserved domain and may have affected ClAPRR2 function in ZXG1555. Subcellular localization analysis showed that ClAPRR2 was expressed in the ZXG1555 cell membrane but was located in the nucleus and cell membrane of COS. Nucleotide polymorphisms and deletions were also detected in the promoter region between parental lines and caused cis-acting element variations. Luciferase activity suggested that promoter variations may not be the main factor in the regulation of ClAPRR2 expression.

Phytoalexin deficient 4 is associated with the lesion mimic trait in watermelon clalm mutant (Citrullus lanatus)

In watermelon (Citrullus lanatus), lesion mimic is a rare, valuable trait that can be used by breeders for selection at early growth stages. In this study, we tested a seven-generation family to determine the inheritance and genetic basis of this trait. As revealed by analysis of the lesion mimic mutant clalm, this trait is controlled by a single dominant gene. Whole genome resequencing-bulked segregant analysis demonstrated that this gene is located on chromosome 4 from 3,760,000 bp to 7,440,000 bp, a region corresponding to a physical distance of 3.68 Mb encompassing approximately 72 annotated genes. There are 6 genes with non-synonymous mutation SNP sites. The predicted target gene, ClCG04G001930, encodes a Phytoalexin deficient 4 (PAD4), a protein that plays an important regulatory role in leaf senescence in many plant species. According to quantitative real-time PCR analysis, the expression level of ClCG04G001930 was significantly higher in the clalm mutant than in normal watermelon. Twenty-five SNPs were identified in the ClCG04G001930 gene of F2 individuals of the clalm mutant. Overexpression the ClCG04G001930 gene, designated as ClPAD4, yielded transgenic lines whose leaves gradually developed chlorotic lesions over 3 weeks. RNA interference of the ClPAD4 yielded transgenic lines whose cotyledon prone to diseased over 2 weeks. Our results suggest that ClPAD4 might be the candidate gene responsible for lesion mimic in the clalm mutant. Our findings may serve as a foundation for elucidating the mechanism underlying the molecular metabolism of programmed cell death and should be useful for marker-assisted selection breeding in watermelon.

Identification of Candidate Genes for Green Rind Color in Watermelon

The color of the rind is one of the most crucial agronomic characteristics of watermelon (Citrullus lanatus L.). Its genetic analysis was conducted to provide the identification of genes regulating rind color and improving the quality of watermelon appearance. In this study, a mapping population of 505 F2 plants, derived from a cross between green (CG058) and light-green (CG265) rinds, along with a high-density genetic linkage (average 0.9 cM distance between bin markers), was used to map and identify possible candidate genes. The green rind trait was determined to be regulated by a single Mendelian locus and was precisely located within a 110 kb genomic site on chromosome nine (Chr 9). In the respective region, two potential genes, Cla97C09G175170 and Cla97C09G175180, were substantially downregulated in the light-green rind in comparison to the green rind. Previous studies revealed that Cla97C09G175170, encoding a two-component response regulator-like protein (APRR2), is possibly involved in the green rind trait in watermelon. Virus-induced gene silencing (VIGS) assay confirmed that ClAPRR2 is a key gene responsible for green rind color. Moreover, qRT-PCR analysis revealed that the transcription levels of multiple key genes in the chlorophyll (Chl) biosynthesis pathway were downregulated in the light-green rind relative to the green rind. The current findings have the potential to clarify the regulatory mechanisms that underlie the color of the watermelon rind. These data would provide valuable insights for the targeted molecular design and development of watermelon rinds.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

High-quality genome assembly and genetic mapping reveal a gene regulating flesh color in watermelon (Citrullus lanatus)

The unique color and type characteristics of watermelon fruits are regulated by many molecular mechanisms. However, it still needs to be combined with more abundant genetic data to fine-tune the positioning. We assembled genomes of two Korean inbred watermelon lines (cv. 242-1 and 159-1) with unique color and fruit-type characteristics and identified 23,921 and 24,451 protein-coding genes in the two genomes, respectively. To obtain more precise results for further study, we resequenced one individual of each parental line and an F₂ population composed of 87 individuals. This identified 1,539 single-nucleotide polymorphisms (SNPs) and 80 InDel markers that provided a high-density genetic linkage map with a total length of 3,036.9 cM. Quantitative trait locus mapping identified 15 QTLs for watermelon fruit quality-related traits, including β-carotene and lycopene content in fruit flesh, fruit shape index, skin thickness, flesh color, and rind color. By investigating the mapping intervals, we identified 33 candidate genes containing variants in the coding sequence. Among them, Cla97C01G008760 was annotated as a phytoene synthase with a single-nucleotide variant (A → G) in the first exon at 9,539,129 bp of chromosome 1 that resulted in the conversion of a lysine to glutamic acid, indicating that this gene might regulate flesh color changes at the protein level. These findings not only prove the importance of a phytoene synthase gene in pigmentation but also explain an important reason for the color change of watermelon flesh.

Primary mapping of quantitative trait loci regulating multivariate horticultural phenotypes of watermelon (Citrullus lanatus L.)

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F_2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32-25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Genetic mapping and molecular characterization of the delayed green gene dg in watermelon (Citrullus lanatus)

Leaf color mutants are common in higher plants that can be used as markers in crop breeding and are important tools in understanding regulatory mechanisms of chlorophyll biosynthesis and chloroplast development. Genetic analysis was performed by evaluating F₁, F₂ and BC₁ populations derived from two parental lines (Charleston gray with green leaf color and Houlv with delayed green leaf color), suggesting that a single recessive gene controls the delayed green leaf color. In this study, the delayed green mutant showed a conditional pale green leaf color at the early leaf development but turned to green as the leaf development progressed. Delayed green leaf plants showed reduced pigment content, photosynthetic, chlorophyll fluorescence parameters, and impaired chloroplast development compared with green leaf plants. The delayed green (dg) locus was mapped to 7.48 Mb on chromosome 3 through bulk segregant analysis approach, and the gene controlling delayed green leaf color was narrowed to 53.54 kb between SNP130 and SNP135 markers containing three candidate genes. Sequence alignment of the three genes indicated that there was a single SNP mutation (G/A) in the coding region of ClCG03G010030 in the Houlv parent, which causes an amino acid change from Arginine to Lysine. The ClCG03G010030 gene encoded FtsH extracellular protease protein family is involved in early delayed green leaf development. The expression level of ClCG03G010030 was significantly reduced in delayed green leaf plants than in green leaf plants. These results indicated that the ClCG03G010030 might control watermelon green leaf color and the single SNP variation in ClCG03G010030 may result in early delayed green leaf color development during evolutionary process.

Fine mapping and identifying candidate gene of Y underlying yellow peel in Cucurbita pepo

As a conspicuous trait, peel color is one of the most important characteristics that affects commodity quality and consumer preferences. The locus Y underlying yellow peel in Cucurbita pepo (zucchini) was first reported in 1922; however, its molecular mechanism is still unknown. In this study, a genetic analysis revealed that yellow peel is controlled by a single dominant genetic factor. Furthermore, Y was mapped in a ~170 kb region on chromosome 10 by bulked segregated analysis (BSA) and fine mapping in F₂ and BC₁ segregating populations. The candidate region harbors fifteen annotated genes, among which Cp4.1LG10g11560 (CpCHLH) is regarded as a promising candidate gene. CpCHLH encodes a magnesium chelatase H subunit involved in chlorophyll biosynthesis, and its mutation can result in a reduction in chlorophyll content and yellow phenotype. Interestingly, a large fragment (~15 kb) duplication containing incomplete CpCHLH was inserted in the candidate interval, resulting in two reformed CpCHLH proteins in the yellow parental line. It is most likely that the reformed CpCHLH proteins act as a malfunctional competitor of the normal CpCHLH protein to interrupt the formation of chlorophyll. Overall, the isolation of Y will shed light on the molecular mechanism of the peel color regulation of zucchini and lay a foundation for breeding.

Genome-wide association analysis reveals a novel QTL CsPC1 for pericarp color in cucumber

BACKGROUND

Cucumber is an important melon crop in the world, with different pericarp colors. However, the candidate genes and the underlying genetic mechanism for such an important trait in cucumber are unknown. In this study, a locus controlling pericarp color was found on chromosome 3 of cucumber genome.

RESULTS

In this study, the light green inbred line G35 and the dark green inbred line Q51 were crossed to produce one F₂ population. Consequently, we identified a major locus CsPC1 (Pericarp color 1). Next, we mapped the CsPC1 locus to a 94-kb region chromosome 3 which contains 15 genes. Among these genes, Csa3G912920, which encodes a GATA transcription factor, was expressed at a higher level in the pericarp of the NIL-1334 line (with light-green pericarp) than in that of the NIL-1325 line (with dark-green pericarp). This study provides a new allele for the improvement of cucumber pericarp color.

CONCLUSION

A major QTL that controls pericarp color in cucumber, CsPC1, was identified in a 94-kb region that harbors the strong candidate gene CsGATA1.

Differences of rhizospheric and endophytic bacteria are recruited by different watermelon phenotypes relating to rind colors formation

To elucidate the biological mechanism of yellow rind formation on watermelon, the characteristics of soil bacterial community structure in rhizosphere and endophytic bacteria in stem of yellow rind watermelon were analyzed. Based on high-throughput sequencing technology, plant stem and rhizosphere soil samples, which collected from yellow and green rind watermelons were used in this paper, respectively. The structural characteristics of the endophytic bacteria in stems and soil bacterial communities in rhizospheres of yellow and green rind watermelons were comparative studied. Firstly, significant different proportions of some dominant bacteria and abundances could be detected between yellow and rind watermelons. Meanwhile, although different abundances of endophytic bacteria could be found, but no significant differences were observed between yellow and green rind watermelons. Moreover, Gemmatimonadota, Myxococcota, WPS-2, norank_f_Gemmatimonadaceae and Bradyrhizobium were the soil dominant bacterial genera in rhizosphere of green rind watermelon. All above results suggest that differences of rhizospheric and endophytic bacteria are exactly recruited as "workers" by different watermelon phenotypes relating to rind color formations.

QTL mapping identifies novel loci and putative candidate genes for rind traits in watermelon

Rind color and pattern are important external fruit quality traits of watermelon influencing the consumer attention and market acceptance. Varying degrees of green and yellow colors in different combinations forming distinct patterns have been observed on watermelon rind. In the current experiment, an attempt has been made to map the QTLs/genes for four important rind traits, namely, stripe color, stripe pattern, interstripe color, and interstripe pattern. The experiment consisted of two mapping populations namely F₂ (Pop I) and BC₁F₂ (Pop II) derived from two parental lines, viz., BIL-53 (characterized by medium green marbled rind pattern and yellowish white blotchy interstripe pattern) and IIHR-140-152 (characterized by dark green solid stripes and greenish wavy interstripe pattern). Linkage mapping identified consistent QTLs across populations on chromosome 9. Comparative genomic analysis of these regions revealed two genes, namely, Cla97C09G175170 and Cla97C09G175150 as potential candidates for stripe and interstripe color. Sequence analysis of Cla97C09G175170 gene in parents along with reference genotypes, viz., 97,103 and Charleston gray suggests a 3 bp deletion on 11th exon to be associated with stripe color polymorphism in watermelon.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03112-7.

Molecular research progress and improvement approach of fruit quality traits in cucumber

Recent molecular studies revealed new opportunities to improve cucumber fruit quality. However, the fruit color and spine traits molecular basis remain vague despite the vast sources of genetic diversity. Cucumber is agriculturally, economically and nutritionally important vegetable crop. China produces three-fourths of the world's total cucumber production. Cucumber fruit quality depends on a number of traits such as the fruit color (peel and flesh color), spine (density, size and color), fruit shape, fruit size, defects, texture, firmness, taste, maturity stage and nutritional composition. Fruit color and spine traits determine critical quality attributes and have been the interest of researchers at the molecular level. Evaluating the molecular mechanisms of fruit quality traits is important to improve production and quality of cucumber varieties. Genes and qualitative trait locus (QTL) that are responsible for cucumber fruit color and fruit spine have been identified. The purpose of this paper is to reveal the molecular research progress of fruit color and spines as key quality traits of cucumber. The markers and genes identified so far could help for marker-assisted selection of the fruit color and spine trait in cucumber breeding and its associated nutritional improvement. Based on the previous studies, peel color and spine density as examples, we proposed a comprehensive approach for cucumber fruit quality traits improvement. Moreover, the markers and genes can be useful to facilitate cloning-mediated genetic breeding in cucumber. However, in the era of climate change, increased human population and high-quality demand of consumers, studies on molecular mechanisms of cucumber fruit quality traits are limited.

A mutation in the intron splice acceptor site of a GA3ox gene confers dwarf architecture in watermelon (Citrullus lanatus L.)

Dwarf architecture is an important trait associated with plant yield, lodging resistance and labor cost. Here, we aimed to identify a gene causing dwarfism in watermelon. The ‘w106’ (dwarf) and ‘Charleston Gray’ (vine) were used as parents to construct F₁ and F₂ progeny. Dwarf architecture of ‘w106’ was mainly caused by longitudinal cell length reduction and was controlled by a single recessive gene. Whole-genome sequencing of two parents and two bulk DNAs of F₂ population localized this gene to a 2.63-Mb region on chromosome 9; this was further narrowed to a 541-kb region. Within this region, Cla015407, encoding a gibberellin 3β-hydroxylase (GA3ox), was the candidate gene. Cla015407 had a SNP mutation (G → A) in the splice acceptor site of the intron, leading to altered splicing event and generating two splicing isoforms in dwarf plants. One splicing isoform retained the intron sequences, while the other had a 13-bp deletion in the second exon of GA3ox transcript, both resulting in truncated proteins and loss of the functional Fe2OG dioxygenase domain in dwarf plants. RNA-Seq analysis indicated that expression of Cla015407 and other GA biosynthetic and metabolic genes were mostly up-regulated in the shoots of dwarf plants compared with vine plants in F₂ population. Measurement of endogenous GA levels indicated that bioactive GA₄ was significantly decreased in the shoots of dwarf plants. Moreover, the dwarf phenotype can be rescued by exogenous applications of GA₃ or GA₄₊₇, with the latter having a more distinct effect than the former. Subcellular localization analyses of GA3ox proteins from two parents revealed their subcellular targeting in nucleus and cytosol. Here, a GA3ox gene controlling dwarf architecture was identified, and loss function of GA3ox leads to GA₄ reduction and dwarfism phenotype in watermelon.

Mapping a Partial Andromonoecy Locus in Citrullus lanatus Using BSA-Seq and GWAS Approaches

The sexual expression of watermelon plants is the result of the distribution and occurrence of male, female, bisexual and hermaphrodite flowers on the main and secondary stems. Plants can be monoecious (producing male and female flowers), andromonoecious (producing male and hermaphrodite flowers), or partially andromonoecious (producing male, female, bisexual, and hermaphrodite flowers) within the same plant. Sex determination of individual floral buds and the distribution of the different flower types on the plant, are both controlled by ethylene. A single missense mutation in the ethylene biosynthesis gene CitACS4, is able to promote the conversion of female into hermaphrodite flowers, and therefore of monoecy (genotype MM) into partial andromonoecy (genotype Mm) or andromonoecy (genotype mm). We phenotyped and genotyped, for the M/m locus, a panel of 207 C. lanatus accessions, including five inbreds and hybrids, and found several accessions that were repeatedly phenotyped as PA (partially andromonoecious) in several locations and different years, despite being MM. A cosegregation analysis between a SNV in CitACS4 and the PA phenotype, demonstrated that the occurrence of bisexual and hermaphrodite flowers in a PA line is not dependent on CitACS4, but conferred by an unlinked recessive gene which we called pa. Two different approaches were performed to map the pa gene in the genome of C. lanatus: bulk segregant analysis sequencing (BSA-seq) and genome wide association analysis studies (GWAS). The BSA-seq study was performed using two contrasting bulks, the monoecious M-bulk and the partially andromonoecious PA-bulk, each one generated by pooling DNA from 20 F2 plants. For GWAS, 122 accessions from USDA gene bank, already re-sequenced by genotyping by sequencing (GBS), were used. The combination of the two approaches indicates that pa maps onto a genomic region expanding across 32.24-36.44 Mb in chromosome 1 of watermelon. Fine mapping narrowed down the pa locus to a 867 Kb genomic region containing 101 genes. A number of candidate genes were selected, not only for their function in ethylene biosynthesis and signalling as well as their role in flower development and sex determination, but also by the impact of the SNPs and indels differentially detected in the two sequenced bulks.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

Comparative Transcriptome Analysis Provides Insights Into Yellow Rind Formation and Preliminary Mapping of the Clyr (Yellow Rind) Gene in Watermelon

As an important appearance trait, the rind color of watermelon fruit affects the commodity value and further determines consumption choices. In this study, a comparative transcriptome analysis was conducted to elucidate the genes and pathways involved in the formation of yellow rind fruit in watermelon using a yellow rind inbred line WT4 and a green rind inbred line WM102. A total of 2,362 differentially expressed genes (DEGs) between WT4 and WM102 at three different stages (0, 7, and 14 DAP) were identified and 9,770 DEGs were obtained by comparing the expression level at 7 DAP and 14 DAP with the former stages of WT4. The function enrichment of DEGs revealed a number of pathways and terms in biological processes, cellular components, and molecular functions that were related to plant pigment metabolism, suggesting that there may be a group of common core genes regulating rind color formation. In addition, next-generation sequencing aided bulked-segregant analysis (BSA-seq) of the yellow rind pool and green rind pool selected from an F₂ population revealed that the yellow rind gene (Clyr) was mapped on the top end of chromosome 4. Based on the BSA-seq analysis result, Clyr was further confined to a region of 91.42 kb by linkage analysis using 1,106 F₂ plants. These results will aid in identifying the key genes and pathways associated with yellow rind formation and elucidating the molecular mechanism of rind color formation in watermelon.

PyBSASeq: a simple and effective algorithm for bulked segregant analysis with whole-genome sequencing data

BACKGROUND

Bulked segregant analysis (BSA), coupled with next-generation sequencing, allows the rapid identification of both qualitative and quantitative trait loci (QTL), and this technique is referred to as BSA-Seq here. The current SNP index method and G-statistic method for BSA-Seq data analysis require relatively high sequencing coverage to detect significant single nucleotide polymorphism (SNP)-trait associations, which leads to high sequencing cost.

RESULTS

We developed a simple and effective algorithm for BSA-Seq data analysis and implemented it in Python; the program was named PyBSASeq. Using PyBSASeq, the significant SNPs (sSNPs), SNPs likely associated with the trait, were identified via Fisher's exact test, and then the ratio of the sSNPs to total SNPs in a chromosomal interval was used to detect the genomic regions that condition the trait of interest. The results obtained this way are similar to those generated via the current methods, but with more than five times higher sensitivity. This approach was termed the significant SNP method here.

CONCLUSIONS

The significant SNP method allows the detection of SNP-trait associations at much lower sequencing coverage than the current methods, leading to ~ 80% lower sequencing cost and making BSA-Seq more accessible to the research community and more applicable to the species with a large genome.

Genetic mapping and development of molecular markers for a candidate gene locus controlling rind color in watermelon

ClCG08G017810 (ClCGMenG) encoding a 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein is associated with formation of dark green versus light green rind color in watermelon. Rind color is an important agronomic trait in watermelon [Citrullus lanatus (Thunb.) Matsum. and Nakai], but the underlying molecular mechanism for this trait is not fully known. In the present study, we identified a single locus on chromosome 8 accounting for watermelon rind color (dark green vs. light green). Genetic analysis of F₁, F₂, and BC₁ populations derived from two parental lines (9904 with dark green rind and Handel with light green rind) revealed that the watermelon rind color (dark green vs. light green) is controlled by a single locus, and dark green is dominant to light green rind. Initial mapping revealed a region of interest spanning 2.07 Mb on chromosome 8. Genetic mapping with CAPS and SNP markers narrowed down the candidate region to 31.4 kb. Gene annotation of the corresponding region in the reference genome revealed the ClCG08G017810 gene sequence encoding the 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein. The sequence alignment of the candidate gene with the two parental lines suggested a nonsynonymous SNP mutation in the coding region of ClCG08G017810, converting an arginine (R) to glycine (G). The SNP might be associated with rind color of 103 watermelon germplasm lines investigated in this study. The qRT-PCR analysis revealed higher expression of ClCG08G017810 in dark green rind than in light green rind. Therefore, ClCG08G017810 is a candidate gene associated with watermelon rind color. The present study facilitates marker-assisted selection useful for the development of cultivars with desirable rind color.

Mapping Cucumber Vein Yellowing Virus Resistance in Cucumber (Cucumis sativus L.) by Using BSA-seq Analysis

Cucumber vein yellowing virus (CVYV) causes severe yield losses in cucurbit crops across Mediterranean countries. The control of this virus is based on cultural practices to prevent the presence of its vector (Bemisia tabaci) and breeding for natural resistance, which requires the identification of the loci involved and the development of molecular markers for linkage analysis. In this work, we mapped a monogenic locus for resistance to CVYV in cucumber by using a Bulked Segregant Analysis (BSA) strategy coupled with whole-genome resequencing. We phenotyped 135 F₃ families from a segregating population between a susceptible pickling cucumber and a resistant Long Dutch type cucumber for CVYV resistance. Phenotypic analysis determined the monogenic and incomplete dominance inheritance of the resistance. We named the locus CsCvy-1. For mapping this locus, 15 resistant and 15 susceptible homozygous F₂ individuals were selected for whole genome resequencing. By using a customized bioinformatics pipeline, we identified a unique region in chromosome 5 associated to resistance to CVYV, explaining more than 80% of the variability. The resequencing data provided us with additional SNP markers to decrease the interval of CsCvy-1 to 625 kb, containing 24 annotated genes. Markers flanking CsCvy-1 in a 5.3 cM interval were developed for marker-assisted selection (MAS) in breeding programs and will be useful for the identification of the target gene in future studies.