Fine mapping and discovery of candidate genes for seed size in watermelon by genome survey sequencing

Genome-wide analysis of Citrus medica ABC transporters reveals the regulation of fruit development by CmABCB19 and CmABCC10

ATP-binding cassette (ABC) transporters are vital for plant growth and development as they facilitate the transport of essential molecules. Despite the family's significance, limited information exists about its functional distinctions in Citrus medica. Our study identified 119 genes encoding ABC transporter proteins in the C. medica genome. Through an evolutionary tree and qPCR analysis, two ABC genes, CmABCB19 and CmABCC10, were implicated in C. medica fruit development, showing upregulation in normal fruits compared to malformed fruits. CmABCB19 was found to localize to the plasma membrane of Nicotiana tabacum, exhibiting indole-3-acetic acid (IAA) efflux activity in the yeast mutant strain yap1. CmABCC10, a tonoplast-localized transporter, exhibited efflux of diosmin, nobiletin, and naringin, with rutin influx in strain ycf1. Transgenic expression of CmABCB19 and CmABCC10 in Arabidopsis thaliana induced alterations in auxin and flavonoid content, impacting silique and seed size. This effect was attributed to the modulation of structural genes in the auxin biosynthesis (YUC5/9, CYP79B2, CYP83B1, SUR1) and flavonoid biosynthesis (4CL2/3, CHS, CHI, FLS1/3) pathways. In summary, the functional characterization of CmABCB19 and CmABCC10 illuminates auxin and flavonoid transport, offering insights into their interplay with biosynthetic pathways and providing a foundation for understanding the transporter's role in fruit development.

Natural Allelic Variations of Bch10G006400 Controlling Seed Size in Chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How)

Seeds are the most important reproductive organs of higher plants, the beginning and end of a plant's lifecycle. They are very important to plant growth and development, and also an important factor affecting yield. In this study, genetic analysis and BSA-seq of the F2 population crossed with the large-seeded material 'J16' and small-seeded material 'FJ5' were carried out, and the seed size locus was initially located within the 1.31 Mb region on chr10. In addition, 2281 F2 plants were used to further reduce the candidate interval to 48.8 Kb. This region contains only one gene encoding the N-acetyltransferase (NAT) protein (Bch10G006400). Transcriptome and expression analysis revealed that the gene was significantly more highly expressed in 'J16' than in 'FJ5'. Variation analysis of Bch10G006400 among parents and 50 chieh-qua germplasms revealed that as well as a nonsynonymous mutation (SNP_314) between parents, two mutations (SNP_400 and InDel_551) were detected in other materials. Combining these three mutations completely distinguished the seed size of the chieh-qua. GO and KEGG enrichment analyses revealed that DGEs played the most important roles in carbohydrate metabolism and plant hormone signal transduction, respectively. The results of this study provide important information for molecular marker-assisted breeding and help to reveal the molecular mechanism of seed size.

A nonsynonymous mutation in BhLS, encoding an acyl-CoA N-acyltransferase leads to fruit and seed size variation in wax gourd (Benincasa hispida)

Wax gourd (Benincasa hispida (Thunb.) Cogn., 2n = 2x = 24) is an economically important vegetable crop cultivated widely in many tropical and subtropical regions, including China, India, and Japan. Both fruit and seeds are prized agronomic attributes in wax gourd breeding and production. However, the genetic mechanisms underlying these traits remain largely unexplored. In this study, we observed a strong correlation between fruit size and seed size variation in our mapping population, indicating genetic control by a single gene, BhLS, with large size being dominant over small. Through bulk segregant analysis sequencing and fine mapping with a large F2 population, we precisely located the BhLS gene within a 47.098-kb physical interval on Chromosome 10. Within this interval, only one gene, Bhi10M000649, was identified, showing homology to Arabidopsis HOOKLESS1. A nonsynonymous mutation (G to C) in the second exon of Bhi10M000649 was found to be significantly associated with both fruit and seed size variation in wax gourd. These findings collectively highlight the pleiotropic effect of the BhLS gene in regulating fruit and seed size in wax gourd. Our results offer molecular insights into the variation of fruit and seed size in wax gourd and establish a fundamental framework for breeding wax gourd cultivars with desired traits.

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.

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Fine mapping of the major gene BhHLS1 controlling seed size in wax gourd (Benincasa hispida)

Introduction/Background

The seed size of wax gourds is an important agronomic trait; however, the associated genes have not yet been reported.

Methods

In this study, we used a high-density genetic map constructed based on F8 recombinant inbred line populations derived from a cross between MY-1 (large seed) and GX-71 (small seed) strains to detect quantitative trait locis (QTLs) for seed-size-related traits in wax gourd over a two-year period.

Results

Two stable QTLs (qSL10 and qSW10) for seed length (SL) and seed width (SW) on chromosome 10 were repeatedly detected over two years (2021-2022). qSL10 had a phenotypic variation rate of 75.30% and 80.80% in 2021 and 2022, respectively. Whereas, qSW10 had a phenotypic variation rate of 66.60% and 73.80% in 2021 and 2022, respectively. Further, a single nucleotide polymorphism mutation was found to cause early termination of Bch10G006400 (BhHLS1) translation in GX-71 through sequencing analysis of candidate genes. Based on gene functional annotation and quantitative real-time PCR analyses, BhHLS1 encoded a probable N-acetyltransferase HLS1-like protein and its expression level was significantly different between parents. Therefore, BhHLS1 is a major candidate gene associated with a one-factor polymorphism regulating the SL and SW of wax gourds. Finally, based on variation in the BhHLS1 sequence, a cleaved amplified polymorphic sequence marker was developed for the molecular marker-assisted breeding of wax gourds.

Discussion

Overall, this study is of great significance for the genetic improvement of seed size, verification of gene functions, and cultivation of specific germplasm resources for wax gourds.

Integrated analysis of biparental and natural populations reveals CRIB domain-containing protein underlying seed coat crack trait in watermelon

The scc locus of the watermelon seed coat crack trait was fine mapped on chromosome 3. Cla97C03G056110 (annotated as CRIB domain-containing protein) was regarded as the most likely candidate gene Seed coat crack (scc) is a special characteristic of watermelon compared with other cucurbit crops. However, information regarding the genetic basis of this trait is limited. We conducted a genetic analysis of six generations derived from PI 192938 (scc) and Cream of Saskatchewan (COS) (non-scc) parental lines and found that the scc trait was regulated by a single recessive gene through two years. Bulk segregant analysis sequencing (BSA-seq) and initial mapping placed the scc locus into an 808.8 kb region on chromosome 3. Evaluation of another 1152 F₂ plants narrowed the scc locus to a 277.11 kb region containing 37 candidate genes. Due to the lack of molecular markers in the fine-mapping interval, we extracted the genome sequence variations in this 277.11 kb region with in silico BSA among seventeen re-sequenced lines (6 scc and 11 non-scc) and finally delimited the scc locus to an 8.34 kb region with only one candidate gene Cla97C03G056110 (CRIB domain-containing protein). Three single nucleotide polymorphism loci in the promoter region of Cla97C03G056110 altered cis-acting elements that were highly correlated with the nature watermelon panel. The expression of Cla97C03G056110 in seed coat tissue was higher in non-scc than in scc lines and was specifically expressed in seed coat compared with fruit flesh.

Identification of the egusi seed trait locus (eg) and its suppressor gene associated with the thin seed coat trait in watermelon

Egusi watermelon has a unique egusi seed type, which could be useful for breeding both edible seeds and edible flesh in watermelon. However, the genetic basis of the unique egusi seed type is not clear. In the present study, we first reported that at least two genes with inhibitory epistasis were responsible for the thin seed coat (unique egusi seed type) in watermelon. Inheritance analysis of five populations, including F₂, BC, and BCF₂, suggested that the thin seed coat trait was controlled by a suppressor gene together with the egusi seed locus (eg) in egusi watermelon. Based on high-throughput sequencing technology, two quantitative trait loci located on chromosome 1 and chromosome 6 were identified for the thin seed coat trait in watermelon. One of the loci, the eg locus on chromosome 6, was finely mapped to a genomic region of 15.7 kb, which contained only one candidate gene. Comparative transcriptome analysis highlighted differentially expressed genes involved in cellulose and lignin synthesis between watermelon genotypes varying in the thickness of the seed coat and provided several potential candidate genes for the thin seed coat trait. Taken together, our data suggest that at least two genes are complementarily involved in the thin seed coat trait and will be useful for cloning novel genes. The results presented here provide a new reference for uncovering egusi seed genetic mechanisms and valuable information for marker-assisted selection in seed coat breeding.

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.

ATP-binding cassette transporters in nonmodel plants

Knowledge about plant ATP-binding cassette (ABC) proteins is of great value for sustainable agriculture, economic yield, and the generation of high-quality products, especially under unfavorable growth conditions. We have learned much about ABC proteins in model organisms, notably Arabidopsis thaliana; however, the importance of research dedicated to these transporters extends far beyond Arabidopsis biology. Recent progress in genomic and transcriptomic approaches for nonmodel and noncanonical model plants allows us to look at ABC transporters from a wider perspective and consider chemodiversity and functionally driven adaptation as distinctive mechanisms during their evolution. Here, by considering several representatives from agriculturally important families and recent progress in functional characterization of nonArabidopsis ABC proteins, we aim to bring attention to understanding the evolutionary background, distribution among lineages and possible mechanisms underlying the adaptation of this versatile transport system for plant needs. Increasing the knowledge of ABC proteins in nonmodel plants will facilitate breeding and development of new varieties based on, for example, genetic variations of endogenous genes and/or genome editing, representing an alternative to transgenic approaches.

Genome-wide association analysis provides molecular insights into the natural variation of watermelon seed size

Seed-consumption watermelon tend to have larger-sized seeds, while flesh-consumed watermelon often require relatively smaller seed. Therefore, the seed size of watermelon has received extensive attention from consumers and breeders. However, the study on the natural variation and genetic mechanism of watermelon seed size is not clear enough. In the present study, 100 seed weight, seed hilum length, seed length, seed width, and seed thickness in 197 watermelon accessions were examined. Furthermore, association analysis was conducted between seed size traits and high-quality SNP data. The results revealed that there was a strong correlation between the five seed traits. And seed enlargement was an important feature during watermelon seed size domestication. Meanwhile, the seed consumption biological species C. mucosospermu and C. lanatus edible seed watermelon had a significantly bigger seed size than other species's. Eleven non-repeating significant SNPs above the threshold line were obtained by GWAS analysis. Four of them on chromosome 5 were considered to be closely associated with seed size traits, i.e. S5: 32250307, S5: 32250454, S5: 32256177, S5: 32260870, which could be used as potential molecular markers for the breeding of watermelon cultivars with target seed size. In addition, combined with gene annotation information and previous reports, five genes near the four significant SNPs may regulate seed size. And qRT-PCR analysis showed that two genes Cla97C05G104360 and Cla97C05G104380, which may be involved in abscisic acid metabolism, may play an important role in regulating the seed size of watermelon. Our findings provide molecular insights into natural variation in watermelon seed size, and gives valuable information of molecular marker-assisted breeding.

Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants

Simple sequence repeats (SSRs) are one of the most important genetic markers and widely exist in most species. Here, we identified 249,822 SSRs from 3,951,919 genes in 112 plants. Then, we conducted a comprehensive analysis of these SSRs and constructed a plant SSR database (PSSRD). Interestingly, more SSRs were found in lower plants than in higher plants, showing that lower plants needed to adapt to early extreme environments. Four specific enriched functional terms in the lower plant Chlamydomonas reinhardtii were detected when it was compared with seven other higher plants. In addition, Guanylate_cyc existed in more genes of lower plants than of higher plants. In our PSSRD, we constructed an interactive plotting function in the chart interface, and users can easily view the detailed information of SSRs. All SSR information, including sequences, primers, and annotations, can be downloaded from our database. Moreover, we developed Web SSR Finder and Batch SSR Finder tools, which can be easily used for identifying SSRs. Our database was developed using PHP, HTML, JavaScript, and MySQL, which are freely available at http://www.pssrd.info/ . We conducted an analysis of the Myb gene families and flowering genes as two applications of the PSSRD. Further analysis indicated that whole-genome duplication and whole-genome triplication played a major role in the expansion of the Myb gene families. These SSR markers in our database will greatly facilitate comparative genomics and functional genomics studies in the future.

CRISPR/Cas9-mediated mutagenesis of ClBG1 decreased seed size and promoted seed germination in watermelon

Abscisic acid (ABA) is a critical regulator of seed development and germination. β-glucosidases (BGs) have been suggested to be contributors to increased ABA content because they catalyze the hydrolysis of ABA-glucose ester to release free ABA. However, whether BGs are involved in seed development is unclear. In this study, a candidate gene, ClBG1, in watermelon was selected for targeted mutagenesis via the CRISPR/Cas9 system. Seed size and weight were significantly reduced in the Clbg1-mutant watermelon lines, which was mainly attributed to decreased cell number resulting from decreased ABA levels. A transcriptome analysis showed that the expression of 1015 and 1429 unique genes was changed 10 and 18 days after pollination (DAP), respectively. Cytoskeleton- and cell cycle-related genes were enriched in the differentially expressed genes of wild type and Clbg1-mutant lines during seed development. Moreover, the expression of genes in the major signaling pathways of seed size control was also changed. In addition, seed germination was promoted in the Clbg1-mutant lines due to decreased ABA content. These results indicate that ClBG1 may be critical for watermelon seed size regulation and germination mainly through the modulation of ABA content and thereby the transcriptional regulation of cytoskeleton-, cell cycle- and signaling-related genes. Our results lay a foundation for dissecting the molecular mechanisms of controlling watermelon seed size, a key agricultural trait of significant economic importance.

Ethylene-responsive factor 4 is associated with the desirable rind hardness trait conferring cracking resistance in fresh fruits of watermelon

Fruit rind plays a pivotal role in alleviating water loss and disease and particularly in cracking resistance as well as the transportability, storability and shelf-life quality of the fruit. High susceptibility to cracking due to low rind hardness is largely responsible for severe annual yield losses of fresh fruits such as watermelon in the field and during the postharvest process. However, the candidate gene controlling the rind hardness phenotype remains unclear to date. Herein, we report, for the first time, an ethylene-responsive transcription factor 4 (ClERF4) associated with variation in rind hardness via a combinatory genetic map with bulk segregant analysis (BSA). Strikingly, our fine-mapping approach revealed an InDel of 11 bp and a neighbouring SNP in the ClERF4 gene on chromosome 10, conferring cracking resistance in F₂ populations with variable rind hardness. Furthermore, the concomitant kompetitive/competitive allele-specific PCR (KASP) genotyping data sets of 104 germplasm accessions strongly supported candidate ClERF4 as a causative gene associated with fruit rind hardness variability. In conclusion, our results provide new insight into the underlying mechanism controlling rind hardness, a desirable trait in fresh fruit. Moreover, the findings will further enable the molecular improvement of fruit cracking resistance in watermelon via precisely targeting the causative gene relevant to rind hardness, ClERF4.

Quantitative Trait Loci for Seed Size Variation in Cucurbits - A Review

Cucurbits (Cucurbitaceae family) include many economically important fruit vegetable crops such as watermelon, pumpkin/squash, cucumber, and melon. Seed size (SS) is an important trait in cucurbits breeding, which is controlled by quantitative trait loci (QTL). Recent advances have deciphered several signaling pathways underlying seed size variation in model plants such as Arabidopsis and rice, but little is known on the genetic basis of SS variation in cucurbits. Here we conducted literature review on seed size QTL identified in watermelon, pumpkin/squash, cucumber and melon, and inferred 14, 9 and 13 consensus SS QTL based on their physical positions in respective draft genomes. Among them, four from watermelon (ClSS2.2, ClSS6.1, ClSS6.2, and ClSS8.2), two from cucumber (CsSS4.1 and CsSS5.1), and one from melon (CmSS11.1) were major-effect, stable QTL for seed size and weight. Whole genome sequence alignment revealed that these major-effect QTL were located in syntenic regions across different genomes suggesting possible structural and functional conservation of some important genes for seed size control in cucurbit crops. Annotation of genes in the four watermelon consensus SS QTL regions identified genes that are known to play important roles in seed size control including members of the zinc finger protein and the E3 ubiquitin-protein ligase families. The present work highlights the utility of comparative analysis in understanding the genetic basis of seed size variation, which may help future mapping and cloning of seed size QTL in cucurbits.

Discovery of the Genomic Region and Candidate Genes of the Scarlet Red Flesh Color (Yscr ) Locus in Watermelon (Citrullus Lanatus L.)

The flesh color of watermelon (Citrullus lanatus) is an important fruit quality trait that helps to determine fruit attractiveness and is potentially beneficial to human health. Previous inheritance analyses determined that a single dominant gene, Yscr , produces the scarlet red flesh color rather than the coral red flesh color in watermelon. However, no genomic region or gene-based molecular markers for the locus Yscr have been reported thus far. In the present study, two high-density genetic maps and whole-genome variation detection aided by genome resequencing were first map the flesh color locus Yscr to a small region on chromosome 6 based on two independent populations derived from two scarlet red-fleshed lines and two coral red-fleshed lines. Two major quantitative trait loci located in the same genomic regions were identified in the F2 and BC1P2 populations and explained 90.36% and 75.1% of the phenotypic variation in flesh color, respectively. Based on the genetic variation in the two parental lines, newly developed PCR-based markers narrowed the Yscr region to 40 Kb. Of the five putative genes in this region, four encoded glycine-rich cell wall structural proteins, which implied that a new regulatory mechanism might occur between scarlet red- and coral red-fleshed in watermelon. Moreover, the genotypes of two newly developed InDel markers (InDel27_fc6 and InDel28_fc6) were completely consistent with the phenotypes in the F2 and BC1P2 populations and all 56 scarlet red-fleshed watermelon accessions. The results presented here provide valuable information for marker-assisted selection of flesh color breeding and the functional validation of candidate genes in watermelon.