Next Article in Journal
Variation Analysis of Starch Properties in Tartary Buckwheat and Construction of Near-Infrared Models for Rapid Non-Destructive Detection
Previous Article in Journal
Ethnobotany around the Virovitica Area in NW Slavonia (Continental Croatia)—Record of Rare Edible Use of Fungus Sarccoscypha coccinea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 15
10.3390/plants13152154
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Huarong Fang
†,
Shishi Huang
†,
Ruirui Li
,
Peng Wang
,
Qingwei Jiang
,
Chuan Zhong
,
Yanjuan Yang
and
Wenjin Yu
*
College of Agriculture, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2024, 13(15), 2154; https://doi.org/10.3390/plants13152154
Submission received: 2 July 2024 / Revised: 22 July 2024 / Accepted: 25 July 2024 / Published: 3 August 2024
(This article belongs to the Section Plant Genetics, Genomics and Biotechnology)
Browse Figures
Review Reports Versions Notes

Abstract

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

1. Introduction

Bottle gourd (Lagenaria siceraria L.), an annual trailing herb, is a popular melon crop with economic importance. Bottle gourd is adapted to warm climates and widely cultivated in tropical-to-temperate regions worldwide [1]. In China, it is primarily grown south of the Changjiang River, particularly in Zhejiang and Guangxi, with an annual planting area of about 133.3 km2 [2]. Bottle gourd has diverse applications: young fruits are consumed as vegetables for their fever-reducing and detoxifying properties, while mature fruits are processed into containers and musical instruments, in addition to being used in handicrafts. Additionally, due to its robust root system and strong grafting compatibility and resistance, bottle gourd seedlings are commonly used as rootstocks for cucurbit crops like watermelon, melon, and cucumber [3,4]. As a horticultural crop, bottle gourd’s commercial value is largely derived from its fruit, with fruit size being a critical trait influencing yield, utility, and consumer preference. Therefore, fruit size is a vital focus in bottle gourd breeding and has the attention of breeders.
The development of fleshy fruit involves complex physiological and biochemical changes. Following fertilization, fruit growth initiates with cell division, progresses through concurrent cell division and expansion to establish fruit set, and concludes with cell expansion. Variations in fruit size are influenced by both the number and size of cells [5]. Numerous factors regulate fruit size, such as cell wall metabolism, the cell cycle, and plant hormone signaling pathways. Cell wall metabolism involves enzymes like glycosyltransferase [6], β-galactosidase [7], and cellulose synthase [8]. The cell cycle is controlled by cell-cycle proteins [9] and cyclin-dependent kinases (CDKs) [10], which influence cell number and fruit size. Plant hormones, including auxin, cytokinin (CK), and gibberellin (GA), interact to impact fruit development and size, with regulatory factors linked to hormone synthesis being crucial. The auxin early-response gene family contains auxin/indoleacetic acid (Aux/IAA), Gretchen Hagen3 (GH3), and small auxin up RNA (SAUR) [11]. Aux/IAA interacts with TIR1/AFB, a component of E3 ubiquitin ligase, allowing the release or inhibition of auxin response factors (ARFs). In tomato, silencing SlAux/IAA17 increases cell size and results in larger fruits [12]. ARF9 negatively regulates cell division during early fruit development and reduces fruit size [13]. ARF106, highly expressed during the cell division and expansion stages, positively regulates apple fruit size by modulating auxin signaling. GH3 binds to amino acids to remove free auxin, and its low expression is associated with high auxin concentration during rapid fruit expansion [11]. MdSAUR36 negatively regulates pericarp cell division, affecting apple fruit size [14]. SAUR041 is a candidate regulator for grape fruit size, and is essential for cell expansion during fruit ripening [15]. In loquat, EjSAUR22 responds to auxins and regulates cell size and fruit expansion [16]. Cytokinin, a primary hormone in plant growth and development, promotes the expansion and growth of kiwi fruit cells [17]. Cytokinin oxidase/dehydrogenase 5 (ZjCKX5) inhibits jujube fruit development, and its overexpression leads to smaller fruits. The transcription factors ZjWRKY23 and ZjWRKY40 target ZjCKX5, downregulating its expression and thereby increasing fruit size in red jujube [18]. The overexpression of AtCKX2 in tomatoes reduces endogenous cytokinin levels in fruit tissue, resulting in thinner peel and smaller fruits [19]. GRAS24, which is involved in GA and auxin signaling, affects tomato fruit development by inhibiting cell division and expansion [20]. The peak period of GA content coincides with the cell division and expansion phases in fruit, highlighting GA’s crucial role in regulating fruit size [21]. Ethylene, a well-known inhibitor of certain tissue and cell growth, is synthesized by 1-aminocyclopropane-1-carboxylate synthase (ACS2). Mutants of acs2 in cucumber and melon produce less ethylene, leading to reduced cell division and ultimately affecting fruit development and size [22,23]. S1DREB3 negatively regulates abscisic acid (ABA) response, and its overexpression alters tomato fruit size [24].
Several transcription factors have been identified as key regulators of fruit size, including members of the ethylene response factor (ERF), AP2, MADS-box, NAC, WOX, YABBY, bHLH, and MYB families. The APETALA2/ethylene responsive factor (AP2/ERF) named PavRAV2 directly represses the expression of PavKLUH, controlling sweet cherry fruit size [25]. In apple, the overexpression of microRNA172 inhibits AP2 transcription, leading to cell shrinkage and a significant reduction in fruit size [26]. Similarly, silencing MdMADS8 or MdMADS9 significantly reduces cell and fruit size [27]. In tomato, the overexpression of the SlNAC transcription factor activated by AP3/PI 1 (SlNAP1) negatively regulates fruit size and weight [28]. Analysis of the SlLAM1 gene, a member of the WOX family, shows that its loss of function can lead to changes in fruit size [29]. Mutation in YABBY causes its depressed expression and promotes polyventricle formation, thereby increasing fruit size [30]. RSL4, a bHLH transcription factor, transmits auxin signals to stimulate cell growth [31]. Three R2R3-MYB transcription factors—SlFSB1, ZmMYB40, and ZmMYB95—interact to regulate the differentiation and expansion of fruit cells [32]. In grape, VvCEB1 is specifically expressed in fruit tissues during the expansion phase and influences the expression of genes related to cell expansion, including those involved in auxin metabolism and signal transduction [33]. In citrus, CsMYB77 negatively regulates fruit ripening and size by modulating abscisic acid and auxin signaling pathways [34].
Fruit size is typically a quantitative trait measured by indicators such as length (L), width (W), L/W ratio, circumference, volume, and weight [35]. Functional genes regulating fruit size have been extensively studied in tomato, including CNR/FW2.2, KLUH/FW3.3, SUN, OVATE, LOCULE NUMBER (LC), WUSCHEL (WUS), FASCIATED (FAS), CLAVATA (CLV), CRC, GRAS2, and so on. Fruit Weight 2.2 (FW2.2) and Fruit Weight 3.2 (FW3.2) were among the first cloned genes involved in regulating cell division and cell number. FW2.2, a member of the cell number regulator (CNR) family, affects fruit growth via intercellular transport mechanisms [36,37]. FW3.2 encodes the KLUH gene of the CYP78A subfamily of P450 enzymes, promoting an increase in cell number in fruit peel tissue and influencing fruit size [38]. SUN and OVATE are pivotal in fruit elongation because they affect cell division [39]. LC (coding for SlWUS) and FAS (coding for SlCLV3) loci are key factors in fruit size variation, regulating size by influencing the number of locules. ENO, a member of the AP2/ERF superfamily, works synergistically with mutations in SlWUS and SlCLV3 to promote cell proliferation by regulating flower meristem activity, thus playing a significant role in increasing fruit size [40]. CRABS CLAW (CRC), belonging to YABBY, is crucial for carpel development; SlCRCa suppresses cell division by regulating related genes and inhibits cell expansion by regulating expansion protein genes and the GA pathway, thereby negatively regulating tomato fruit size [41]. The FAS locus inversions at YABBY and SlCLV3 reduce gene expression and locule number [42]. Silencing GRAS2 inhibits ovary growth and cell expansion, resulting in smaller fruits and lower weight [43]. In addition, cell size regulator (CSR), identified as a key determinant of tomato fruit weight, has been cloned and characterized, and the CSR-D allele increases fruit weight primarily by enlarging pericarp cell size rather than increasing the number of cell layers [44]. SlPZF1, encoding a zinc finger protein family member (C2H2), is preferentially expressed in the exocarp during tomato fruit development, and can control fruit size by influencing pericarp cell size and mediating the cell-cycle regulatory factors interacting with it [45].
Many advances have been made in identifying genes that regulate fruit size in melon crops. Weng et al. [46] identified 12 quantitative trait loci (QTLs) for fruit size (FS1.1, FS1.2, FS2.1, FS2.2, FS3.1, FS3.2, FS3.3, FS4.1, FS5.1, FS6.1, FS6.2, and FS7.1). In cucumber, the CsSUN at the FS1.2 locus is a homolog of the tomato fruit-shape gene SUN [47]. CsACS is a candidate gene affecting slender cucumber fruit development [48]. CsFUL1 negatively regulates fruit length by influencing the cell division and expansion mediated by CsSUP and regulating the expression of CsPIN1 and CsPIN7 [49]. SF2, encoding histone deacetylase, targets genes related to hormone synthesis and signal response pathways, and influences cell division and fruit length by regulating the balance of ethylene, cytokinin, and polyamine [50]. Short-fruit 1 (sf1) is a single recessive gene that likely regulates fruit length through various hormone biosynthesis and signal transduction pathways [51]. In melon, genes such as SUN, OVATE, LC, FAS, CNR/FW2.2, and SlKLUH/FW3.3 have also been identified to control fruit size [52]. Compared to other Cucurbitaceae crops, research on bottle gourd fruit shape has only emerged in recent years. Xu et al. [53] assembled a high-quality reference genome ‘ZAAS_Lsic_2.0′ for bottle gourd and located a dominant QTL controlling fruit shape on chromosome 6. Zhang et al. [54] conducted transcriptome analysis on two gourd species with diverse fruit sizes, revealing significant differences in the expression of genes related to cell wall metabolism, the cell cycle, and phytohormones, with ERF family transcription factors being the most abundant. However, to date, no specific genes related to fruit size in bottle gourd have been cloned.
In this study, bulked segregant analysis (BSA) combined with transcriptome sequencing (RNA-seq) was utilized to identify key genes regulating bottle gourd fruit size. A small-fruited gourd (H28) was used as the female parent and a large-fruited gourd (H17) as the male parent to measure fruit size at different growth stages. This allowed for the analysis of dynamic changes in fruit size and differences in cell number and size between the parental pericarps through microscopic observation. Then, we initially localized an interval regulating gourd fruit size, expanded the F2 population, and developed single-nucleotide polymorphism (SNP) and Insertion–Deletion (InDel) markers within the candidate interval to refine the mapping range. Through BSA-based mapping and RNA-seq, combined with gene annotation, six genes were finally predicted as candidate genes. This is not only helpful for identifying genes regulating fruit size in bottle gourd, but also lays the foundation for elucidating the genetic mechanism underlying this trait. Furthermore, it provides valuable theoretical and practical insight for the molecular-assisted breeding of bottle gourd, accelerating the development of new cultivars with optimized fruit size.

2. Results

2.1. Phenotypic Evaluation of Fruit Size in Bottle Gourd

The fruit size of the parent plants in this study exhibited visible variations (Figure 1), which were analyzed by measuring the transverse diameter (maximum width of the fruit) and longitudinal diameter (length of the fruit from top to bottom) of the parent and F1 fruits. The results demonstrated significant differences in fruit dimensions between the female H28 and male H17 plants, with H28 producing smaller fruits than H17 (Table 1). Specifically, H28 fruits had a width of 48.47 ± 3.73 mm and a length of 67.99 ± 3.41 mm, categorizing them as small fruits. In contrast, H17 fruits had a width of 163.18 ± 9.63 mm and a length of 295.76 ± 49.23 mm, categorizing them as large fruits. The F1 fruits measured 121.32 ± 7.60 mm in width and 176.02 ± 32.46 mm in length, with mean values between the parents favoring the large-fruited H17.

2.2. Dynamic Changes in Fruit Development Process in Bottle Gourd

The fruit widths and lengths of H28, H17, and their F1 progeny were measured every three days from 3 days before pollination (−3 DAP) to 30 days after pollination (DAP). Measurements were taken over twelve periods, with ten fruits from different plants measured at each period to calculate the average values. These averages were used to plot the parental growth curves and analyze the dynamic changes in the growth processes of the parental fruits (Figure 2). As the fruits developed, both the fruit widths and lengths of H28, H17, and F1 gradually increased, with the differences between the parents becoming increasingly significant. Starting from −3 DAP, the fruit sizes of H28 and H17 exhibited noticeable differences, with H17 consistently having larger widths and lengths than H28. The size of F1 was between the two parents, but gradually tended towards H17. The most substantial growth in H28′s fruit width and length occurred at 0 DAP-9 DAP, during which the growth rate was at its peak. The average width of H28 fruits increased from 8.39 mm to 42.16 mm, while the average length grew from 16.62 mm to 62.10 mm. From 9 DAP to 15 DAP, the growth rate slowed down, and the fruit dimensions stabilized by 15 DAP. At maturity (30 DAP), the average fruit width and length of H28 were 48.04 mm and 67.54 mm, respectively. Meanwhile, the fruit size of H17 exhibited rapid growth from 0 DAP to 15 DAP, with the average width increasing from 12.00 mm to 146.00 mm and the average length from 30.59 mm to 230.99 mm. From 15 DAP to 24 DAP, the growth rate slowed, and the fruit width and length stabilized after 24 DAP. At maturity (30 DAP), the average width and length of H17 fruits were 168.37 mm and 280.89 mm, respectively. Consequently, the difference in fruit width and length between H28 and H17 reached its maximum at 30 DAP.

2.3. Paraffin Section Observation of Fruit Skin in Bottle Gourd

The growth change curves of the gourd fruits revealed significant differences in fruit width and length between H28 and H17. To further investigate these differences at the cellular level, bottle gourd flesh samples were collected at −3, 0, 3, 6, 12, 18, 24, and 30 DAP. These samples were prepared as paraffin slices in transverse and longitudinal sections to analyze the cell number and area at eight different developmental stages. As the gourd fruits developed and expanded, the cell area in both the transverse and longitudinal sections gradually increased (Figure 3A, Supplementary Figure S1). In the transverse and longitudinal sections of the parental fruit skin, the cell number in the large-fruited H17 was significantly higher than in the small-fruited H28 from −3 DAP onward, while the cell area in H17 was significantly smaller than in H28 from 0 DAP (Figure 3B, Supplementary Figure S2). These results indicate that both cell number and cell area are critical factors contributing to the differences in fruit size between H28 and H17, with cell number having an earlier impact on fruit size.

2.4. Mapping of Candidate Genes Related to Fruit Size in Bottle Gourd

Through BSA analysis, a total of 43.98 Gb of clean data was obtained from the parent and extreme pools, with each sample achieving a Q30 percentage above 91% and a GC content over 33%. The comparison of clean reads from the parent and F2 samples to the reference genome ‘ZAAS_Lsic_2.0′ showed an average alignment efficiency of 98.86%, an average coverage depth of 30×, and a genome coverage of 99.29% (with at least one base covered). These results confirm the sequencing data’s high quality and suitability for subsequent mutation detection and correlation analysis. A total of 771,950 SNPs were obtained from the four pools, including 294,636 high-quality SNPs used to calculate the SNP index between the two extreme pools. Preliminary BSA-seq mapping results, using the Euclidean Distance (ED) and ΔSNP-index algorithms, located genes regulating bottle gourd fruit size within the 13,950,000 bp to 19,960,000 bp interval on Chr9, spanning 6.01 Mb (Figure 4A). This region contained 642 genes, including 180 non-synonymous genes.
To further narrow the candidate range, four KASP markers were developed every 1–2 Mb within the initial region, and 500 F2 plants were genotyped. Using phenotypic analysis, the fruit size interval was refined to between markers KS16.5 (16,536,349 bp) and KS18.5 (18,589,323 bp), covering a length of 2.05 Mb. Subsequently, the F2 population was expanded to 3000 plants, and several InDel markers (with a base number difference ≥3 bp) were further developed to screen recombinant plants through genotypic–phenotypic analysis. Finally, the candidate interval for fruit size was positioned between markers IS17.74 (17,747,353 bp) and IS18.18 (18,185,825 bp), with five and two recombinant plants on the left and right sides, respectively (Figure 4B). The refined region spans 0.44 Mb and contains 44 genes, 10 of which exhibit non-synonymous mutations in both parents.

2.5. Comparison of H28 and H17 Fruit Transcriptomes

To determine the reason for the fruit size difference between H28 and H17, we collected the fruits at the period of rapid fruit expansion (6 DAP) and compared their transcriptomes. A total of 22.74 GB of clean data was obtained (an average of 5.68 Gb per sample), with each sample achieving a Q30 percentage over 93% (Table 2), indicating high-quality sequencing and suitability for further analysis. The clean reads from each sample were compared with the ‘ZAAS_Lsic_2.0’ reference genome, identifying 1,166 differentially expressed genes (DEGs).
Gene expression analysis using DESeq identified DEGs between H28 and H17 based on two criteria: an absolute value of log2FoldChange >1 and a significance p-value < 0.05. The results revealed a total of 2433 DEGs, with 1632 genes up-regulated and 801 genes down-regulated in the small-fruited H28 (Figure 5A,B). To elucidate the main functions of these DEGs, Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed. The GO analysis showed that the DEGs were enriched in cellular components and biological processes, with a significant number involved in biological processes such as cell wall development and plant post-embryonic morphogenesis (Figure 5C, D). The KEGG analysis indicated that DEGs were enriched in metabolic pathways, including sugar metabolism, phytohormone signaling, photosynthesis, and MAPK signaling (Figure 5E).

2.6. Candidate Gene Prediction of Fruit Size

Based on the resequencing results, the sequences of 44 genes within the 0.44 Mb candidate interval were initially analyzed. Ten genes (HG_GLEAN_10001518, HG_GLEAN_10001525, HG_GLEAN_10001529, HG_GLEAN_10001543, HG_GLEAN_10001546, HG_GLEAN_10001548, HG_GLEAN_10001549, HG_GLEAN_10001550, HG_GLEAN_10001556, and HG_GLEAN_10001558) were found to have non-synonymous mutations in both parents. Through RNA-seq, the expression of all candidate genes within the mapped interval was analyzed, revealing that four genes (HG_GLEAN_10001522, HG_GLEAN_10001544, HG_GLEAN_10001548, and HG_GLEAN_10001561) were differentially expressed between the parents. In total, 13 genes had either non-synonymous mutations or differential expression (Supplementary Table S1), with HG_GLEAN_10001548 exhibiting both. With the use of functional annotation, six genes were finally predicted as candidate regulators of fruit size in bottle gourd (Table 3).

3. Discussion

Fruit size in bottle gourd is an important quality trait and a key indicator of yield and economic value. Understanding the dynamic changes in fruit development at the physiological level is essential for studying fruit size differences in bottle gourd. In this study, the size of large-fruited H17 and small-fruited H28 was measured and analyzed at different developmental stages. The results indicated that the growth pattern of bottle gourd fruits follows a single ‘S’ curve: initially increasing and then flattening. Fruit size differences between the parents were evident at the ovary stage (−3 DAP). The fruit width and length continued to increase, with the differences peaking at 30 DAP, consistent with Yan et al.’s findings on fruit size variation in wax gourd [55]. In addition, many studies have shown that cell division and expansion determine the number and size of fruit cells, influencing overall fruit size. Studies in blueberries [56] and loquat [57] suggest that cell size affects final fruit size, while cell number is a more crucial factor in horticultural crops like apple [58], pear [59], sweet cherry [60], and plum [61]. The cytological observation of paraffin sections revealed that large-fruited H17 had significantly more cells than small-fruited H28 at the ovary stage (−3 DAP), consistent with the initial size differences observed. Significant differences in cell area between the parents were also observed from 0 DAP. Previous studies indicate that post-fertilization fruit growth starts with cell division, followed by continuous cell division and expansion to promote fruit formation. Our study confirms these findings. Importantly, it was observed that H17 fruit enlargement is influenced by both cell number and cell area, with cell number being the more critical factor. Therefore, cell division is likely the primary contributor to fruit size differences in bottle gourd, suggesting it is a key target for future breeding efforts.
In horticultural plant research, exploring the genes and molecular mechanisms regulating fruit size is a hot topic. However, little is known about the genes influencing bottle gourd fruit size. In this study, BSA-seq was employed to map the candidate genes within the range of 17,747,353 bp to 18,185,825 bp on Chr9, covering a physical distance of 0.44 Mb. This interval contained 44 genes, 10 of which had non-synonymous mutations. Subsequently, RNA-seq identified four DEGs within this interval. With the addition of functional annotation, six genes were predicted to regulate fruit size in bottle gourd. Among these, HG_GLEAN_10001548, which has both non-synonymous mutations and upregulated expression, encodes a protein with a pentapeptide repeat (PPR) sequence. PPR proteins are crucial for mitochondrial function and nuclear development, acting as essential RNA-binding proteins in plants and influencing the expression of organelle mRNA transcripts [62]. ClaPPRs have been implicated in regulating fruit development and ripening in watermelon [63]. The QTL fw3.2 in tomato, which controls fruit weight, has been linked with seven putative genes, among which ORF4 encodes a protein that is highly consistent with Arabidopsis PNM1, belonging to the PPR family [64,65]. Similarly, the Ca12g10030 (PPR) site in pepper is associated with fruit weight, with the loss of PPR function leading to fruit development defects and seed sterility [66,67,68]. These findings emphasize the pivotal role of PPR proteins in controlling fruit size and weight. Additionally, HG_GLEAN_10001518 encodes the transcription factor DIVARICATA subtype X1, belonging to the MYB family, that is known to regulate the uneven distribution of auxin signals and participate in ABA signal transduction [69]. HG_GLEAN_10001525 encodes the cell division control protein 6 homology B-like protein, which is involved in DNA replication and controlled by cyclin-dependent kinase (CDK) phosphorylation [70]. HG_GLEAN_10001544 is annotated as zinc transporter 6, which is a member of the ZIP protein family and is involved in the uptake, transport, and distribution of metal ions, influencing fruit growth and development [71,72,73,74]. HG_GLEAN_10001556 and HG_GLEAN_10001558 both encode UDP-glycosyltransferase 87A1-like (UGT) proteins, which are involved in hormone regulation and secondary metabolism, impacting cell wall metabolism and cell size increase. The overexpression of GSA1, which encodes UGT, can increase grain size in rice, while Arabidopsis UGT75D1 negatively regulates epidermal cell growth in cotyledons [75,76]. Cytokinins primarily exist as glycosides in plants, and the overexpression of rice UGT Os6 in Arabidopsis thaliana significantly increases cytokinin glycoside content. Although these predicted genes in the mapped interval are promising, further studies are required to pinpoint the key gene regulating fruit size in bottle gourd.
In this study, GO analysis revealed that DEGs are mainly enriched in processes such as cell wall development and plant postembryonic morphogenesis. KEGG analysis indicated enrichment in pathways related to sugar metabolism, phytohormone signaling, and MAPK signaling, aligning with previous research findings [77]. It has been shown that cell wall development significantly influences cell division and expansion during fruit development. Representative cell wall proteins and modifications can regulate cell size during plant growth and fruit development [78,79]. Since the cell wall is primarily composed of polysaccharides, the sugar metabolism pathway likely affects cell wall biosynthesis, thereby influencing cell division and expansion [80]. Hormones also play a critical role in regulating fruit size and development, often acting in concert. Several functional genes and transcriptional regulators, including MYB and bZIP, which are related to the candidate genes in this study, have been implicated in the hormonal regulation of fruit size. MYB is involved in auxin biosynthesis [81]. For instance, the exogenous application of GA3 to grape berries upregulated vrax2 expression (a member of the R2R3-MYB family), and the overexpression of vrax2 increased fruit size [82]. Similarly, bZIP controls cell elongation by regulating the biosynthesis and transduction of auxin and GA [83,84]. Therefore, future research on bottle gourd fruit size should focus on the key pathways related to cell wall synthesis and metabolism, as well as plant hormones, to identify genes or transcription factors controlling fruit size and explore the underlying regulatory mechanisms.

4. Conclusions

We used BSA-seq combined with RNA-seq to identify key genes regulating the fruit size of bottle gourd. The fine mapping of the dominant locus revealed that it was on Chr9 in a 0.44 Mb interval containing 44 genes. Six genes within the localized region were predicted to be candidate regulators associated with bottle gourd fruit size through the analysis of sequence variation and transcriptional differences, together with gene annotation. Our findings will aid the cloning of genes controlling bottle gourd fruit size, and provide information about the potential molecular mechanisms of the candidate genes regulating fruit size from a transcriptomic perspective.

5. Materials and Methods

5.1. Plant Materials and Phenotypic Evaluation of Fruit Size

Bottle gourd self-inbred lines H28 and H17, preserved by the Vegetable Research Group at the Agricultural College of Guangxi University, were chosen as the female (P1) and male (P2) parents, respectively. H28 produces small fruits with an average width of 48.47 mm and a length of 67.99 mm, whereas H17 produces large fruits with an average width of 163.18 mm and a length of 295.76 mm. Both parent lines exhibit a cucurbit shape, but with significant differences in fruit size (Figure 1). H28 and H17 were crossed to obtain the F1 generation, which was subsequently self-crossed to conduct an F2 segregating population. All plants were cultivated in the experimental field of Guangxi University in 2021.
Fruits at full maturity (40 DAP) were selected, and their transverse diameter (maximum width of the fruit) and longitudinal diameter (length of the fruit from top to bottom) were measured using an electronic vernier caliper with an accuracy of 0.01 mm. These measurements were used to evaluate the mature fruit size of H28 and H17, with three fruits measured per plant to calculate the average value. Additionally, parental fruits were collected at twelve developmental stages: −3, 0, 3, 6, 9, 12, 15, 18, 21, 24, 27, and 30 days after pollination (DAP). The width and length of the parent fruits were measured at each stage to compare dynamic size changes during fruit development, with three fruits measured to obtain an average in each period.

5.2. Paraffin Section and Cytological Observation

The parent fruits were harvested at eight representative developmental stages (−3, 0, 3, 6, 12, 18, 24, and 30 DAP). The lower and middle parts of the bottle gourd fruits were selected, the skin was scraped off, and the flesh tissue was cut into approximately 2 cm3 cubes for transverse and longitudinal paraffin sectioning. The paraffin sectioning process involved fixing and dehydrating the pre-treated samples, embedding them in paraffin, and sectioning them using a microtome (RM2016 Leica, Wetzlar, Germany). The prepared sections were then dewaxed and stained. Observations were made using a light microscope (Nikon Eclipse E100, Tokyo, Japan), and panoramic scanning software (Slide Viewer 2.5) was used to select appropriate views and take pictures. Three views per period were randomly selected to count the cell numbers, and thirty cells per view with complete and clear boundaries were selected to calculate the cell area. The data were collated and analyzed using ImageJ 2023, Excel 2021, and SPSS 18.

5.3. BSA-Seq Mapping Approach

Thirty plants with the largest fruits and thirty plants with the smallest fruits were selected from five hundred F2 plants. Young leaves were collected from these sixty plants to construct extreme mixed pools. Association analysis was conducted on the two F2 pools and the parental pools using the ‘ZAAS_Lsic_2.0’ genome as the reference. Pooled DNA samples were prepared for library construction and sequenced on the Illumina HiSeq™PE150 platform (San Diego, CA, USA). The raw reads obtained from high-throughput sequencing were analyzed and converted to sequencing reads after base calling. These reads were then filtered to obtain clean reads, ensuring quality for subsequent analysis. The clean reads were compared with the reference genome for mutation detection. Finally, the Euclidean Distance (ED) and ΔSNP-index methods identified the regions associated with the target traits [85,86].

5.4. Fine Mapping

To narrow the preliminary region, Kompetitive Allele-Specific PCR (KASP) markers were designed for each 1–2 Mb distance within the candidate interval based on BSA-seq data. PCR amplification was prepared following the manufacturer’s instructions (LGC Genomics, Shanghai, China). The PCR reaction volume was 3 μL, consisting of 1.0 μL of DNA (8–15 ng μL−1), 1.5 μL of 2× master mix, and 0.5 μL of primer mix. The amplification used landing PCR with the following conditions: heat treatment at 95 °C for 15 min; denaturation at 95 °C for 20 s; annealing and extension between 65 and 55 °C for 1 min, 10 cycles, reducing 1.0 °C each cycle; denaturation at 95 °C for 20 s; annealing and extension at 57 °C for 1 min, 26 cycles; then holding at 4 °C in the dark. After amplification, fluorescence scanning and genotyping were performed. Four pairs of KASP markers were designed to identify recombinant plants from the 500 F2 plants. Subsequently, the population was expanded to 3000 plants, and new InDel markers within the mapped interval were developed to screen recombinant plants and determine their genotypes. Finally, the most likely candidate region was inferred using genotype–phenotype analysis. Primers are listed in Supplementary Table S2.

5.5. RNA Sequencing (RNA-Seq) Analysis

Fruit flesh from both parents at the rapid fruit expansion stage (6 DAP) was selected and sent to Paisonol Company (Shanghai, China) for RNA-seq and bioinformatics analysis. Total RNA was extracted from the tissue samples, and its purity, quantity, and integrity were assessed using a NanoDrop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, Waltham, MA, USA). Library construction and high-throughput sequencing were performed on the Illumina HiSeq2500 platform (San Diego, CA, USA). Clean data were filtered from raw data using SOAPnuke v1.5.6 and compared with the reference genome ‘ZAAS_Lsic_2.0’ using HISAT2 v2.0.4 and Bowtie2 v2.2.5. DESeq2 software was used to analyze the differentially expressed genes (DEGs) related to fruit size and enrich the GO [87] and KEGG [88] pathways in the parents.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13152154/s1, Figure S1. Longitudinal section of H17 and H28 fruit at different developmental stages. (A) H17: -3 DAP, (B) H28: -3 DAP, (C) H17: 0 DAP, (D) H28: 0 DAP, (E) H17: 3 DAP, (F) H28: 3 DAP, (G) H17: 6 DAP, (H) H28: 6 DAP, (I) H17: 12 DAP, (J) H28: 12 DAP, (K) H17: 18 DAP, (L) H28: 18 DAP, (M) H17: 24 DAP, (N) H28: 24 DAP, (O) H17: 30 DAP, (P) H28: 30 DAP; Bar = 50 μm. Figure S2. Comparison of the average cell area and number in longitudinal sections of H17 and H28 fruit at different periods. (A) Statistical diagram of cell number. (B) Statistical diagram of cell area. Table S1. Nonsynonymous mutations, expression mode and functional annotations of all genes in the BSA-seq mapped region. Table S2. Primer sequences for polymorphic InDel markers.

Author Contributions

Conceptualization, H.F. and S.H.; writing—original draft preparation, H.F.; writing—review and editing, P.W., C.Z. and W.Y.; methodology, S.H.; investigation, S.H., R.L. and Q.J.; data curation, H.F. and S.H.; project administration, W.Y.; and funding acquisition, W.Y. and Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Guangxi Natural Science Foundation (2023GXNSFAA026508) (2024GXNSFAA010492).

Data Availability Statement

The data presented in this study are available in this article and in the Supplementary Materials.

Acknowledgments

We are grateful to Guangxi University.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Erickson, D.L.; Smith, B.D.; Clarke, A.C.; Sandweiss, D.H.; Tuross, N. An Asian origin for a 10,000-year-old domesticated plant in the Americas. Proc. Natl. Acad. Sci. USA 2005, 102, 18315–18320. [Google Scholar] [CrossRef] [PubMed]
  2. Li, Y.W.; Shi, J.S.; Wang, B.G.; Wang, Y.; Wu, X.Y.; Lu, Z.F.; Li, G.J.; Wu, X.H. Phenotypic identification and genetic diversity analysis of local germplasm resources of bottle gourd in Zhejiang province. J. Plant Genet. Resour. 2020, 21, 1135–1147. [Google Scholar] [CrossRef]
  3. Pang, Z.L.; Mo, X.P.; Luo, H.; Sun, X.W.; Deng, D.C.; Yan, Q.P.; Sun, P.; Yin, B.Y.; Yang, L.H.; Dong, J.N.; et al. Breeding of gourd type rootstock cultivar ‘Xuefengqiangzhen No.1’. Chin. Melon 2020, 33, 66–68. [Google Scholar] [CrossRef]
  4. Tan, M.M.; Zhang, Y.; Yan, L.X.; Deng, H.J.; Li, H.L.; He, Z.Q.; Xu, Z.F. Effect of rootstock grafting on growth characteristics of Muskmelon with thin skin. Agric. Technol. 2021, 41, 64–66. [Google Scholar] [CrossRef]
  5. He, C.Y.; Wang, L.; Yan, L.X.; Li, Q.R.; Yong, B.; Zhu, W.W. Evolutionary developmental mechanisms underlying the origin and diversification of the fruits. Sci. China Life Sci. 2019, 49, 301–319. [Google Scholar] [CrossRef]
  6. Scheible, W.R.; Pauly, M. Glycosyltransferases and cell wall biosynthesis: Novel players and insights. Curr. Opin. Plant Biol. 2004, 7, 285–295. [Google Scholar] [CrossRef] [PubMed]
  7. Roach, M.J.; Mokshina, N.Y.; Badhan, A.; Snegireva, A.V.; Hobson, N.; Deyholos, M.K.; Gorshkova, T.A. Development of cellulosic secondary walls in flax fibers requires β-galactosidase. Plant Physiol. 2011, 156, 1351–1363. [Google Scholar] [CrossRef] [PubMed]
  8. Hu, H.; Zhang, R.; Tao, Z.; Li, X.; Li, Y.; Huang, J.; Li, X.; Han, X.; Feng, S.; Zhang, G.; et al. Cellulose synthase mutants distinctively affect cell growth and cell wall integrity for plant biomass production in Arabidopsis. Plant Cell Physiol. 2018, 59, 1144–1157. [Google Scholar] [CrossRef]
  9. Mathieu-Rivet, E.; Gevaudant, F.; Sicard, A.; Salar, S.; Do, P.T.; Mouras, A.; Fernie, A.R.; Gibon, Y.; Rothan, C.; Chevalier, C.; et al. Functional analysis of the anaphase promoting complex activator CCS52A highlights the crucial role of endo-reduplication for fruit growth in tomato. Plant. J. 2010, 62, 727–741. [Google Scholar] [CrossRef]
  10. Jiao, Q.Q.; Cui, M.; Chen, X.S.; Zhang, Z.Y.; Wang, N.; Jiang, Y.M.; Mao, Z.Q.; Shu, J.; Shen, X.; Jiang, S.H.; et al. MdCKS1 expression involved in fruit size in apple cultivar ‘Fuji’ and ‘Ralls’. Russ. J. Plant Physiol. 2021, 68, 1059–1068. [Google Scholar] [CrossRef]
  11. Devoghalaere, F.; Doucen, T.; Guitton, B.; Keeling, J.; Payne, W.; Ling, T.J.; Ross, J.J.; Hallett, I.C.; Gunaseelan, K.; Dayatilake, G.A.; et al. A genomics approach to understanding the role of auxin in apple (Malus × domestica) fruit size control. BMC Plant Biol. 2012, 12, 7. [Google Scholar] [CrossRef] [PubMed]
  12. Bu, H.; Sun, X.; Yue, P.; Qiao, J.; Sun, J.; Wang, A.; Yuan, H.; Yu, W. The MdAux/IAA2 transcription repressor regulates cell and fruit size in apple fruit. Int. J. Mol. Sci. 2022, 23, 9454. [Google Scholar] [CrossRef] [PubMed]
  13. De Jong, M.; Wolters-Arts, M.; Schimmel, B.C.J.; Stultiens, C.L.M.; de Groot, P.F.M.; Powers, S.J.; Tikunov, Y.M.; Bovy, A.G.; Mariani, C.; Vriezen, W.H.; et al. Solanum lycopersicum AUXIN RESPONSE FACTOR 9 regulates cell division activity during early tomato fruit development. J. Exp. Bot. 2015, 66, 3405–3416. [Google Scholar] [CrossRef] [PubMed]
  14. Tian, Z.; Wu, B.; Liu, J.; Zhang, L.; Wu, T.; Wang, Y.; Han, Z.; Zhang, X. Genetic variations in MdSAUR36 participate in the negative regulation of mesocarp cell division and fruit size in Malus species. Mol. Breed. 2024, 44, 1. [Google Scholar] [CrossRef] [PubMed]
  15. Li, M.; Chen, R.; Gu, H.; Cheng, D.; Guo, X.; Shi, C.; Li, L.; Xu, G.; Gu, S.; Wu, Z.; et al. Grape small auxin upregulated RNA (SAUR) 041 is a candidate regulator of berry size in grape. Int. J. Mol. Sci. 2021, 22, 11818. [Google Scholar] [CrossRef] [PubMed]
  16. Peng, Z.; Li, W.; Gan, X.; Zhao, C.; Paudel, D.; Su, W.; Lv, J.; Lin, S.; Liu, Z.; Yang, X. Genome-wide analysis of SAUR gene family identifies a candidate associated with fruit size in Loquat (Eriobotrya japonica Lindl.). Int. J. Mol. Sci. 2022, 23, 13271. [Google Scholar] [CrossRef] [PubMed]
  17. Nardozza, S.; Cooney, J.; Boldingh, H.L.; Hewitt, K.G.; Trower, T.; Jones, D.; Thrimawithana, A.H.; Allan, A.C.; Richardson, A.C. Phytohormone and transcriptomic analysis reveals endogenous cytokinins affect kiwifruit growth under restricted carbon supply. Metabolites 2020, 10, 23. [Google Scholar] [CrossRef] [PubMed]
  18. Zhao, X.; Zhao, Z.; Cheng, S.; Wang, L.; Luo, Z.; Ai, C.; Liu, Z.; Liu, P.; Wang, L.; Wang, J.; et al. ZjWRKY23 and ZjWRKY40 promote fruit size enlargement by targeting and downregulating cytokinin oxidase/dehydrogenase 5 expression in Chinese Jujube. J. Agric. Food. Chem. 2023, 71, 18046–18058. [Google Scholar] [CrossRef]
  19. Gan, L.; Song, M.; Wang, X.; Yang, N.; Li, H.; Liu, X.; Li, Y. Cytokinins is involved in regulation of tomato pericarp thickness and fruit size. Hortic. Res. 2022, 9, uhab041. [Google Scholar] [CrossRef]
  20. Huang, W.; Peng, S.; Xian, Z.; Lin, D.; Hu, G.; Yang, L.; Ren, M.; Li, Z. Overexpression of a tomato miR171 target gene SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis. Plant Biotechnol. J. 2017, 15, 472–488. [Google Scholar] [CrossRef]
  21. He, H.; Yamamuro, C. Interplays between auxin and GA signaling coordinate early fruit development. Hortic. Res. 2022, 9, uhab078. [Google Scholar] [CrossRef]
  22. Ishiki, Y.; Oda, A.; Yaegashi, Y.; Orihara, Y.; Arai, T.; Hirabayashi, T.; Nakagawa, H.; Sato, T. Cloning of an auxin-responsive 1-aminocyclopropane-1-carboxylate synthase gene (CMe-ACS2) from melon and the expression of ACS genes in etiolated melon seedlings and melon fruits. Plant Sci. 2000, 159, 173–181. [Google Scholar] [CrossRef] [PubMed]
  23. Xin, T.; Zhang, Z.; Li, S.; Zhang, S.; Li, Q.; Zhang, Z.H.; Huang, S.; Yang, X. Genetic regulation of ethylene dosage for cucumber fruit elongation. Plant. Cell 2019, 31, 1063–1076. [Google Scholar] [CrossRef] [PubMed]
  24. Gupta, A.; Upadhyay, R.K.; Prabhakar, R.; Tiwari, N.; Garg, R.; Sane, V.A.; Sane, A.P. SlDREB3, a negative regulator of ABA responses, controls seed germination, fruit size and the onset of ripening in tomato. Plant Sci. 2022, 319, 111249. [Google Scholar] [CrossRef] [PubMed]
  25. Qi, X.; Liu, L.; Liu, C.; Song, L.; Dong, Y.; Chen, L.; Li, M. Sweet cherry AP2/ERF transcription factor, PavRAV2, negatively modulates fruit size by directly repressing PavKLUH expression. Physiol. Plant 2023, 175, e14065. [Google Scholar] [CrossRef] [PubMed]
  26. Yao, J.L.; Xu, J.; Cornille, A.; Tomes, S.; Karunairetnam, S.; Luo, Z.; Bassett, H.; Whitworth, C.; Rees-George, J.; Ranatunga, C.; et al. A microRNA allele that emerged prior to apple domestication may underlie fruit size evolution. Plant. J. 2015, 84, 417–427. [Google Scholar] [CrossRef] [PubMed]
  27. Ireland, H.S.; Yao, J.L.; Tomes, S.; Sutherland, P.W.; Nieuwenhuizen, N.; Gunaseelan, K.; Winz, R.A.; David, K.M.; Schaffer, R.J. Apple SEPALLATA1/2-like genes control fruit flesh development and ripening. Plant. J. 2013, 73, 1044–1056. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, J.; Zheng, C.; Shao, X.; Hu, Z.; Li, J.; Wang, P.; Wang, A.; Yu, J.; Shi, K. Transcriptomic and genetic approaches reveal an essential role of the NAC transcription factor SlNAP1 in the growth and defense response of tomato. Hortic. Res. 2020, 7, 209. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, C.; Zhao, B.; He, L.; Zhou, S.; Liu, Y.; Zhao, W.; Guo, S.; Wang, R.; Bai, Q.; Li, Y.; et al. The WOX family transcriptional regulator SlLAM1 controls compound leaf and floral organ development in Solanum lycopersicum. J. Exp. Bot. 2021, 72, 1822–1835. [Google Scholar] [CrossRef]
  30. Cong, B.; Barrero, L.S.; Tanksley, S.D. Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication. Nat. Genet. 2008, 40, 800–804. [Google Scholar] [CrossRef]
  31. Yi, K.; Menand, B.; Bell, E.; Dolan, L. A basic helix-loop-helix transcription factor controls cell growth and size in root hairs. Nat. Genet. 2010, 42, 108–264. [Google Scholar] [CrossRef]
  32. Machemer-Noonan; Katja, M. Interplay between Unusual MYB Transcription Factors and Their Role in Cell Size Regulation in Plants; The Ohio State University: Columbus, OH, USA, 2014. [Google Scholar]
  33. Nicolas, P.; Lecourieux, D.; Gomes, E.; Delrot, S.; Lecourieux, F. The grape berry-specific basic helixloophelix transcription factor VvCEB1 affects cell size. J. Exp. Bot. 2013, 64, 991–1003. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, L.; Xu, Y.; Li, Y.; Zheng, S.; Zhao, Z.; Chen, M.; Yang, H.; Yi, H.; Wu, J. Transcription factor CsMYB77 negatively regulates fruit ripening and fruit size in citrus. Plant Physiol. 2024, 194, 867–883. [Google Scholar] [CrossRef]
  35. Qiao, S.; Xu, Y.; Hu, Q.; Dong, W.; He, S.; Qi, X.; Sun, Y. Transcriptome analysis of sponge gourd (Luffa cylindrica) reveals candidate genes associated with fruit size. Agronomy 2022, 12, 1810. [Google Scholar] [CrossRef]
  36. Frary, A.; Nesbitt, T.C.; Frary, A.; Grandillo, S.; van der Knaap, E.; Cong, B.; Liu, J.P.; Meller, J.; Elber, R.; Alpert, K.B.; et al. Fw2.2: A quantitative trait locus key to the evolution of tomato fruit size. Science 2000, 289, 85–88. [Google Scholar] [CrossRef] [PubMed]
  37. Beauchet, A.; Gevaudant, F.; Gonzalez, N.; Chevalier, C. In search of the still unknown function of FW2.2/CELL NUMBER REGULATOR, a major regulator of fruit size in tomato. J. Exp. Bot. 2021, 72, 5300–5311. [Google Scholar] [CrossRef] [PubMed]
  38. Chakrabarti, M.; Zhang, N.; Sauvage, C.; Munos, S.; Blanca, J.; Canizares, J.; Diez, M.J.; Schneider, R.; Mazourek, M.; McClead, J.; et al. A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc. Natl. Acad. Sci. USA 2013, 110, 17125–17130. [Google Scholar] [CrossRef]
  39. Wang, Y.; Clevenger, J.P.; Illa-Berenguer, E.; Meulia, T.; van der Knaap, E.; Sun, L. A comparison of sun, ovate, fs8.1 and auxin application on tomato fruit shape and gene expression. Plant Cell Physiol. 2019, 60, 1067–1081. [Google Scholar] [CrossRef]
  40. Yuste-Lisbona, F.J.; Fernandez-Lozano, A.; Pineda, B.; Bretones, S.; Ortiz-Atienza, A.; Garcia-Sogo, B.; Mueller, N.A.; Angosto, T.; Capel, J.; Moreno, V.; et al. ENO regulates tomato fruit size through the floral meristem development network. Proc. Natl. Acad. Sci. USA 2020, 117, 8187–8195. [Google Scholar] [CrossRef]
  41. Yang, T.; He, Y.; Niu, S.; Zhang, Y. A YABBY gene CRABS CLAW a (CRCa) negatively regulates flower and fruit sizes in tomato. Plant Sci. 2022, 320, 111285. [Google Scholar] [CrossRef]
  42. Xu, C.; Liberatore, K.L.; MacAlister, C.A.; Huang, Z.; Chu, Y.; Jiang, K.; Brooks, C.; Ogawa-Ohnishi, M.; Xiong, G.; Pauly, M.; et al. A cascade of arabinosyltransferases controls shoot meristem size in tomato. Nat. Genet. 2015, 47, 784. [Google Scholar] [CrossRef]
  43. Li, M.; Wang, X.; Li, C.; Li, H.; Zhang, J.; Ye, Z. Silencing GRAS2 reduces fruit weight in tomato. J. Integr. Plant Biol. 2018, 60, 498–513. [Google Scholar] [CrossRef]
  44. Mu, Q.I.; Huang, Z.; Chakrabarti, M.; Illa-Berenguer, E.; Liu, X.; Wang, Y.; Ramos, A.; van der Knaap, E. Fruit weight is controlled by Cell Size Regulator encoding a novel protein that is expressed in maturing tomato fruits. PLoS Genet. 2017, 13, e1006930. [Google Scholar] [CrossRef]
  45. Zhao, F.; Zhang, J.; Weng, L.; Li, M.; Wang, Q.; Xiao, H. Fruit size control by a zinc finger protein regulating pericarp cell size in tomato. Mol. Hortic. 2021, 1, 6. [Google Scholar] [CrossRef]
  46. Weng, Y.; Colle, M.; Wang, Y.; Yang, L.; Rubinstein, M.; Sherman, A.; Ophir, R.; Grumet, R. QTL mapping in multiple populations and development stages reveals dynamic quantitative trait loci for fruit size in cucumbers of different market classes. Theor. Appl. Genet. 2015, 128, 1747–1763. [Google Scholar] [CrossRef]
  47. Pan, Y.; Liang, X.; Gao, M.; Liu, H.; Meng, H.; Weng, Y.; Cheng, Z. Round fruit shape in WI7239 cucumber is controlled by two interacting quantitative trait loci with one putatively encoding a tomato SUN homolog. Theor. Appl. Genet. 2017, 130, 573–586. [Google Scholar] [CrossRef]
  48. Tan, J.; Tao, Q.; Niu, H.; Zhang, Z.; Li, D.; Gong, Z.; Weng, Y.; Li, Z. A novel allele of monoecious (m) locus is responsible for elongated fruit shape and perfect flowers in cucumber (Cucumis sativus L.). Theor. Appl. Genet. 2015, 128, 2483–2493. [Google Scholar] [CrossRef]
  49. Zhao, J.Y. Molecular Mechanism of CsFUL1 and CsSUP Regulating Fruit Elongation in Cucumber. Ph.D. Thesis, China Agricultural University, Beijing, China, 2017. [Google Scholar] [CrossRef]
  50. Zhang, Z. Identification and Functional Analysis of SF2 in Controlling of Fruit Length in Cucumis sativus L. Ph.D. Thesis, Northwest Agriculture and Forestry University, Xianyang, China, 2019. [Google Scholar] [CrossRef]
  51. Wang, L.; Cao, C.; Zheng, S.; Zhang, H.; Liu, P.; Ge, Q.; Li, J.; Ren, Z. Transcriptomic analysis of short-fruit 1 (sf1) reveals new insights into the variation of fruit-related traits in Cucumis sativus. Sci. Rep. 2017, 7, 2950. [Google Scholar] [CrossRef]
  52. Monforte, A.J.; Diaz, A.; Cano-Delgado, A.; van der Knaap, E. The genetic basis of fruit morphology in horticultural crops: Lessons from tomato and melon. J. Exp. Bot. 2014, 65, 4625–4637. [Google Scholar] [CrossRef] [PubMed]
  53. Xu, P.; Wang, Y.; Sun, F.; Wu, R.; Du, H.; Wang, Y.; Jiang, L.; Wu, X.; Wu, X.; Yang, L.; et al. Long-read genome assembly and genetic architecture of fruit shape in the bottle gourd. Plant. J. 2021, 107, 956–968. [Google Scholar] [CrossRef]
  54. Zhang, H.; Tan, J.; Zhang, M.; Huang, S.; Chen, X. Comparative transcriptomic analysis of two bottle gourd accessions differing in fruit size. Genes 2020, 11, 359. [Google Scholar] [CrossRef] [PubMed]
  55. Yan, J.; Cheng, J.; Xie, D.; Wang, Y.; Wang, M.; Yang, S.; Jiang, B.; Chen, L.; Cai, J.; Liu, W. A nonsynonymous mutation in BhLS, encoding an acyl-CoA N-acyltransferase leads to fruit and seed size variation in wax gourd (Benincasa hispida). Theor. Appl. Genet. 2024, 137, 100. [Google Scholar] [CrossRef]
  56. CanoMedrano, R.; Darnell, R.L. Cell number and cell size in parthenocarpic vs. pollinated blueberry (Vaccinium ashei) fruits. Ann. Bot. 1997, 80, 419–425. [Google Scholar] [CrossRef]
  57. Su, W.; Shao, Z.; Wang, M.; Gan, X.; Yang, X.; Lin, S. EjBZR1 represses fruit enlargement by binding to the EjCYP90 promoter in loquat. Hortic. Res. 2021, 8, 152. [Google Scholar] [CrossRef]
  58. Harada, T.; Kurahashi, W.; Yanai, M.; Wakasa, Y.; Satoh, T. Involvement of cell proliferation and cell enlargement in increasing the fruit size of Malus species. Sci. Hortic. 2005, 105, 447–456. [Google Scholar] [CrossRef]
  59. Zhang, C.; Tanabe, K.; Wang, S.; Tamura, F.; Yoshida, A.; Matsumoto, K. The impact of cell division and cell enlargement on the evolution of fruit size in Pyrus pyrifolia. Ann. Bot. 2006, 98, 537–543. [Google Scholar] [CrossRef] [PubMed]
  60. Olmstead, J.W.; Lezzoni, A.F.; Whiting, M.D. Genotypic differences in sweet cherry fruit size are primarily a function of cell number. J. Am. Soc. Hortic. Sci. 2007, 132, 697–703. [Google Scholar] [CrossRef]
  61. Cerri, M.; Rosati, A.; Famiani, F.; Reale, L. Fruit size in different plum species (genus Prunus L.) is determined by post-bloom developmental processes and not by ovary characteristics at anthesis. Sci. Hortic. 2019, 255, 1–7. [Google Scholar] [CrossRef]
  62. Qi, W.; Yang, Y.; Feng, X.; Zhang, M.; Song, R. Mitochondrial function and maize kernel development requires Dek2, a pentatricopeptide repeat protein involved in nad1 mRNA splicing. Genetics 2017, 205, 239. [Google Scholar] [CrossRef]
  63. Subburaj, S.; Tu, L.; Lee, K.; Park, G.S.; Lee, H.; Chun, J.P.; Lim, Y.P.; Park, M.W.; McGregor, C.; Lee, G.J. A genome-wide analysis of the pentatricopeptide repeat (PPR) gene family and PPR-derived markers for flesh color in watermelon (Citrullus lanatus). Genes 2020, 11, 1125. [Google Scholar] [CrossRef]
  64. Zhang, N.; Brewer, M.T.; van der Knaap, E. Fine mapping of fw3.2 controlling fruit weight in tomato. Theor. Appl. Genet. 2012, 125, 273–284. [Google Scholar] [CrossRef] [PubMed]
  65. Nimmakayala, P.; Abburi, V.L.; Saminathan, T.; Alaparthi, S.B.; Almeida, A.; Davenport, B.; Nadimi, M.; Davidson, J.; Tonapi, K.; Yadav, L.; et al. Genome-wide diversity and association mapping for capsaicinoids and fruit weight in Capsicum annuum L. Sci. Rep. 2016, 6, 38081. [Google Scholar] [CrossRef] [PubMed]
  66. Cushing, D.A.; Forsthoefel, N.R.; Gestaut, D.R.; Vernon, D.M. Arabidopsis emb 175 and other ppr knockout mutants reveal essential roles for pentatricopeptide repeat (PPR) proteins in plant embryogenesis. Planta 2005, 221, 424–436. [Google Scholar] [CrossRef] [PubMed]
  67. Liu, Y.; Xiu, Z.; Meeley, R.; Tan, B. Empty pericarp5 encodes a pentatricopeptide repeat protein that is required for mitochondrial RNA editing and seed development in maize. Plant. Cell 2013, 25, 868–883. [Google Scholar] [CrossRef]
  68. Li, X.; Zhang, Y.; Hou, M.; Sun, F.; Shen, Y.; Xiu, Z.; Wang, X.; Chen, Z.; Sun, S.S.M.; Small, I.; et al. Small kernel 1 encodes a pentatricopeptide repeat protein required for mitochondrial nad7 transcript editing and seed development in maize (Zea mays) and rice (Oryza sativa). Plant. J. 2014, 79, 797–809. [Google Scholar] [CrossRef] [PubMed]
  69. Fang, Q.; Wang, Q.; Mao, H.; Xu, J.; Wang, Y.; Hu, H.; He, S.; Tu, J.; Cheng, C.; Tian, G.; et al. AtDIV2, an R-R-type MYB transcription factor of Arabidopsis, negatively regulates salt stress by modulating ABA signaling. Plant Cell Rep. 2018, 37, 1499–1511. [Google Scholar] [CrossRef] [PubMed]
  70. Davis, A.J.; Yan, Z.; Martinez, B.; Mumby, M.C. Protein phosphatase 2A is targeted to cell division control protein 6 by a calcium-binding regulatory subunit. J. Biol. Chem. 2008, 283, 16104–16114. [Google Scholar] [CrossRef] [PubMed]
  71. Grotz, N.; Guerinot, M.L. Molecular aspects of Cu, Fe and Zn homeostasis in plants. Biochim. Biophys. Acta 2006, 1763, 595–608. [Google Scholar] [CrossRef]
  72. Sasaki, A.; Yamaji, N.; Mitani-Ueno, N.; Kashino, M.; Ma, J.F. A node-localized transporter OsZIP3 is responsible for the preferential distribution of Zn to developing tissues in rice. Plant. J. 2015, 84, 374–384. [Google Scholar] [CrossRef]
  73. Lee, S.; Kim, S.A.; Lee, J.; Guerinot, M.L.; An, G. Zinc deficiency-inducible OsZIP8 encodes a plasma membrane-localized Zinc transporter in Rice. Mol. Cells 2010, 29, 551–558. [Google Scholar] [CrossRef]
  74. Vert, G.; Grotz, N.; Dedaldechamp, F.; Gaymard, F.; Guerinot, M.L.; Briat, J.F.; Curie, C. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant. Cell 2002, 14, 1223–1233. [Google Scholar] [CrossRef]
  75. Dong, N.Q.; Sun, Y.; Guo, T.; Shi, C.L.; Zhang, Y.M.; Kan, Y.; Xiang, Y.H.; Zhang, H.; Yang, Y.B.; Li, Y.C.; et al. UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice. Nat. Commun. 2020, 11, 2629. [Google Scholar] [CrossRef]
  76. Zhang, G.Z.; Jin, S.H.; Jiang, X.Y.; Dong, R.R.; Li, P.; Li, Y.J.; Hou, B.K. Ectopic expression of UGT75D1, a glycosyltransferase preferring indole-3-butyric acid, modulates cotyledon development and stress tolerance in seed germination of Arabidopsis thaliana. Plant Mol. Biol. 2016, 90, 77–93. [Google Scholar] [CrossRef] [PubMed]
  77. Li, P.; Lei, K.; Li, Y.; He, X.; Wang, S.; Liu, R.; Ji, L.; Hou, B. Identification and characterization of the first cytokinin glycosyltransferase from rice. Rice 2019, 12, 19. [Google Scholar] [CrossRef]
  78. Huang, M.; Zhu, X.; Bai, H.; Wang, C.; Gou, N.; Zhang, Y.; Chen, C.; Yin, M.; Wang, L.; Wuyun, T. Comparative anatomical and transcriptomics reveal the larger cell size as a major contributor to larger fruit size in Apricot. Int. J. Mol. Sci. 2023, 24, 8748. [Google Scholar] [CrossRef]
  79. Marowa, P.; Ding, A.; Kong, Y. Expansins: Roles in plant growth and potential applications in crop improvement. Plant Cell Rep. 2016, 35, 949–965. [Google Scholar] [CrossRef]
  80. Di Marzo, M.; Viana, V.E.; Banfi, C.; Cassina, V.; Corti, R.; Herrera-Ubaldo, H.; Babolin, N.; Guazzotti, A.; Kiegle, E.; Gregis, V.; et al. Cell wall modifications by α-XYLOSIDASE1 are required for control of seed and fruit size in Arabidopsis. J. Exp. Bot. 2022, 73, 1499–1515. [Google Scholar] [CrossRef]
  81. Kwon, Y.; Kim, J.H.; Nguyen, H.N.; Jikumaru, Y.; Kamiya, Y.; Hong, S.W.; Lee, H. A novel Arabidopsis MYB-like transcription factor, MYBH, regulates hypocotyl elongation by enhancing auxin accumulation. J. Exp. Bot. 2013, 64, 3911–3922. [Google Scholar] [CrossRef]
  82. Li, W.F.; Zhou, Q.; Ma, Z.H.; Zuo, C.W.; Chu, M.Y.; Mao, J.; Chen, B.H. Regulatory mechanism of GA (3) application on grape (Vitis vinifera L.) berry size. Plant Physiol. Biochem. 2024, 210, 108543. [Google Scholar] [CrossRef]
  83. Fukazawa, J.; Sakai, T.; Ishida, S.; Yamaguchi, I.; Kamiya, Y.; Takahashi, Y. Repression of shoot growth, a bZIP transcriptional activator, regulates cell elongation by controlling the level of gibberellins. Plant. Cell 2000, 12, 901–915. [Google Scholar] [CrossRef]
  84. Weiste, C.; Droge-Laser, W. The Arabidopsis transcription factor bZIP11 activates auxin-mediated transcription by recruiting the histone acetylation machinery. Nat. Commun. 2014, 5, 3883. [Google Scholar] [CrossRef]
  85. Hill, J.T.; Demarest, B.L.; Bisgrove, B.W.; Gorsi, B.; Su, Y.C.; Yost, H.J. MMAPPR: Mutation Mapping Analysis Pipeline for Pooled RNA-seq. Genome Res. 2013, 23, 687–697. [Google Scholar] [CrossRef]
  86. Fekih, R.; Takagi, H.; Tamiru, M.; Abe, A.; Natsume, S.; Yaegashi, H.; Sharma, S.; Sharma, S.; Kanzaki, H.; Matsumura, H.; et al. MutMap+: Genetic Mapping and Mutant Identification without Crossing in Rice. PLoS ONE 2013, 8, e68529. [Google Scholar] [CrossRef]
  87. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  88. Kanehisa, M.; Goto, S.; Kawashima, S.; Okuno, Y.; Hattori, M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32, 277–280. [Google Scholar] [CrossRef] [PubMed]
Plants 13 02154 g001
Figure 1. Bottle gourd fruits: male H28, female H17, their F1, and some F2. (A) Mature fruit of parent and F1. (B) Ovary of parent and F1. (C) Mature fruit of some F2.
Figure 1. Bottle gourd fruits: male H28, female H17, their F1, and some F2. (A) Mature fruit of parent and F1. (B) Ovary of parent and F1. (C) Mature fruit of some F2.
Plants 13 02154 g001
Plants 13 02154 g002
Figure 2. Dynamic changes in the fruit growth process in bottle gourd. (A) The growth curve of fruit width of parent and their F1 from −3 DAP to 30 DAP. (B) The growth curve of fruit length of parent and their F1 from −3 DAP to 30 DAP.
Figure 2. Dynamic changes in the fruit growth process in bottle gourd. (A) The growth curve of fruit width of parent and their F1 from −3 DAP to 30 DAP. (B) The growth curve of fruit length of parent and their F1 from −3 DAP to 30 DAP.
Plants 13 02154 g002
Plants 13 02154 g003
Figure 3. Microscopic observation and comparison of the average cell area and number in transverse sections of H17 and H28 fruit at different periods. (A) Paraffin section of H17 and H28 fruit from −3 DAP to 30 DAP; bar = 50 μm. (B) Statistical diagram of cell number. (C) Statistical diagram of cell area. *, p < 0.05; **, p < 0.01.
Figure 3. Microscopic observation and comparison of the average cell area and number in transverse sections of H17 and H28 fruit at different periods. (A) Paraffin section of H17 and H28 fruit from −3 DAP to 30 DAP; bar = 50 μm. (B) Statistical diagram of cell number. (C) Statistical diagram of cell area. *, p < 0.05; **, p < 0.01.
Plants 13 02154 g003
Plants 13 02154 g004
Figure 4. Mapping of fruit size in bottle gourd. (A) Chromosome distribution in terms of ΔSNP-index correlation value. The colored dots represent the calculated ΔSNP-index values, the black lines represent the fitted ΔSNP-index values, and the red line represents the 99% threshold line. (B) Candidate interval of fruit-size genes in bottle gourd. Black represents the genotype of small-fruit H28, white represents the genotype of large-fruit H17, and gray represents the F1 genotype.
Figure 4. Mapping of fruit size in bottle gourd. (A) Chromosome distribution in terms of ΔSNP-index correlation value. The colored dots represent the calculated ΔSNP-index values, the black lines represent the fitted ΔSNP-index values, and the red line represents the 99% threshold line. (B) Candidate interval of fruit-size genes in bottle gourd. Black represents the genotype of small-fruit H28, white represents the genotype of large-fruit H17, and gray represents the F1 genotype.
Plants 13 02154 g004
Plants 13 02154 g005
Figure 5. Transcriptomic analysis of fruit size of H28 and H17. (A) Number of DEGs. (B) Heat map showing the expression pattern of genes in blue and red tissues. Horizontal represents genes and each column represents a sample, with red representing highly expressed genes and green representing low-expressed genes. (C) Rich factor of GO enrichment analysis of DEGs. (D) Histogram of GO enrichment analysis of DEGs. (E) Histogram of KEGG enrichment analysis of DEGs.
Figure 5. Transcriptomic analysis of fruit size of H28 and H17. (A) Number of DEGs. (B) Heat map showing the expression pattern of genes in blue and red tissues. Horizontal represents genes and each column represents a sample, with red representing highly expressed genes and green representing low-expressed genes. (C) Rich factor of GO enrichment analysis of DEGs. (D) Histogram of GO enrichment analysis of DEGs. (E) Histogram of KEGG enrichment analysis of DEGs.
Plants 13 02154 g005
Table 1. Analysis of width and length of bottle gourd fruit in parents and F1.
Table 1. Analysis of width and length of bottle gourd fruit in parents and F1.
PopulationFruit Width (mm)Fruit Length (mm)
P1 (H28)48.47 ± 3.7367.99 ± 3.41
P2 (H17)163.18 ± 9.63295.76 ± 49.23
F1 (H28 × H17)121.32 ± 7.60176.02 ± 32.46
Table 2. RNA-Seq reads obtained from H28 and H17 fruits.
Table 2. RNA-Seq reads obtained from H28 and H17 fruits.
SampleClean Reads (Strip)Clean Bases (bp)Q30 (%)
H286d139,275,1505,930,547,65094.13
H286d238,048,2305,745,282,73094.46
H176d139,540,4905,970,613,99093.97
H176d233,727,9205,092,915,92094.05
Table 3. Functional annotation and relative expression mode of candidate genes.
Table 3. Functional annotation and relative expression mode of candidate genes.
Gene IDGene Function AnnotationNon-Synonymous MutationRelative Expression Mode
HG_GLEAN_10001518Transcription factor DIVARICATA subtype X1, MYB familyYESNodiff
HG_GLEAN_10001525Cell division control protein 6 homology B-like proteinYESNodiff
HG_GLEAN_10001544Zinc transporter 6, ZIP protein familyNOUpregulation
HG_GLEAN_10001548Pentapeptide-repeat (PPR) sequence proteinYESUpregulation
HG_GLEAN_10001556UDP-glycosyltransferase 87A1-like (UGT) proteinYESNodiff
HG_GLEAN_10001558UDP-glycosyltransferase 87A1-like (UGT) proteinYESNodiff
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fang, H.; Huang, S.; Li, R.; Wang, P.; Jiang, Q.; Zhong, C.; Yang, Y.; Yu, W. Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.). Plants 2024, 13, 2154. https://doi.org/10.3390/plants13152154

AMA Style

Fang H, Huang S, Li R, Wang P, Jiang Q, Zhong C, Yang Y, Yu W. Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.). Plants. 2024; 13(15):2154. https://doi.org/10.3390/plants13152154

Chicago/Turabian Style

Fang, Huarong, Shishi Huang, Ruirui Li, Peng Wang, Qingwei Jiang, Chuan Zhong, Yanjuan Yang, and Wenjin Yu. 2024. "Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.)" Plants 13, no. 15: 2154. https://doi.org/10.3390/plants13152154

APA Style

Fang, H., Huang, S., Li, R., Wang, P., Jiang, Q., Zhong, C., Yang, Y., & Yu, W. (2024). Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.). Plants, 13(15), 2154. https://doi.org/10.3390/plants13152154

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop