Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement.

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops.

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Structural variation of GL1 gene determines the trichome formation in Brassica juncea.

KEY MESSAGE: A 448 kb region on chromosome B02 was delimited to be associated with trichome trait in Brassica juncea, in which the BjuVB02G54610 gene with a structural variation of 3 kb structure variation (SV) encoding a MYB transcription factor was predicted as the possible candidate gene. Mustards (Brassica juncea) are allopolyploid crops in the worldwide, and trichomes are essential quality attributes that significantly influence its taste and palpability in vegetable-use cultivars. As important accessory tissues from specialized epidermal cells, trichomes also play an important role in mitigating biotic and abiotic stresses. In this study, we constructed a F2 segregating population using YJ27 with intensive trichome leaves and 03B0307 with glabrous leaves as parents. By bulked segregant analysis (BSA-seq), we obtained a 2.1 Mb candidate region on B02 chromosome associated with the trichome or glabrous trait formation. Then, we used 13 Kompetitive Allele Specific PCR (KASP) markers for fine mapping and finally narrowed down the candidate region to about 448 kb in length. Interestingly, among the region, there was a 3 kb sequence deletion that located on the BjuVB02G54610 gene in the F2 individuals with trichome leaves. Genotyping results of F2 populations confirmed this deletion (R2 = 81.44%) as a major QTL. Natural population re-sequencing analysis and genotyping results further validated the key role of the 3 kb structure variation (SV) of insertion/deletion type in trichome development in B. juncea. Our findings provide important information on the formation of trichomes and potential target gene for breeding vegetable mustards.

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon.

Drought stress has detrimental effects on crop productivity worldwide. A strong root system is crucial for maintaining water and nutrients uptake under drought stress. Wild watermelons possess resilient roots with excellent drought adaptability. However, the genetic factors controlling this trait remain uninvestigated. In this study, we conducted a bulk segregant analysis (BSA) on an F2 population consisting of two watermelon genotypes, wild and domesticated, which differ in their lateral root development under drought conditions. We identified two quantitative trait loci (qNLR_Dr. Chr01 and qNLR_Dr. Chr02) associated with the lateral root response to drought. Furthermore, we determined that a small region (0.93 Mb in qNLR_Dr. Chr01) is closely linked to drought adaptation through quantitative trait loci (QTL) validation and fine mapping. Transcriptome analysis of the parent roots under drought stress revealed unique effects on numerous genes in the sensitive genotype but not in the tolerant genotype. By integrating BSA, fine mapping, and the transcriptome, we identified six genes, namely L-Ascorbate Oxidase (AO), Cellulose Synthase-Interactive Protein 1 (CSI1), Late Embryogenesis Abundant Protein (LEA), Zinc-Finger Homeodomain Protein 2 (ZHD2), Pericycle Factor Type-A 5 (PFA5), and bZIP transcription factor 53-like (bZIP53-like), that might be involved in the drought adaptation. Our findings provide valuable QTLs and genes for marker-assisted selection in improving water-use efficiency and drought tolerance in watermelon. They also lay the groundwork for the genetic manipulation of drought-adapting genes in watermelon and other Cucurbitacea species.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops.

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

A natural mutation of the NST1 gene arrests secondary cell wall biosynthesis in the seed coat of a hull-less pumpkin accession.

Hull-less pumpkins (Cucurbita pepo L.) are naturally occurring novel variants known as oilseed or naked-seeded pumpkins, and are characterized by the absence of a normal lignified seed coat. Due to a specialized seed coat structure, these variants serve as a good model for studying seed coat formation and simplify the processing of pumpkin seeds. However, causal genes for this hull-less trait still remain unknown. Here, by bulked segregant analysis and fine mapping, we found that mutation of a single gene, NAC SECONDARY WALL THICKENING PROMOTING FACTOR 1 (NST1), accounts for the hull-less trait. A 14-bp sequence insertion in the CpNST1 gene causes premature termination of CpNST1 translation, leading to lack of secondary cell wall (SCW) biosynthesis in hull-less seed coats. In situ hybridization analysis provided further evidence for the role of CpNST1 in pumpkin seed coat SCW biosynthesis. Interestingly, through secondary cell wall compositional analysis, we found that the main SCW components differed among cell layers in the seed coat. RNA-seq analysis indicated an upstream role of CpNST1 in the SCW biosynthesis network. Collectively, our findings provide mechanistic insight into seed coat SCW biosynthesis, and a target gene for breeders to introduce this hull-less trait for commercial exploitation.

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

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

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 F2 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.

Characterization of watermelon anther and its programmed cell death-associated events during dehiscence under cold stress.

KEY MESSAGE: The 'neglected' thermophile fruit crop of watermelon was first used as a model crop to study the PCD associated with anther dehiscence in cold-exposed condition during anther development. Anther dehiscence ensures normal pollen release and successful fertilization at fruit-setting stages in flowering plants. However, most researches pertinent to anther dehiscence are centered on model plant and/or major field crops under optimal growth condition. Due to anther indehiscence in cold condition, crop plants of thermophile tropical or subtropical fruit crops fail to accomplish timely pollination and fertilization, resulting in a great yield loss annually. Herein, we developed an ideal model crop for studying the programmed cell death (PCD) associated with anther dehiscence under low-temperature stress using the S-shaped spiral anther in watermelon as instead. Our results revealed that, including the tapetal cell layers, both cells of the interlocular septum and the stomium were blocked in PCD associated with anther dehiscence at 15 °C. Likewise, TUNEL assays visualized the evidence that low temperature at 15 °C interferes with not only the PCD of tapetal cells, but also the PCD of interlocular septum and stomium. Furthermore, the expressions of genes correlated with PCD of tapetum and stomium were significantly inhibited at 15 °C, suggesting that low temperature affects anther dehiscence by inhibiting PCD of sporophytic tissue-related gene expressions. The findings of the current research provide mechanistic insights into anther indehiscence leading to poor fruit-setting for thermophile fruit crop such as watermelon under adverse cold condition in flowering.