Identification of suitable reference genes for gene expression normalization in qRT-PCR analysis in watermelon

GRAS family member LATERAL SUPPRESSOR regulates the initiation and morphogenesis of watermelon lateral organs

The lateral organs of watermelon (Citrullus lanatus), including lobed leaves, branches, flowers, and tendrils, together determine plant architecture and yield. However, the genetic controls underlying lateral organ initiation and morphogenesis remain unclear. Here, we found that knocking out the homologous gene of shoot branching regulator LATERAL SUPPRESSOR in watermelon (ClLs) repressed the initiation of branches, flowers, and tendrils and led to developing round leaves, indicating that ClLs undergoes functional expansion compared with its homologs in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and tomato (Solanum lycopersicum). Using ClLs as the bait to screen against the cDNA library of watermelon, we identified several ClLs-interacting candidate proteins, including TENDRIL (ClTEN), PINOID (ClPID), and APETALA1 (ClAP1). Protein-protein interaction assays further demonstrated that ClLs could directly interact with ClTEN, ClPID, and ClAP1. The mRNA in situ hybridization assay revealed that the transcriptional patterns of ClLs overlapped with those of ClTEN, ClPID, and ClAP1 in the axillary meristems and leaf primordia. Mutants of ClTEN, ClPID, and ClAP1 generated by the CRISPR/Cas9 gene editing system lacked tendrils, developed round leaves, and displayed floral diapause, respectively, and all these phenotypes could be observed in ClLs knockout lines. Our findings indicate that ClLs acts as lateral organ identity protein by forming complexes with ClTEN, ClPID, and ClAP1, providing several gene targets for transforming the architecture of watermelon.

Comprehensive Characterization of B-Box Zinc Finger Genes in Citrullus lanatus and Their Response to Hormone and Abiotic Stresses

Plant B-BOX (BBX) zinc finger transcription factors play crucial roles in growth and development and the stress response. Although the BBX family has been characterized in various plants, systematic analysis in watermelon is still lacking. In this study, 25 watermelon ClBBX genes were identified. ClBBXs were grouped into five clades (Clade I, II, III, IV, and V) based on their conserved domains and phylogenetic relationships. Most of the ClBBXs (84%) might be localized in the nuclei or cytoplasm. The classification of ClBBXs was consistent with their gene structures. They were unevenly distributed in nine chromosomes except for Chr4 and Chr10, with the largest number of six members in Chr2. Segmental duplications were the major factor in ClBBX family expansion. Some BBXs of watermelon and Arabidopsis evolved from a common ancestor. In total, 254 hormonal and stress-responsive cis elements were discovered in ClBBX promoters. ClBBXs were differentially expressed in tissues, and the expression levels of ClBBX15 and 16 were higher in aboveground tissues than in roots, while the patterns of ClBBX21a, 21b, 21c, 28 and 30b were the opposite. With salicylic acid, methyl jasmonate and salt stress conditions, 17, 18 and 18 ClBBXs exhibited significant expression changes, respectively. In addition, many ClBBXs, including ClBBX29b, 30a and 30b, were also responsive to cold and osmotic stress. In summary, the simultaneous response of multiple ClBBXs to hormonal or abiotic stress suggests that they may have functional interactions in the stress hormone network. Clarifying the roles of key ClBBXs in transcriptional regulation and mediating protein interactions will be an important task. Our comprehensive characterization of the watermelon ClBBX family provides vital clues for the in-depth analysis of their biological functions in stress and hormone signaling pathways.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Comparative bioinformatics analysis and abiotic stress responses of expansin proteins in Cucurbitaceae members: watermelon and melon

Watermelon and melon are members of the Cucurbitaceae family including economically significant crops in the world. The expansin protein family, which is one of the members of the cell wall, breaks down the non-covalent bonds between cell wall polysaccharides, causing pressure-dependent cell expansion. Comparative bioinformatics and molecular characterization analysis of the expansin protein family were carried out in the watermelon (Citrullus lanatus) and melon (Cucumis melo) plants in the study. Gene expression levels of expansin family members were analyzed in leaf and root tissues of watermelon and melon under ABA, drought, heat, cold, and salt stress conditions by quantitative real-time PCR analysis. After comprehensive searches, 40 expansin proteins (22 ClaEXPA, 14 ClaEXPLA, and 4 ClaEXPB) in watermelon and 43 expansin proteins (19 CmEXPA, 15 CmEXPLA, 3 CmEXPB, and 6 CmEXPLB) in melon were identified. The greatest orthologous genes were identified with soybean expansin genes for watermelon and melon. However, the latest divergence time between orthologous genes was determined with poplar expansin genes for watermelon and melon expansin genes. ClaEXPA-04, ClaEXPA-09, ClaEXPB-01, ClaEXPB-03, and ClaEXPLA-13 genes in watermelon and CmEXPA-12, CmEXPA-10, and CmEXPLA-01 genes in melon can be involved in tissue development and abiotic stress response of the plant. The current study combining bioinformatics and experimental analysis can provide a detailed characterization of the expansin superfamily which has roles in growth and reaction to the stress of the plant. The study ensures detailed data for future studies examining gene functions including the roles in plant growth and stress conditions.

Genetic mapping and molecular characterization of the delayed green gene dg in watermelon (Citrullus lanatus)

Leaf color mutants are common in higher plants that can be used as markers in crop breeding and are important tools in understanding regulatory mechanisms of chlorophyll biosynthesis and chloroplast development. Genetic analysis was performed by evaluating F₁, F₂ and BC₁ populations derived from two parental lines (Charleston gray with green leaf color and Houlv with delayed green leaf color), suggesting that a single recessive gene controls the delayed green leaf color. In this study, the delayed green mutant showed a conditional pale green leaf color at the early leaf development but turned to green as the leaf development progressed. Delayed green leaf plants showed reduced pigment content, photosynthetic, chlorophyll fluorescence parameters, and impaired chloroplast development compared with green leaf plants. The delayed green (dg) locus was mapped to 7.48 Mb on chromosome 3 through bulk segregant analysis approach, and the gene controlling delayed green leaf color was narrowed to 53.54 kb between SNP130 and SNP135 markers containing three candidate genes. Sequence alignment of the three genes indicated that there was a single SNP mutation (G/A) in the coding region of ClCG03G010030 in the Houlv parent, which causes an amino acid change from Arginine to Lysine. The ClCG03G010030 gene encoded FtsH extracellular protease protein family is involved in early delayed green leaf development. The expression level of ClCG03G010030 was significantly reduced in delayed green leaf plants than in green leaf plants. These results indicated that the ClCG03G010030 might control watermelon green leaf color and the single SNP variation in ClCG03G010030 may result in early delayed green leaf color development during evolutionary process.

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

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

Micronutrients Affect Expression of Induced Resistance Genes in Hydroponically Grown Watermelon against Fusarium oxysporum f. sp. niveum and Meloidogyne incognita

The soil-borne pathogens, particularly Fusarium oxysporum f. sp. niveum (FON) and southern root-knot nematode (RKN, Meloidogyne incognita) are the major threats to watermelon production in the southeastern United States. The role of soil micronutrients on induced resistance (IR) to plant diseases is well-documented in soil-based media. However, soil-based media do not allow us to determine the contribution of individual micronutrients in the induction of IR. In this manuscript, we utilized hydroponics-medium to assess the effect of controlled application of micronutrients, including iron (Fe), manganese (Mn), and zinc (Zn) on the expression of important IR genes (PR1, PR5, and NPR1 from salicylic acid (SA) pathway, and VSP, PDF, and LOX genes from jasmonic acid (JA) pathway) in watermelon seedlings upon inoculation with either FON or RKN or both. A subset of micronutrient-treated plants was inoculated (on the eighth day of micronutrient application) with FON and RKN (single or mixed inoculation). The expression of the IR genes in treated and control samples was evaluated using qRT-PCR. Although, significant phenotypic differences were not observed with respect to the severity of wilt symptoms or RKN galling with any of the micronutrient treatments within the 30-day experimental period, differences in the induction of IR genes were considerably noticeable. However, the level of gene expression varied with sampling period, type and concentration of micronutrients applied, and pathogen inoculation. In the absence of pathogens, micronutrient applications on the seventh day, in general, downregulated the expression of the majority of the IR genes. However, pathogen inoculation preferentially either up- or down-regulated the expression levels of the IR genes at three days post-inoculation depending on the type and concentration of micronutrients. The results demonstrated here indicate that micronutrients in watermelon may potentially make watermelon plants susceptible to infection by FON and RKN. However, upon infection the IR genes are significantly up-regulated that they may potentially aid the prevention of further infection via SA- and JA-pathways. This is the first demonstration of the impact of micronutrients affecting IR in watermelon against FON and RKN infection.

Selection and Validation of Reference Genes for qRT-PCR Gene Expression Analysis in Kengyilia melanthera

Kengyilia is a newly established genus. Most species in this genus survive in hash environment, which might be an indicator of an acquirement of stress resistance genes and the potential for molecular breeding in Triticeae species. Quantitative real-time PCR (qRT-PCR) is a widely used technique with varied sensitivity heavily dependent on the optimal level of the reference genes. K. melanthera is a typical psammophyte species which has high drought resistance. The reference genes of K. melanthera are not yet reported. This study aims to evaluate the expression stability of 14 candidate reference genes (EF1A, GAPDH, ACT1, UBI, TUBB3, TIPRL, CACS, PPP2R1B, TUBA1A, EIF4A1, CYPA3, TCTP, ABCG11L, and FBXO6L) under five treatments (drought, heat, cold, salt, and ABA) and find the most stable and suitable one even upon stressed conditions. The software NormFinder, GeNorm, BestKeeper, and RefFinder were used for data analysis. In general, the genes CACS and PPP2R1B are concluded to have the best overall performance under the various treatments. With the ABA treatment, TCTP and TIPRL show the best stability. CACS and TCTP, as well as TIPRL and CYPA3, were most stable under the treatments of cold and salt, respectively. CACS and FBXO6L were ranked the highest with the heat treatment and drought treatment, respectively. Finally, the Catalase-1 (CAT1) gene was used to verify the reliability of the above reference genes. Accordingly, CAT1’s expression pattern remained unchanged after normalization with stable reference genes. This study provides beneficial information about the stability and reliability of potential reference genes for qRT-PCR in K. melanthera.

Identification and Characterization Roles of Phytoene Synthase (PSY) Genes in Watermelon Development

Phytoene synthase (PSY) plays an essential role in carotenoid biosynthesis. In this study, three ClPSY genes were identified through the watermelon genome, and their full-length cDNA sequences were cloned. The deduced proteins of the three ClPSY genes were ranged from 355 to 421 amino acid residues. Phylogenetic analysis suggested that the ClPSYs are highly conserved with bottle gourd compared to other cucurbit crops PSY proteins. Variation in ClPSY1 expression in watermelon with different flesh colors was observed; ClPSY1 was most highly expressed in fruit flesh and associated with the flesh color formation. ClPSY1 expression was much lower in the white-fleshed variety than the colored fruits. Gene expression analysis of ClPSY genes in root, stem, leaf, flower, ovary and flesh of watermelon plants showed that the levels of ClPSY2 transcripts found in leaves was higher than other tissues; ClPSY3 was dominantly expressed in roots. Functional complementation assays of the three ClPSY genes suggested that all of them could encode functional enzymes to synthesize the phytoene from Geranylgeranyl Pyrophosphate (GGPP). Some of the homologous genes clustered together in the phylogenetic tree and located in the synteny chromosome region seemed to have similar expression profiles among different cucurbit crops. The findings provide a foundation for watermelon flesh color breeding with regard to carotenoid synthesis and also provide an insight for the further research of watermelon flesh color formation.

Selection of suitable candidate genes for mRNA expression normalization in bulbil development of Pinellia ternata

Pinellia ternata (Thunb.) Breit. (Abbreviated as P. ternata). It is a commonly prescribed Chinese traditional medicinal herb for the treatment of phlegm, cough, and morning sick. Bulbil reproduction is one of the main reproductive methods of P. ternata. The accurate quantification of gene expression patterns associated with bulbil development might be helpful to explore the molecular mechanism involved in P. ternata reproduction. Quantitative real-time PCR was the most preferred method for expression profile and function analysis of mRNA. However, the reference genes in different tissues of P. ternata in different periods of bulbil development have not been studied in detail. In present study, the expression stability of eight candidate reference genes were determined with programs: geNorm, NormFinder, BestKeeper, and refFinder. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was identified as the top- rated reference gene in all samples of P. ternata, while different combinations of reference gene proved to be the most stable depending on development stage and tissue type. Furthermore, the reliability of GAPDH expression was verified by six P. ternata related genes in hormone and nutrient biosynthesis pathways, and the expression profiles of these genes were agreed with the results of RNA-seq digital gene expression analysis. These results can contribute to studies of gene expression patterns and functional analysis of P. ternata involved in bulbil development.

Transcriptome Profiling to Dissect the Role of Genome Duplication on Graft Compatibility Mechanisms in Watermelon

Watermelon (Citrullus lanatus) is a popular crop worldwide. Compared to diploid seeded watermelon, triploid seedless watermelon cultivars are in great demand. Grafting in triploid and tetraploid watermelon produces few seedlings. To learn more about how genome duplication affects graft compatibility, we compared the transcriptomes of tetraploid and diploid watermelons grafted on squash rootstock using a splicing technique. WGCNA was used to compare the expression of differentially expressed genes (DEGs) between diploid and tetraploid watermelon grafted seedlings at 0, 3, and 15 days after grafting (DAG). Only four gene networks/modules correlated significantly with phenotypic characteristics. We found 11 genes implicated in hormone, AOX, and starch metabolism in these modules based on intramodular significance and RT-qPCR. Among these genes, two were linked with IAA (r² = 0.81), one with ZR (r² = 0.85) and one with POD (r² = 0.74). In the MElightsteelblue1 module, Cla97C11G224830 gene was linked with CAT (r² = 0.81). Two genes from the MEivory module, Cla97C07G139710 and Cla97C04G077300, were highly linked with SOD (r² = 0.72). Cla97C01G023850 and Cla97C01G006680 from the MEdarkolivegreen module were associated with sugars and starch (r² = 0.87). Tetraploid grafted seedlings had higher survival rates and hormone, AOX, sugar, and starch levels than diploids. We believe that compatibility is a complicated issue that requires further molecular research. We found that genome duplication dramatically altered gene expression in the grafted plants' IAA and ZR signal transduction pathways and AOX biosynthesis pathways, regulating hormone levels and improving plant survival.

Grafting Watermelon Onto Pumpkin Increases Chilling Tolerance by Up Regulating Arginine Decarboxylase to Increase Putrescine Biosynthesis

Low temperature is a major environmental factor that severely impairs plant growth and productivity. Watermelon (Citrullus lanatus) is a chilling-sensitive crop. Grafting of watermelon onto pumpkin rootstock is an effective technique to increase the chilling tolerance of watermelon when exposure to short-time chilling stress. However, the mechanism by which pumpkin rootstock increases chilling tolerance remains poorly understood. Under 10°C/5°C (day/night) chilling stress treatment, pumpkin-grafted watermelon seedlings showed higher chilling tolerance than self-grafted watermelon plants with significantly reduced lipid peroxidation and chilling injury (CI) index. Physiological analysis revealed that pumpkin rootstock grafting led to the notable accumulation of putrescine in watermelon seedlings under chilling conditions. Pre-treat foliar with 1 mM D-arginine (inhibitor of arginine decarboxylase, ADC) increased the electrolyte leakage (EL) of pumpkin-grafted watermelon leaves under chilling stress. This result can be ascribed to the decrease in transcript levels of ADC, ornithine decarboxylase, spermidine synthase, and polyamine oxidase genes involved in the synthesis and metabolism of polyamines. Transcriptome analysis showed that pumpkin rootstock improved chilling tolerance in watermelon seedlings by regulating differential gene expression under chilling stress. Pumpkin-grafted seedling reduced the number and expression level of differential genes in watermelon scion under chilling stress. It specifically increased the up-regulated expression of ADC (Cla97C11G210580), a key gene in the polyamine metabolism pathway, and ultimately promoted the accumulation of putrescine. In conclusion, pumpkin rootstock grafting increased the chilling tolerance of watermelon through transcription adjustments, up regulating the expression level of ADC, and promoting the synthesis of putrescine, which ultimately improved the chilling tolerance of pumpkin-grafted watermelon plants.

Pistachio hull polysaccharide protects Nile tilapia against LPS-induced excessive inflammatory responses and oxidative stress, possibly via TLR2 and Nrf2 signaling pathways

Polysaccharides are polymeric carbohydrates found in living organisms, which have several physiological functions. In the present study, Nile tilapia (Oreochromis niloticus) were fed diets containing three levels (0%, 0.2%, and 0.6%) of Pistacia vera hull polysaccharide (PHP) for 45 days and then injected with PBS or bacterial lipopolysaccharide (LPS). Before the LPS challenge, Nile tilapia fed 0.2% and 0.6% PHP showed significantly increased mean final weight and weight gain compared to those received 0% PHP. The specific growth rate and feed conversion ratio were significantly improved in the treatment fed 0.6% PHP compared to the remaining groups. After LPS challenge, the activities of liver antioxidant enzymes, including superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase exhibited the highest values in the 0.6% PHP group. Malondialdehyde (MDA) levels were significantly augmented in the model (fed 0% PHP diet and injected with LPS) and 0.2% PHP groups compared to the control. However, MDA showed decreased levels in the 0.6% PHP group. LPS induced higher mRNA and/or protein levels of Toll-like receptor 2 (TLR2), nuclear factor kappa B (NF-κB), myeloid differentiation primary response protein 88 (Myd88), tumor necrosis factor α (TNF-α), interleukin 1β (IL-1β), and interferon γ (IFN-γ) in Nile tilapia liver. However, PHP administration significantly upregulated the expression of interleukin 10 (IL-10), nuclear erythroid 2-related factor 2 (Nrf2), SOD, and CAT, but markedly suppressed TLR2, NF-κB, Myd88, and pro-inflammatory cytokine expression and/or production in the liver. The findings of the current study indicated that PHP has positive effects on growth performance, immune gene-related expression, and antioxidative activities. We can conclude that PHP can attenuate LPS-induced oxidative stress and inflammatory responses in vivo, possibly via induction of Nrf2 and blockade of TLR2/Myd88/NF-κB pathways in Nile tilapia.

Comparative genomic analysis of expansin superfamily gene members in zucchini and cucumber and their expression profiles under different abiotic stresses

Zucchini and cucumber belong to the Cucurbitaceae family, a group of economical and nutritious food plants that is consumed worldwide. Expansin superfamily proteins are generally localized in the cell wall of plants and are known to possess an effect on cell wall modification by causing the expansion of this region. Although the whole genome sequences of cucumber and zucchini plants have been resolved, the determination and characterization of expansin superfamily members in these plants using whole genomic data have not been implemented yet. In the current study, a genome-wide analysis of zucchini (Cucurbita pepo) and cucumber (Cucumis sativus) genomes was performed to determine the expansin superfamily genes. In total, 49 and 41 expansin genes were identified in zucchini and cucumber genomes, respectively. All expansin superfamily members were subjected to further bioinformatics analysis including gene and protein structure, ontology of the proteins, phylogenetic relations and conserved motifs, orthologous relations with other plants, targeting miRNAs of those genes and in silico gene expression profiles. In addition, various abiotic stress responses of zucchini and cucumber expansin genes were examined to determine their roles in stress tolerance. CsEXPB-04 and CsEXPA-11 from cucumber and CpEXPA-20 and CpEXPLA-14 from zucchini can be candidate genes for abiotic stress response and tolerance in addition to their roles in the normal developmental processes, which are supported by the gene expression analysis. This work can provide new perspectives for the roles of expansin superfamily genes and offers comprehensive knowledge for future studies investigating the modes of action of expansin proteins.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01108-w.

Transcriptome-wide N6-methyladenosine (m6A) methylation in watermelon under CGMMV infection

BACKGROUND

Cucumber green mottle mosaic virus (CGMMV) causes substantial global losses in cucurbit crops, especially watermelon. N6-methyladenosine (m⁶A) methylation in RNA is one of the most important post-transcriptional modification mechanisms in eukaryotes. It has been shown to have important regulatory functions in some model plants, but there has been no research regarding m⁶A modifications in watermelon.

RESULTS

We measured the global m⁶A level in resistant watermelon after CGMMV infection using a colorimetric method. And the results found that the global m⁶A level significantly decreased in resistant watermelon after CGMMV infection. Specifically, m⁶A libraries were constructed for the resistant watermelon leaves collected 48 h after CGMMV infection and the whole-genome m⁶A-seq were carried out. Numerous m⁶A modified peaks were identified from CGMMV-infected and control (uninfected) samples. The modification distributions and motifs of these m⁶A peaks were highly conserved in watermelon transcripts but the modification was more abundant than in other reported crop plants. In early response to CGMMV infection, 422 differentially methylated genes (DMGs) were identified, most of which were hypomethylated, and probably associated with the increased expression of watermelon m⁶A demethylase gene ClALKBH4B. Gene Ontology (GO) analysis indicated quite a few DMGs were involved in RNA biology and stress responsive pathways. Combined with RNA-seq analysis, there was generally a negative correlation between m⁶A RNA methylation and transcript level in the watermelon transcriptome. Both the m⁶A methylation and transcript levels of 59 modified genes significantly changed in response to CGMMV infection and some were involved in plant immunity.

CONCLUSIONS

Our study represents the first comprehensive characterization of m⁶A patterns in the watermelon transcriptome and helps to clarify the roles and regulatory mechanisms of m⁶A modification in watermelon in early responses to CGMMV.

Selection and validation of reference genes for qRT-PCR analysis during fruit ripening of red pitaya (Hylocereus polyrhizus)

Red pitaya (Hylocereus polyrhizus) is widely cultivated in southern and southwestern China. To provide a basis for studying the molecular mechanisms of the ripening of this fruit, we carried out RNA sequencing (RNA-seq) analysis to identify differentially and stably expressed unigenes. The latter may serve as a resource of potential reference genes for normalization of target gene expression determined using quantitative real-time PCR (qRT-PCR). We selected 11 candidate reference genes from our RNA-seq analysis of red pitaya fruit ripening (ACT7, EF-1α, IF-4α, PTBP, PP2A, EF2, Hsp70, GAPDH, DNAJ, TUB and CYP), as well as β-ACT, which has been used as a reference gene for pitayas in previous studies. We then comprehensively evaluated their expression stability during fruit ripening using four statistical methods (GeNorm, NormFinder, BestKeeper and DeltaCt) and merged the four outputs using the online tool RefFinder for the final ranking. We report that PTBP and DNAJ showed the most stable expression patterns, whereas CYP and ACT7 showed the least stable expression patterns. The relative gene expression of red pitaya sucrose synthase and 4, 5-dihydroxyphenylalanine-extradiol-dioxygenase as determined by quantitative real-time PCR and normalized to PTBP and DNAJ was consistent with the RNA-seq results, suggesting that PTBP and DNAJ are suitable reference genes for studies of red pitaya fruit ripening.

Comparative Transcriptome Analysis in Homo- and Hetero-Grafted Cucurbit Seedlings

Watermelon (Citrullus lanatus) is a valuable horticultural crop with nutritional benefits grown worldwide. It is almost exclusively cultivated as grafted scions onto interspecific squash rootstock (Cucurbita maxima × Cucurbita moschata) to improve the growth and yield and to address the problems of soilborne diseases and abiotic stress factors. This study aimed to examine the effect of grafting (homo- and hetero-grafting) on the transcriptome level of the seedlings. Therefore, we compared homo-grafted watermelon (WW) with non-grafted watermelon control (W), homo-grafted squash (SS) with non-grafted squash control (S), hetero-grafted watermelon onto squash (WS) with SS, and WS with WW. Different numbers of differentially expressed genes (DEGs) were identified in each comparison. In total, 318 significant DEGs were detected between the transcriptomes of hetero-grafts and homo-grafts at 16 h after grafting. Overall, a significantly higher number of downregulated transcripts was detected among the DEGs. Only one gene showing increased expression related to the cytokinin synthesis was common in three out of four comparisons involving WS, SS, and S. The highest number of differentially expressed (DE) transcripts (433) was detected in the comparison between SS and S, followed by the 127 transcripts between WW and W. The study provides a description of the transcriptomic nature of homo- and hetero-grafted early responses, while the results provide a start point for the elucidation of the molecular mechanisms and candidate genes for the functional analyses of hetero-graft and homo-graft systems in Cucurbitaceae and generally in the plants.

The role of watermelon caffeic acid O-methyltransferase (ClCOMT1) in melatonin biosynthesis and abiotic stress tolerance

Melatonin is a pleiotropic signaling molecule that regulates plant growth and responses to various abiotic stresses. The last step of melatonin synthesis in plants can be catalyzed by caffeic acid O-methyltransferase (COMT), a multifunctional enzyme reported to have N-acetylserotonin O-methyltransferase (ASMT) activity; however, the ASMT activity of COMT has not yet been characterized in nonmodel plants such as watermelon (Citrullus lanatus). Here, a total of 16 putative O-methyltransferase (ClOMT) genes were identified in watermelon. Among them, ClOMT03 (Cla97C07G144540) was considered a potential COMT gene (renamed ClCOMT1) based on its high identities (60.00-74.93%) to known COMT genes involved in melatonin biosynthesis, expression in almost all tissues, and upregulation under abiotic stresses. The ClCOMT1 protein was localized in the cytoplasm. Overexpression of ClCOMT1 significantly increased melatonin contents, while ClCOMT1 knockout using the CRISPR/Cas-9 system decreased melatonin contents in watermelon calli. These results suggest that ClCOMT1 plays an essential role in melatonin biosynthesis in watermelon. In addition, ClCOMT1 expression in watermelon was upregulated by cold, drought, and salt stress, accompanied by increases in melatonin contents. Overexpression of ClCOMT1 enhanced transgenic Arabidopsis tolerance against such abiotic stresses, indicating that ClCOMT1 is a positive regulator of plant tolerance to abiotic stresses.

Permissive action of H2O2 mediated ClUGT75 expression for auxin glycosylation and Al3+- tolerance in watermelon

Although Al³⁺-toxicity is one of the limiting factors for crop production in acidic soils, little is known about the Al³⁺-tolerance mechanism in watermelon, a fairly acid-tolerant crop. This work aimed to identify the interaction between the H₂O₂ scavenging pathway and auxin glycosylation relevant to watermelon Al³⁺-tolerance. By analyzing expressions of hormone-related ClUGTs and antioxidant enzyme genes in Al³⁺-tolerant (ZJ) and Al³⁺-sensitive (NBT) cultivars, we identified ClUGT75s (B1, B2, and D1) and ClSOD1-2-ClCAT as crucial components associated with Al³⁺-tolerance. Al³⁺-stress significantly increased H₂O₂ content by 92.7% in NBT and 42.3% in ZJ, accompanied by less Al³⁺-, auxin (IAA and IBA), and MDA contents in ZJ than NBT. These findings coincided with significant ClSOD1-2 expression and stable dismutation activity in NBT than ZJ. Hence, higher H₂O₂ content in the root apex of NBT than ZJ correlated with a significant increase in auxin content and ClSOD1-2 up-regulation. Moreover, Al³⁺-activated ClUGT75D1 and ClUGT75B2 in ZJ coincided with no considerable change in IBA content, suggesting that glycosylation-mediated changes in IBA content might be relevant to Al³⁺-tolerance in watermelon. Furthermore, exogenous H₂O₂ and IBA indicated ClUGT75D1 modulating IBA is likely dependent on H₂O₂ background. We hypothesize that a higher H₂O₂ level in NBT represses ClUGT75, resulting in increased auxin than those in ZJ roots. Thus, excess in both H₂O₂ and auxin aggravated the inhibition of root elongation under Al³⁺-stress. Our findings provide insights on the permissive action of H₂O₂ in the mediation of auxin glycosylation by ClUGT75 in root apex for Al³⁺-tolerance in watermelon.

Positive Interaction between H2O2 and Ca2+ Mediates Melatonin-Induced CBF Pathway and Cold Tolerance in Watermelon (Citrullus lanatus L.)

Cold stress is a major environmental factor that detrimentally affects plant growth and development. Melatonin has been shown to confer plant tolerance to cold stress through activating the C-REPEAT BINDING FACTOR (CBF) pathway; however, the underlying modes that enable this function remain obscure. In this study, we investigated the role of H₂O₂ and Ca²⁺ signaling in the melatonin-induced CBF pathway and cold tolerance in watermelon (Citrullus lanatus L.) through pharmacological, physiological, and genetic approaches. According to the results, melatonin induced H₂O₂ accumulation, which was associated with the upregulation of respiratory burst oxidase homolog D (ClRBOHD) during the early response to cold stress in watermelon. Besides, melatonin and H₂O₂ induced the accumulation of cytoplasmic free Ca²⁺ ([Ca²⁺]_cyt) in response to cold. This was associated with the upregulation of cyclic nucleotide-gated ion channel 2 (ClCNGC2) in watermelon. However, blocking of Ca²⁺ influx channels abolished melatonin- or H₂O₂-induced CBF pathway and cold tolerance. Ca²⁺ also induced ClRBOHD expression and H₂O₂ accumulation in early response to cold stress in watermelon. Inhibition of H₂O₂ production in watermelon by RBOH inhibitor or in Arabidopsis by AtRBOHD knockout compromised melatonin-induced [Ca²⁺]_cyt accumulation and melatonin- or Ca²⁺-induced CBF pathway and cold tolerance. Overall, these findings indicate that melatonin induces RBOHD-dependent H₂O₂ generation in early response to cold stress. Increased H₂O₂ promotes [Ca²⁺]_cyt accumulation, which in turn induces H₂O₂ accumulation via RBOHD, forming a reciprocal positive-regulatory loop that mediates melatonin-induced CBF pathway and subsequent cold tolerance.

Evolutionary gain of oligosaccharide hydrolysis and sugar transport enhanced carbohydrate partitioning in sweet watermelon fruits

How raffinose (Raf) family oligosaccharides, the major translocated sugars in the vascular bundle in cucurbits, are hydrolyzed and subsequently partitioned has not been fully elucidated. By performing reciprocal grafting of watermelon (Citrullus lanatus) fruits to branch stems, we observed that Raf was hydrolyzed in the fruit of cultivar watermelons but was backlogged in the fruit of wild ancestor species. Through a genome-wide association study, the alkaline alpha-galactosidase ClAGA2 was identified as the key factor controlling stachyose and Raf hydrolysis, and it was determined to be specifically expressed in the vascular bundle. Analysis of transgenic plants confirmed that ClAGA2 controls fruit Raf hydrolysis and reduces sugar content in fruits. Two single-nucleotide polymorphisms (SNPs) within the ClAGA2 promoter affect the recruitment of the transcription factor ClNF-YC2 (nuclear transcription factor Y subunit C) to regulate ClAGA2 expression. Moreover, this study demonstrates that C. lanatus Sugars Will Eventually Be Exported Transporter 3 (ClSWEET3) and Tonoplast Sugar Transporter (ClTST2) participate in plasma membrane sugar transport and sugar storage in fruit cell vacuoles, respectively. Knocking out ClAGA2, ClSWEET3, and ClTST2 affected fruit sugar accumulation. Genomic signatures indicate that the selection of ClAGA2, ClSWEET3, and ClTST2 for carbohydrate partitioning led to the derivation of modern sweet watermelon from non-sweet ancestors during domestication.

Identification of Key Gene Networks Associated With Cell Wall Components Leading to Flesh Firmness in Watermelon

Flesh firmness of watermelon is an important quality trait for commercial fruit values, including fruit storability, transportability, and shelf life. To date, knowledge of the gene networks underlying this trait is still limited. Herein, we used weighted genes co-expression network analysis (WGCNA) based on correlation and the association of phenotypic data (cell wall contents) with significantly differentially expressed genes between two materials, a near isogeneic line "HWF" (with high average flesh firmness) and inbred line "203Z" (with low average flesh firmness), to identify the gene networks responsible for changes in fruit flesh firmness. We identified three gene modules harboring 354 genes; these gene modules demonstrated significant correlation with water-soluble pectin, cellulose, hemicellulose, and protopectin. Based on intramodular significance, eight genes involved in cell wall biosynthesis and ethylene pathway are identified as hub genes within these modules. Among these genes, two genes, Cla012351 (Cellulose synthase) and Cla004251 (Pectinesterase), were significantly correlated with cellulose (r 2 = 0.83) and protopectin (r 2 = 0.81); three genes, Cla004120 (ERF1), Cla009966 (Cellulose synthase), and Cla006648 (Galactosyltransferase), had a significant correlation with water-soluble pectin (r 2 = 0.91), cellulose (r 2 = 0.9), and protopectin (r 2 = 0.92); and three genes, Cla007092 (ERF2a), Cla004119 (probable glycosyltransferase), and Cla018816 (Xyloglucan endotransglucosylase/hydrolase), were correlated with hemicellulose (r 2 = 0.85), cellulose (r 2 = 0.8), and protopectin (r 2 = 0.8). This study generated important insights of biosynthesis of a cell wall structure and ethylene signaling transduction pathway, the mechanism controlling the flesh firmness changes in watermelon, which provide a significant source to accelerate future functional analysis in watermelon to facilitate crop improvement.

Identification and selection of reference genes for gene expression analysis by quantitative real-time PCR in Suaeda glauca's response to salinity

Quantitative real-time polymerase chain reaction (qPCR) using a stable reference gene is widely used for gene expression research. Suaeda glauca L. is a succulent halophyte and medicinal plant that is extensively used for phytoremediation and extraction of medicinal compounds. It thrives under high-salt conditions, which promote the accumulation of high-value secondary metabolites. However, a suitable reference gene has not been identified for gene expression standardization in S. glauca under saline conditions. Here, 10 candidate reference genes, ACT7, ACT11, CCD1, TUA5, UPL1, PP2A, DREB1D, V-H⁺-ATPase, MPK6, and PHT4;5, were selected from S. glauca transcriptome data. Five statistical algorithms (ΔCq, geNorm, NormFinder, BestKeeper, and RefFinder) were applied to determine the expression stabilities of these genes in 72 samples at different salt concentrations in different tissues. PP2A and TUA5 were the most stable reference genes in different tissues and salt treatments, whereas DREB1D was the least stable. The two reference genes were sufficient to normalize gene expression across all sample sets. The suitability of identified reference genes was validated with MYB and AP2 in germinating seeds of S. glauca exposed to different NaCl concentrations. Our study provides a foundational framework for standardizing qPCR analyses, enabling accurate gene expression profiling in S. glauca.

Methyl jasmonate mediates melatonin-induced cold tolerance of grafted watermelon plants

Root-shoot communication has a critical role in plant adaptation to environmental stress. Grafting is widely applied to enhance the abiotic stress tolerance of many horticultural crop species; however, the signal transduction mechanism involved in this tolerance remains unknown. Here, we show that pumpkin- or figleaf gourd rootstock-enhanced cold tolerance of watermelon shoots is accompanied by increases in the accumulation of melatonin, methyl jasmonate (MeJA), and hydrogen peroxide (H₂O₂). Increased melatonin levels in leaves were associated with both increased melatonin in rootstocks and MeJA-induced melatonin biosynthesis in leaves of plants under cold stress. Exogenous melatonin increased the accumulation of MeJA and H₂O₂ and enhanced cold tolerance, while inhibition of melatonin accumulation attenuated rootstock-induced MeJA and H₂O₂ accumulation and cold tolerance. MeJA application induced H₂O₂ accumulation and cold tolerance, but inhibition of JA biosynthesis abolished rootstock- or melatonin-induced H₂O₂ accumulation and cold tolerance. Additionally, inhibition of H₂O₂ production attenuated MeJA-induced tolerance to cold stress. Taken together, our results suggest that melatonin is involved in grafting-induced cold tolerance by inducing the accumulation of MeJA and H₂O₂. MeJA subsequently increases melatonin accumulation, forming a self-amplifying feedback loop that leads to increased H₂O₂ accumulation and cold tolerance. This study reveals a novel regulatory mechanism of rootstock-induced cold tolerance.

Identification and characterization of a natural SNP variant in ALTERNATIVE OXIDASE gene associated with cold stress tolerance in watermelon

Alternative oxidase (AOX) is a mitochondrial enzyme encoded by a small nuclear gene family, which contains the two subfamilies, AOX1 and AOX2. In the present study on watermelon (Citrullus lanatus), only one ClAOX gene, belonging to AOX2 subfamily but having a similar gene structure to AtAOX1a, was found in the watermelon genome. The expression analysis suggested that ClAOX had the constitutive expression feature of AOX2 subfamily, but was cold inducible, which is normally considered an AOX1 subfamily feature. Moreover, one single nucleotide polymorphism (SNP) in ClAOX sequence, which led to the change from Lys (N) to Asn (K) in the 96^th amino acids, was found among watermelon subspecies. Ectopic expression of two ClAOX alleles in the Arabidopsis aox1a knock-out mutant indicated that ClAOX^K-expressing plants had stronger cold tolerance than aox1a mutant and ClAOX^N-expressing plants. Our findings suggested watermelon genome contained a single ClAOX that possessed the expression features of both AOX1 and AOX2 subfamilies. A naturally existing SNP in ClAOX differentiated the cold tolerance of transgenic Arabidopsis plants, impling a possibility this gene might be a functional marker for stress-tolerance breeding.

The effects of pomegranate peel extract on the gene expressions of antioxidant enzymes in a rat model of alloxan-induced diabetes

This study was conducted to evaluate the anti-diabetic and antioxidant effects of hydroalcoholic pomegranate peel extract (APE) in alloxan-induced diabetes rat models. We divided 60 rats into the following six equal groups (n = 10): Healthy control; diabetic control (100 mg/kg alloxan); sham + glibenclamide (10 mg/kg); diabetic + glibenclamide (10 mg/kg); sham + APE (200 mg/kg) and diabetic + APE (200 mg/kg). After 8 weeks, kidneys were taken out for biochemical and molecular studies. Following APE treatment, biochemical parameters including malondialdehyde (MDA), and glutathione (GSH), glutathione peroxidase (GPx), catalase (CAT), superoxide dismutase (SOD) significantly induced in the treated group as compared with the control group (p < 0.05). Also, gene expression of GPx (3-fold), CAT (2.6-fold), and SOD (1.5-fold) were increased as compared to controls (p < 0.05). Overall, our results indicated that pomegranate can be used as an antioxidant agent to reduce complications from diseases associated with oxidative stress.

Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp

Quantitative real-time PCR (qRT-PCR) is commonly used to measure gene expression to further explore gene function, while suitable reference genes must be stably expressed under different experimental conditions to obtain accurate and reproducible data for relative quantification. Taxol or paclitaxel is an important anticancer compound mainly identified in Taxus spp. The molecular mechanism of the regulation of taxol biosynthesis is current research goal. However, in the case of Taxus spp., few reports were published on screening suitable reference genes as internal controls for qRT-PCR. Here, eight reference genes were selected as candidate reference genes for further study. Common statistical algorithms geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to analyze the data from samples collected from a cell line of Taxus × media under various experimental conditions and from tissues of Taxus chinensis var. mairei. The expression patterns of TcMYC under salicylic acid treatment differed significantly, with the best and worst reference genes in the cell line. This study screened out suitable reference genes (GAPDH1 and SAND) under different treatments and tissues for the accurate and reliable normalization of the qRT-PCR expression data of Taxus spp. At the same time, this study will aid future research on taxol biosynthesis-related genes expression in Taxus spp., and can also be directly used to other related species.

Identification of key gene networks controlling organic acid and sugar metabolism during watermelon fruit development by integrating metabolic phenotypes and gene expression profiles

The organoleptic qualities of watermelon fruit are defined by the sugar and organic acid contents, which undergo considerable variations during development and maturation. The molecular mechanisms underlying these variations remain unclear. In this study, we used transcriptome profiles to investigate the coexpression patterns of gene networks associated with sugar and organic acid metabolism. We identified 3 gene networks/modules containing 2443 genes highly correlated with sugars and organic acids. Within these modules, based on intramodular significance and Reverse Transcription Quantitative polymerase chain reaction (RT-qPCR), we identified 7 genes involved in the metabolism of sugars and organic acids. Among these genes, Cla97C01G000640, Cla97C05G087120 and Cla97C01G018840 (r² = 0.83 with glucose content) were identified as sugar transporters (SWEET, EDR6 and STP) and Cla97C03G064990 (r² = 0.92 with sucrose content) was identified as a sucrose synthase from information available for other crops. Similarly, Cla97C07G128420, Cla97C03G068240 and Cla97C01G008870, having strong correlations with malic (r² = 0.75) and citric acid (r² = 0.85), were annotated as malate and citrate transporters (ALMT7, CS, and ICDH). The expression profiles of these 7 genes in diverse watermelon genotypes revealed consistent patterns of expression variation in various types of watermelon. These findings add significantly to our existing knowledge of sugar and organic acid metabolism in watermelon.

Acidovorax citrulli Type III Effector AopP Suppresses Plant Immunity by Targeting the Watermelon Transcription Factor WRKY6

Acidovorax citrulli (Ac) is the causal agent of bacterial fruit blotch (BFB), and BFB poses a threat to global watermelon production. Despite its economic importance, the molecular mechanisms underlying Ac pathogenicity and virulence are not well understood, particularly with regard to its type III secreted effectors. We identify a new effector, AopP, in Ac and confirm its secretion and translocation. AopP suppresses reactive oxygen species burst and salicylic acid (SA) content and significantly contributes to virulence. Interestingly, AopP interacts with a watermelon transcription factor, ClWRKY6, in vivo and in vitro. ClWRKY6 shows typical nuclear localization, and AopP and ClWRKY6 co-localize in the nucleus. Ac infection, SA, and the pathogen-associated molecular pattern flg22 ^Ac promote ClWRKY6 production, suggesting that ClWRKY6 is involved in plant immunity and SA signaling. Furthermore, ClWRKY6 positively regulates PTI and SA production when expressed in Nicotiana benthamiana. Importantly, AopP reduces ClWRKY6 mRNA and ClWRKY6 protein levels, suggesting that AopP suppresses plant immunity by targeting ClWRKY6. In summary, we identify a novel effector associated with the virulence mechanism of Ac, which interacts with the transcription factor of the natural host, watermelon. The findings of this study provide insights into the mechanisms of watermelon immune responses and may facilitate molecular breeding for bacterial fruit blotch resistance.

Subcellular distribution of aluminum associated with differential cell ultra-structure, mineral uptake, and antioxidant enzymes in root of two different Al+3-resistance watermelon cultivars

Crop plants, such as watermelon, suffer from severe Aluminum (Al³⁺)-toxicity in acidic soils with their primary root elongation being first arrested. However, the significance of apoplastic or symplastic Al³⁺-toxicity in watermelon root is scarcely reported. In this work, we identified a medium fruit type (ZJ) and a small fruit type (NBT) as Al⁺³-tolerant and sensitive based on their differential primary root elongation rate respectively, and used them to show the effects of symplastic besides apoplastic Al distribution in the watermelon's root. Although the Al content was higher in the root of NBT than ZJ, Al⁺³ allocated in their apoplast, vacuole and plastid fractions were not significantly different between the two cultivars. Thus, only a few proportion of Al⁺³ differentially distributed in the nucleus and mitochondria corresponded to interesting differential morphological and physiological disorders recorded in the root under Al⁺³-stress. The symplastic amount of Al⁺³ substantially induced the energy efficient catalase pathway in ZJ, and the energy consuming ascorbate peroxidase pathway in NBT. These findings coincided with obvious starch granule visibility in the root ultra-structure of ZJ than NBT, suggesting a differential energy was used in supporting the root elongation and nutrient uptake for Al⁺³-tolerance in the two cultivars. This work provides clues that could be further investigated in the identification of genetic components and molecular mechanisms associated with Al⁺³-tolerance in watermelon.

Screening and verification of reference genes for analysis of gene expression in winter rapeseed (Brassica rapa L.) under abiotic stress

Winter rapeseed (Brassica rapa L.) is the main oilseed crop in northern China and can safely overwinter at 35 (i.e., Tianshui, China) to 48 degrees north latitude (i.e., Altai, Heilongjiang, Raohe, and Xinjiang, China). In order to identify stable reference genes to understand the molecular mechanisms of stress tolerance in winter rapeseed, internal reference genes of winter rapeseed under four abiotic stresses were analyzed using GeNorm, NormFinder, BestKeeper, and RefFinder software. The most stable combinations of internal reference genes were β-actin and SAND in cold-stressed leaves, β-actin and EF1a in cold-stressed roots, F-box and SAND in high temperature-stressed leaves, and PP2A and RPL in high temperature-stressed roots, SAND and PP2A in NaCl-stressed leaves, RPL and UBC in NaCl-stressed roots, RPL and PP2A in PEG-stressed leaves, and PP2A and RPL in PEG-stressed roots. Expression profiles of PXG3 were used to verify these results. The stable reference genes identified in this study are useful tools for identifying stress-responsive genes to understand the molecular mechanisms of stress tolerance in winter rapeseed.

Identification and Functional Analysis of AopN, an Acidovorax Citrulli Effector that Induces Programmed Cell Death in Plants

Bacterial fruit blotch (BFB), caused by Acidovorax citrulli, seriously affects watermelon and other cucurbit crops, resulting in significant economic losses. However, the pathogenicity mechanism of A. citrulli is not well understood. Plant pathogenic bacteria often suppress the plant immune response by secreting effector proteins. Thus, identifying A. citrulli effector proteins and determining their functions may improve our understanding of the underlying pathogenetic mechanisms. In this study, a novel effector, AopN, which is localized on the cell membrane of Nicotiana benthamiana, was identified. The functional analysis revealed that AopN significantly inhibited the flg22-induced reactive oxygen species burst. AopN induced a programmed cell death (PCD) response. Unlike its homologous protein, the ability of AopN to induce PCD was dependent on two motifs of unknown functions (including DUP4129 and Cpta_toxin), but was not dependent on LXXLL domain. More importantly, the virulence of the aopN mutant of A. citrulli in N. benthamiana significantly decreased, indicating that it was a core effector. Further analysis revealed that AopN interacted with watermelon ClHIPP and ClLTP, which responds to A. citrulli strain Aac5 infection at the transcription level. Collectively, these findings indicate that AopN suppresses plant immunity and activates the effector-triggered immunity pathway.

Transcriptomic Analysis of Short-Term Salt Stress Response in Watermelon Seedlings

Watermelon (Citrullus lanatus L.) is a widely popular vegetable fruit crop for human consumption. Soil salinity is among the most critical problems for agricultural production, food security, and sustainability. The transcriptomic and the primary molecular mechanisms that underlie the salt-induced responses in watermelon plants remain uncertain. In this study, the photosynthetic efficiency of photosystem II, free amino acids, and transcriptome profiles of watermelon seedlings exposed to short-term salt stress (300 mM NaCl) were analyzed to identify the genes and pathways associated with response to salt stress. We observed that the maximal photochemical efficiency of photosystem II decreased in salt-stressed plants. Most free amino acids in the leaves of salt-stressed plants increased many folds, while the percent distribution of glutamate and glutamine relative to the amino acid pool decreased. Transcriptome analysis revealed 7622 differentially expressed genes (DEGs) under salt stress, of which 4055 were up-regulated. The GO analysis showed that the molecular function term “transcription factor (TF) activity” was enriched. The assembled transcriptome demonstrated up-regulation of 240 and down-regulation of 194 differentially expressed TFs, of which the members of ERF, WRKY, NAC bHLH, and MYB-related families were over-represented. The functional significance of DEGs associated with endocytosis, amino acid metabolism, nitrogen metabolism, photosynthesis, and hormonal pathways in response to salt stress are discussed. The findings from this study provide novel insights into the salt tolerance mechanism in watermelon.

Analysis of differentially expressed genes and pathways associated with male sterility lines in watermelon via bulked segregant RNA-seq

Genic male sterility (GMS) is a common and important trait, which is widely used for the production of hybrid seeds. However, the molecular mechanism of GMS in watermelon remains poorly understood. In this study, we comparatively analyzed the transcriptome profiles of sterile and fertile floral buds using the bulked segregant analysis (BSA) and transcriptome sequencing (RNA-seq). A total of 2507 differentially expressed genes (DEGs) including 593 up-regulated and 1914 down-regulated, were identified to be related to male sterility in watermelon line Se18. Gene ontology (GO) analysis showed that 57 GO terms were significantly enriched, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed plant hormone signal transduction, glycolysis/gluconeogenesis, starch and sucrose metabolism, plant-pathogen interaction, phenylpropanoid biosynthesis pathways were obviously enriched. Furthermore, the efficiency of the RNA-seq analysis was validated by quantitative real-time PCR (qRT-PCR). Among the DEGs, some valuable candidate genes involved in pollen development were identified, such as gene Cla000029, a bHLH transcription factor and homologous to MS1 in Arabidopsis. Moreover, other DEGs including MYB gene Cla012590 (MYB26), Cla017100 (MYB21), etc., also provide useful information for further understanding the function of key genes involved in pollen development. This study provides new insights into the global network of male sterility in watermelon.

MicroRNA-202-3p Targets Brain-Derived Neurotrophic Factor and Is Involved in Depression-Like Behaviors

BACKGROUND

Brain-derived neurotrophic factor (BDNF) and microRNA (miRNA) play crucial roles in the etiology of depression. However, the molecular mechanisms underlying this disease are not fully understood. The primary objective of this study was to investigate the relationship between miR-202-3p and BDNF in a chronic unpredictable mild stress (CUMS) model.

METHODS

Depression model was established with chronic mild unpredictable mild stimulation (CUMS) combined with solitary feeding. The expression levels of miR-202-3p and BDNF in rat hippocampus were measured by qRT-PCR. The novelty inhibition feeding test (NSFT), sucrose preference test (SPT), and forced swimming test (FST) were used to evaluate the functions of miR-202-3p and BDNF. Target gene prediction and screening and luciferase reporter assay were used to verify the target of miR-202-3p. The expression levels of BNDF, CREB1 and p-CREB1 were detected by Western blot.

RESULTS

Upregulation of miR-202-3p was associated with decreased expression of BDNF in the hippocampus of the CUMS model. Antidepressant was observed when LV-BDNF or LV-si-miR-202-3p was injected into the hippocampus. In addition, in the rat hippocampus and cultured nerve cells, the expression levels of BDNF and cyclic AMP response element binding protein 1 (CREB1), which is a target gene of BDNF, were reduced after LV-miR-202-3p injection. Overexpression of miR-202-3p aggravated depressive behavior and decreased the expression levels of BDNF. Luciferase reporter assay also confirmed that BDNF was a target of miR-202-3p.

CONCLUSION

Silencing miR-202-3p can reduce the damage to hippocampal nerve in CUMS rats; the mechanism may be related to the upregulation of BNDF expression. miR-202-3p may be an effective target for the treatment of depression.

Genome and Transcriptome Analyses Provide Insight Into the Omega-3 Long-Chain Polyunsaturated Fatty Acids Biosynthesis of Schizochytrium limacinum SR21

Schizochytrium sp. is the best natural resource for omega-3 long-chain polyunsaturated fatty acids. We report a high-quality genome sequence of Schizochytrium limacinum SR21, which has a 63 Mb genome size, with a contig N50 of 2.67 Mb and 6,838 protein-coding genes. Phylogenomic and comparative genomic analyses revealed that DHA-producing Schizochytrium and Aurantiochytrium strains were highly similar and possessed similar genes. Analysis of the fatty acid synthase (FAS) for LC-PUFAs production results in the annotation of all genes in map00062 and map01212. A gene cluster and 10 ORFs related to PKS pathway were found in the genome. 1,402 differentially expressed genes (DEGs) of the treated groups (0.5 g/L yeast extract) were identified by comparing with the control groups (1.0 g/L yeast extract) at 36 h. A weighted gene coexpression network analysis revealed that 2 of 7 modules correlated highly with the fatty acid and DHA contents. The DEGs and transcription factors were significantly correlated with fatty acid biosynthesis, including MYB, Zinc Finger and ACOX. The results showed that these hub genes are regulated by genes involved in fatty acid biosynthesis pathways. The results providing an important reference for further research on promoting fatty acid and DHA accumulation in S. limacinum SR21.

Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plants under differential biotic stresses

The selection of reliable reference genes (RGs) for normalization under given experimental conditions is necessary to develop an accurate qRT-PCR assay. To the best of our knowledge, only a small number of RGs have been rigorously identified and used in tea plants (Camellia sinensis (L.) O. Kuntze) under abiotic stresses, but no critical RG identification has been performed for tea plants under any biotic stresses till now. In the present study, we measured the mRNA transcriptional levels of ten candidate RGs under five experimental conditions; these genes have been identified as stable RGs in tea plants. By using the ΔCt method, geNorm, NormFinder and BestKeeper, CLATHRIN1 and UBC1, TUA1 and SAND1, or SAND1 and UBC1 were identified as the best combination for normalizing diurnal gene expression in leaves, stems and roots individually; CLATHRIN1 and GAPDH1 were identified as the best combination for jasmonic acid treatment; ACTIN1 and UBC1 were identified as the best combination for Toxoptera aurantii-infested leaves; UBC1 and GAPDH1 were identified as the best combination for Empoasca onukii-infested leaves; and SAND1 and TBP1 were identified as the best combination for Ectropis obliqua regurgitant-treated leaves. Furthermore, our results suggest that if the processing time of the treatment was long, the best RGs for normalization should be recommended according to the stability of the proposed RGs in different time intervals when intragroup differences were compared, which would strongly increase the accuracy and sensitivity of target gene expression in tea plants under biotic stresses. However, when the differences of intergroup were compared, the RGs for normalization should keep consistent across different time points. The results of this study provide a technical guidance for further study of the molecular mechanisms of tea plants under different biotic stresses.

MicroRNA-579-3p Exerts Neuroprotective Effects Against Ischemic Stroke via Anti-Inflammation and Anti-Apoptosis

BACKGROUND/AIMS

Multiple studies have found that microRNAs (miRNAs) are involved in the development of cerebral ischemia. MiR-579-3p can inhibit inflammatory responses and apoptosis, leading to ischemia/reperfusion (I/R) damage. However, the mechanism of how miR-579-3p actions in brain I/R injury remains unclear. This study aimed to investigate the mechanism of the role of miR-579-3p in brain I/R injury.

METHODS

A rat model of cerebral ischemia-reperfusion injury was established by suture method. The effects of miR-579-3p on cerebral infarction size, brain water content, and neurological symptoms were evaluated. Flow cytometry was used to detect apoptosis. ELISA was used to detect the level of inflammatory factors. Western blot was used to detect the expression of P65, NCOA1, Bcl-2 and Bax. The relationship between miR-579-3p and NCOA1 was analyzed by bioinformatics analysis and luciferase assay.

RESULTS

Overexpression of miR-579-3p reduced infarct volume, brain water content and neurological deficits. Overexpression of miR-579-3p inhibited the expression level of the inflammatory cytokines, such as TNF-α, IL-6, COX-2 and iNOS, and increased the expression level of IL-10. MiR-579-3p overexpression inhibited NF-кB activity by reducing NRIP1. In addition, miR-579-3p could reduce the apoptotic rate of cortical neurons. Overexpression of miR-579-3p inhibited the activity of caspase-3, increased the expression level of anti-apoptotic gene Bcl-2 in neurons, and decreased the expression level of apoptotic gene Bax.

CONCLUSION

miR-579-3p can be used to treat brain I/R injury, and its neuroprotective effect may be ascribed to the reduction of inflammation and apoptosis.

Functional Relevance of Citrulline in the Vegetative Tissues of Watermelon During Abiotic Stresses

A non-protein amino acid, citrulline, is a compatible solute involved in the maintenance of cellular osmolarity during abiotic stresses. Despite its significance, a coherent model indicating the role of citrulline during stress conditions has not yet emerged. We have used watermelon, naturally rich in citrulline, as a model to understand its accumulation during drought stress and nitrogen perturbation using transcriptomic and metabolomic analysis. Experiments were performed in the semi-controlled environment, and open field to study the accumulation of drought-induced citrulline in the vegetative tissues of watermelon by monitoring the stress treatments using physiological measurements. The amino acid profiling of leaves and stems in response to drought stress showed up to a 38 and 16-fold increase in citrulline content, respectively. Correlation between amino acids indicated a concomitant activation of a metabolic pathway that included citrulline, its precursor (ornithine), and catabolic product (arginine). Consistent with its accumulation, the gene expression analysis and RNA-Sequencing confirmed activation of citrulline biosynthesis-related genes - Ornithine carbamoyl-transferase (OTC), N-acetylornithine deacetylase (AOD) and Carbamoyl phosphate synthases (CPS), and down-regulation of catabolic genes; Arginosuccinate lyase (ASL) and Arginosuccinate synthases (ASS) in drought-stressed leaf tissues. Based on the relative abundance in the nitrogen-depleted vegetative tissues and down-regulation of genes involved in citrulline biosynthesis, we also demonstrated that the nitrogen status of the plant regulates citrulline. Taken together, these data provide further insights into the metabolic and molecular mechanisms underlying the amino acid metabolism under environmental stress and the significance of non-protein amino acid citrulline in plants.

Identification of watermelon heat shock protein members and tissue-specific gene expression analysis under combined drought and heat stresses

Heat shock protein (Hsp) gene family members in the watermelon genome were identified and characterized by bioinformatics analysis. In addition, expression profiles of genes under combined drought and heat stress conditions were experimentally analyzed. In the watermelon genome, 39 genes belonging to the sHsp family, 101 genes belonging to the Hsp40 family, 23 genes belonging to the Hsp60 family, 12 genes belonging to the Hsp70 family, 6 genes belonging to the Hsp90 family, and 102 genes belonging to the Hsp100 family were found. It was also observed that the proteins in the same cluster in the phylogenetic trees had similar motif patterns. When the estimated 3-dimensional structures of the Hsp proteins were examined, it was determined that the α-helical structure was dominant in almost all families. The most orthologous relationship appeared to be between watermelon, soybean, and poplar in the ClaHsp gene families. For tissue-specific gene expression analysis under combined stress conditions, expression analysis of one representative Hsp gene each from root, stem, leaf, and shoot tissues was performed by real-time PCR. A significant increase was detected usually at 30 min in almost all tissues. This study provides extensive information for watermelon Hsps, and can enhance our knowledge about the relationships between Hsp genes and combined stresses.

Wheat straw increases the defense response and resistance of watermelon monoculture to Fusarium wilt

BACKGROUND

Wheat straw is a rich resource worldwide. Straw return is an effective strategy to alleviate soil-borne diseases on monoculture watermelon. Previous studies focus on soil structure, physical and chemical properties; however, little is known about the molecular responses on host plant.

RESULTS

No significant difference on the population of Fusarium oxysporum f.sp. niveum race 1(Fon1) in rhizosphere soil was found between control (no addition of wheat straw) and the treated groups (addition of 1% (T1) or 2% (T2) wheat straw). RNA-Seq analysis showed that 3419 differentially expressed genes were clustered into 8 profiles. KEGG analysis revealed that phenylpropanoid biosynthesis and plant hormone signal transduction were involved in wheat straw induced response in monoculture watermelon. Genes in lignin biosynthesis were found to be upregulated, and the lignin and auxin contents were higher in T1 and T2 compared to the control. Lignin was also enriched and the Fon1 population decreased in watermelon roots treated with wheat straw. The enzyme activities of phenylalanine ammonia-lyase and peroxidase were increased.

CONCLUSIONS

Our data suggest that the addition of wheat straw enhances the defense response to Fon1 infection in watermelon through increasing lignin and auxin biosynthesis.

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

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

Expression and Localization Patterns of a Small Heat Shock Protein that Interacts with the Helicase Domain of Cucumber Green Mottle Mosaic Virus

Cucumber green mottle mosaic virus (CGMMV), a member of the genus Tobamovirus (family Virgaviridae), is an economically important virus that has detrimental effects on cucurbit crops worldwide. Understanding the interaction between host factors and CGMMV viral proteins will facilitate the design of new strategies for disease control. In this study, a yeast two-hybrid assay revealed that the CGMMV helicase (HEL) domain interacts with a Citrullus lanatus small heat shock protein (sHSP), and we verified this observation by performing in vitro GST pull-down and in vivo coimmunoprecipitation assays. Measurement of the levels of accumulated sHSP transcript revealed that sHSP is upregulated on initial CGMMV infection in both Nicotiana benthamiana and C. lanatus plants, although not in the systemically infected leaves. We also found that the subcellular localization of the sHSP was altered after CGMMV infection. To further validate the role of sHSP in CGMMV infection, we produced and assayed N. benthamiana transgenic plants with up- and down-regulated sHSP expression. Overexpression of sHSP inhibited viral RNA accumulation and retarded disease development, whereas sHSP silencing had no marked effect on CGMMV infection. Therefore, we postulate that the identified sHSP may be one of the factors modulating host defense mechanisms in response to CGMMV infection and that the HEL domain interaction may inhibit this sHSP function to promote viral infection.

Identification of internal control genes for circular RNAs

OBJECTIVE

At present, no studies have established internal control genes for circular RNA (circRNA) analyses. We aimed to identify reference circRNAs for real-time quantitative PCR (RT-qPCR).

RESULTS

After analyzing the RNA-seq data, we obtained 50 circRNAs that were expressed in all samples. We ranked these 50 circRNAs according to their stability and obtained the six most stable circRNAs. We further evaluated the stability of the six circRNAs and three linear control genes (i.e., GAPDH, β-actin and 18S rRNA) in 22 cell lines. Our results indicated that hsa_circ_0000284 (circHIPK3) and hsa_circ_0000471 (circN4BP2L2) were the two most stable genes. After removing linear RNAs or including the cells treated with Adriamycin, NH₄Cl and shikonin, the two most stable genes were hsa_circ_0000471 and hsa_circ_0000284. The amplification efficiency was 100% for hsa_circ_0000471 and 95% for hsa_circ_0000284.

CONCLUSIONS

In conclusion, since the stability of circRNAs is higher than that of linear RNAs, hsa_circ_0000284 and hsa_circ_0000471 may be used as reference genes not only for circRNAs but also for other kinds of RNAs. The findings in the present study fill the gap of lacking reference genes in the detection of circRNAs.

Mesoporous Silica Containers and Programmed Catalytic Hairpin Assembly/Hybridization Chain Reaction Based Electrochemical Sensing Platform for MicroRNA Ultrasensitive Detection with Low Background

In this work, based on mesoporous silica containers (MSNs) with the programmed enzyme-free DNA assembly amplification of catalytic hairpin assembly (CHA) and hybridization chain reaction (HCR), an ultrasensitive electrochemical sensing platform with low background is developed for the detection of microRNA (miRNA). Herein, the electrochemical reporter methylene blue (MB) was sealed in the pores of MSNs by the double-stranded DNA (dsDNA) gate of hairpin DNA H1 and anchor DNA. In the absence of target, neither the CHA nor the HCR process happened, which enabled a low background. After target was added, DNA H1 was displaced from the MSNs surface and participated in the CHA process with the assistance of hairpin DNA H2, which accelerated the release of MB from the MSNs pore. Meanwhile, the CHA products H1-H2 were hybridized with the capture probes (SH-CP) on the electrode surface, which further initiated the HCR process. The released MB from the MSNs will effectively intercalate into long dsDNA polymers of HCR products, resulting in a significant electrochemical response. Taking miRNA-21 as the model target, the proposed sensing platform achieves a satisfactory detection limit down to 0.037 fM, which is lower than that of electrochemical assay with amplification methods. In addition, the strategy shows good selectivity against other miRNAs and is capable in practical analytes. Benefitting from the features of being label-free and enzyme-free and having low background, high sensitivity, and selectivity, this strategy shows great potential in bioanalysis and clinical diagnostics.

Selection of reference genes for flowering pathway analysis in the masting plants, Celmisia lyallii and Chionochloa pallens, under variable environmental conditions

Mast flowering is characterised by mass synchronised flowering at irregular intervals over a wide geographical area. An understanding of the molecular drivers of mast flowering requires expression analysis of key developmentally regulated gene(s). Reverse transcription-quantitative PCR is the gold standard technique used to assess expression of target gene(s) and to validate high-throughput sequencing data. Selection and validation of appropriate reference gene(s), used as normalisation factors in transcript abundance analysis, is an essential step to avoid ambiguous expression results. Eight candidate reference genes were assessed to select the best internal normalisation factors in naturally growing masting plants Chionochloa pallens and Celmisia lyallii. Statistical packages geNorm, Normfinder, BestKeeper, ΔCt and RefFinder were used to determine the expression stability in plants translocated to different altitudes and sampled across the season. GAPDH and PP2a in Celmisia and ExP and THP in Chionochloa were found to be the best pairs of reference genes for normalisation of the gene expression data. Our study revealed environmentally-induced changes in reference gene expression, information that will be utilised as we investigate flowering phenology of masting plants under global climatic change.

Comparative Analysis of Calcium-Dependent Protein Kinase in Cucurbitaceae and Expression Studies in Watermelon

Both the calcium-dependent protein kinases (CDPKs) and CDPK-related kinases (CRKs) play numerous roles in plant growth, development, and stress response. Despite genome-wide identification of both families in Cucumis, comparative evolutionary and functional analysis of both CDPKs and CRKs in Cucurbitaceae remain unclear. In this study, we identified 128 CDPK and 56 CRK genes in total in six Cucurbitaceae species (C. lanatus, C. sativus, C. moschata, C. maxima, C. pepo, and L. siceraria). Dot plot analysis indicated that self-duplication of conserved domains contributed to the structural variations of two CDPKs (CpCDPK19 and CpCDPK27) in C. pepo. Using watermelon genome as reference, an integrated map containing 25 loci (16 CDPK and nine CRK loci) was obtained, 16 of which (12 CDPK and four CRK) were shared by all seven Cucurbitaceae species. Combined with exon-intron organizations, topological analyses indicated an ancient origination of groups CDPK IV and CRK. Moreover, the evolutionary scenario of seven modern Cucurbitaceae species could also be reflected on the phylogenetic trees. Expression patterns of ClCDPKs and ClCRKs were studied under different abiotic stresses. Some valuable genes were uncovered for future gene function exploration. For instance, both ClCDPK6 and its ortholog CsCDPK14 in cucumber could be induced by salinity, while ClCDPK6 and ClCDPK16, as well as their orthologs in Cucumis, maintained high expression levels in male flowers. Collectively, these results provide insights into the evolutionary history of two gene families in Cucurbitaceae, and indicate a subset of candidate genes for functional characterizations in the future.

Comprehensive Analysis of TIFY Transcription Factors and Their Expression Profiles under Jasmonic Acid and Abiotic Stresses in Watermelon

The TIFY gene family is plant-specific and encodes proteins involved in the regulation of multiple biological processes. Here, we identified 15 TIFY genes in the watermelon genome, which were divided into four subfamilies (eight JAZs, four ZMLs, two TIFYs, and one PPD) in the phylogenetic tree. The ClTIFY genes were unevenly located on eight chromosomes, and three segmental duplication events and one tandem duplication event were identified, suggesting that gene duplication plays a vital role in the expansion of the TIFY gene family in watermelon. Further analysis of the protein architectures, conserved domains, and gene structures provided additional clues for understanding the putative functions of the TIFY family members. Analysis of qRT-PCR and RNA-seq data revealed that the detected ClTIFY genes had preferential expression in specific tissues. qRT-PCR analysis revealed that nine selected TIFY genes were responsive to jasmonic acid (JA) and abiotic stresses including salt and drought. JA activated eight genes and suppressed one gene, among which ClJAZ1 and ClJAZ7 were the most significantly induced. Salt and drought stress activated nearly all the detected genes to different degrees. These results lay a foundation for further functional characterization of TIFY family genes in Citrullus lanatus.

A Single Nucleotide Deletion in an ABC Transporter Gene Leads to a Dwarf Phenotype in Watermelon

Dwarf habit is one of the most important traits in crop plant architecture, as it can increase plant density and improved land utilization, especially for protected cultivation, as well as increasing lodging resistance and economic yield. At least four dwarf genes have been identified in watermelon, but none of them has been cloned. In the current study, the Cldw-1 gene was primary-mapped onto watermelon chromosome 9 by next-generation sequencing-aided bulked-segregant analysis (BSA-seq) of F₂ plants derived from a cross between a normal-height line, WT4, and a dwarf line, WM102, in watermelon. The candidate region identified by BSA-seq was subsequently validated and confirmed by linkage analysis using 30 simple sequence repeat (SSR) markers in an F₂ population of 124 plants. The Cldw-1 gene was further fine-mapped by chromosome walking in a large F₂ population of 1,053 plants and was delimited into a candidate region of 107.00 kb. Six genes were predicted to be in the candidate region, and only one gene, Cla010337, was identified to have two single nucleotide polymorphisms (SNPs) and a single nucleotide deletion in the exons in the dwarf line, WM102. A derived cleaved amplified polymorphic sequence (dCAPS) marker was developed from the single nucleotide deletion, co-segregated with the dwarf trait in both the F₂ population and a germplasm collection of 165 accessions. Cla010337 encoded an ATP-binding cassette transporter (ABC transporter) protein, and the expression levels of Cla010337 were significantly reduced in all the tissues tested in the dwarf line, WM102. The results of this study will be useful in achieving a better understanding of the molecular mechanism of the dwarf plant trait in watermelon and for the development of marker-assisted selection (MAS) for new dwarf cultivars.

A point mutation resulting in a 13 bp deletion in the coding sequence of Cldf leads to a GA-deficient dwarf phenotype in watermelon

The dwarf architecture is an important and valuable agronomic trait in watermelon breeding and has the potential to increase fruit yield and reduce labor cost in crop cultivation. However, the molecular basis for dwarfism in watermelon remains largely unknown. In this study, a recessive dwarf allele (designated as Cldf (Citrullus lanatus dwarfism)) was fine mapped in a 32.88 kb region on chromosome 09 using F₂ segregation populations derived from reciprocal crossing of a normal line M08 and a dwarf line N21. Gene annotation of the corresponding region revealed that the Cla015407 gene encoding a gibberellin 3β-hydroxylase functions as the best possible candidate gene for Cldf. Sequence analysis showed that the fourth polymorphism site (a G to A point mutation) at the 3' AG splice receptor site of the intron leads to a 13 bp deletion in the coding sequence of Cldf in dwarf line N21 and thus results in a truncated protein lacking the conserved domain for binding 2-oxoglutarate. In addition, the dwarf phenotype of Cldf could be rescued by exogenous GA₃ application. Phylogenetic analysis suggested that the small multigene family GA3ox (GA3 oxidase) in cucurbit species may originate from three ancient lineages in Cucurbitaceae. All these data support the conclusion that Cldf is a GA-deficient mutant, which together with the cosegregated marker can be used for breeding new dwarf cultivars.