Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Long-read bitter gourd (Momordica charantia) genome and the genomic architecture of nonclassic domestication

The genetic architecture of quantitative traits is determined by both Mendelian and polygenic factors, yet classic examples of plant domestication focused on selective sweep of newly mutated Mendelian genes. Here we report the chromosome-level genome assembly and the genomic investigation of a nonclassic domestication example, bitter gourd (Momordica charantia), an important Asian vegetable and medicinal plant of the family Cucurbitaceae. Population resequencing revealed the divergence between wild and South Asian cultivars about 6,000 y ago, followed by the separation of the Southeast Asian cultivars about 800 y ago, with the latter exhibiting more extreme trait divergence from wild progenitors and stronger signs of selection on fruit traits. Unlike some crops where the largest phenotypic changes and traces of selection happened between wild and cultivar groups, in bitter gourd large differences exist between two regional cultivar groups, likely reflecting the distinct consumer preferences in different countries. Despite breeding efforts toward increasing female flower proportion, a gynoecy locus exhibits complex patterns of balanced polymorphism among haplogroups, with potential signs of selective sweep within haplogroups likely reflecting artificial selection and introgression from cultivars back to wild accessions. Our study highlights the importance to investigate such nonclassic example of domestication showing signs of balancing selection and polygenic trait architecture in addition to classic selective sweep in Mendelian factors.

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.

Molecular Mapping and Candidate Gene Analysis for GA3 Responsive Short Internode in Watermelon (Citrullus lanatus)

Plants with shorter internodes are suitable for high-density planting, lodging resistance and the preservation of land resources by improving yield per unit area. In this study, we identified a locus controlling the short internode trait in watermelon using Zhengzhouzigua (long internode) and Duan125 (short internode) as mapping parents. Genetic analysis indicated that F1 plants were consistent with long internode plants, which indicates that the long internode was dominant over the short internode. The observed F2 and BC1 individuals fitted the expected phenotypic segregation ratios of 3:1 and 1:1, respectively. The locus was mapped on chromosome 9 using a bulked segregant analysis approach. The region was narrowed down to 8.525 kb having only one putative gene, Cla015407, flanking by CAPS90 and CAPS91 markers, which encodes gibberellin 3β-hydroxylase (GA 3β-hydroxylase). The sequence alignment of the candidate gene between both parents revealed a 13 bp deletion in the short internode parent, which resulted in a truncated protein. Before GA3 application, significantly lower GA3 content and shorter cell length were obtained in the short internode plants. However, the highest GA3 content and significant increase in cell length were observed in the short internode plants after exogenous GA3 application. In the short internode plants, the expression level of the Cla015407 was threefold lower than the long internode plants in the stem tissue. In general, our results suggested that Cla015407 might be the candidate gene responsible for the short internode phenotype in watermelon and the phenotype is responsive to exogenous GA3 application.