Search

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32–25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Fusarium wilt is one of the most destructive and less controllable diseases in melon, which is usually caused byfusarium oxysporum. In this study, transcriptome sequencing and Yeast Two-Hybrid (Y2H) methods were used for quantification of differentially expressed genes (DEGs) involved infusarium oxysporum(f. sp.melonisrace 1) stress-induced mechanisms in contrasted melon varieties (M4-45 “susceptible” and MR-1 “resistant”). The interaction factors ofFom-2resistance genes were also explored in response to the plant-pathogen infection mechanism. Transcriptomic analysis exhibited total 1,904 new genes; however, candidate DEGs analysis revealed a total of 144 specific genes (50 upregulated and 94 downregulated) for M4-45 variety and 104 specific genes (71 upregulated and 33 downregulated) for MR-1 variety, respectively. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway depicted some candidate DEGs, including Phenylalanine metabolism, phenylpropane biosynthesis, plants-pathogen interaction, and signal transduction of plant hormones, which were mainly involved in disease resistance metabolic pathways. The weighted gene co-expression network analysis (WGCNA) analysis revealed a strong correlation module and exhibited the disease resistance-related genes encoding course proteins, transcription factors, protein kinase, benzene propane biosynthesis path, plants-pathogen interaction pathway, and glutathione S-transferase. Meanwhile, the resistance-related specific genes expression was relatively abundant in MR-1 compared to the M4-45, and cell wall-associated receptor kinases (MELO3C008452andMELO3C008453), heat shock protein (Cucumis_melo_newGene_172), defensin-like protein (Cucumis_melo_newGene_5490), and disease resistance response protein (MELO3C016325), activator response protein (MELO3C021623), leucine-rich repeat receptor protein kinase (MELO3C024412), lactyl glutathione ligase (Cucumis_melo_newGene_36), and unknown protein (MELO3C007588) were persisted by exhibiting the upregulated expressions. At the transcription level, the interaction factors between the candidate genes in response to thefusarium oxysporuminduced stress, and Y2H screening signified the main contribution of MYB transcription factors (MELO3C009678andMELO3C014597), BZIP (MELO3C011839andMELO3C019349), unknown proteins, and key enzymes in the ubiquitination process (4XM334FK014). The candidate genes were further verified in exogenously treated melon plants withf. oxysporum(Fom-2, Race 1), Abscisic acid (ABA), Methyl Jasmonite (MeJA), and Salicylic acid (SA), using the fluorescence quantitative polymerase chain reaction (qRT-PCR) analysis. The overall expression results indicated that the SA signal pathway is involved in effective regulation of theFom-2gene activity.

Melon is an important Cucurbitaceae crop. Field observations had shown that the green stigmas of melon are more attractive to pollinators than yellow stigmas. In this study, F2 and F2:3 populations obtained by crossing MR-1 (green stigma) and M4-7 (yellow stigma) were used for genetic analysis and mapping. A genetic map of 1,802.49 cm was constructed with 116 cleaved amplified polymorphism sequence (CAPS) markers. Two stable quantitative trait loci (QTLs) linked to the trait of stigma color were identified on chromosomes 2 (SC2.1) and 8 (SC8.1), respectively. An expanded F2 population was used to narrow down the confidence regions of SC2.1 and SC8.1. As a result, SC2.1 was further mapped to a 3.6 cm region between CAPS markers S2M3 and S2B1-3, explaining 9.40% phenotypic variation. SC8.1 was mapped to a 3.7-cm region between CAPS markers S8E7 and S8H-1, explaining 25.92% phenotypic variation. This study broadens our understanding of the mechanisms of stigma color regulation and will be of benefit to the breeding of melon.

Fruit pedicel (FP) is an important determinant of premium fruit quality that directly affects commercial market value. However, in-depth molecular and genetic basis of pedicel-related traits has not been identified in watermelon. Herein, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic regions controlling FP traits based on newly derived whole-genome single nucleotide polymorphism based cleaved amplified polymorphism sequence (SNP-CAPS) markers. Next-generation sequencing based whole-genome re-sequencing of two watermelon parent lines revealed 98.30 and 98.40% of average coverage, 4,989,869 SNP variants, and 182,949 CAPS loci pairs across the reference genome, respectively. A total of 221 sets of codominant markers exhibited 46.42% polymorphism rate and were effectively genotyped within 100-F2:3 derived mapping population. The developed linkage map covered a total of 2,630.49 cM genetic length with averaged 11.90 cM, and depicted a valid marker-trait association. In total, 6 QTLs (qFPL4.1, qFPW4.1, qFPD2.1, qFPD2.2, qFPD8.1, qFPD10.1) were mapped with five major effects and one minor effect between the whole genome adjacent markers positioned over distinct chromosomes (02, 04, 08, 10), based on the ICIM-ADD mapping approach. These significant QTLs were similarly mapped in delimited flanking regions of 675.10, 751.38, 859.24, 948.39, and 947.51 kb, which collectively explained 8.64–13.60% PVE, respectively. A highly significant and positive correlation was found among the observed variables. To our knowledge, we first time reported the mapped QTLs/genes affecting FP traits of watermelon, and our illustrated outcomes will deliver the potential insights for fine genetic mapping as well as functional gene analysis through MAS-based breeding approaches.

The melon fruit surface groove (fsg) not only affects peel structure and causes stress-induced fruit cracking but also fits consumers’ requirements in different regions. In this study, genetic inheritance analysis of three F2 populations derived from six parental lines revealed that the fsg trait is controlled by a simple recessive inherited gene. Through bulked segregant analysis sequencing (BSA-seq), the Cmfsg locus was detected in an 8.96 Mb interval on chromosome 11 and then initially mapped to a region of approximately 1.15 Mb. Further fine mapping with a large F2 population including 1,200 plants narrowed this region to 207 kb containing 11 genes. A genome-wide association study (GWAS) with 187 melon accessions also produced the same chromosome region for the Cmfsg locus. Due to the rare molecular markers and lack of mutations in the coding and promoter regions of the 11 candidate genes in the fine-mapped interval, we conducted in silico BSA to explore the natural melon panel to predict candidate genes for the Cmfsg locus. A 1.07 kb segment upstream of MELO3C019694.2 (annotated as the AGAMOUS MADS-box transcription factor) exhibited a correlation with the grooved and non-grooved accessions among the F2 individuals, and a natural panel consisted of 17 melon accessions. The expression level of MELO3C019694.2 in the pericarp was higher in grooved lines than in non-grooved lines and was specifically expressed in fruit compared with other tissues (female flower, male flower, root, and leaf). This work provides fundamental information for further research on melon fsg trait formation and molecular markers for melon breeding.

The three-seeded pod number is an important trait that positively influences soybean yield. Soybean variety with increased three-seeded pod number contributes to the seed number/plant and higher yield. The candidate genes of the three-seeded pod may be the key for improving soybean yield. In this study, identification and validation of candidate genes for three-seeded pod has been carried out. First, a total of 36 quantitative trait locus (QTL) were detected from the investigation of recombinant inbred lines including 147 individuals derived from a cross between Charleston and Dongning 594 cultivars. Five consensus QTLs were integrated. Second, an introgressed line CSSL-182 carrying the target segment for the trait from the donor parent was selected to verify the consensus QTL based on its phenotype. Third, a secondary group was constructed by backcrossing with CSSL-182, and two QTLs were confirmed. There were a total of 162 genes in the two QTLs. The mining of candidate genes resulted in the annotation of eight genes with functions related to pod and seed sets. Finally, haplotype analysis and quantitative reverse transcriptase real-time PCR were carried to verify the candidate genes. Four of these genes had different haplotypes in the resource group, and the differences in the phenotype were highly significant. Moreover, the differences in the expression of the four genes during pod and seed development were also significant. These four genes were probably related to the development process underlying the three-seeded pod in soybean. Herein, we discuss the past and present studies related to the three-seeded pod trait in soybean.

Watermelon fruit texture and quality are determined by flesh firmness. As a quality trait, flesh firmness is controlled by multigenes. Defining the key regulatory factors of watermelon flesh firmness is of great significance for watermelon genetic breeding. In this study, the hard-flesh egusi seed watermelon PI186490 was used as the male parent, the soft-flesh cultivated watermelon W1-1 was used as the female parent, and 175 F2 generations were obtained from selfing F1. Primary mapping of the major genes controlling center flesh firmness was achieved by bulked-segregant analysis (BSA)-Seq analysis and molecular marker technology. Finally, major genes were delimited in the physical interval between 6,210,787 and 7,742,559 bp on chromosome 2 and between 207,553 and 403,137 bp on chromosome 8. The content of each cell wall component and hormone was measured, and comparative transcriptome analysis was performed during fruit development in watermelon. The protopectin, cellulose, hemicellulose, indole-3-acetic acid (IAA) and abscisic acid (ABA) contents were measured, and paraffin sections were made during the three fruit developmental stages. The results revealed that protopectin, celluloses, and hemicelluloses exhibited similar trends for flesh firmness, while the IAA and ABA concentrations continued to decrease with fruit ripening. Paraffin sections showed that PI186490 cells were more numerous, were more tightly packed, had clearer cell wall edges and had thicker cell walls than W1-1 cells at every developmental stage. Comparative transcriptome analysis was conducted on RNA samples of flesh during fruit development and ripening in W1-1 and PI186490. The results from the localization interval transcriptome analysis showed that Cla016033 (DUF579 family member), which may influence the cell wall component contents to adjust the flesh firmness in watermelon fruit, was different in W1-1 and PI186490 and that Cla012507 (MADS-box transcription factor) may be involved in the regulation of fruit ripening and affect the hardness of watermelon fruit.

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

Melon is an important agricultural and economic vegetable crop worldwide. The genetic male sterility mutant (ms-5) has a recessive nuclear gene that controls the male sterility germplasm. Male sterility could reduce the cost of F1 seed production in melon, but heterozygous fertile plants should be removed before pollination. In this study, bulked segregant analysis combined with specific length amplified fragment sequencing was applied to map the single nuclear male sterility recessive gene. A 30-kb candidate region on chromosome 9 located on scaffold 000048 and spanning 2,522,791 to 2,555,104 bp was identified and further confirmed by cleavage amplified polymorphic sequence markers based on parental line resequencing data and classical mapping of 252 F2 individuals. Gene prediction indicated that six annotated genes are present in the 30-kb candidate region. Quantitative RT-PCR revealed significant differences in the expression level of the LOC103498166 ABORTED MICROSPORES (AMS) gene in male-sterile lines (ms-5) and male-fertile (HM1-1) lines during the 2-mm (tetrad) and 5-mm (the first pollen mitosis) periods, and negative regulation of the AMS candidate gene transcription factor was also detected. Sequencing and cluster analysis of the AMS transcription factor revealed five single-nucleotide polymorphisms between the parental lines. The data presented herein suggest that the AMS transcription factor is a possible candidate gene for single nuclear male sterility in melon. The results of this study will help breeders to identify male-sterile and -fertile plants at seeding as marker-assisted selection methods, which would reduce the cost of seed production and improve the use of male-sterile lines in melon.