Assessment of Suitable Reference Genes for Quantitative Gene Expression Studies in Melon Fruits

The ethylene-responsive transcription factor ERF024 is a novel regulator of climacteric fruit ripening in melon

Fruit ripening is an essential developmental stage in Angiosperms triggered by hormonal signals such as ethylene, a major player in climacteric ripening. Melon is a unique crop showing both climacteric and non-climacteric cultivars, offering an ideal model for dissecting the genetic mechanisms underpinning this process. The major quantitative trait locus ETHQV8.1 was previously identified as a key regulator of melon fruit ripening. Here, we narrowed down ETHQV8.1 to a precise genomic region containing a single gene, the transcription factor CmERF024. Functional validation using CRISPR/Cas9 knock-out plants unequivocally identified CmERF024 as the causal gene governing ETHQV8.1. The erf024 mutants exhibited suppression of ethylene production, leading to a significant delay and attenuation of fruit ripening. Integrative multi-omic analyses encompassing RNA-seq, DAP-seq, and DNase-seq revealed the association of CmERF024 with chromatin accessibility and gene expression dynamics throughout fruit ripening. Our data suggest CmERF024 as a novel regulator of climacteric fruit ripening in melon.

The maturation profile triggers differential expression of sugar metabolism genes in melon fruits

Melons are commercially important crops that requires specific quality attributes for successful commercialization, including accumulation of sugars, particularly sucrose. This trait can be influenced by various factors, such as the type of ripening. Cucumis melo L. is an ideal species for studying sugar metabolism because it has both climacteric and non-climacteric cultivars. Thus, this study aimed to examine the gene expression of sucrose metabolism candidates using RT-qPCR, in conjunction with postharvest physiological analyzes and high-performance liquid chromatography-based sugar quantification, in the melon cultivars 'Gaúcho' (climacteric) and 'Eldorado' (non-climacteric). The results showed that sucrose synthase 1 played a role in the synthesis and accumulation of sucrose in both cultivars, whereas sucrose synthase 2 was more highly expressed in 'Gaúcho', contributing to lower hexose content. Invertase inhibitor 1 was more highly expressed in 'Eldorado' and may be involved in sugar-induced maturation. Neutral α-galactosidase had distinct functions, playing a role in substrate synthesis for the growth of young 'Eldorado' fruits, whereas in mature 'Gaúcho' fruits it participated in the metabolism of raffinose family oligosaccharides for sucrose accumulation. The expression of trehalose-6-phosphate synthase genes indicated a greater involvement of these enzymes in the sugar regulation in 'Gaúcho' melons. These findings shed light on the intraspecific differences related to fruit quality attributes in different types of maturation and contribute to a deeper understanding of the underlying molecular mechanisms involved in the metabolism of sugars in melons, which can inform breeding programs aimed at improving fruit quality attributes in this crop.

Genome-wide identification of the OMT gene family in Cucumis melo L. and expression analysis under abiotic and biotic stress

Background

O-methyltransferase (OMT)-mediated O-methylation is a frequent modification that occurs during natural product biosynthesis, and it increases the diversity and stability of secondary metabolites. However, detailed genome-wide identification and expression analyses of OMT gene family members have not been performed in melons. In this study, we aimed to perform the genome-wide identification of OMT gene family members in melon to identify and clarify their actions during stress.

Methods

Genome-wide identification of OMT gene family members was performed using data from the melon genome database. The Cucumis melo OMT genes (CmOMTs) were then compared with the genes from two representative monocotyledons and three representative dicotyledons. The basic information, cis-regulatory elements in the promoter, predicted 3-D-structures, and GO enrichment results of the 21 CmOMTs were analyzed.

Results

In our study, 21 CmOMTs (named CmOMT1-21) were obtained by analyzing the melon genome. These genes were located on six chromosomes and divided into three groups composed of nine, six, and six CmOMTs based on phylogenetic analysis. Gene structure and motif descriptions were similar within the same classes. Each CmOMT gene contains at least one cis-acting element associated with hormone transport regulation. Analysis of cis-acting elements illustrated the potential role of CmOMTs in developmental regulation and adaptations to various abiotic and biotic stresses. The RNA-seq and quantitative real-time PCR (qRT-PCR) results indicated that NaCl stress significantly induced CmOMT6/9/14/18 and chilling and high temperature and humidity (HTH) stresses significantly upregulated CmOMT14/18. Furthermore, the expression pattern of CmOMT18 may be associated with Fusarium oxysporum f. sp. melonis race 1.2 (FOM1.2) and powdery mildew resistance. Our study tentatively explored the biological functions of CmOMT genes in various stress regulation pathways and provided a conceptual basis for further detailed studies of the molecular mechanisms.

Mutagenesis of the melon Prv gene by CRISPR/Cas9 breaks papaya ringspot virus resistance and generates an autoimmune allele with constitutive defense responses

The majority of plant disease resistance (R) genes encode nucleotide binding-leucine-rich repeat (NLR) proteins. In melon, two closely linked NLR genes, Fom-1 and Prv, were mapped and identified as candidate genes that control resistance to Fusarium oxysporum f.sp. melonis races 0 and 2, and to papaya ringspot virus (PRSV), respectively. In this study, we validated the function of Prv and showed that it is essential for providing resistance against PRSV infection. We generated CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9] mutants using Agrobacterium-mediated transformation of a PRSV-resistant melon genotype, and the T1 progeny proved susceptible to PRSV, showing strong disease symptoms and viral spread upon infection. Three alleles having 144, 154, and ~3 kb deletions, respectively, were obtained, all of which caused loss of resistance. Interestingly, one of the Prv mutant alleles, prvΔ154, encoding a truncated product, caused an extreme dwarf phenotype, accompanied by leaf lesions, high salicylic acid levels, and defense gene expression. The autoimmune phenotype observed at 25 °C proved to be temperature dependent, being suppressed at 32 °C. This is a first report on the successful application of CRISPR/Cas9 to confirm R gene function in melon. Such validation opens up new opportunities for molecular breeding of disease resistance in this important vegetable crop.

Identification of suitable reference genes for quantitative reverse transcription PCR in Luffa (Luffa cylindrica)

Reverse transcription real-time quantitative PCR is widely used to quantify gene expression. Reference genes are usually used as internal controls to measure the target gene expression level. To date, there is no consensus on the use of systematically validated reference genes in different tissues of Luffa. This study evaluated the expression stability of 11 candidate reference genes in different tissues using five algorithms (BestKeeper, comparative delta-Ct method, GeNorm, NormFinder, and RefFinder). Protein phosphatase 2A was the most stable gene, while alpha Tubulin was the least stable. The relative expression of ethylene-related genes in different tissues was also analyzed to reveal their role in sex determination. This study provides the basis for using suitable reference genes to evaluate targeted gene expression.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01182-8.

Expression Pattern of Entire Cytochrome P450 Genes and Response of Defenses in a Metabolic-Herbicide-Resistant Biotype of Polypogon fugax

Enhanced herbicide metabolism mediated by cytochrome P450s has been proposed as one of the major mechanisms of resistance to fenoxaprop-P-ethyl in a metabolic-herbicide-resistant biotype of Asia minor bluegrass (Polypogon fugax Nees ex Steud.). Upon pre-treatment with the P450 inhibitor piperonyl butoxide, a remarkable reduction in metabolic rates of the phytotoxic fenoxaprop-P has been observed in the resistant plants, implying that constitutive and/or fenoxaprop-P-ethyl-induced up-regulation of specific P450 isoforms are involved in the fenoxaprop-P-ethyl resistance. However, which P450 gene(s) were responsible for the metabolic resistance is still unknown. In this present study, based on the abundant gene resources of P. fugax established previously, a total of 48 putative P450 genes were isolated from the metabolic-herbicide-resistant plants and used for gene expression analysis. The most suitable reference genes for accurate normalization of real-time quantitative PCR data were first identified in P. fugax and recognized as actin (ACT), 18S rRNA (18S), and ribulose-1,5-bisphosphate carboxylase oxygenase (RUBP) under fenoxaprop-P-ethyl stress, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and elongation factor 1α (EF1α) under mesosulfuron-methyl stress, and ACT, EF1α, eukaryotic initiation factor 4a (EIF4A), and 25S rRNA (25S) at different growth stages. Expression analysis of the putative P450 genes revealed that six genes, respectively, annotated as CYP709B1, CYP71A1-4, CYP711A1, CYP78A9, P450-11, and P450-39 were up-regulated more than 10-fold in the resistant plants by fenoxaprop-P-ethyl treatment, and all of them exhibited constitutively and/or herbicide-induced higher transcript levels in the fenoxaprop-P-ethyl-resistant than in the susceptible plants. Three genes, respectively, annotated as CYPRO4, CYP313A4, and CYP51H11 constantly up-regulated in the resistant than in the susceptible plants after fenoxaprop-P-ethyl treatment. Up-regulated expressions of these specific P450 genes were consistent with the higher P450 contents determined in the resistant plants. These results will help to elucidate the mechanisms for P450-mediated metabolic-herbicide resistance in P. fugax as well as other grass weed species.

Comprehensive evaluation of candidate reference genes for quantitative real-time PCR-based analysis in Caucasian clover

The forage species Caucasian clover (Trifolium ambiguum M. Bieb.), a groundcover plant, is resistant to both cold and drought. However, reference genes for qRT-PCR-based analysis of Caucasian clover are lacking. In this study, 12 reference genes were selected on the basis of transcriptomic data. These genes were used to determine the most stably expressed genes in various organs of Caucasian clover under cold, salt and drought stress for qRT-PCR-based analysis. Reference gene stability was analyzed by geNorm, NormFinder, BestKeeper, the ∆Ct method and RefFinder. Under salt stress, RCD1 and PPIL3 were the most stable reference genes in the leaves, and NLI1 and RCD1 were the most stable references genes in the roots. Under low-temperature stress, APA and EFTu-GTP were the most stable reference genes in the leaves, and the RCD1 and NLI2 genes were highly stable in the roots. Under 10% PEG-6000 stress, NLI1 and NLI2 were highly stable in the leaves, and RCD1 and PPIL3 were the most stable in the roots. Overall, RCD1 and NLI2 were the most stable reference genes in organs under normal conditions and across all samples. The most and least stable reference genes were validated by assessing their appropriateness for normalization via WRKY genes.

Comparative transcriptome profiling reveals that brassinosteroid-mediated lignification plays an important role in garlic adaption to salt stress

Garlic (Allium sativum L.) is an economically important vegetable crop which is used worldwide for culinary and medicinal purposes. Soil salinity constrains the yield components of garlic. Understanding the responsive mechanism of garlic to salinity is crucial to improve its tolerance. To address this problem, two garlic cultivars differing in salt tolerance were used to investigate the long-term adaptive responses to salt stress at phenotype and transcriptome levels. Phenotypic analysis showed four-week salt stress significantly decreased the yield components of salt-sensitive cultivar. Transcriptomes of garlics were de novo assembled and mined for transcriptional activities regulated by salt stress. The results showed that photosynthesis, energy allocation, and secondary metabolism were commonly enriched in both sensitive and tolerant genotypes. Moreover, distinct responsive patterns were also observed between the two genotypes. Compared with the salt-tolerant genotype, most transcripts encoding enzymes in the phenylpropanoid biosynthesis pathway were coordinately down regulated in the salt-sensitive genotype, resulting in alternation of the content and composition of lignin. Meanwhile, transcripts encoding the enzymes in the brassinosteroid (BR) biosynthesis pathway were also systematically down regulated in the salt-sensitive genotypes. Taken together, these results suggested that BR-mediated lignin accumulation possibly plays an important role in garlic adaption to salt stress. These findings expand the understanding of responsive mechanism of garlic to salt stress.

Characterization, identification and expression profiling of genome-wide R-genes in melon and their putative roles in bacterial fruit blotch resistance

Background

Bacterial fruit blotch (BFB), a disease caused by Acidovorax citrulli, results in significant economic losses in melon. The causal QTLs and genes for resistance to this disease have yet to be identified. Resistance (R)-genes play vital roles in resistance to plant diseases. Since the complete genome sequence of melon is available and genome-wide identification of R-genes has been performed for this important crop, comprehensive expression profiling may lead to the identification of putative candidate genes that function in the response to BFB.

Results

We identified melon accessions that are resistant and susceptible to BFB through repeated bioassays and characterized all 70 R-genes in melon, including their gene structures, chromosomal locations, domain organizations, motif distributions, and syntenic relationships. Several disease resistance-related domains were identified, including NBS, TIR, LRR, CC, RLK, and DUF domains, and the genes were categorized based on the domains of their encoded proteins. In addition, we profiled the expression patterns of the genes in melon accessions with contrasting levels of BFB resistance at 12 h, 1 d, 3 d, and 6 d after inoculation with A. citrulli. Six R-genes exhibited consistent expression patterns (MELO3C023441, MELO3C016529, MELO3C022157, MELO3C022146, MELO3C025518, and MELO3C004303), with higher expression levels in the resistant vs. susceptible accession.

Conclusion

We identified six putative candidate R-genes against BFB in melon. Upon functional validation, these genes could be targeted for manipulation via breeding and biotechnological approaches to improve BFB resistance in melon in the future.

Transcriptome profiling of non-climacteric 'yellow' melon during ripening: insights on sugar metabolism

BACKGROUND

The non-climacteric 'Yellow' melon (Cucumis melo, inodorus group) is an economically important crop and its quality is mainly determined by the sugar content. Thus, knowledge of sugar metabolism and its related pathways can contribute to the development of new field management and post-harvest practices, making it possible to deliver better quality fruits to consumers.

RESULTS

The RNA-seq associated with RT-qPCR analyses of four maturation stages were performed to identify important enzymes and pathways that are involved in the ripening profile of non-climacteric 'Yellow' melon fruit focusing on sugar metabolism. We identified 895 genes 10 days after pollination (DAP)-biased and 909 genes 40 DAP-biased. The KEGG pathway enrichment analysis of these differentially expressed (DE) genes revealed that 'hormone signal transduction', 'carbon metabolism', 'sucrose metabolism', 'protein processing in endoplasmic reticulum' and 'spliceosome' were the most differentially regulated processes occurring during melon development. In the sucrose metabolism, five DE genes are up-regulated and 12 are down-regulated during fruit ripening.

CONCLUSIONS

The results demonstrated important enzymes in the sugar pathway that are responsible for the sucrose content and maturation profile in non-climacteric 'Yellow' melon. New DE genes were first detected for melon in this study such as invertase inhibitor LIKE 3 (CmINH3), trehalose phosphate phosphatase (CmTPP1) and trehalose phosphate synthases (CmTPS5, CmTPS7, CmTPS9). Furthermore, the results of the protein-protein network interaction demonstrated general characteristics of the transcriptome of young and full-ripe melon and provide new perspectives for the understanding of ripening.

Development of Molecular Marker Linked with Bacterial Fruit Blotch Resistance in Melon (Cucumis melo L.)

Bacterial fruit blotch (BFB) causes losses in melon marketable yield. However, until now, there has been no information about the genetic loci responsible for resistance to the disease or their pattern of inheritance. We determined the inheritance pattern of BFB resistance from a segregating population of 491 F₂ individuals raised by crossing BFB-resistant (PI 353814) and susceptible (PI 614596) parental accessions. All F₁ plants were resistant to Acidovorax citrulli strain KACC18782, and F₂ plants segregated with a 3:1 ratio for resistant and susceptible phenotypes, respectively, in a seedling bioassay experiment, indicating that BFB resistance is controlled by a monogenic dominant gene. In an investigation of 57 putative disease-resistance related genes across the melon genome, only the MELO3C022157 gene (encoding TIR-NBS-LRR domain), showing polymorphism between resistant and susceptible parents, revealed as a good candidate for further investigation. Cloning, sequencing and quantitative RT-PCR expression of the polymorphic gene MELO3C022157 located on chromosome 9 revealed multiple insertion/deletions (InDels) and single nucleotide polymorphisms (SNPs), of which the SNP A²⁰³⁵T in the second exon of the gene caused loss of the LRR domain and truncated protein in the susceptible accession. The InDel marker MB157-2, based on the large (504 bp) insertion in the first intron of the susceptible accession, was able to distinguish resistant and susceptible accessions among 491 F₂ and 22 landraces/inbred accessions with 98.17% and 100% detection accuracy, respectively. This novel PCR-based, co-dominant InDel marker represents a practical tool for marker-assisted breeding aimed at developing BFB-resistant melon accessions.

Reference genes for accurate evaluation of expression levels in Trichophyton interdigitale grown under different carbon sources, pH levels and phosphate levels

Tinea pedis is a type of dermatophytosis caused by anthropophilic keratinolytic fungi such as Trichophyton interdigitale. Quantitative reverse transcription PCR (RT-qPCR) is a reliable and reproducible technique for measuring changes in target gene expression across various biological conditions. A crucial aspect of accurate normalization is the choice of appropriate internal controls. To identify reference genes for accurate evaluation of expression levels in T. interdigitale, the transcription levels of eight candidate reference genes (adp-rf, β-act, ef1-α, gapdh, psm1, sdha, rpl2 and ubc) and one target gene (Tri m4) were analysed by RT-qPCR after growing the dermatophyte under different environmental conditions. The results obtained from expression stability evaluations with NormFinder, geNorm, BestKeeper, and RefFinder software demonstrated that adp-rf and psm1 were the most stable internal control genes across all experimental conditions. The present study constitutes the first report of the identification and validation of reference genes for RT-qPCR normalization for T. interdigitale grown under different environmental conditions resembling the conditions encountered by fungi during invasion of skin.

Selection of the optimal reference genes for expression analyses in different materials of Eriobotrya japonica

BACKGROUND

Loquat (Eriobotrya japonica) is a subtropical tree bearing fruit that ripens during late spring and early summer, which is the off-season for fruit production. The specific flowering habit of loquat, which starts in fall and ends in winter, has attracted an increasing number of researchers who believe that it may represent an ideal model for studying flowering shift adaptations to climate change in Rosaceae. These studies require an understanding of gene expression patterns within the fruit and other tissues of this plant. Although ACTINs (ACTs) have previously been used as reference genes (RGs) for gene expression studies in loquats, a comprehensive analysis of whether these RGs are optimal for normalizing RT-qPCR data has not been performed.

RESULTS

In this study, 11 candidate RGs (RIBOSOMAL-LIKE PROTEIN4 (RPL4), RIBOSOMAL-LIKE PROTEIN18 (RPL18), Histone H3.3 (HIS3), Alpha-tubulin-3 (TUA3), S-Adenosyl Methionine Decarboxylase (SAMDC), TIP41-like Family Protein (TIP41), (UDP)-glucose Pyrophosphorylase (UGPase), 18S ribosomal RNA (18S), Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH), Plasma Intrinsic Protein 2 (PIP2) and ACTIN(ACT)) were assessed to determine their expression stability in 23 samples from different tissues or organs of loquat. Integrated expression stability evaluations using five computational statistical methods (GeNorm, NormFinder, ΔCt, BestKeeper, and RefFinder) suggested that a RG set, including RPL4, RPL18, HIS3 and TUA3, was the most stable one across all of the tested loquat samples. The expression pattern of EjCDKB1;2 in the tested loquat tissues normalized to the selected RG set demonstrated its reliability.

CONCLUSIONS

This study reveals the reliable RGs for accurate normalization of gene expression in loquat. In addition, our findings demonstrate an efficient system for identifying the most effective RGs for different organs, which may be applied to related rosaceous crops.

Genome-wide identification of new reference genes for RT-qPCR normalization in CGMMV-infected Lagenaria siceraria

Lagenaria siceraria is an economically important cucurbitaceous crop, but suitable reference genes (RGs) to use when the plants are infected by cucumber green mottle mosaic virus (CGMMV) have not been determined. Sixteen candidate RGs of both leaf and fruit and 18 candidate RGs mostly from separate RNA-Seq datasets of bottle gourd leaf or fruit were screened and assessed by RT-qPCR. The expression stability of these genes was determined and ranked using geNorm, NormFinder, BestKeeper and RefFinder. Comprehensive analysis resulted in the selection of LsCYP, LsH3, and LsTBP as the optimal RGs for bottle gourd leaves, and LsP4H, LsADP, and LsTBP for fruits. LsWD, LsGAPDH, and LsH3 were optimal for use in both leaves and fruits under the infection of CGMMV. Isopentenyl transferase (IPT) and DNA-directed RNA polymerase (DdRP) were used to validate the applicability of the most stable identified RGs from bottle gourd in response to CGMMV. All the candidate RGs performed in RT-qPCR consistently with the data from the transcriptome database. The results demonstrated that LsWD, LsGAPDH and LsH3 were the most suitable internal RGs for the leaf, and LsH3, LsGAPDH, LsP4H and LsCYP for the fruit.

Identification of Optimal Reference Genes for Expression Analysis in Radish (Raphanus sativus L.) and Its Relatives Based on Expression Stability

Radish (Raphanus sativus) is an important cruciferous root crop with a close relationship to Chinese cabbage (Brassica rapa). RT-qPCR is used extensively to evaluate the expression levels of target genes, and accurate measurement of target gene expression with this method is determined by the valid reference genes used for data nomalization in different experimental conditions. Screening for appropriate reference genes with stable expression based on RT-qPCR data is important for gene expression and functional analysis research in radish and its relatives. However, many researches have thought that almost no single reference gene is widely suitable for all experimental conditions, and few researchers have paid attention to the validation of reference genes in radish gene expression analysis. In the present study, 12 candidate reference genes were selected for analysis. Their expression in 28 samples, including 20 radish samples from different organs and conditions, four Chinese cabbage organs and four organs of their distant hybrid, was assessed by RT-qPCR and then five software tools-ΔCt, geNorm, NormFinder, BestKeeper and RefFinder-were used to compare their expression stability. The results showed that the most suitable reference genes were different in different organs and conditions. GAPDH, DSS1, and UP2 were optimal reference genes for gene expression analysis in all organs and conditions in radish. UPR, GSNOR1, and ACTIN2/7 were the most stable reference genes in different radish organs. UP2 and GAPDH were suitable reference genes for radish pistil development studies. RPII, UBC9, and GAPDH had the most stable expression in radish under various stresses. DSS1, UP2, and TEF2 were the optimal reference genes for Chinese cabbage organs, whereas TUA was optimal for the distant hybrid. UP2, and TEF2 were appropriate reference genes for all of the samples together. The optimal reference genes we identified, UP2, GAPDH, UPR, and GSNOR1 were verified by normalizing the expression patterns of YAB3, RPL, and FUL. These results will provide important information for selecting target reference genes in different research contexts and improve the accuracy and precision of gene expression analysis for radish, Chinese cabbage and their distant hybrid.

Selection of Reliable Reference Genes for Gene Expression Analysis under Abiotic Stresses in the Desert Biomass Willow, Salix psammophila

Salix psammophila is a desert shrub willow that has extraordinary adaptation to abiotic stresses and plays an important role in maintaining local ecosystems. Moreover, S. psammophila is regarded as a promising biomass feedstock because of its high biomass yields and short rotation coppice cycle. However, few suitable reference genes (RGs) for quantitative real-time polymerase chain reaction (qRT-PCR) constrain the study on normalization of gene expression in S. psammophila until now. Here, we investigated the expression stabilities of 14 candidate RGs across tissue types and under four abiotic stress treatments, including heat, cold, salt, and drought treatments. After calculation of PCR efficiencies, three different software, NormFinder, geNorm, and BestKeeper were employed to analyze systematically the qRT-PCR data, and the outputs were merged by RankAggreg software. The optimal RGs selected for gene expression analysis were EF1α (Elongation factor-1 alpha) and OTU (OTU-like cysteine protease family protein) for different tissue types, UBC (Ubiquitin-conjugating enzyme E2) and LTA4H (Leukotriene A-4 hydrolase homolog) for heat treatment, HIS (Histone superfamily protein H3) and ARF2 (ADP-ribosylation factor 2) for cold treatment, OTU and ACT7 (Actin 7) for salt treatment, UBC and LTA4H for drought treatment. The expression of UBC, ARF2, and VHAC (V-type proton ATPase subunit C) varied the least across tissue types and under abiotic stresses. Furthermore, the relative genes expression profiles of one tissue-specific gene WOX1a (WUSCHEL-related homeobox 1a), and four stress-inducible genes, including Hsf-A2 (Heat shock transcription factors A2), CBF3 (C-repeat binding factor 3), HKT1 (High-Affinity K⁺ Transporter 1), and GST (Glutathione S-transferase), were conducted to confirm the validity of the RGs in this study. These results provided an important RGs application guideline for gene expression characterization in S. psammophila.