Identification and Application of BhAPRR2 Controlling Peel Colour in Wax Gourd (Benincasa hispida)

Identification of Candidate Genes for Green Rind Color in Watermelon

The color of the rind is one of the most crucial agronomic characteristics of watermelon (Citrullus lanatus L.). Its genetic analysis was conducted to provide the identification of genes regulating rind color and improving the quality of watermelon appearance. In this study, a mapping population of 505 F2 plants, derived from a cross between green (CG058) and light-green (CG265) rinds, along with a high-density genetic linkage (average 0.9 cM distance between bin markers), was used to map and identify possible candidate genes. The green rind trait was determined to be regulated by a single Mendelian locus and was precisely located within a 110 kb genomic site on chromosome nine (Chr 9). In the respective region, two potential genes, Cla97C09G175170 and Cla97C09G175180, were substantially downregulated in the light-green rind in comparison to the green rind. Previous studies revealed that Cla97C09G175170, encoding a two-component response regulator-like protein (APRR2), is possibly involved in the green rind trait in watermelon. Virus-induced gene silencing (VIGS) assay confirmed that ClAPRR2 is a key gene responsible for green rind color. Moreover, qRT-PCR analysis revealed that the transcription levels of multiple key genes in the chlorophyll (Chl) biosynthesis pathway were downregulated in the light-green rind relative to the green rind. The current findings have the potential to clarify the regulatory mechanisms that underlie the color of the watermelon rind. These data would provide valuable insights for the targeted molecular design and development of watermelon rinds.

Map-based cloning revealed BhAPRR2 gene regulating the black peel formation of mature fruit in wax gourd (Benincasa hispida)

KEY MESSAGE

Map-based cloning revealed BhAPRR2, encoding a two-component response-regulating protein that regulates the black peel formation of mature fruit in wax gourd. Wax gourd is an economically significant vegetable crop, and peel color is a crucial agronomic trait that influences its commercial value. Although genes controlling light green or white peel have been cloned in wax gourd, the genetic basis and molecular mechanism underlying black peel remain unclear. Here, we confirmed that the peel color of wax gourd is a qualitative trait governed by single gene, with black being dominant over green. Through bulked segregant analysis sequencing (BSA-seq) and map-based cloning, we identified Bh.pf3chr5g483 as the candidate gene. This gene encodes a two-component response-regulating protein and is homologous to APRR2, referred to as BhAPRR2. Compared to P170, the BhAPRR2 in YD1 exhibits multiple mutations in both its coding and promoter regions. Notably, the mutations in the coding region do not affect its nuclear localization or transcriptional activation activity. However, the mutations in the promoter region substantially increase its expression in the peel of YD1, potentially contributing to the black peel phenotype observed in this variety. Furthermore, we developed an insertion/deletion (InDel) marker based on a 93-base pair (bp) insertion/deletion mutation in the promoter region of BhAPRR2, which achieved up to 95.8% phenotypic accuracy in a natural population comprising 165 wax gourd germplasms. In summary, our findings suggest that mutations in the promoter region of BhAPRR2 may contribute to the development of black peel in wax gourd. This discovery provides new insights into the molecular and genetic mechanisms underlying peel color diversity and offers a valuable molecular marker for wax gourd breeding efforts.

Fine mapping and functional validation of the candidate gene BhGA2ox3 for fruit pedicel length in wax gourd (Benincasa hispida)

KEY MESSAGE

The gene regulating fruit pedicel length in wax gourd was finely mapped to a 211 kb region on chromosome 8. The major gene, Bch08G017310 (BhGA2ox3), was identified through forward genetics. Fruit pedicel length (FPL) is a crucial trait in wax gourd (Benincasa hispida) that affects fruit development and cultivation management. However, the key regulatory genes and mechanisms of FPL in wax gourds remain poorly understood. In this study, we constructed an F2 population using wax gourd plants with long fruit pedicels (GF-7-1-1) and short fruit pedicels (YSB-1-1-2) as parents. Through BSA-seq, we initially localised the FPL candidate gene to an 8.4 Mb region on chromosome 8, which was further narrowed down to a 1.1 Mb region via linkage analysis. A large F2 population of 2163 individuals was used to screen for recombinants, and the locus was ultimately narrowed to within a 211 kb (62,299,856-62,511,174 bp) region. Sequence and expression analyses showed that Bch08G017310 (named BhGA2ox3) is a strong candidate gene for FPL in wax gourds. It encodes gibberellin (GA) 2-beta-dioxygenase, a member of the GA 2-oxidase (GA2ox) family. Cytology showed that GA treatment significantly elongated the fruit pedicels and enlarged the cells in the plants with short fruit pedicels. Ectopic expression of BhGA2ox3 showed that BhGA2ox3 overexpression in Arabidopsis thaliana resulted in significantly shorter fruit pedicels. This study lays a theoretical foundation for the regulatory mechanism of FPL in wax gourds and molecular breeding.

A point mutation in the zinc-finger transcription factor CqLOL1 controls the green flesh color in chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How)

Introduction

Flesh color is an essential trait in chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How); however, the inheritance and molecular basis of green flesh trait remain unclear.

Methods

In the present study, two F2 populations, derived from 1742 (white flesh) × FJ3211 (green flesh) and J16 (white flesh) × FJ5 (green flesh), were used to identify the green flesh (Cqgf) locus.

Results

Genetic analysis revealed that the presence of green flesh was a quantitative trait that closely followed a normal distribution. Combining the results from QTL mapping and BSA-seq analysis, the Cqgf locus was preliminarily determined to be located on chromosome 05 and was narrowed down to a 2.55-Mb interval by linkage analysis. A large J16 × FJ5 F2 population comprising 3,180 individuals was subsequently used to screen the recombinants, and the Cqgf locus was fine-mapped to a region of 329.70 kb that harbors six genes. One of the candidate genes, Bch05G003700, the zinc-finger transcription factor LOL1 (lsd one like 1 protein; CqLOL1), was the strongest candidate gene for the Cqgf locus according to sequence variation and expression analysis. Additionally, a point mutation (A > C) in CqLOL1 resulted in the substitution of threonine (T) with proline (P) in the amino acid sequence, showing a complete relationship linked with flesh color in a panel of 45 germplasms.

Discussion

The study suggests that CqLOL1 promotes the accumulation of chlorophyll content in chieh-qua and lead to green flesh. Our findings establish a theoretical and technical foundation for breeding different flesh color lines and elucidating the underlying mechanisms of flesh color in chieh-qua.

Mapping and transcriptomic profiling reveal that the KNAT6 gene is involved in the dark green peel colour of mature pumpkin fruit (Cucurbita maxima L.)

KEY MESSAGE

We identified a 580 bp deletion of CmaKNAT6 coding region influences peel colour of mature Cucurbita maxima fruit. Peel colour is an important agronomic characteristic affecting commodity quality in Cucurbit plants. Genetic mapping of fruit peel colour promotes molecular breeding and provides an important basis for understanding the regulatory mechanism in Cucurbit plants. In the present study, the Cucurbita maxima inbred line '9-6' which has a grey peel colour and 'U3-3-44' which has a dark green peel colour in the mature fruit stage, were used as plant materials. At 5-40 days after pollination (DAP), the contents of chlorophyll a, chlorophyll b, total chlorophyll and carotenoids in the 'U3-3-44' peels were significantly greater than those in the '9-6' peels. In the epicarp of the '9-6' mature fruit, the presence of nonpigmented cell layers and few chloroplasts in each cell in the pigmented layers were observed. Six generations derived by crossing '9-6' and 'U3-3-44' were constructed, and the dark green peel was found to be controlled by a single dominant locus, which was named CmaMg (mature green peel). Through bulked-segregant analysis sequencing (BSA-seq) and insertion-deletion (InDel) markers, CmaMg was mapped to a region of approximately 449.51 kb on chromosome 11 using 177 F2 individuals. Additionally, 1703 F2 plants were used for fine mapping to compress the candidate interval to a region of 32.34 kb. Five coding genes were in this region, and CmaCh11G000900 was identified as a promising candidate gene according to the reported function, sequence alignment, and expression analyses. CmaCh11G000900 (CmaKNAT6) encodes the homeobox protein knotted-1-like 6 and contains 4 conserved domains. CmaKNAT6 of '9-6' had a 580 bp deletion, leading to premature transcriptional termination. The expression of CmaKNAT6 tended to increase sharply during the early fruit development stage but decrease gradually during the late period of fruit development. Allelic diversity analysis of pumpkin germplasm resources indicated that the 580 bp deletion in the of CmaKNAT6 coding region was associated with peel colour. Subcellular localization analysis indicated that CmaKNAT6 is a nuclear protein. Transcriptomic analysis of the inbred lines '9-6' and 'U3-3-44' indicated that genes involved in chlorophyll biosynthesis were more enriched in 'U3-3-44' than in '9-6'. Additionally, the expression of transcription factor genes that positively regulate chlorophyll synthesis and light signal transduction pathways was upregulated in 'U3-3-44'. These results lay a foundation for further studies on the genetic mechanism underlying peel colour and for optimizing peel colour-based breeding strategies for C. maxima.

Map-based cloning reveals Cpgp gene encoding an APRR2 protein to regulate the green fruit peel formation in Cucurbita pepo

Fruit peel color is a major factor that influences fruit quality and customers' demand. However, the molecular mechanisms underlying the green fruit peel color trait of Cucurbita pepo L. remain unknown. Two parental lines, RP16 and RP38, were used to study the fruit peel color trait in C. pepo. The parental line RP16 shows white peel color, whereas RP38 exhibits green peel color. 384 F2 populations were used to identify the inheritance pattern associated with green fruit and white fruit peel in Cucurbita pepo L. 293 F2 individuals were white, and 91 F2 individuals were green, resulting in a ratio of 3:1. Hence, white peel is dominant over the green fruit peel trait, and a single recessive green peel gene (Cpgp) controls the green fruit peel. The fruit chlorophyll (Chll) content decreases as fruit matures in the RP16 line. In contrast, Chll increases during the fruit growing periods on fruit peels of the RP38 line. The BSA-sequence analysis revealed the Cpgp locus on Chr5, within a 2.3 Mb region. Subsequent fine-mapping analysis, using 699 F2 plants, narrowed down this region to 23.90 kb on the same chromosome. Within this region, two annotated genes, namely Cp4.1LG05g02070 and Cp4.1LG05g02060, are present. These genes are predicted to encode a two-component Arabidopsis Pseudo-Response Regulator 2-like protein (APRR2), which may be involved in green pigmentation processes in plants. Consequently, sequence alignment and gene expression analyses at various fruit development stages supported that Cp4.1LG05g02070 may be the primary candidate gene responsible for regulating the green fruit peel color trait in Cucurbita pepo L. This study may provide a basis for further study on the basic mechanisms that control the fruit peel colors in Cucurbita spp.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01492-7.

Quantitative Trait Loci Sequencing and Genetic Mapping Reveal Two Main Regulatory Genes for Stem Color in Wax Gourds

Stem color is an important agronomic trait of wax gourds. However, its regulatory genes have not been identified. In this study, 105 inbred lines constructed from two parents (GX-71 and MY-1) were sequenced and quantitative trait loci sequencing was used to mine the genes that regulate stem color in wax gourds. The results identified two quantitative trait loci related to stem color, qSC5 and qSC12, located on Chr05 (11,134,567-16,459,268) and Chr12 (74,618,168-75,712,335), respectively. The qSC5 had a phenotypic variation rate of 36.9% and a maximum limit of detection of 16.9. And the qSC12 had a phenotypic variation rate of 20.9%, and a maximum limit of detection of 11.2. Bch05G003950 (named BchAPRR2) and Bch12G020400 were identified as candidate genes involved in stem color regulation in wax gourds. The chlorophyll content and expression of BchAPRR2 and Bch12G020400 were significantly higher in green-stemmed wax gourds than in white-stemmed ones. Therefore, BchAPRR2 and Bch12G020400 were considered the main and secondary regulatory genes for wax gourd stem color, respectively. Finally, InDel markers closely linked to BchAPRR2 were developed to validate the prediction of wax gourd stem color traits in 55 germplasm lines, with an accuracy of 81.8%. These findings lay the foundation for exploring the genetic regulation of wax gourd stem color and future research on wax gourd breeding.

Fine mapping and transcriptome profiling reveal CpAPRR2 to modulate immature fruit rind color formation in zucchini (Cucurbita pepo)

KEY MESSAGE

A large fragment deletion of CpAPRR2, encoding a two-component response regulator-like protein, which influences immature white rind color formation in zucchini (Cucurbita pepo). Fruit rind color is an important agronomic trait that affects commodity quality and consumer choice in zucchini (Cucurbita pepo). However, the molecular mechanism controlling rind color is unclear. We characterized two zucchini inbred lines: '19' (dark green rind) and '113' (white rind). Genetic analysis revealed white immature fruit rind color to be controlled by a dominant locus (CpW). Combining bulked segregant analysis sequencing (BSA-seq) and Kompetitive Allele-Specific PCR (KASP) markers, we mapped the CpW locus to a 100.4 kb region on chromosome 5 and then narrow down the candidate region to 37.5 kb using linkage analysis of 532 BC1 and 1613 F2 individuals, including 6 coding genes. Among them, Cp4.1LG05g02070 (CpAPRR2), encoding a two-component response regulator-like protein, was regarded to be a promising candidate gene. The expression level of CpAPRR2 in dark green rind was significantly higher than that in white rind and was induced by light. A deletion of 2227 bp at the 5' end of CpAPRR2 in '113' might explain the white phenotype. Further analysis of allelic diversity in zucchini germplasm resources revealed rind color to be associated with the deletion of CpAPRR2. Subcellular localization analysis indicated that CpAPRR2 was a nuclear protein. Transcriptome analysis using near-isogenic lines with dark green (DG) and white (W) rind indicated that genes involved in photosynthesis and porphyrin metabolism pathways were enriched in DG compared with W. Additionally, chlorophyll synthesis-related genes were upregulated in DG. These results identify mechanisms of zucchini rind color and provide genetic resources for breeding.

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

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

Potential locus W and candidate gene McPRR2 associated with pericarp pigment accumulation in bitter gourd (Momordica charantia L.) revealed via BSA-seq analysis

Pericarp color is a prominent agronomic trait that exerts a significant impact on consumer and breeder preferences. Genetic analysis has revealed that the pericarp color of bitter gourd is a quantitative trait. However, the underlying mechanism for this trait in bitter gourd remains largely unknown. In the present study, we employed bulked segregant analysis (BSA) to identify the candidate genes responsible for bitter gourd pericarp color (specifically, dark green versus white) within F2 segregation populations resulting from the crossing of B07 (dark green pericarp) and A06 (white pericarp). Through genomic variation, genetic mapping, and expression analysis, we identified a candidate gene named McPRR2, which was a homolog of Arabidopsis pseudo response regulator 2 (APRR2) encoded by LOC111023472. Sequence alignment of the candidate gene between the two parental lines revealed a 15-bp nucleotide insertion in the coding region of LOC111023472, leading to a premature stop codon and potentially causing a loss-of-function mutation. qRT-PCR analysis demonstrated that the expression of McPRR2 was significantly higher in B07 compared to A06, and it was primarily expressed in the immature fruit pericarp. Moreover, overexpression of McPRR2 in tomato could enhance the green color of immature fruit pericarp by increasing the chlorophyll content. Consequently, McPRR2 emerged as a strong candidate gene regulating the bitter gourd pericarp color by influencing chlorophyll accumulation. Finally, we developed a molecular marker linked to pericarp color, enabling the identification of genotypes in breeding populations. These findings provided valuable insights into the genetic improvement of bitter gourd pericarp color.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

Genome-wide identification of bZIP gene family and expression analysis of BhbZIP58 under heat stress in wax gourd

BACKGROUND

The basic leucine zipper (bZIP) transcription factor family is one of the most abundant and evolutionarily conserved gene families in plants. It assumes crucial functions in the life cycle of plants, including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. Although the genome of wax gourd has been published, little is known about the functions, evolutionary background, and gene expression patterns of the bZIP gene family, which limits its utilization.

RESULTS

A total of 61 bZIP genes (BhbZIPs) were identified from wax gourd (Benincasa hispida) genome and divided into 12 subgroups. Whole-genome duplication (WGD) and dispersed duplication (DSD) were the main driving forces of bZIP gene family expansion in wax gourd, and this family may have undergone intense purifying selection pressure during the evolutionary process. We selected BhbZIP58, only one in the member of subgroup B, to study its expression patterns under different stresses, including heat, salt, drought, cold stress, and ABA treatment. Surprisingly, BhbZIP58 had a dramatic response under heat stress. BhbZIP58 showed the highest expression level in the root compared with leaves, stem, stamen, pistil, and ovary. In addition, BhbZIP58 protein was located in the nucleus and had transcriptional activation activity. Overexpression of BhbZIP58 in Arabidopsis enhanced their heat tolerance.

CONCLUSIONS

In this study, bZIP gene family is systematically bioinformatically in wax gourd for the first time. Particularly, BhbZIP58 may have an important role in heat stress. It will facilitate further research on the bZIP gene family regarding their evolutionary history and biological functions.

McAPRR2: The Key Regulator of Domesticated Pericarp Color in Bitter Gourd

Pericarp color is a crucial commercial trait influencing consumer preferences for bitter gourds. However, until now, the gene responsible for this trait has remained unidentified. In this study, we identified a gene (McAPRR2) controlling pericarp color via a genome-wide association study (GWAS) utilizing the resequencing data of 106 bitter gourd accessions. McAPRR2 exhibits three primary haplotypes: Hap1 is a wild type with a green pericarp, Hap2 is a SA (South Asian) and SEA (Southeast Asia) type with a green pericarp, and Hap3 is primarily a SEA type with a light green pericarp. The McAPRR2 haplotype is significantly correlated with both pericarp color and ecological type. Importantly, McAPRR2 with the light green pericarp demonstrated premature termination due to a 15 bp sequence insertion. The phylogenetic tree clustered according to pericarp color and ecological type, using SNPs located in the McAPRR2 gene and its promoter. High πwild/SEA and πSA/SEA values indicate high nucleotide diversity between wild and SEA types and between SA and SEA types in the McAPRR2 gene. The haplotypes, phylogenetic tree, and nucleotide diversity of McAPRR2 suggest that McAPRR2 has undergone domestication selection. This study identifies McAPRR2 as the key gene determining pericarp color in bitter gourds and introduces a novel insight that McAPRR2 is subject to domestication selection.

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

Introduction/Background

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

Methods

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

Results

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

Discussion

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

Map-based cloning of the APRR2 gene controlling green stigma in bitter gourd (Momordica charantia)

Bitter gourd is an economically important vegetable and medicinal crop distinguished by its bitter fruits. Its stigma color is widely used to assess the distinctiveness, uniformity, and stability of bitter gourd varieties. However, limited researches have been dedicated to genetic basis of its stigma color. In this study, we employed bulked segregant analysis (BSA) sequencing to identify a single dominant locus McSTC1 located on pseudochromosome 6 through genetic mapping of an F₂ population (n =241) derived from the cross between green and yellow stigma parental lines. An F₂-derived F₃ segregation population (n = 847) was further adopted for fine mapping, which delimited the McSTC1 locus to a 13.87 kb region containing one predicted gene McAPRR2 (Mc06g1638), a homolog of the Arabidopsis two-component response regulator-like gene AtAPRR2. Sequence alignment analysis of McAPRR2 revealed that a 15 bp insertion at exon 9 results in a truncated GLK domain of its encoded protein, which existed in 19 bitter gourd varieties with yellow stigma. A genome-wide synteny search of the bitter gourd McAPRR2 genes in Cucurbitaceae family revealed its close relationship with other cucurbits APRR2 genes that are corresponding to white or light green fruit skin. Our findings provide insights into the molecular marker-assisted breeding of bitter gourd stigma color and the mechanism of gene regulation for stigma color.

Comparative transcriptome analysis of molecular mechanisms underlying adventitious root developments in Huangshan Bitter tea (Camellia gymnogyna Chang) under red light quality

As the formation of adventitious roots (AR) is an important component of in vitro regeneration of tea plants, the propagation and preservation of Huangshan Bitter tea (Camellia gymnogyna Chang) cuttings have been hindered due to its lower rooting rate. As light is a crucial environmental factor that affects AR formation, this study aimed to investigate the special role of red light (RL) in the formation of AR in Huangshan Bitter tea plants, which has not been well understood. Huangshan Bitter tea plants were induced with white light (control, WL) and red light (660 nm, RL) qualities 36 days after induced treatment (DAI) to investigate dynamic AR formation and development, anatomical observation, hormones content change, and weighted gene co-expression network analysis (WGCNA) of the transcriptome. Results showed that RL promoted the rooting rate and root characteristics compared to WL. Anatomical observations demonstrated that root primordium was induced earlier by RL at the 4 DAI. RL positively affected IAA, ZT and GA₃ content and negatively influenced ABA from the 4 to 16 DAI. RNA-seq and analysis of differential expression genes (DEGs) exhibited extensive variation in gene expression profiles between RL and WL. Meanwhile, the results of WGCNA and correlation analysis identified three highly correlated modules and hub genes mainly participated in 'response to hormone', 'cellular glucan metabolic progress', and 'response to auxin'. Furthermore, the proportion of transcription factors (TFs) such as ethylene response factor (ERF), myeloblastosis (MYB), basic helix-loop-helix (bHLH), and WRKYGQK (WRKY) were the top four in DEGs. These results suggested that the AR-promoting potential of red light was due to complex hormone interactions in tea plants by regulating the expression of related genes. This study provided an important reference to shorten breeding cycles and accelerate superiority in tea plant propagation and preservation.

Forward genetic studies reveal LsAPRR2 as a key gene in regulating the green color of pericarp in bottle gourd (Lagenaria siceraria)

The fruit peel color is an important factor that affects its quality. However, genes involved in regulating pericarp color in bottle gourd (Lagenaria siceraria) have not been explored to date. Genetic analysis of color traits in bottle gourd peel through a genetic population of six generations demonstrated that the green color of peels is inherited as a single gene dominant trait. Combined phenotype-genotype analysis of recombinant plants using BSA-seq mapped the candidate gene to a 22.645 Kb interval at the head end of chromosome 1. We observed that the final interval contained only one gene, LsAPRR2 (HG_GLEAN_10010973). Sequence and spatiotemporal expression analyses of LsAPRR2 unraveled two nonsynonymous mutations (A→G) and (G→C) in the parental CDS sequences. Further, LsAPRR2 expression was higher in all green-skinned bottle gourds (H16) at various stages of fruit development than in white-skinned bottle gourds (H06). Cloning and sequence comparison of the two parental LsAPRR2 promoter regions indicated 11 bases insertion and 8 SNPs mutations in the region -991~-1033, upstream of the start codon in white bottle gourd. Proof of GUS reporting system, Genetic variation in this fragment significantly reduced the expression of LsAPRR2 in the pericarp of white bottle gourd. In addition, we developed a tightly linked (accuracy 93.88%) InDel marker for the promoter variant segment. Overall, the current study provides a theoretical basis for comprehensive elucidation of the regulatory mechanisms underlying the determination of bottle gourd pericarp color. This would further help in the directed molecular design breeding of bottle gourd pericarp.

Fine Mapping and Identification of SmAPRR2 Regulating Rind Color in Eggplant (Solanum melongena L.)

Rind color is an economically important agronomic trait in eggplant that impacts consumer preferences. In this study, bulked segregant analysis and competitive allele-specific PCR were employed to identify the candidate gene for eggplant rind color through constructing a 2794 F₂ population generated from a cross between "BL01" (green pericarp) and "B1" (white pericarp). Genetic analysis of rind color revealed that a single dominant gene controls green color of eggplant peel. Pigment content measurement and cytological observations demonstrated that chlorophyll content and chloroplast number in BL01 were higher than in B1. A candidate gene (EGP19168.1) was fine-mapped to a 20.36 Kb interval on chromosome 8, which was predicted to encode the two-component response regulator-like protein Arabidopsis pseudo-response regulator2 (APRR2). Subsequently, allelic sequence analysis revealed that a SNP deletion (ACT→AT) in white-skinned eggplant led to a premature termination codon. Genotypic validation of 113 breeding lines using the Indel marker closely linked to SmAPRR2 could predict the skin color (green/white) trait with an accuracy of 92.9%. This study will be valuable for molecular marker-assisted selection in eggplant breeding and provides theoretical foundation for analyzing the formation mechanism of eggplant peel color.

Fine mapping and identifying candidate gene of Y underlying yellow peel in Cucurbita pepo

As a conspicuous trait, peel color is one of the most important characteristics that affects commodity quality and consumer preferences. The locus Y underlying yellow peel in Cucurbita pepo (zucchini) was first reported in 1922; however, its molecular mechanism is still unknown. In this study, a genetic analysis revealed that yellow peel is controlled by a single dominant genetic factor. Furthermore, Y was mapped in a ~170 kb region on chromosome 10 by bulked segregated analysis (BSA) and fine mapping in F₂ and BC₁ segregating populations. The candidate region harbors fifteen annotated genes, among which Cp4.1LG10g11560 (CpCHLH) is regarded as a promising candidate gene. CpCHLH encodes a magnesium chelatase H subunit involved in chlorophyll biosynthesis, and its mutation can result in a reduction in chlorophyll content and yellow phenotype. Interestingly, a large fragment (~15 kb) duplication containing incomplete CpCHLH was inserted in the candidate interval, resulting in two reformed CpCHLH proteins in the yellow parental line. It is most likely that the reformed CpCHLH proteins act as a malfunctional competitor of the normal CpCHLH protein to interrupt the formation of chlorophyll. Overall, the isolation of Y will shed light on the molecular mechanism of the peel color regulation of zucchini and lay a foundation for breeding.

A transcriptome analysis of Benincasa hispida revealed the pathways and genes involved in response to Phytophthora melonis infection

Wilt disease caused by Phytophthora melonis infection is one of the most serious threats to Benincasa hispida production. However, the mechanism of the response of B. hispida to a P. melonis infection remains largely unknown. In the present study, two B. hispida cultivars with different degrees of resistance to P. melonis were identified: B488 (a moderately resistant cultivar) and B214 (a moderately susceptible cultivar). RNA-seq was performed on P. melonis-infected B488 and B214 12 hours post infection (hpi). Compared with the control, 680 and 988 DEGs were respectively detected in B488 and B214. A KEGG pathway analysis combined with a cluster analysis revealed that phenylpropanoid biosynthesis, plant-pathogen interaction, the MAPK signaling pathway-plant, and plant hormone signal transduction were the most relevant pathways during the response of both B488 and B214 to P. melonis infection, as well as the differentially expressed genes in the two cultivars. In addition, a cluster analysis of transcription factor genes in DEGs identified four genes upregulated in B488 but not in B214 at 6 hpi and 12 hpi, which was confirmed by qRT-PCR. These were candidate genes for elucidating the mechanism of the B. hispida response to P. melonis infection and laying the foundation for the improvement of B. hispida.

Genome-Wide Identification of the SAUR Gene Family in Wax Gourd (Benincasa hispida) and Functional Characterization of BhSAUR60 during Fruit Development

The wax gourd (Benincasa hispida) is an important vegetable crop whose fruits contain nutrients and metabolites. Small auxin upregulated RNA (SAUR) genes constitute the largest early auxin-responsive gene family and regulate various biological processes in plants, although this gene family has not been studied in the wax gourd. Here, we performed genome-wide identification of the SAUR gene family in wax gourds and analyzed their syntenic and phylogenetic relationships, gene structures, conserved motifs, cis-acting elements, and expression patterns. A total of 68 SAUR (BhSAUR) genes were identified, which were distributed on nine chromosomes with 41 genes in two clusters. More than half of the BhSAUR genes were derived from tandem duplication events. The BhSAUR proteins were classified into seven subfamilies. BhSAUR gene promoters contained cis-acting elements involved in plant hormone and environmental signal responses. Further expression profiles showed that BhSAUR genes displayed different expression patterns. BhSAUR60 was highly expressed in fruits, and overexpression led to longer fruits in Arabidopsis. In addition, the plants with overexpression displayed longer floral organs and wavy stems. In conclusion, our results provide a systematic analysis of the wax gourd SAUR gene family and facilitate the functional study of BhSAUR60 during wax gourd fruit development.

Development and Application of InDel Markers Linked to Fruit-Shape and Peel-Colour Genes in Wax Gourd

The wax gourd is commonly grown in many countries because of its high nutritional and economic value. While the genes for the fruit shape and peel colour of wax gourd have been reported, the InDel markers linked to these genes remain undeveloped. In this study, the InDel markers linked to fruit-shape (Bch02G016830) and peel-colour (Bch05G003950) genes were developed from resequenced data. We used 120 inbred lines, 536 isolated populations, and 4 commercial hybrids to evaluate the validity and application value of the InDel markers. The accuracy rates of nine pairs of fruit-shape InDel markers (GX1-GX9) were 84.16–91.66% in 120 inbred lines. The accuracy rates of 27 pairs of peel-colour InDel markers (PS1-PS27) within approximately 3.0 Mb upstream and 3.0 Mb downstream of the peel-colour gene were 100% and those of 6 pairs of peel-colour InDel markers (PS28-PS33) within 3.0–20 Mb upstream and downstream of the peel-colour gene were 55.83–90% in 120 inbred lines. The purity of four commercial hybrids determined using GX1, GX2, PS13, and PS14 was highly consistent with the field results for purity determination. Our results provide important information for genetic linkage map construction, molecular-marker-assisted selective breeding, and purity determination of wax gourd hybrids.

Fine Mapping and Functional Analysis of Major Regulatory Genes of Soluble Solids Content in Wax Gourd (Benincasa hispida)

Soluble solids content (SSC) is an important quality trait of wax gourd, but reports about its regulatory genes are scarce. In this study, the SSC regulatory gene BhSSC2.1 in wax gourd was mined via quantitative trait locus (QTL) mapping based on high-density genetic mapping containing 12 linkage groups (LG) and bulked segregant analysis (BSA)-seq. QTL mapping and BSA-seq revealed for the first time that the SSC QTL (107.658–108.176 cM) of wax gourd was on Chr2 (LG2). The interpretable phenotypic variation rate and maximum LOD were 16.033% and 6.454, respectively. The QTL interval contained 13 genes. Real-time fluorescence quantitative expression analysis, functional annotation, and sequence analysis suggested that Bch02G016960, named BhSSC2.1, was a candidate regulatory gene of the SSC in wax gourd. Functional annotation of this gene showed that it codes for a NADP-dependent malic enzyme. According to BhSSC2.1 sequence variation, an InDel marker was developed for molecular marker-assisted breeding of wax gourd. This study will lay the foundation for future studies regarding breeding and understanding genetic mechanisms of wax gourd.

Genome-wide association analysis reveals a novel QTL CsPC1 for pericarp color in cucumber

BACKGROUND

Cucumber is an important melon crop in the world, with different pericarp colors. However, the candidate genes and the underlying genetic mechanism for such an important trait in cucumber are unknown. In this study, a locus controlling pericarp color was found on chromosome 3 of cucumber genome.

RESULTS

In this study, the light green inbred line G35 and the dark green inbred line Q51 were crossed to produce one F₂ population. Consequently, we identified a major locus CsPC1 (Pericarp color 1). Next, we mapped the CsPC1 locus to a 94-kb region chromosome 3 which contains 15 genes. Among these genes, Csa3G912920, which encodes a GATA transcription factor, was expressed at a higher level in the pericarp of the NIL-1334 line (with light-green pericarp) than in that of the NIL-1325 line (with dark-green pericarp). This study provides a new allele for the improvement of cucumber pericarp color.

CONCLUSION

A major QTL that controls pericarp color in cucumber, CsPC1, was identified in a 94-kb region that harbors the strong candidate gene CsGATA1.