QTL mapping in multiple populations and development stages reveals dynamic quantitative trait loci for fruit size in cucumbers of different market classes

Genetic and molecular regulation of fruit development in cucumber

Fruit development can be generally classified into a set of biologically sequential stages including fruit initiation, growth, and ripening. Cucumber, a globally important vegetable crop, displays two important features during fruit development: parthenocarpy at fruit initiation and prematurity at harvest for consumption. Therefore, fruit growth plays essential role for cucumber yield and quality formation, and has become the research hot spot in cucumber fruit development. Here, we describe recent advances in molecular mechanisms underlying fruit growth in cucumber, include key players and regulatory networks controlling fruit length variation, fruit neck elongation, and locule development. We also provide insights into future directions for scientific research and breeding strategies in cucumber.

The mutation of CsSUN, an IQD family protein, is responsible for the short and fat fruit (sff) in cucumber (Cucumis sativus L.)

The fruit shape of cucumber is an important agronomic trait, and mining regulatory genes, especially dominant ones, is vital for cucumber breeding. In this study, we identified a short and fat fruit mutant, named sff, from an EMS mutagenized population. Compared to the CCMC (WT), sff (MT) exhibited reduced fruit length and increased dimeter. Segregation analysis revealed that the sff phenotype is controlled by a semi-dominant single gene with dosage effects. Through map-based cloning, the SFF locus was narrowed down to a 52.6 kb interval with two SNPs (G651A and C1072T) in the second and third exons of CsaV3_1G039870, which encodes an IQD family protein, CsSUN. The G651A within the IQ domain of CsSUN was identified as the unique SNP among 114 cucumber accessions, and it was the primary cause of the functional alteration in CsSUN. By generating CsSUN knockout lines in cucumber, we confirmed that CsSUN was responsible for sff mutant phenotype. The CsSUN is localized to the plasma membrane. CsSUN exhibited the highest expression in the fruit with lower expression in sff compared to WT. Histological observations suggest that the sff mutant phenotype is due to increased transverse cell division and inhibited longitudinal cell division. Transcriptome analysis revealed that CsSUN significantly affected the expression of genes related to cell division, expansion, and auxin signal transduction. This study unveils CsSUN's crucial role in shaping cucumber fruit and offers novel insights for cucumber breeding.

Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.)

Fruit size is a crucial agronomic trait in bottle gourd, impacting both yield and utility. Despite its significance, the regulatory mechanism governing fruit size in bottle gourd remains largely unknown. In this study, we used bottle gourd (small-fruited H28 and large-fruited H17) parent plants to measure the width and length of fruits at various developmental stages, revealing a single 'S' growth curve for fruit expansion. Paraffin section observations indicated that both cell number and size significantly influence bottle gourd fruit size. Through bulked segregant analysis and combined genotype-phenotype analysis, the candidate interval regulating fruit size was pinpointed to 17,747,353 bp-18,185,825 bp on chromosome 9, encompassing 0.44 Mb and including 44 genes. Parental fruits in the rapid expansion stage were subjected to RNA-seq, highlighting that differentially expressed genes were mainly enriched in pathways related to cell wall biosynthesis, sugar metabolism, and hormone signaling. Transcriptome and resequencing analysis, combined with gene function annotation, identified six genes within the localized region as potential regulators of fruit size. This study not only maps the candidate interval of genes influencing fruit size in bottle gourd through forward genetics, but also offers new insights into the potential molecular mechanisms underlying this trait through transcriptome analysis.

CucumberAI: Cucumber Fruit Morphology Identification System Based on Artificial Intelligence

Cucumber is an important vegetable crop that has high nutritional and economic value and is thus favored by consumers worldwide. Exploring an accurate and fast technique for measuring the morphological traits of cucumber fruit could be helpful for improving its breeding efficiency and further refining the development models for pepo fruits. At present, several sets of measurement schemes and standards have been proposed and applied for the characterization of cucumber fruits; however, these manual methods are time-consuming and inefficient. Therefore, in this paper, we propose a cucumber fruit morphological trait identification framework and software called CucumberAI, which combines image processing techniques with deep learning models to efficiently identify up to 51 cucumber features, including 32 newly defined parameters. The proposed tool introduces an algorithm for performing cucumber contour extraction and fruit segmentation based on image processing techniques. The identification framework comprises 6 deep learning models that combine fruit feature recognition rules with MobileNetV2 to construct a decision tree for fruit shape recognition. Additionally, the framework employs U-Net segmentation models for fruit stripe and endocarp segmentation, a MobileNetV2 model for carpel classification, a ResNet50 model for stripe classification and a YOLOv5 model for tumor identification. The relationships between the image-based manual and algorithmic traits are highly correlated, and validation tests were conducted to perform correlation analyses of fruit surface smoothness and roughness, and a fruit appearance cluster analysis was also performed. In brief, CucumberAI offers an efficient approach for extracting and analyzing cucumber phenotypes and provides valuable information for future cucumber genetic improvements.

Identification of Key Genes of Fruit Shape Variation in Jujube with Integrating Elliptic Fourier Descriptors and Transcriptome

Jujube (Ziziphus jujuba) exhibits a rich diversity in fruit shape, with natural occurrences of gourd-like, flattened, and other special shapes. Despite the ongoing research into fruit shape, studies integrating elliptical Fourier descriptors (EFDs) with both Short Time-series Expression Miner (STEM) and weighted gene co-expression network analysis (WGCNA) for gene discovery remain scarce. In this study, six cultivars of jujube fruits with distinct shapes were selected, and samples were collected from the fruit set period to the white mature stage across five time points for shape analysis and transcriptome studies. By combining EFDs with WGCNA and STEM, the study aimed to identify the critical periods and key genes involved in the formation of jujube fruit shape. The findings indicated that the D25 (25 days after flowering) is crucial for the development of jujube fruit shape. Moreover, ZjAGL80, ZjABI3, and eight other genes have been implicated to regulate the shape development of jujubes at different periods of fruit development, through seed development and fruit development pathway. In this research, EFDs were employed to precisely delineate the shape of jujube fruits. This approach, in conjunction with transcriptome, enhanced the precision of gene identification, and offered an innovative methodology for fruit shape analysis. This integration facilitates the advancement of research into the morphological characteristics of plant fruits, underpinning the development of a refined framework for the genetic underpinnings of fruit shape variation.

Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L)

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20 cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087 bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

Identification of major QTLs for drought tolerance in soybean, together with a novel candidate gene, GmUAA6

Drought tolerance is a complex trait in soybean that is controlled by polygenetic quantitative trait loci (QTLs). In this study, wilting score, days-to-wilting, leaf relative water content, and leaf relative conductivity were used to identify QTLs associated with drought tolerance in recombinant inbred lines derived from a cross between a drought-sensitive variety, Lin, and a drought-tolerant variety, Meng. A total of 33 drought-tolerance QTLs were detected. Of these 17 were major QTLs. In addition, 15 were novel drought-tolerance QTLs. The most predominant QTL was on chromosome 11. This was detected in at least three environments. The overlapped mapping interval of the four measured traits was 0.2 cM in genetic distance (about 220 kb in physical length). Glyma.11g143500 (designated as GmUAA6), which encodes a UDP-N-acetylglucosamine transporter, was identified as the most likely candidate gene. The allele of GmUAA6 from Lin (GmUAA6Lin) was associated with improved soybean drought tolerance. Overexpression of GmUAA6Lin in Arabidopsis and soybean hairy roots enhanced drought tolerance. Furthermore, a 3-bp insertion/deletion (InDel) in the coding sequence of GmUAA6 explained up to 49.9% of the phenotypic variation in drought tolerance-related traits, suggesting that this InDel might be used in future marker-assisted selection of drought-tolerant lines in soybean breeding programs.

Morphological Analyses and QTL Mapping of Mottled Leaf in Zucchini (Cucurbita pepo L.)

The mottled leaf is one of the agronomic traits of zucchini and can be applied as a marker trait in aggregation breeding. However, the genetic mechanism responsible for mottled leaf has yet to be elucidated. In the present study, we used two inbred lines (line '19': silver mottled leaf; line '113': normal leaf) as parents for the physiological and genetic analysis of mottled leaf. The synthesis and net photosynthetic rate of chlorophyll were not significantly affected in the mottled areas of leaves. However, we detected a large space between the palisade parenchyma in the leaf mottle area of line '19', which may have caused the mottled leaf phenotype. Light also plays an important role in the formation of mottled leaf, and receiving light during the early stages of leaf development is a necessary factor. Genetic analysis has previously demonstrated that mottled leaf is a quantitative trait that is controlled by multiple genes. Based on the strategy of quantitative trait locus sequencing (QTL-seq), two QTLs were identified on chromosomes 1 and 17, named CpML1.1 and CpML17.1, respectively. Two major loci were identified using R/qtl software version 1.66 under greenhouse conditions in April 2019 (2019A) and April 2020 (2020A) and under open cultivation conditions in May 2020 (2020M). The major QTL, CpML1.1, was located in a 925.2-kb interval on chromosome 1 and explained 10.51%-24.15% of the phenotypic variation. The CpML17.1 was located in a 719.7-kb interval on chromosome 17 and explained 16.25%-38.68% of the phenotypic variation. Based on gene annotation, gene sequence alignment, and qRT-PCR analysis, the Cp4.1LG01g23790 at the CpML1.1 locus encoding a protein of the TPX2 family (target protein of Xklp2) may be a candidate gene for mottled leaf in zucchini. Our findings may provide a theoretical basis for the formation of mottled leaf and provide a foundation for the fine mapping of genes associated with mottled leaf. Molecular markers closely linked to mottled leaf can be used in molecular-assisted selection for the zucchini mottled leaf breeding.

QTL mapping and transcriptomic analysis of fruit length in cucumber

A total of 151 recombinant inbred lines (RILs) were derived from the cross between 'Cucumis sativus L. hardwickii' (HW) and a cultivated Northern Chinese inbred line 'XinTaiMiCi' (XTMC). We used resequencing to construct the genetic map and analyze the genetic background of RIL population, and combined with the phenotypes of RIL population and the analysis of RNA-seq data, we located the major loci controlling the fruit length of cucumber and related analysis. A genetic map containing 600 bin markers was constructed via re-sequencing. Based on the phenotype data collected in two different seasons (spring 2021 and autumn 2022), the major quantitative trait loci (QTLs) controlling cucumber fruit length were located and their transcriptomic analysis carried out. The results revealed three QTLs (Fl2.1, Fl4.1, and Fl6.1) detected repeatedly in the two seasons, of which Fl4.1 was the dominant QTL. From the functional annotation of corresponding genes there, we discovered the gene Csa4G337340 encoding an auxin efflux carrier family protein. The expression of that gene was significantly lower in XTMC and the long-fruit RIL lines than in HW and the short-fruit RIL lines; hence, we speculated the gene could be negatively correlated with the fruit length of cucumber. Transcriptomic analysis showed that 259 differentially expressed genes (DEGs) were enriched in the plant hormone signal transduction pathway. In addition, among those DEGs, 509 transcription factors were detected, these distributed in several transcription factor gene families, such as bHLH, AP2/ErF -ERF, C2H2, and NAC. Therefore, we concluded that the major gene controlling the fruit length of cucumber is located in the interval of Fl4.1, whose gene Csa4G337340 may be involved in the negative regulation of fruit length. Further, genes related to plant hormone signal transduction and several transcription factors were also found involved in the regulation of cucumber fruit length. Our results provide a reference for the fine mapping of major genes and analyzing the mechanism of cucumber fruit length.

QTL mapping and stability analysis of trichome density in zucchini (Cucurbita pepo L.)

Trichomes provide an excellent model for studying cell differentiation and proliferation. The aboveground tissues of plants with long dense trichomes (LDTs) can cause skin itching in people working in a zucchini field, in which management, pollination, and fruit harvesting are difficult. In this study, an F2 population was constructed with the LDT inbred line "16" and the sparse micro trichome (SMT) inbred line "63" for QTL analysis of type I and II trichome density. Two QTLs were identified on chromosomes 3 and 15 using the QTL-seq method. Additionally, 191 InDel markers were developed on 20 chromosomes, a genetic map was constructed for QTL mapping, and three QTLs were identified on chromosomes 3, 6, and 15. Two QTLs, CpTD3.1 and CpTD15.1, were identified in both QTL-seq and genetic map-based QTL analyses, and CpTD15.1 was the major-effect QTL. The stability of CpTD3.1 and CpTD15.1 was confirmed using data from F2 plants under different environmental conditions. The major-effect QTL CpTD15.1 was located between markers chr15-4991349 and chr15-5766791, with a physical distance of 775.44 kb, and explained 12.71%-29.37% of the phenotypic variation observed in the three environments. CpTD3.1 was located between markers chr3-218350 and chr3-2891236, in a region with a physical distance of 2,672.89 kb, and explained 5.00%-10.64% of the phenotypic variation observed in the three environments. The functional annotations of the genes within the CpTD15.1 region were predicted, and five genes encoding transcription factors regulating trichome development were selected. Cp4.1LG15g04400 encoded zinc finger protein (ZFP) and harbored nonsynonymous SNPs in the conserved ring finger domain between the two parental lines. There were significant differences in Cp4.1LG15g04400 expression between "16" and "63", and a similar pattern was found between germplasm resources of LDT lines and SMT lines. It was presumed that Cp4.1LG15g04400 might regulate trichome density in zucchini. These results lay a foundation for better understanding the density of multicellular nonglandular trichomes and the regulatory mechanism of trichome density in zucchini.

Molecular and genetic regulations of fleshy fruit shape and lessons from Arabidopsis and rice

Fleshy fruit shape is an important external quality trait influencing the usage of fruits and consumer preference. Thus, modification of fruit shape has become one of the major objectives for crop improvement. However, the underlying mechanisms of fruit shape regulation are poorly understood. In this review we summarize recent progress in the genetic basis of fleshy fruit shape regulation using tomato, cucumber, and peach as examples. Comparative analyses suggest that the OFP-TRM (OVATE Family Protein - TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathways are probably conserved in regulating fruit shape by primarily modulating cell division patterns across fleshy fruit species. Interestingly, cucumber homologs of FRUITFULL (FUL1), CRABS CLAW (CRC) and 1-aminocyclopropane-1-carboxylate synthase 2 (ACS2) were found to regulate fruit elongation. We also outline the recent progress in fruit shape regulation mediated by OFP-TRM and IQD pathways in Arabidopsis and rice, and propose that the OFP-TRM pathway and IQD pathway coordinate regulate fruit shape through integration of phytohormones, including brassinosteroids, gibberellic acids, and auxin, and microtubule organization. In addition, functional redundancy and divergence of the members of each of the OFP, TRM, and IQD families are also shown. This review provides a general overview of current knowledge in fruit shape regulation and discusses the possible mechanisms that need to be addressed in future studies.

QTL Mapping for Ovary- and Fruit-Related Traits in Cucumis sativus-C. hystrix Introgression Line IL52

IL52 is a valuable introgression line obtained from interspecific hybridization between cultivated cucumber (Cucumis sativus L., 2n = 14) and the wild relative species C. hystrix Chakr. (2n = 24). IL52 exhibits high resistance to a number of diseases, including downy mildew, powdery mildew, and angular leaf spot. However, the ovary- and fruit-related traits of IL52 have not been thoroughly investigated. Here, we conducted quantitative trait loci (QTL) mapping for 11 traits related to ovary size, fruit size, and flowering time using a previously developed 155 F_7:8 RIL population derived from a cross between CCMC and IL52. In total, 27 QTL associated with the 11 traits were detected, distributed on seven chromosomes. These QTL explained 3.61% to 43.98% of the phenotypic variance. Notably, we identified a major-effect QTL (qOHN4.1) on chromosome 4 associated with the ovary hypanthium neck width and further delimited it into a 114-kb candidate region harboring 13 candidate genes. Furthermore, the QTL qOHN4.1 is co-localized with the QTL detected for ovary length, mature fruit length, and fruit neck length, all residing within the consensus QTL FS4.1, suggesting a plausible pleiotropic effect.

Genetic Mapping and QTL Analysis of Fruit Traits in Melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important horticultural cash crop and its quality traits directly affect consumer choice and market price. These traits are controlled by genetic as well as environmental factors. In this study, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic loci controlling quality traits of melons (i.e., exocarp and pericarp firmness and soluble solid content) based on newly derived whole-genome single nucleotide polymorphism-based cleaved amplified polymorphic sequence (SNP-CAPS) markers. Specifically, SNPs of two melon varieties, M4-5 and M1-15, as revealed by whole-genome sequencing, were converted to the CAPS markers, which were used to construct a genetic linkage map comprising 12 chromosomes with a total length of 1414.88 cM, in the F₂ population of M4-5 and M1-15. The six identified QTLs included: SSC6.1 and SSC11.1 related to soluble solid content; EF12.1 associated with exocarp firmness; and EPF3.1, EPF3.2 and EPF7.1 related to edible pericarp firmness. These genes were located on five chromosomes (3, 6, 7, 11, and 12) in the flanking regions of the CAPS markers. Moreover, the newly developed CAPS markers will be useful in guiding genetic engineering and molecular breeding in melon.

CsTRM5 regulates fruit shape via mediating cell division direction and cell expansion in cucumber

Fruit shape and size are important appearance and yield traits in cucumber, but the underlying genes and their regulatory mechanisms remain poorly understood. Here we identified a mutant with spherical fruits from an Ethyl Methane Sulfonate (EMS)-mutagenized library, named the qiu mutant. Compared with the cylindrical fruit shape in 32X (wild type), the fruit shape in qiu was round due to reduced fruit length and increased fruit diameter. MutMap analysis narrowed the candidate gene in the 6.47 MB range on Chr2, harboring the FS2.1 locus reported previously. A single-nucleotide polymorphism (SNP) (11359603) causing a truncated protein of CsaV3_2G013800, the homolog of tomato fruit shape gene SlTRM5, may underlie the fruit shape variation in the qiu mutant. Knockout of CsTRM5 by the CRISPR-Cas9 system confirmed that CsaV3_2G013800/CsTRM5 was the causal gene responsible for qiu. Sectioning analysis showed that the spherical fruit in qiu resulted mainly from increased and reduced cell division along the transverse and longitudinal directions, respectively. Meanwhile, the repressed cell expansion contributed to the decreased fruit length in qiu. Transcriptome profiling showed that the expression levels of cell-wall-related genes and abscisic acid (ABA) pathway genes were significantly upregulated in qiu. Hormone measurements indicated that ABA content was greatly increased in the qiu mutant. Exogenous ABA application reduced fruit elongation by inhibiting cell expansion in cucumber. Taken together, these data suggest that CsTRM5 regulates fruit shape by affecting cell division direction and cell expansion, and that ABA participates in the CsTRM5-mediated cell expansion during fruit elongation in cucumber.

Natural variation in CRABS CLAW contributes to fruit length divergence in cucumber

Fruit length is a key domestication trait that affects crop yield and appearance. Cucumber (Cucumis sativus) fruits vary from 5 to 60 cm in length. Despite the identification of several regulators and multiple quantitative trait loci (QTLs) underlying fruit length, the natural variation, and molecular mechanisms underlying differences in fruit length are poorly understood. Through map-based cloning, we identified a nonsynonymous polymorphism (G to A) in CRABS CLAW (CsCRC) as underlying the major-effect fruit size/shape QTL FS5.2 in cucumber. The short-fruit allele CsCRCA is a rare allele that has only been found in round-fruited semi-wild Xishuangbanna cucumbers. A near-isogenic line (NIL) homozygous for CsCRCA exhibited a 34∼39% reduction in fruit length. Introducing CsCRCG into this NIL rescued the short-fruit phenotype, and knockdown of CsCRCG resulted in shorter fruit and smaller cells. In natural cucumber populations, CsCRCG expression was positively correlated with fruit length. Further, CsCRCG, but not CsCRCA, targets the downstream auxin-responsive protein gene CsARP1 to regulate its expression. Knockout of CsARP1 produced shorter fruit with smaller cells. Hence, our work suggests that CsCRCG positively regulates fruit elongation through transcriptional activation of CsARP1 and thus enhances cell expansion. Using different CsCRC alleles provides a strategy to manipulate fruit length in cucumber 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.

Morphological and Genetic Diversity of Cucumber (Cucumis sativus L.) Fruit Development

Cucumber (Cucumis sativus L.) fruits, which are eaten at an immature stage of development, can vary extensively in morphological features such as size, shape, waxiness, spines, warts, and flesh thickness. Different types of cucumbers that vary in these morphological traits are preferred throughout the world. Numerous studies in recent years have added greatly to our understanding of cucumber fruit development and have identified a variety of genetic factors leading to extensive diversity. Candidate genes influencing floral organ establishment, cell division and cell cycle regulation, hormone biosynthesis and response, sugar transport, trichome development, and cutin, wax, and pigment biosynthesis have all been identified as factors influencing cucumber fruit morphology. The identified genes demonstrate complex interplay between structural genes, transcription factors, and hormone signaling. Identification of genetic factors controlling these traits will facilitate breeding for desired characteristics to increase productivity, improve shipping, handling, and storage traits, and enhance consumer-desired qualities. The following review examines our current understanding of developmental and genetic factors driving diversity of cucumber fruit morphology.

Sigma factor binding protein 1 (CsSIB1) is a putative candidate of the major-effect QTL dm5.3 for downy mildew resistance in cucumber (Cucumis sativus)

The dm5.3 major-effect QTL in cucumber encodes a homolog of Arabidopsis sigma factor binding protein 1 (CsSIB1). CsSIB1 positively regulates defense responses against downy mildew in cucumber through the salicylic acid (SA) biosynthesis/signaling pathway. Downy mildew (DM) caused by the oomycete pathogen Pseudoperonospora cubensis is an important disease of cucumber and other cucurbits. Our knowledge on molecular mechanisms of DM resistance is still limited. In this study, we reported identification and functional characterization of the candidate gene for the major-effect QTL, dm5.3 for DM resistance originated from PI 197088. The dm5.3 QTL was Modelized through marker-assisted development of near isogenic lines (NILs). NIL-derived segregating populations were used for fine mapping which narrowed the dm5.3 locus down to a 144 kb region. Based on multiple lines of evidence, we show that CsSIB1 (CsGy5G027140) that encodes the VQ motif-containing sigma factor binding protein 1 as the most likely candidate for dm5.3. Local association analysis identified a haplotype consisting of 7 SNPs inside the coding and promoter region of CsSIB1 that was associated with DM resistance. Expression of CsSIB1 was up-regulated with P. cubensis infection. Transcriptome profiling of NILs in response to P. cubensis inoculation revealed key players and associated gene networks in which increased expression of CsSIB1 antagonistically promoted salicylic acid (SA) but suppressed jasmonic acid (JA) biosynthesis/signaling pathways. Our work provides novel insights into the function of CsSIB1/dm5.3 as a disease resistance (R) gene. The roles of sigma factor binding protein genes in pathogen defense in cucumber were also discussed.

Relationship between Key Environmental Factors and the Architecture of Fruit Shape and Size in Near-Isogenic Lines of Cucumber (Cucumis sativus L.)

Fruit shape and size are complex traits influenced by numerous factors, especially genetics and environment factors. To explore the mechanism of fruit shape and size development in cucumber, a pair of near-isogenic lines (NIL) Ln35 and Ln37 were used. The fruit length and diameter, cell length and diameter, and related gene expression were measured. Both the fruit length, diameter, and cell length and diameter showed sigmate curves in the two lines. The cell length and diameter were significantly positively correlated with fruit length and diameter both in two lines. The expression of CsACS2 and CsLNG showed significant positive correlations with fruit length and diameter increment in Ln35, and there was no correlation in Ln37. Furthermore, there were significant positive correlations between fruit size and thermal effectiveness (TE), as well as between fruit size and photosynthetic active radiation (PAR), both in two lines. Two models using logistic regression were formulated to assess the relationships among fruit length and diameter in Ln35 and Ln37, respectively, based on thermal effectiveness and photosynthetic active radiation (TEP). The coefficient R² values of the models were 0.977 and 0.976 in Ln35, and 0.987 and 0.981 in Ln37, respectively. The root mean square error (RMSE) was 12.012 mm and 4.338 mm in Ln35, and 5.17 mm and 7.082 mm in Ln37, respectively, which illustrated the accurate and efficient of these models. These biologically interpreted parameters will provide precision management for monitoring fruit growth and forecasting the time of harvesting under different temperatures and light conditions.

Phenotypic Characterization and Fine Mapping of a Major-Effect Fruit Shape QTL FS5.2 in Cucumber, Cucumis sativus L., with Near-Isogenic Line-Derived Segregating Populations

Cucumber (Cucumis sativus L.) fruit size/shape (FS) is an important yield and quality trait that is quantitatively inherited. Many quantitative trait loci (QTLs) for fruit size/shape have been identified, but very few have been fine-mapped or cloned. In this study, through marker-assisted foreground and background selections, we developed near-isogenic lines (NILs) for a major-effect fruit size/shape QTL FS5.2 in cucumber. Morphological and microscopic characterization of NILs suggests that the allele of fs5.2 from the semi-wild Xishuangbanna (XIS) cucumber (C. s. var. xishuangbannesis) reduces fruit elongation but promotes radial growth resulting in shorter but wider fruit, which seems to be due to reduced cell length, but increased cellular layers. Consistent with this, the NIL carrying the homozygous XIS allele (fs5.2) had lower auxin/IAA contents in both the ovary and the developing fruit. Fine genetic mapping with NIL-derived segregating populations placed FS5.2 into a 95.5 kb region with 15 predicted genes, and a homolog of the Arabidopsis CRABS CLAW (CsCRC) appeared to be the most possible candidate for FS5.2. Transcriptome profiling of NIL fruits at anthesis identified differentially expressed genes enriched in the auxin biosynthesis and signaling pathways, as well as genes involved in cell cycle, division, and cell wall processes. We conclude that the major-effect QTL FS5.2 controls cucumber fruit size/shape through regulating auxin-mediated cell division and expansion for the lateral and longitudinal fruit growth, respectively. The gibberellic acid (GA) signaling pathway also plays a role in FS5.2-mediated fruit elongation.

The CsHEC1-CsOVATE module contributes to fruit neck length variation via modulating auxin biosynthesis in cucumber

Fruit neck is the proximal portion of the fruit with undesirable taste that has detrimental effects on fruit shape and commercial value in cucumber. Despite the dramatic variations in fruit neck length of cucumber germplasms, the genes and regulatory mechanisms underlying fruit neck elongation remain mysterious. In this study, we found that Cucumis sativus HECATE1 (CsHEC1) was highly expressed in fruit neck. Knockout of CsHEC1 resulted in shortened fruit neck and decreased auxin accumulation, whereas overexpression of CsHEC1 displayed the opposite effects, suggesting that CsHEC1 positively regulated fruit neck length by modulating local auxin level. Further analysis showed that CsHEC1 directly bound to the promoter of the auxin biosynthesis gene YUCCA4 (CsYUC4) and activated its expression. Enhanced expression of CsYUC4 resulted in elongated fruit neck and elevated auxin content. Moreover, knockout of CsOVATE resulted in longer fruit neck and higher auxin. Genetic and biochemical data showed that CsOVATE physically interacted with CsHEC1 to antagonize its function by attenuating the CsHEC1-mediated CsYUC4 transcriptional activation. In cucumber germplasms, the expression of CsHEC1 and CsYUC4 positively correlated with fruit neck length, while that of CsOVATE showed a negative correlation. Together, our results revealed a CsHEC1-CsOVATE regulatory module that confers fruit neck length variation via CsYUC4-mediated auxin biosynthesis in cucumber.

Ectopic Expression of CsSUN in Tomato Results in Elongated Fruit Shape via Regulation of Longitudinal Cell Division

Fruit shape, an important agronomic trait of cucumber (Cucumis sativus L.), is tightly controlled by a series of genes such as CsSUN, a homologue of SlSUN that is responsible for the tomato (Solanum lycopersicum) fruit shape via the modulation of cell division. However, the direct genetic evidence about the CsSUN-mediated regulation of fruit shape is still scarce, limiting our mechanistic understanding of the biological functions of CsSUN. Here, we introduced CsSUN into the round-fruited tomato inbred line ‘SN1′ (wild type, WT) via the Agrobacterium tumefaciens-mediated method. The high and constitutive expression of CsSUN was revealed by real-time PCR in all the tested tissues of the transgenic plants, especially in the fruits and ovaries. Phenotypic analyses showed that the ectopic expression of CsSUN increased fruit length while it decreased fruit diameter, thus leading to the enhanced fruit shape index in the transgenic tomato lines relative to the WT. Additionally, the reduction in the seed size and seed-setting rate and the stimulation of seed germination were observed in the CsSUN-expressed tomato. A histological survey demonstrated that the elongated fruits were mainly derived from the significant increasing of the longitudinal cell number, which compensated for the negative effects of decreased cell area in the central columellae. These observations are different from action mode of SlSUN, thus shedding new insights into the SUN-mediated regulation of fruit shape.

Research Progress on the Leaf Morphology, Fruit Development and Plant Architecture of the Cucumber

Cucumber (Cucumis sativus L.) is an annual climbing herb that belongs to the Cucurbitaceae family and is one of the most important economic crops in the world. The breeding of cucumber varieties with excellent agronomic characteristics has gained more attention in recent years. The size and shape of the leaves or fruit and the plant architecture are important agronomic traits that influence crop management and productivity, thus determining the crop yields and consumer preferences. The growth of the plant is precisely regulated by both environmental stimuli and internal signals. Although significant progress has been made in understanding the plant morphological regulation of Arabidopsis, rice, and maize, our understanding of the control mechanisms of the growth and development of cucumber is still limited. This paper reviews the regulation of phytohormones in plant growth and expounds the latest progress in research regarding the genetic regulation pathways in leaf development, fruit size and shape, branching, and plant type in cucumber, so as to provide a theoretical basis for improving cucumber productivity and cultivation efficiency.

Serial-Omics and Molecular Function Study Provide Novel Insight into Cucumber Variety Improvement

Cucumbers are rich in vitamins and minerals. The cucumber has recently become one of China’s main vegetable crops. More specifically, the adjustment of the Chinese agricultural industry’s structure and rapid economic development have resulted in increases in the planting area allocated to Chinese cucumber varieties and in the number of Chinese cucumber varieties. After complete sequencing of the “Chinese long” genome, the transcriptome, proteome, and metabolome were obtained. Cucumber has a small genome and short growing cycle, and these traits are conducive to the application of molecular breeding techniques for improving fruit quality. Here, we review the developments and applications of molecular markers and genetic maps for cucumber breeding and introduce the functions of gene families from the perspective of genomics, including fruit development and quality, hormone response, resistance to abiotic stress, epitomizing the development of other omics, and relationships among functions.

Functional copy number variation of CsSHINE1 is associated with fruit skin netting intensity in cucumber, Cucumis sativus

Fruit skin netting in cucumber (Cucumis sativus) is associated with important fruit quality attributes. Two simply inherited genes H (Heavy netting) and Rs (Russet skin) control skin netting, but their molecular basis is unknown. Here, we reported map-based cloning and functional characterization of the candidate gene for the Rs locus that encodes CsSHINE1 (CsSHN1), an AP2 domain containing ethylene-responsive transcription factor protein. Comparative phenotypic analysis in near-isogenic lines revealed that fruit with netted skin had different epidermal structures from that with smooth skin including thicker cuticles, smaller, palisade-shaped epidermal and sub-epidermal cells with heavily suberized and lignified cell walls, higher peroxidase activities, which suggests multiple functions of CsSHN1 in regulating fruit skin netting and epidermal cell patterning. Among three representative cucumber inbred lines, three haplotypes at three polymorphic sites were identified inside CsSHN1: a functional copy in Gy14 (wild type) with light fruit skin netting, a copy number variant with two tandemly arrayed functional copies in WI7120 with heavy skin netting, and a loss-of-function copy in 9930 with smooth skin. The expression level of CsSHN1 in fruit exocarp of three lines was positively correlated with the skin netting intensity. Comparative analysis between cucumber and melon revealed conserved and divergent genetic mechanisms underlying fruit skin netting/reticulation that may reflect the different selection histories in the two crops. A discussion was made on genetic basis of fruit skin netting in the context of natural and artificial selections of fruit quality-related epidermal features during cucumber breeding.

Promoter variation in a homeobox gene, CpDll, is associated with deeply lobed leaf in Cucurbita pepo L

CpDll, encoding an HD-Zip I transcription factor, positively regulates formation of deeply lobed leaf shape in zucchini, Cucurbita pepo, which is associated with sequence variation in its promoter region. Leaf shape is an important horticultural trait in zucchini (Cucurbita pepo L.). Deeply lobed leaves have potential advantages for high-density planting and hybrid production. However, little is known about the molecular basis of deeply lobed leaf formation in this important vegetable crop. Here, we conducted QTL analysis and fine mapping of the deeply lobed leaf (CpDll) locus using recombinant inbred lines and large F₂ populations developed from crosses between the deeply lobed leaf HM-S2, and entire leaf Jin-GL parental lines. We show that CpDll exhibited incomplete dominance for the deeply lobed leaf shape in HM-S2. Map-based cloning provided evidence that CpCll encodes a type I homeodomain (HD)- and Leu zipper (Zip) element-containing transcription factor. Sequence analysis between HM-S2 and Jin-GL revealed no sequence variations in the coding sequences, whereas a number of variations were identified in the promoter region between them. DUAL-LUC assays revealed significantly stronger promoter activity in HM-S2 than that in Jin-GL. There was also significantly higher expression of CpDll in the leaf base of deeply lobed leaves of HM-S2 compared with entire leaf Jin-GL. Comparative analysis of CpDll gene homologs in nine cucurbit crop species (family Cucurbitaceae) revealed conservation in both structure and function of this gene in regulation of deeply lobed leaf formation. Our work provides new insights into the molecular basis of leaf lobe formation in pumpkin/squash and other cucurbit crops. This work also facilitates marker-assisted selection for leaf shape in zucchini breeding.

CmFSI8/CmOFP13 encoding an OVATE family protein controls fruit shape in melon

Fruit shape is an important quality and yield trait in melon (Cucumis melo). Although some quantitative trait loci for fruit shape have been reported in in this species, the genes responsible and the underlying mechanisms remain poorly understood. Here, we identified and characterized a gene controlling fruit shape from two melon inbred lines, B8 with long-horn fruit and HP22 with flat-round fruit. Genetic analysis suggested that the shape was controlled by a single and incompletely dominant locus, which we designate as CmFSI8/CmOFP13. This gene was finely mapped to a 53.7-kb interval on chromosome 8 based on bulked-segregant analysis sequencing and map-based cloning strategies. CmFSI8/CmOFP13 encodes an OVATE family protein (OFP) and is orthologous to AtOFP1 and SlOFP20. The transcription level of CmFSI8/CmOFP13 in the ovary of HP22 was significantly higher than that in B8, and sequence analysis showed that a 12.5-kb genomic variation with a retrotransposon insertion identified in the promoter was responsible for elevating the expression, and this ultimately caused the differences in fruit shape. Ectopic overexpression of CmFSI8/CmOFP13 in Arabidopsis led to multiple phenotypic changes, including kidney-shaped leaves and shortened siliques. Taken together, our results demonstrate the involvement of an OFP in regulating fruit shape in melon, and our improved understanding of the molecular mechanisms will enable us to better manipulate fruit shape in breeding.

Mapping of genetic loci controlling fruit linked morphological traits of melon using developed CAPS markers

BACKGROUND

Fruit morphology traits are important commercial traits that directly affect the market value. However, studying the genetic basis of these traits in un-explored botanical groups is a fundamental objective for crop genetic improvement through marker-assisted breeding.

METHODS AND RESULTS

In this study, a quantitative trait loci (QTLs) mapping strategy was used for dissecting the genomic regions of fruit linked morphological traits by single nucleotide polymorphism (SNP) based cleaved amplified polymorphism sequence (CAPS) molecular markers. Next-generation sequencing was done for the genomic sequencing of two contrasted melon lines (climacteric and non-climacteric), which revealed 97% and 96% of average coverage over the reference melon genome database, respectively. A total of 57.51% non-synonymous SNPs and 42.49% synonymous SNPs were found, which produced 149 sets of codominant markers with a 24% polymorphism rate. Total 138-F₂ derived plant populations were genotyped for linkage mapping and composite interval mapping based QTL mapping exposed 6 genetic loci, positioned over distinct chromosomes (02, 04, 08, 09, and 12) between the flanking intervals of CAPS markers, which explained an unlinked polygenic architecture in genome. Three minor QTLs of fruit weight (FWt2.1, FWt4.1, FWt9.1), one major QTL of fruit firmness (FrFir8.1), one major QTL of fruit length (FL12.1), and one major QTL of fruit shape (FS12.1) were determined and collectively explained the phenotypic variance from 5.64 to 15.64%. Fruit phenotypic correlation exhibited the significant relationship and principal component analysis also identified the potential variability. Multiple sequence alignments also indicated the significant base-mutations in the detected genetic loci, respectively.

CONCLUSION

In short, our illustrated genetic loci are expected to provide the reference insights for fine QTL mapping and candidate gene(s) mining through molecular genetic breeding approaches aimed at developing the new varieties.

Genetic and Molecular Regulation Mechanisms in the Formation and Development of Vegetable Fruit Shape

Vegetable crops have a long history of cultivation worldwide and rich germplasm resources. With its continuous development and progress, molecular biology technology has been applied to various fields of vegetable crop research. Fruit is an important organ in vegetable crops, and fruit shape can affect the yield and commercialization of vegetables. In nature, fruits show differences in size and shape. Based on fruit shape diversity, the growth direction and coordination mechanism of fruits remain unclear. In this review, we discuss the latest research on fruit shape. In addition, we compare the current theories on the molecular mechanisms that regulate fruit growth, size, and shape in different vegetable families.

Fine mapping and identification of the candidate gene BFS for fruit shape in wax gourd (Benincasa hispida)

Non-synonymous mutations in the BFS gene, which encodes the IQD protein, are responsible for the shape of wax gourd fruits. Fruit shape is an important agronomic trait in wax gourds. Therefore, in this study, we employed bulked segregant analysis (BSA) to identify a candidate gene for fruit shape in wax gourds within F₂ populations derived by crossing GX-71 (long cylindrical fruit, fruit shape index = 4.56) and MY-1 (round fruit, fruit shape index = 1.06) genotypes. According to BSA, the candidate gene is located in the 17.18 Mb region on chromosome 2. Meanwhile, kompetitive allele-specific PCR (KASP) markers were used to reduce it to a 19.6 Kb region. Only one gene was present within the corresponding region of the reference genome, namely Bch02G016830 (designated BFS). Subsequently, BFS was sequenced in six wax gourd varieties with different fruit shapes. Sequence analysis revealed two non-synonymous mutations in the round wax gourd and one non-synonymous mutation in the cylindrical wax gourd. Quantitative real‑time PCR (qRT-PCR) analysis further showed that the expression of BFS in round fruits was significantly higher than in long cylindrical fruits at the ovary formation stage. Therefore, BFS is a candidate gene for determination wax gourd shape. The predicted protein encoded by the BFS gene belongs to the IQ67-domain protein family, which have the structural characteristics of scaffold proteins and coordinate Ca²⁺ CaM signaling from the membrane to the nucleus. Ultimately, two derived cleaved amplified polymorphic sequence (dCAPS) markers were developed to facilitate marker-assisted selection for wax gourds breeding.

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.

Genomic Prediction and the Practical Breeding of 12 Quantitative-Inherited Traits in Cucumber (Cucumis sativus L.)

Genomic prediction is an effective way for predicting complex traits, and it is becoming more essential in horticultural crop breeding. In this study, we applied genomic prediction in the breeding of cucumber plants. Eighty-one cucumber inbred lines were genotyped and 16,662 markers were identified to represent the genetic background of cucumber. Two populations, namely, diallel cross population and North Carolina II population, having 268 combinations in total were constructed from 81 inbred lines. Twelve cucumber commercial traits of these two populations in autumn 2018, spring 2019, and spring 2020 were collected for model training. General combining ability (GCA) models under five-fold cross-validation and cross-population validation were applied to model validation. Finally, the GCA performance of 81 inbred lines was estimated. Our results showed that the predictive ability for 12 traits ranged from 0.38 to 0.95 under the cross-validation strategy and ranged from -0.38 to 0.88 under the cross-population strategy. Besides, GCA models containing non-additive effects had significantly better performance than the pure additive GCA model for most of the investigated traits. Furthermore, there were a relatively higher proportion of additive-by-additive genetic variance components estimated by the full GCA model, especially for lower heritability traits, but the proportion of dominant genetic variance components was relatively small and stable. Our findings concluded that a genomic prediction protocol based on the GCA model theoretical framework can be applied to cucumber breeding, and it can also provide a reference for the single-cross breeding system of other crops.

Long-read genome assembly and genetic architecture of fruit shape in the bottle gourd

The bottle gourd (Lagenaria siceraria, Cucurbitaceae) is an important horticultural crop exhibiting tremendous diversity in fruit shape. The genetic architecture of fruit shape variation in this species remains unknown. We assembled a long-read-based, high-quality reference genome (ZAAS_Lsic_2.0) with a contig N50 value over 390-fold greater than the existing reference genomes. We then focused on dissection of fruit shape using a one-step geometric morphometrics-based functional mapping approach. We identified 11 quantitative trait loci (QTLs) responsible for fruit shape (fsQTLs), reconstructed their visible effects and revealed syntenic relationships of bottle gourd fsQTLs with 12 fsQTLs previously reported in cucumber, melon or watermelon. Homologs of several well-known and newly identified fruit shape genes, including SUN, OFP, AP2 and auxin transporters, were comapped with bottle gourd QTLs.

CsKTN1 for a katanin p60 subunit is associated with the regulation of fruit elongation in cucumber (Cucumis sativus L.)

We identified a short fruit3 (sf3) mutant in cucumber. Map-based cloning revealed that CsKTN1 gene encodes a katanin p60 subunit, which is associated with the regulation of fruit elongation. Fruit length is an important horticultural trait for both fruit yield and quality of cucumber (Cucumis sativus L.). Knowledge on the molecular regulation of fruit elongation in cucumber is very limited. In this study, we identified and characterized a cucumber short fruit3 (sf3) mutant. Histological examination indicated that the shorter fruit in the mutant was due to reduced cell numbers. Genetic analysis revealed that the phenotype of the sf3 mutant was controlled by a single gene with semi-dominant inheritance. By map-based cloning and Arabidopsis genetic transformation, we showed that Sf3 was a homolog of KTN1 (CsKTN1) encoding a katanin p60 subunit. A non-synonymous mutation in the fifth exon of CsKTN1 resulted in an amino acid substitution from Serine in the wild type to Phenylalanine in the sf3 mutant. CsKTN1 expressed in all tissues of both the wild type and the sf3 mutant. However, there was no significant difference in CsKTN1 expression levels between the wild type and the sf3 mutant. The hormone quantitation and RNA-seq analysis suggested that auxin and gibberellin contents are decreased in sf3 by changing the expression levels of genes related with auxin and gibberellin metabolism and signaling. This work helps understand the function of the katanin and the molecular mechanisms of fruit growth regulation in cucumber.

Identification of quantitative trait loci for powdery mildew resistance in highly resistant cucumber (Cucumis sativus L.) using ddRAD-seq analysis

Powdery mildew, caused by Podosphaera xanthii (syn. Sphaerotheca fuliginea ex Fr. Poll.), is one of the most economically important foliar diseases in cucumber (Cucumis sativus L.). Cucumber parental line 'Kyuri Chukanbohon Nou 5 Go', developed from weedy cucumber line CS-PMR1, is highly resistant to powdery mildew and is promising breeding material. We performed quantitative trait locus (QTL) analysis using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in a population from a cross between 'Kyuri Chukanbohon Nou 5 Go' and the Japanese native cultivar 'Kaga-aonaga-fushinari', which is susceptible to powdery mildew. The resistance of the population and its parents was evaluated using leaf disc assays and image analysis. We detected one major QTL on Chr. 5 that was effective at both 20°C and 25°C and one minor QTL on Chr. 1 effective at 20°C. We detected two additional QTLs in subpopulation: one on Chr. 3 effective at 20°C and one on Chr. 5 effective at both 20°C and 25°C in a position different from the major QTL. The resistance alleles at all four QTLs were contributed by 'Kyuri Chukanbohon Nou 5 Go'. The results of this study can be used to develop practical DNA markers tightly linked to genes for powdery mildew resistance.

Quantitative trait loci for horticulturally important traits defining the Sikkim cucumber, Cucumis sativus var. sikkimensis

QTL mapping identified simply inherited genes and quantitative trait loci underlying morphologically characteristic traits of the Sikkim cucumber, which reveals their genetic basis during crop evolution. The data suggest the Sikkim cucumber as an ecotype of cultivated cucumber not worthy of formal taxonomic recognition. The Sikkim cucumber, Cucumis sativus var. sikkimensis, is featured with some morphological traits like black spine, brown fruit with fine and heavy netting, as well as large hollow in mature fruit. Despite its establishment as a botanical variety ~ 150 years ago, and its wide use as an important source of disease resistances in cucumber breeding, little is known about its taxonomic status and genetic basis of those characteristic traits. Here we reported QTL mapping with segregating populations derived from two Sikkim-type inbred lines, WI7088D and WI7120, and identification of 48 QTL underlying phenotypic variation for 18 horticulturally important traits. We found that the fruit spine and skin colors in the two populations were controlled by the previously cloned pleiotropic B (black spine) locus. The fruit netting in WI7088D and WI7120 was controlled by the well-known H (Heavy netting) and a novel Rs (Russet skin) locus, which was delimited to a 271-kb region on Chr5 and ~ 736-kb region on Chr1, respectively. A single major-effect QTL was detected for flowering time in each population (ft1.1 for WI7088D and ft6.2 for WI7120). Fifteen, six and five QTL were identified for fruit size, hollow size and flesh thickness variation in the two populations, respectively. No major structural changes were found between the Sikkim and cultivated cucumbers. Except for the rare allele at the Rs locus, there seem no private QTL/alleles identified from this study supporting the Sikkim cucumber as an ecotype of C. sativus, not worthy of formal taxonomic recognition.

A durum wheat adult plant stripe rust resistance QTL and its relationship with the bread wheat Yr80 locus

KEY MESSAGE

A stripe rust resistance QTL in durum wheat maps near the bread wheat Yr80 locus with the latter reduced to 15 candidate genes. Some wheat adult plant resistance (APR) genes provide partial resistance in the later stages of plant development to rust diseases and are an important component in protecting wheat crops from these fungal pathogens. These genes provide protection in both bread wheat and durum wheat. Here, we have mapped APR to wheat stripe rust, caused by the fungal pathogen Puccinia striiformis f. sp. tritici, in a cross between durum cultivars Stewart and Bansi. Two resistance QTLs derived from the Stewart parent were identified in multi-generational field trials. One QTL is located on chromosome 1BL and maps to the previously identified Yr29/Lr46/Sr58/Pm39 multi-pathogen APR locus. The second locus, located on chromosome 3BL, maps near the recently described bread wheat APR gene, Yr80. Fine mapping in durum and bread wheat families shows that the durum 3BL locus and Yr80 are closely located, with the later APR gene reduced to 15 candidate genes present in the Chinese Spring genome sequence. Distorted segregation of the durum 3BL region was observed with the Stewart locus preferentially transmitted through pollen when compared with the equivalent Bansi region.

Identification of fruit size associated quantitative trait loci featuring SLAF based high-density linkage map of goji berry (Lycium spp.)

BACKGROUND

Goji (Lycium spp., 2n = 24) is a fruit bearing woody plant popular as a superfood for extensive medicinal and nutritional advantages. Fruit size associated attributes are important for evaluating small-fruited goji berry and plant architecture. The domestication traits are regulated quantitatively in crop plants but few studies have attempted on genomic regions corresponding to fruit traits.

RESULTS

In this study, we established high-resolution map using specific locus amplified fragment (SLAF) sequencing for de novo SNPs detection based on 305 F₁ individuals derived from L. chinense and L. barbarum and performed quantitative trait loci (QTL) analysis of fruit size related traits in goji berry. The genetic map contained 3495 SLAF markers on 12 LGs, spanning 1649.03 cM with 0.47 cM average interval. Female and male parents and F₁ individuals` sequencing depth was 111.85-fold and 168.72-fold and 35.80-fold, respectively. The phenotype data were collected for 2 successive years (2018-2019); however, two-year mean data were combined in an extra year (1819). Total 117 QTLs were detected corresponding to multiple traits, of which 78 QTLs in 2 individual years and 36 QTLs in extra year. Six Promising QTLs (qFW10-6.1, qFL10-2.1, qLL10-2.1, qLD10-2.1, qLD12-4.1, qLA10-2.1) were discovered influencing fruit weight, fruit length and leaf related attributes covering an interval ranged from 27.32-71.59 cM on LG10 with peak LOD of 10.48 and 14.6% PVE. Three QTLs targeting fruit sweetness (qFS3-1, qFS5-2) and fruit firmness (qFF10-1) were also identified. Strikingly, various traits QTLs were overlapped on LG10, in particular, qFL10-2.1 was co-located with qLL10-2.1, qLD10-2.1 and qLA10-2.1 among stable QTLs, harbored tightly linked markers, while qLL10-1 was one major QTL with 14.21 highest LOD and 19.3% variance. As LG10 harbored important traits QTLs, we might speculate that it could be hotspot region regulating fruit size and plant architectures.

CONCLUSIONS

This report highlights the extremely saturated linkage map using SLAF-seq and novel loci contributing fruit size-related attributes in goji berry. Our results will shed light on domestication traits and further strengthen molecular and genetic underpinnings of goji berry; moreover, these findings would better facilitate to assemble the reference genome, determining potential candidate genes and marker-assisted breeding.

Identification of QTLs for grain yield and other traits in tropical maize under Striga infestation

Striga is an important biotic factor limiting maize production in sub-Saharan Africa and can cause yield losses as high as 100%. Marker-assisted selection (MAS) approaches hold a great potential for improving Striga resistance but requires identification and use of markers associated with Striga resistance for adequate genetic gains from selection. However, there is no report on the discovery of quantitative trait loci (QTL) for resistance to Striga in maize under artificial field infestation. In the present study, 198 BC₁S₁ families obtained from a cross involving TZEEI 29 (Striga resistant inbred line) and TZEEI 23 (Striga susceptible inbred line) plus the two parental lines were screened under artificial Striga-infested conditions at two Striga-endemic locations in Nigeria in 2018, to identify QTL associated with Striga resistance indicator traits, including grain yield, ears per plant, Striga damage and number of emerged Striga plants. Genetic map was constructed using 1,386 DArTseq markers distributed across the 10 maize chromosomes, covering 2076 cM of the total genome with a mean spacing of 0.11 cM between the markers. Using composite interval mapping (CIM), fourteen QTL were identified for key Striga resistance/tolerance indicator traits: 3 QTL for grain yield, 4 for ears per plant and 7 for Striga damage at 10 weeks after planting (WAP), across environments. Putative candidate genes which encode major transcription factor families WRKY, bHLH, AP2-EREBPs, MYB, and bZIP involved in plant defense signaling were detected for Striga resistance/tolerance indicator traits. The QTL detected in the present study would be useful for rapid transfer of Striga resistance/tolerance genes into Striga susceptible but high yielding maize genotypes using MAS approaches after validation. Further studies on validation of the QTL in different genetic backgrounds and in different environments would help verify their reproducibility and effective use in breeding for Striga resistance/tolerance.

Localization of quantitative trait loci for cucumber fruit shape by a population of chromosome segment substitution lines

Cucumber fruit shape, a significant agronomic trait, is controlled by quantitative trait loci (QTLs). Feasibility of chromosome segment substitution lines (CSSLs) is well demonstrated to map QTLs, especially the minor-effect ones. To detect and identify QTLs with CSSLs can provide new insights into the underlying mechanisms regarding cucumber fruit shape. In the present study, 71 CSSLs were built from a population of backcross progeny (BC4F2) by using RNS7 (a round-fruit cucumber) as the recurrent parent and CNS21 (a long-stick-fruit cucumber) as the donor parent in order to globally detect QTLs for cucumber fruit shape. With the aid of 114 InDel markers covering the whole cucumber genome, 21 QTLs were detected for fruit shape-related traits including ovary length, ovary diameter, ovary shape index, immature fruit length, immature fruit diameter, immature fruit shape index, mature fruit length, mature fruit diameter and mature fruit shape index, and 4 QTLs for other traits including fruit ground and flesh color, and seed size were detected as well. Together our results provide important resources for the subsequent theoretical and applied researches on cucumber fruit shape and other traits.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

The major-effect quantitative trait locus Fnl7.1 encodes a late embryogenesis abundant protein associated with fruit neck length in cucumber

Fruit neck length (FNL) is an important quality trait in cucumber because it directly affects its market value. However, its genetic basis remains largely unknown. We identified a candidate gene for FNL in cucumber using a next-generation sequencing-based bulked segregant analysis in F₂ populations, derived from a cross between Jin5-508 (long necked) and YN (short necked). A quantitative trait locus (QTL) on chromosome 7, Fnl7.1, was identified through a genome-wide comparison of single nucleotide polymorphisms between long and short FNL F₂ pools, and it was confirmed by traditional QTL mapping in multiple environments. Fine genetic mapping, sequences alignment and gene expression analysis revealed that CsFnl7.1 was the most likely candidate Fnl7.1 locus, which encodes a late embryogenesis abundant protein. The increased expression of CsFnl7.1 in long-necked Jin5-508 may be attributed to mutations in the promoter region upstream of the gene body. The function of CsFnl7.1 in FNL control was confirmed by its overexpression in transgenic cucumbers. CsFnl7.1 regulates fruit neck development by modulating cell expansion. Probably, this is achieved through the direct protein-protein interactions between CsFnl7.1 and a dynamin-related protein CsDRP6 and a germin-like protein CsGLP1. Geographical distribution differences of the FNL phenotype were found among the different cucumber types. The East Asian and Eurasian cucumber accessions were highly enriched with the long-necked and short-necked phenotypes, respectively. A further phylogenetic analysis revealed that the Fnl7.1 locus might have originated from India. Thus, these data support that the CsFnl7.1 has an important role in increasing cucumber FNL.

Cucumber Fruit Size and Shape Variations Explored from the Aspects of Morphology, Histology, and Endogenous Hormones

Fruit size and shape are important qualities and yield traits in cucumber (Cucumis sativus L.), but the factors that influence fruit size and shape remain to be explored. In this study, we investigated the dynamic changes of fruit size and shape from the aspects of morphology, cellular levels and endogenous hormones for nine typical cucumber inbred lines. The results show that fruit length had a strong positive correlation to the cell number in the longitudinal section of fruit throughout the four stages of 0, 6, 12, and 30 DAA (days after anthesis). However, the significant negative correlations were found between fruit length and the fruit cell size at 12 and 30 DAA. Furthermore, fruit diameter was positively correlated to the cell number in the cross section at all the investigated fruit growth stages. The indole-3-acetic acid (IAA) content showed significant positive correlations to the fruit length at all fruit growth stages of −6, −3, 0, 3, 6, 9 and 12 DAA, but IAA content and fruit diameter showed significant negative correlations for all the stages except for at −6 DAA. The trans-zeatin riboside (tZR), zeatin (ZT), gibberellic acid (GA₃) and jasmonic acid (JA) content had a positive or negative correlation with fruit length or diameter only at certain stages. Neither fruit length nor diameter had significant correlations to abscisic acid (ABA) content. These results indicate that variations in fruit size and shape of different cucumber inbred lines mainly result from the differences in fruit cell number and endogenous IAA content. The present work is the first to propose cucumber fruit size and shape changes from the combined aspects of morphology, cellular levels, and endogenous hormones.

Origin and domestication of Cucurbitaceae crops: insights from phylogenies, genomics and archaeology

Some of the World's most valuable crops, including watermelon, honey melon, cucumber, squash, zucchini and pumpkin, belong to the family Cucurbitaceae. We review insights on their domestication from new phylogenies, archaeology and genomic studies. Ancestral state estimation on the most complete Cucurbitaceae phylogeny to date suggests that an annual life cycle may have contributed to domestication. Domestication started c. 11 000 years ago in the New World and Asia, and apparently more recently in Africa. Some cucurbit crops were domesticated only once, others multiple times (e.g. melon from different Asian and African populations). Most wild cucurbit fruits are bitter and nonpalatable to humans, and nonbitterness of the pulp apparently was a trait favoured early during domestication, with genomic data showing how bitterness loss was achieved convergently. The genetic pathways underlying lycopene accumulation, red or orange pulp colour, and fruit size and shape are only just beginning to be understood. The study of cucurbit domestication in recent years has benefitted from the increasing integration of archaeological and genomic data with insights from herbarium collections, the most efficient way to understand species' natural geographic ranges and climate adaptations.

Gene regulatory network controlling carpel number variation in cucumber

The WUSCHEL-CLAVATA3 pathway genes play an essential role in shoot apical meristem maintenance and floral organ development, and under intense selection during crop domestication. The carpel number is an important fruit trait that affects fruit shape, size and internal quality in cucumber, but the molecular mechanism remains elusive. Here, we found that CsCLV3 expression was negatively correlated with carpel number in cucumber cultivars. CsCLV3-RNAi led to increased number of petals and carpels, whereas overexpression of CsWUS resulted in more sepals, petals and carpels, suggesting that CsCLV3 and CsWUS function as a negative and a positive regulator for carpel number variation, respectively. Biochemical analyses indicated that CsWUS directly bound to the promoter of CsCLV3 and activated its expression. Overexpression of CsFUL1^A , a FRUITFULL-like MADS-box gene, resulted in more petals and carpels. CsFUL1^A can directly bind to the CsWUS promoter to stimulate its expression. Furthermore, we found that auxin participated in carpel number variation in cucumber through interaction of CsARF14 with CsWUS. Therefore, we have identified a gene regulatory pathway involving CsCLV3, CsWUS, CsFUL1^A and CsARF14 in determining carpel number variation in an important vegetable crop - cucumber.

The MADS-Box Gene CsSHP Participates in Fruit Maturation and Floral Organ Development in Cucumber

Cucumber is an important vegetable crop bearing fleshy pepo fruit harvested immature. Fruits left unpicked in time during summer production, as well as unfavorable environmental conditions during post-harvest shelf, will cause cucumber fruits to turn yellow and ripen, and thus impair the market value. Identification of maturity-related genes is of great agricultural and economic importance for cucumber production. Here, we isolated and characterized a MADS-box gene, Cucumis sativus SHATTERPROOF (CsSHP) in cucumber. Expression analysis indicated that CsSHP was specifically enriched in reproductive organs including stamens and carpels. Ectopic expression of CsSHP was unable to rescue the indehiscence silique phenotype of shp1 shp2 mutant plant in Arabidopsis. Instead, overexpression of CsSHP resulted in early flowering, precocious phenotypes, and capelloid organs in wild-type Arabidopsis. Biochemical analysis indicated that CsSHP directly interacted with cucumber SEPALLATA (SEP) proteins. CsSHP expression increased significantly during the yellowing stage of cucumber ripening, and was induced by exogenous application of abscisic acid (ABA). Therefore, CsSHP may participate in fruit maturation through the ABA pathway and floral organ specification via interaction with CsSEPs to form protein complex in cucumber.

Map-based cloning of qBWT-c12 discovered brassinosteroid-mediated control of organ size in cotton

Assuring fiber yield stability is the primary objective for cotton breeders since the world population is on the rise, and the demand for cotton fiber is increasing every year. Thus, enhancing average cotton boll weight (BWT) could improve seed cotton production, and ultimately to increase cotton fiber yield. This study accomplished the map-based cloning of a novel boll weight regulating locus, qBWT-c12, in cotton. Bulk segregation analysis detected linked markers, aided in the detection of a stable BWT regulating locus, qBWT-c12, on Chr12 in a novel boll size mutant, BS41. Progeny evaluation confined the qBWT-c12 to a 0.89 cM interval between the AD-A12_07 and AD-FM_44 markers in recombinant derived F₃ and F₄ populations. Homology mapping detected a 40 bp insertion-deletion (InDel) site in the AD-FM_44 clone sequence situated +341 downstream of GhBRH1_A12, which showed complete linkage to the BWT phenotype. The suppressed expression of GhBRH1_A12 suggested its putative involvement during early boll development events in BS41. Although brassinosteroid (BR) biosynthesis and signaling pathway genes were up regulated in different tissues, but the organ growth was suppressed leading to dwarf plants, smaller leaves, and de-morphed smaller bolls in BS41. Thus, a disruption in the BR signal cascade is anticipated and could be related to lower GhBRH1_A12 expression in BS41.This study firstly reported the genetic dissection of boll size regulation of G. barbadense in G. hirsutum background using map-based cloning of a BWT regulating locus, qBWT-c12. Moreover, it also emphasized the putative role GhBRH1_A12 in regulating BR homeostasis and its potential to modulate plant growth and boll development in cotton.

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.

Quantitative Trait Loci Mapping and Candidate Gene Analysis of Low Temperature Tolerance in Cucumber Seedlings

Cucumber (Cucumis sativus L.) is an economically important vegetable crop worldwide, but it is sensitive to low temperatures. Cucumber seedlings exposed to long-term low temperature stress (LT), i.e., below 20°C during the day, and 8°C at night, exhibit leaf yellowing, accelerated senescence, and reduced yield, therefore posing a threat to cucumber production. Studying the underlying mechanisms involved in LT tolerance in cucumber seedlings, and developing germplasm with improved LT-tolerance could provide fundamental solutions to the problem. In this study, an F₂ population was generated from two parental lines, "CG104" (LT-tolerant inbred line) and "CG37" (LT-sensitive inbred line), to identify loci that are responsible for LT tolerance in cucumber seedlings. Replicated phenotypic analysis of the F₂-derived F₃ family using a low-temperature injury index (LTII) suggested that the LT tolerance of cucumber seedlings is controlled by multiple genes. A genetic map of 990.8 cM was constructed, with an average interval between markers of 5.2 cM. One quantitative trait loci (QTL) named qLTT5.1 on chromosome 5, and two QTLs named qLTT6.1 and qLTT6.2 on chromosome 6 were detected. Among them, qLTT6.2 accounted for 26.8 and 24.1% of the phenotypic variation in two different experiments. Single-nucleotide polymorphism (SNP) variations within the region of qLTT6.2 were analyzed using two contrasting in silico bulks generated from the cucumber core germplasm. Result showed that 214 SNPs were distributed within the 42-kb interval, containing three candidate genes. Real-time quantitative reverse transcription PCR and sequence analysis suggested that two genes Csa6G445210, an auxin response factor, and Csa6G445230, an ethylene-responsive transmembrane protein, might be candidate genes responsible for LT tolerance in cucumber seedlings. This study furthers the understanding of the molecular mechanism underlying LT tolerance in cucumber seedlings, and provides new markers for molecular breeding.

Mapping quantitative trait loci for yield-related traits and predicting candidate genes for grain weight in maize

Quantitative trait loci (QTLs) mapped in different genetic populations are of great significance for marker-assisted breeding. In this study, an F_2:3 population were developed from the crossing of two maize inbred lines SG-5 and SG-7 and applied to QTL mapping for seven yield-related traits. The seven traits included 100-kernel weight, ear length, ear diameter, cob diameter, kernel row number, ear weight, and grain weight per plant. Based on an ultra-high density linkage map, a total of thirty-three QTLs were detected for the seven studied traits with composite interval mapping (CIM) method, and fifty-four QTLs were indentified with genome-wide composite interval mapping (GCIM) methods. For these QTLs, Fourteen were both detected by CIM and GCIM methods. Besides, eight of the thirty QTLs detected by CIM were identical to those previously mapped using a F₂ population (generating from the same cross as the mapping population in this study), and fifteen were identical to the reported QTLs in other recent studies. For the fifty-four QTLs detected by GCIM, five of them were consistent with the QTLs mapped in the F₂ population of SG-5 × SG-7, and twenty one had been reported in other recent studies. The stable QTLs associated with grain weight were located on maize chromosomes 2, 5, 7, and 9. In addition, differentially expressed genes (DEGs) between SG-5 and SG-7 were obtained from the transcriptomic profiling of grain at different developmental stages and overlaid onto the stable QTLs intervals to predict candidate genes for grain weight in maize. In the physical intervals of confirmed QTLs qKW-7, qEW-9, qEW-10, qGWP-6, qGWP-8, qGWP-10, qGWP-11 and qGWP-12, there were 213 DEGs in total. Finally, eight genes were predicted as candidate genes for grain size/weight. In summary, the stable QTLs would be reliable and the candidate genes predicted would be benefit for maker assisted breeding or cloning.