Selection of reference genes for gene expression studies in zucchini (Cucurbita pepo) using qPCR

Genome-Wide Identification of the WRKY Gene Family and Functional Characterization of CpWRKY5 in Cucurbita pepo

The WRKY gene family is crucial for regulating plant growth and development. However, the WRKY gene is rarely studied in naked kernel formation in hull-less Cucurbita pepo L. (HLCP), a natural mutant that lacks the seed coat. In this research, 76 WRKY genes were identified through bioinformatics-based methods in C. pepo, and their phylogenetics, conserved motifs, synteny, collinearity, and temporal expression during seed coat development were analyzed. The results showed that 76 CpWRKYs were identified and categorized into three main groups (I-III), with Group II further divided into five subgroups (IIa-IIe). Moreover, 31 segmental duplication events were identified in 49 CpWRKY genes. A synteny analysis revealed that C. pepo shared more collinear regions with cucumber than with melon. Furthermore, quantitative RT-PCR (qRT-PCR) results indicated the differential expression of CpWRKYs across different varieties, with notable variations in seed coat development between HLCP and CP being attributed to differences in CpWRKY5 expression. To investigate this further, CpWRKY5-overexpression tobacco plants were generated, resulting in increased lignin content and an upregulation of related genes, as confirmed by qRT-PCR. This study offers valuable insights for future functional investigations of CpWRKY genes and presents novel information for understanding the regulation mechanism of lignin synthesis.

Exogenous Application of dsRNA for Protection against Tomato Leaf Curl New Delhi Virus

Tomato leaf curl New Delhi virus (ToLCNDV) is an emerging plant pathogen, fast spreading in Asian and Mediterranean regions, and is considered the most harmful geminivirus of cucurbits in the Mediterranean. ToLCNDV infects several plant and crop species from a range of families, including Solanaceae, Cucurbitaceae, Fabaceae, Malvaceae and Euphorbiaceae. Up to now, protection from ToLCNDV infection has been achieved mainly by RNAi-mediated transgenic resistance, and non-transgenic fast-developing approaches are an urgent need. Plant protection by the delivery of dsRNAs homologous to a pathogen target sequence is an RNA interference-based biotechnological approach that avoids cultivating transgenic plants and has been already shown effective against RNA viruses and viroids. However, the efficacy of this approach against DNA viruses, particularly Geminiviridae family, is still under study. Here, the protection induced by exogenous application of a chimeric dsRNA targeting all the coding regions of the ToLCNDV DNA-A was evaluated in zucchini, an important crop strongly affected by this virus. A reduction in the number of infected plants and a delay in symptoms appearance, associated with a tendency of reduction in the viral titer, was observed in the plants treated with the chimeric dsRNA, indicating that the treatment is effective against geminiviruses but requires further optimization. Limits of RNAi-based vaccinations against geminiviruses and possible causes are discussed.

A Genome-Wide Identification and Comparative Analysis of the Heavy-Metal-Associated Gene Family in Cucurbitaceae Species and Their Role in Cucurbita pepo under Arsenic Stress

The heavy-metal-associated (HMA) proteins are a class of PB1-type ATPases related to the intracellular transport and detoxification of metals. However, due to a lack of information regarding the HMA gene family in the Cucurbitaceae family, a comprehensive genome-wide analysis of the HMA family was performed in ten Cucurbitaceae species: Citrullus amarus, Citrullus colocynthis, Citrullus lanatus, Citrullus mucosospermus, Cucumis melo, Cucumis sativus, Cucurbita maxima, Cucurbita moschata, Cucurbita pepo, and Legenaria siceraria. We identified 103 Cucurbit HMA proteins with various members, ranging from 8 (Legenaria siceraria) to 14 (Cucurbita pepo) across species. The phylogenetic and structural analysis confirmed that the Cucurbitaceae HMA protein family could be further classified into two major clades: Zn/Co/Cd/Pb and Cu/Ag. The GO-annotation-based subcellular localization analysis predicted that all HMA gene family members were localized on membranes. Moreover, the analysis of conserved motifs and gene structure (intron/exon) revealed the functional divergence between clades. The interspecies microsynteny analysis demonstrated that maximum orthologous genes were found between species of the Citrullus genera. Finally, nine candidate HMA genes were selected, and their expression analysis was carried out via qRT-PCR in root, leaf, flower, and fruit tissues of C. pepo under arsenic stress. The expression pattern of the CpeHMA genes showed a distinct pattern of expression in root and shoot tissues, with a remarkable expression of CpeHMA6 and CpeHMA3 genes from the Cu/Ag clade. Overall, this study provides insights into the functional analysis of the HMA gene family in Cucurbitaceae species and lays down the basic knowledge to explore the role and mechanism of the HMA gene family to cope with arsenic stress conditions.

In Vivo Antiphytoviral and Aphid Repellency Activity of Essential Oils and Hydrosols from Mentha suaveolens and Foeniculum vulgare to Control Zucchini Yellow Mosaic Virus and Its Vector Aphis gossypii

In recent years, natural compounds have gained attention in many fields due to their wide-range biological activity. In particular, essential oils and their associated hydrosols are being screened to control plant pests, exerting antiviral, antimycotic and antiparasitic actions. They are more quickly and cheaply produced and are generally considered safer for the environment and non-target organisms than conventional pesticides. In this study, we report the evaluation of the biological activity of two essential oils and their corresponding hydrosols obtained from Mentha suaveolens and Foeniculum vulgare in the control of zucchini yellow mosaic virus and its vector, Aphis gossypii, in Cucurbita pepo plants. The control of the virus was ascertained with treatments applied either concurrently with or after virus infection; choice tests were performed to verify repellency activity against the aphid vector. The results indicated that treatments could decrease virus titer as measured using real-time RT-PCR, while the experiments on the vector showed that the compounds effectively repelled aphids. The extracts were also chemically characterized using gas chromatography-mass spectrometry. Mentha suaveolens and Foeniculum vulgare hydrosol extracts mainly comprised fenchone and decanenitrile, respectively, while essential oils analysis returned a more complex composition, as expected.

In vivo Antiphytoviral Activity of Essential Oils and Hydrosols From Origanum vulgare, Thymus vulgaris, and Rosmarinus officinalis to Control Zucchini Yellow Mosaic Virus and Tomato Leaf Curl New Delhi Virus in Cucurbita pepo L

In the last decades, the interest in biological activity of natural compounds has been growing. In plant protection, essential oils have been reported to exhibit antiviral, antimycotic, and antiparasitic activities, and are regarded as promising for the formulation of safe antimicrobial agents. Attention has also been focused on hydrosols, the by-products of hydro-distillation of essential oils. Their production is easy, fast, and cheap, and they seem to arise less concern for human health than essential oils. Plant viruses represent a major concern for agricultural crops since no treatment compound is available for virus control. This work was aimed at evaluating the antiphytoviral effectiveness of treatments with three essential oils and corresponding hydrosols extracted from Origanum vulgare, Thymus vulgaris, and Rosmarinus officinalis on Cucurbita pepo plants infected by zucchini yellow mosaic virus or tomato leaf curl New Delhi virus. Treatments were applied either concurrently or after virus inoculation to ascertain an inhibition or curative activity, respectively. Symptoms were observed and samplings were performed weekly. Virus titer and expression levels of phenylalanine ammonia lyase gene (PAL) were measured on treated and untreated infected plants by real-time PCR. PAL gene plays an important role in plant defense response as it is involved in tolerance/resistance to phytopathogens. Results indicated that treatments were effective against tomato leaf curl New Delhi virus whether applied simultaneously with the inoculation or after. A major inhibition was observed with O. vulgare essential oil and hydrosol, resulting in 10^–4-fold decrease of virus titer 3 weeks after treatment. Curative activity gave maximum results with all three essential oils and T. vulgaris and R. officinalis hydrosols, recording from 10^–2-fold decrease to virus not detected 4 weeks after treatment. An induction of PAL gene expression was recorded at 12 d.p.i. and then was restored to the levels of untreated control. This allows to hypothesize an early plant defense response to virus infection, possibly boosted by treatments. Plant extracts’ composition was characterized by gas chromatography-mass spectrometry. Phenols were largely main components of O. vulgare and T. vulgaris extracts (carvacrol and thymol, respectively), while extracts from R. officinalis were based on monoterpene hydrocarbons (essential oil) and oxygenated monoterpenes (hydrosol).

Gene Expression in Zucchini Fruit Development

The study of fruit development in zucchini via gene expression has proven to be applicable in breeding programs. Phenotypic and transcriptomic studies of fruit set and parthenocarpy have been previously developed and some relevant genes have been reported. From these studies, three genotypes (MUCU-16, Whitaker, and Cavili) and six genes (CpAUX22, CpIAA4, CpIAMT-1, CpPIN5, CpCYCD6-1, and CpEXPLB1) were selected. The expression of these genes was analyzed in each genotype under three different treatments (pollination, auxin-treatment and non-treatment) during one week post anthesis. Also, a phenotyping analysis was conducted. The different nature of the samples and the genes selected allowed associations between different fruit traits and fruit development stages. There was a rapid response of CpAUX22 and CpIAA4 to the auxin treatment. Also, these genes and the CpIAMT-1 became more overexpressed in pollinated samples over time. The CpPIN5 gene increased its expression over time in all genotypes while CpCYCD6-1 was overexpressed in the early stages of fruit development in all samples. The CpEXPLB1 was highly up-regulated in non-treated samples, suggesting a relationship with fruit abortion. The overexpression of CpAUX22 and the non-overexpression of CpEXPLB1 in early stages may be associated with fruit growth in zucchini.

Stable reference gene selection for quantitative real-time PCR normalization in passion fruit (Passiflora edulis Sims.)

BACKGROUND

Passiflora edulis is a tropical fruit with high nutrient and medicinal values that is widely planted in southern China. However, the molecular biology of P. edulis has not been well studied. There are few reports regarding the choice of reference genes for gene expression studies of passion fruit.

METHODS AND RESULTS

By using three algorithms, implemented in geNorm, NormFinder and BestKeeper, we have selected ten candidate reference genes to explore their transcriptional expression stability in various tissues and under cold stress conditions. EF1 and HIS were stably expressed in five tissues. Ts and OTU were stably in vegetative organs. 50 S and Liom were stably in reproductive organs. The transcriptional abundance of EF1 and UBQ was stable in cold-treated and recovery treated leaf samples of P. edulis. In all samples, EF1 and Ts exhibited the highest expression stability. Evaluation of selected genes using simple statistical methods (ANOVA and post hoc analysis). Overall, EF1 emerged as the optimum reference gene for qRT-PCR normalize in P. edulis. In addition, the qRT-PCR analysis revealed that expression of ICE1 increases with the duration of cold treatment.

CONCLUSIONS

In this study, we successfully screened stable reference genes from 10 candidates in P. edulis and verified the results by analyzing the expression level of ICE1. The results provide reliable and effective reference genes for future research on gene expression analysis in P. edulis, and lay a foundation for follow-up research on functional genes in P. edulis.

Tolerance Screening for Phytophthora capsici and Fusarium solani f. sp. cucurbitae in Cucurbita spp. and Gene Expression Responses in Mutant Families

Cucurbita species can be affected by soil-borne pathogens, such as Phytopthora capsici and Fusarium solani f. sp. cucurbitae (Fsc). Diverse commercial and conserved lines of Cucurbita spp. were tested. C. pepo subsp. pepo genotypes showed the highest susceptibility to both pathogens. The tolerance to P. capsici and Fsc was then screened in a zucchini mutant population. Two M3 mutant lines (Cp107 and Cp116) with a high occurrence of tolerant individuals to Phytophthora capsici were obtained from a screening of 160 M2 mutant lines. The M3 lines presented higher tolerance than the background MUCU-16. Furthermore, in the inoculated samples, both mutants overexpressed CpDEF and expressed more CpPAL and CpChiIV than the susceptible control. It has been previously shown that this expression pattern could be associated with tolerance in the P. capsici - Cucurbita spp. pathosystem. The M3 lines obtained could be applied in breeding programs, as they are likely to be compatible with the highly susceptible C. pepo subsp. pepo genotype.

Genome-wide identification and expression analysis of the VQ gene family in Cucurbita pepo L

VQ protein is a plant specific protein, which plays an important role in plant growth and development and biological and abiotic stress response. This study aimed to systematically analyze for the first time the VQ of Cucurbita pepo and understand their expression patterns in response to different stimuli. Herein, 44 VQ genes were identified, which were divided into eight groups (I-VIII) based on phylogenetic analysis. Two genes (CpVQ1 and CpVQ2) could not be located on the chromosome, whereas the remaining CpVQ genes were randomly distributed on the chromosomes, except for chromosomes 15 and 18. Noteworthy, the main event driving the expansion of the VQ gene family was chromosome fragment duplication. Based on qRT-PCR analysis, VQ genes are expressed in different tissues, and VQ genes are differentially regulated under a variety of abiotic stresses and powdery mildew stress, indicating that they play an important role in plant stress response and other aspects. This report presents the first systematic analysis of VQ genes from C. pepo and provides a solid foundation for further research of the specific functions of VQ proteins.

Gene Expression in Cucurbita spp. Root and Crown during Phytophthora capsici Infection

Phytophtora capsici causes major diseases in cucurbit crops worldwide. In this study, we inoculated this pathogen into Cucurbita pepo subsp. pepo susceptible MUCU-16 and C. moschata tolerant M63. The gene expression of plant pathogenesis-related proteins chitinase (CpChiIV), lignin-forming peroxidase (CpLPOX), and defensin (CpDEF) and hormone-related enzymes salicylic acid (CpPAL) and ethylene (CpACO) was analyzed for two weeks post-inoculation in root and crown tissues. Differentially expressed genes were found between genotypes, tissues, days post-inoculation, and inoculated/non-inoculated samples. After inoculation, CpPAL and CpChiIV (crown) were downregulated in MUCU-16, while CpLPOX and CpDEF were upregulated in M63. In inoculated samples, higher expression changes were presented on days 10-14 than on day 3 for CpACO, CpLPOX, and CpDEF genes. Overexpression was higher for CpDEF compared to the other tested genes, indicating good suitability as a marker of biotic stress. The overexpression of CpDEF was higher in crown than in roots for both inoculated genotypes. The basal expression of CpPAL and CpDEF was higher in MUCU-16, but after inoculation, CpPAL and CpDEF gene expression were higher in M63. These changes suggest an association between CpDEF upregulation and tolerance, and between CpPAL downregulation and susceptibility.

CmRCC1 Gene From Pumpkin Confers Cold Tolerance in Tobacco by Modulating Root Architecture and Photosynthetic Activity

Low-temperature stress is the main limiting factor of cucurbit crop cultivation as it affects crop yield and quality. The identification of genes involved in cold tolerance is a crucial aspect of pumpkin rootstock breeding. Here, we examined the function of a pumpkin Regulator of Chromosome Condensation 1 (CmRCC1) gene in the root development and cold stress responses of tobacco (Nicotiana benthamiana). CmRCC1 expression was differentially induced in pumpkin root, stem, and leaf under cold stress. Transient transformation showed that CmRCC1 is located in the nucleus. CmRCC1 overexpression in tobacco increased the gravitropic set-point angle in lateral roots, as well as root diameter and volume. The expression of auxin polar transport factors, PIN1 and PIN3, decreased and increased in CmRCC1-overexpressed plants, respectively. Yeast two-hybrid verification and luciferase complementation imaging assay showed that CmRCC1 interacts with CmLAZY1. Furthermore, the decreases in maximum quantum yield of PS II, the effective quantum yield of PS II, and electron transfer rate and the increases in quantum yield of nonregulated energy dissipation and malondialdehyde content were compromised in transgenic plants compared with wild-type plants under cold stress. The results suggest that CmRCC1 plays an important role in the regulation of root architecture and positively modulates cold tolerance.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Revealing the Complexity of Sweepovirus-Deltasatellite-Plant Host Interactions: Expanded Natural and Experimental Helper Virus Range and Effect Dependence on Virus-Host Combination

Sweepoviruses are begomoviruses (genus Begomovirus, family Geminiviridae) with ssDNA genomes infecting sweet potato and other species of the family Convolvulaceae. Deltasatellites (genus Deltasatellite, family Tolecusatellitidae) are small-size non-coding DNA satellites associated with begomoviruses. In this study, the genetic diversity of deltasatellites associated with sweepoviruses infecting Ipomoea indica plants was analyzed by further sampling the populations where the deltasatellite sweet potato leaf curl deltasatellite 1 (SPLCD1) was initially found, expanding the search to other geographical areas in southern continental Spain and the Canary Islands. The sweepoviruses present in the samples coinfected with deltasatellites were also fully characterized by sequencing in order to define the range of viruses that could act as helper viruses in nature. Additionally, experiments were performed to assess the ability of a number of geminivirids (the monopartite tomato leaf deformation virus and the bipartite NW begomovirus Sida golden yellow vein virus, the bipartite OW begomovirus tomato leaf curl New Delhi virus, and the curtovirus beet curly top virus) to transreplicate SPLCD1 in their natural plant hosts or the experimental host Nicotiana benthamiana. The results show that SPLCD1 can be transreplicated by all the geminivirids assayed in N. benthamiana and by tomato leaf curl New Delhi virus in zucchini. The presence of SPLCD1 did not affect the symptomatology caused by the helper viruses, and its effect on viral DNA accumulation depended on the helper virus-host plant combination.

Selection and validation reference genes for qRT-PCR normalization in different cultivars during fruit ripening and softening of peach (Prunus persica)

Quantitative real-time PCR (qRT-PCR) has been emerged as an effective method to explore the gene function and regulatory mechanisms. However, selecting appropriate reference gene (s) is a prerequisite for obtaining accurate qRT-PCR results. Peach is one of important fruit in Rosaceae and is widely cultivated worldwide. In this study, to explore reliable reference gene (s) in peach with different types during fruit ripening and softening (S1-S4), nine candidate reference genes (EF-1α, GAPDH, TBP, UBC, eIF-4α, TUB-A, TUB-B, ACTIN, and HIS) were selected from the whole-genome data. Then, the expression levels of the nine selected genes were detected using qRT-PCR in three peach types, including 'Hakuho' (melting type), 'Xiacui' (stony hard type), 'Fantasia' and 'NJC108' (non-melting type) cultivars were detected using qRT-PCR. Four software (geNorm, NormFinder, BestKeeper and RefFinder) were applied to evaluate the expression stability of these candidate reference genes. Gene expression was characterized in different peach types during fruit ripening and softening stages. The overall performance of each candidate in all samples was evaluated. The Actin gene (ACTIN) was a suitable reference gene and displayed excellent stability in 'Total' set, 'Hakuho' samples, S3 and S4 fruit developmental stages. Ubiquitin C gene (UBC) showed the best stability in most independent samples, including 'Fantasia', 'NJC108', S2 sets. Elongation factor-1α gene (EF-1α) was the most unstable gene across the set of all samples, 'NJC108' and S2 sets, while showed the highest stability in 'Xiacui' samples. The stability of candidate reference genes was further verified by analyzing the relative expression level of ethylene synthase gene of Prunus persica (PpACS1) in fruit ripening and softening periods of 'Hakuho'. Taken together, the results from this study provide a basis for future research on the mining of important functional genes, expression patterns and regulatory mechanisms in peach.

Reference gene selection for qRT-PCR analyses of luffa (Luffa cylindrica) plants under abiotic stress conditions

Selecting suitable internal reference genes is an important prerequisite for the application of quantitative real-time PCR (qRT-PCR). However, no systematic studies have been conducted on reference genes in luffa. In this study, seven reference genes were selected, and their expression levels in luffa plants exposed to various simulated abiotic stresses [i.e., cold, drought, heat, salt, H₂O₂, and abscisic acid (ABA) treatments] were analyzed by qRT-PCR. The stability of the reference gene expression levels was validated using the geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The results indicated that EF-1α was the most stably expressed and suitable reference gene overall and for the heat, cold, and ABA treatments. Additionally, UBQ expression was stable following the salt treatment, whereas TUB was identified as a suitable reference gene for H₂O₂ and drought treatments. The reliability of the selected reference genes was verified by analyzing the expression of copper/zinc superoxide dismutase (Cu/Zn-SOD) gene in luffa. When the most unstable reference genes were used for data normalizations, the resulting expression patterns had obvious biases when compared with the expression patterns for the most ideal reference genes used alone or combined. These results will be conducive to more accurate quantification of gene expression levels in luffa.

Zucchini Plants Alter Gene Expression and Emission of (E)-β-Caryophyllene Following Aphis gossypii Infestation

Zucchini (Cucurbita pepo L.) is widely cultivated in temperate regions. One of the major production challenges is the damage caused by Aphis gossypii (Homoptera: Aphididae), a polyphagous aphid, which can negatively affect its host plant, both directly by feeding and indirectly by vectoring viruses. To gain insights into the transcriptome events that occur during the zucchini-aphid interaction and to understand the early-to-late defense response through gene expression profiles, we performed RNA-sequencing (RNA-Seq) on zucchini leaves challenged by A. gossypii (24, 48, and 96 h post-infestation; hpi). Data analysis indicated a complex and dynamic pattern of gene expression and a transient transcriptional reconfiguration that involved more than 700 differentially expressed genes (DEGs), including a large number of defense-related genes. The down-regulation of key genes of plant immunity, such as leucine-rich repeat (LRR) protein kinases, transcription factors, and genes associated with direct (i.e., protease inhibitors, cysteine peptidases, etc.) and indirect (i.e., terpene synthase) defense responses, suggests the aphid ability to manipulate plant immune responses. We also investigated the emission of volatile organic compounds (VOCs) from infested plants and observed a reduced emission of (E)-β-caryophyllene at 48 hpi, likely the result of aphid effectors, which reflects the down-regulation of two genes involved in the biosynthesis of terpenoids. We showed that (E)-β-caryophyllene emission was modified by the duration of plant infestation and by aphid density and that this molecule highly attracts Aphidius colemani, a parasitic wasp of A. gossypii. With our results we contributed to the identification of genes involved in cucurbit plant interactions with phloem feeders. Our findings may also help pave the way toward developing tolerant zucchini varieties and to identify molecules for sustainable management of harmful insect populations.

Selection of suitable reference genes for quantitative real time PCR in different Tulasnella isolates and orchid-fungus symbiotic germination system

Under natural conditions, mycorrhizal symbiosis accompanies nearly the entire life cycle of orchids from seed germination through to flowering and fruiting. Tulasnella-like orchid mycorrhizal fungi are the most common mycorrhizal fungi found in association with orchid species. Presently suitable reference genes have not been systematically selected for the quantification of gene expression via Real-Time Quantitative Reverse Transcription PCR (RT-qPCR). We evaluated 12 candidate Tulasnella genes in nine different Tulasnella isolates and in the Dendrobium-fungal symbiotic germination associations followed by statistical analysis using the programs Bestkeeper, geNorm, and Normfinder to analyze the expression stability of the individual genes. The results showed that the EF2, UBC, and PP2A genes had the highest rankings with relatively stable expression levels across the different genotypes and during the symbiotic seed germination process by the three programs, and may be suitable for RT-qPCR normalization. Furthermore, the gene encoding C-5 Sterol desaturase (C5SD) was selected to verify the reliability of EF2, UBC, and PP2A expression during the Tulasnella-Dendrobium symbiotic seed germination process. This study is the first systematic exploration of optimal reference genes for gene expression studies during the colonization of orchid seeds by the mycorrhizal fungus Tulasnella.

Expression Analysis of XTH in Stem Swelling of Stem Mustard and Selection of Reference Genes

Accurate analysis of gene expression requires selection of appropriate reference genes. In this study, we report analysis of eight candidate reference genes (ACTIN, UBQ, EF-1α, UBC, IF-4α, TUB, PP2A, and HIS), which were screened from the genome and transcriptome data in Brassica juncea. Four statistical analysis softwares geNorm, NormFinder, BestKeeper, and RefFinder were used to test the reliability and stability of gene expression of the reference genes. To further validate the stability of reference genes, the expression levels of two CYCD3 genes (BjuB045330 and BjuA003219) were studied. In addition, all genes in the xyloglucan endotransglucosylase/hydrolase (XTH) family were identified in B. juncea and their patterns at different periods of stem enlargement were analyzed. Results indicated that UBC and TUB genes showed stable levels of expression and are recommended for future research. In addition, XTH genes were involved in regulation of stem enlargement expression. These results provide new insights for future research aiming at exploring important functional genes, their expression patterns and regulatory mechanisms for mustard development.

Reference genes for qRT-PCR normalisation in different tissues, developmental stages, and stress conditions of Hypericum perforatum

Hypericum perforatum L. is a widely known medicinal herb used mostly as a remedy for depression because it contains high levels of naphthodianthrones, phloroglucinols, alkaloids, and some other secondary metabolites. Quantitative real-time PCR (qRT-PCR) is an optimized method for the efficient and reliable quantification of gene expression studies. In general, reference genes are used in qRT-PCR analysis because of their known or suspected housekeeping roles. However, their expression level cannot be assumed to remain stable under all possible experimental conditions. Thus, the identification of high quality reference genes is essential for the interpretation of qRT-PCR data. In this study, we investigated the expression of 14 candidate genes, including nine housekeeping genes (HKGs) (ACT2, ACT3, ACT7, CYP1, EF1-α, GAPDH, TUB-α, TUB-β, and UBC2) and five potential candidate genes (GSA, PKS1, PP2A, RPL13, and SAND). Three programs-GeNorm, NormFinder, and BestKeeper-were applied to evaluate the gene expression stability across four different plant tissues, four developmental stages and a set of abiotic stress and hormonal treatments. Integrating all of the algorithms and evaluations revealed that ACT2 and TUB-β were the most stable combination in different developmental stages samples and all of the experimental samples. ACT2, TUB-β, and EF1-α were identified as the three most applicable reference genes in different tissues and stress-treated samples. The majority of the conventional HKGs performed better than the potential reference genes. The obtained results will aid in improving the credibility of the standardization and quantification of transcription levels in future expression studies on H. perforatum.

Selection of reference genes for quantitative real-time PCR analysis in cucumber (Cucumis sativus L.), pumpkin (Cucurbita moschata Duch.) and cucumber-pumpkin grafted plants

Background

Quantitative real-time PCR (qRT-PCR) is a commonly used high-throughput technique to measure mRNA transcript levels. The accuracy of this evaluation of gene expression depends on the use of optimal reference genes. Cucumber-pumpkin grafted plants, made by grafting a cucumber scion onto pumpkin rootstock, are superior to either parent plant, as grafting conveys many advantages. However, although many reliable reference genes have been identified in both cucumber and pumpkin, none have been obtained for cucumber-pumpkin grafted plants.

Methods

In this work, 12 candidate reference genes, including eight traditional genes and four novel genes identified from our transcriptome data, were selected to assess their expression stability. Their expression levels in 25 samples, including three cucumber and three pumpkin samples from different organs, and 19 cucumber-pumpkin grafted samples from different organs, conditions, and varieties, were analyzed by qRT-PCR, and the stability of their expression was assessed by the comparative ΔCt method, geNorm, NormFinder, BestKeeper, and RefFinder.

Results

The results showed that the most suitable reference gene varied dependent on the organs, conditions, and varieties. CACS and 40SRPS8 were the most stable reference genes for all samples in our research. TIP41 and CACS showed the most stable expression in different cucumber organs, TIP41 and PP2A were the optimal reference genes in pumpkin organs, and CACS and 40SRPS8 were the most stable genes in all grafted cucumber samples. However, the optimal reference gene varied under different conditions. CACS and 40SRPS8 were the best combination of genes in different organs of cucumber-pumpkin grafted plants, TUA and RPL36Aa were the most stable in the graft union under cold stress, LEA26 and ARF showed the most stable expression in the graft union during the healing process, and TIP41 and PP2A were the most stable across different varieties of cucumber-pumpkin grafted plants. The use of LEA26, ARF and LEA26+ARF as reference genes were further verified by analyzing the expression levels of csaCYCD3;1, csaRUL, cmoRUL, and cmoPIN in the graft union at different time points after grafting.

Discussion

This work is the first report of appropriate reference genes in grafted cucumber plants and provides useful information for the study of gene expression and molecular mechanisms in cucumber-pumpkin grafted plants.

Lateral Root Initiation in the Parental Root Meristem of Cucurbits: Old Players in a New Position

While in most higher plants, including the model system Arabidopsis thaliana, the formation of lateral root primordia is induced in the elongation zone of the parental root, in seven plant families, including Cucurbitaceae, an alternative root branching mechanism is established such that lateral roots are initiated directly in the apical meristem of the parental root. In Arabidopsis, the transcription factor GATA23 and MEMBRANE-ASSOCIATED KINASE REGULATOR4 (MAKR4) are involved in the gene regulatory network of lateral root initiation. Among all marker genes examined, these are the earliest known marker genes up-regulated by auxin during lateral root initiation. In this study, putative functional orthologs of Arabidopsis GATA23 and MAKR4 were identified in cucumber (Cucumis sativus) and squash (Cucurbita pepo). Both cucurbits contained 26 genes encoding GATA family transcription factors and only one MAKR4 gene. Phylogenetic and transcriptional analysis of up-regulation by auxin led to the identification of GATA23 putative functional orthologs in Cucurbitaceae - CpGATA24 and CsGATA24. In squash, CpMAKR4 was up-regulated by naphthylacetic acid (NAA) and, similar to MAKR4 in Arabidopsis, indole-3-butyric acid (IBA). A detailed analysis of the expression pattern of CpGATA24 and CpMAKR4 in squash roots from founder cell specification until emergence of lateral root primordia was carried out using promoter-fluorescent reporter gene fusions and confocal microscopy. Their expression was induced in the protoxylem, and then expanded to founder cells in the pericycle. Thus, while the overall expression pattern of these genes was significantly different from that in Arabidopsis, in founder cells their expression was induced in the same order as in Arabidopsis. Altogether, these findings suggest that in Cucurbitaceae the putative functional orthologs of GATA23 and MAKR4 might play a role in founder cell specification and primordium positioning during lateral root initiation. The role of the protoxylem in auxin transport as a trigger of founder cells specification and lateral root initiation is discussed.

Selection of Reference Genes for qRT-PCR Analysis in Lentinula edodes after Hot-Air Drying

Volatile sulfur compounds gradually develop in Lentinula edodes after hot-air drying, and many genes are involved in the generation of these sulfur compounds. The expression stability of reference genes may vary in a particular experimental treatment when analyzing their expressions by quantitative real-time polymerase chain reaction (qRT-PCR). In this study, the expression profile of 17 candidate genes was assessed in L. edodes under treatment at 50 °C for 0, 1, 2, and 3 h, and the expression stability of each reference gene was analyzed by three statistical algorithms, including geNorm, NormFinder, and BestKeeper. Results indicated that the two optimal reference genes for mycelium and fruiting body were CAC and DAHP as well as CAC and NUP, respectively. Additionally, CAC and DAHP were found to be the two most stable reference genes across the mycelium and fruiting body set. Our results will provide a genetic foundation for further research on the metabolism genes of sulfur compounds in L. edodes.

The Effect of Methyl Jasmonate and Temperature on the Transient Expression of Recombinant Proteins in Cucurbita pepo L

The aim of this study is to assess the effect of methyl jasmonate (MeJA) and temperature on the valuable pharmaceuticals expression in a virus-mediated transient expression system, and so the Zuchini Yellow Mosaic Virus (ZYMV) based vector was used for transferring the GFP reporter gene and recombinant tissue plasminogen activator (rtPA) gene (K2S) to cucurbit (Cucurbita pepo L.). MeJA, temperature and time (days after inoculation), were evaluated as a factorial experiment in a completely randomized design (CRD). At first, the effect of all treatment combinations on GFP expression was assessed. At this step, the ELISA test was used to select the optimum treatment combination. ELISA method revealed the significant difference between applied treatments. The optimized treatment significantly increased the expression of rtPA compared to the control. The Real-Time PCR reaction for both GFP and rtPA genes showed no significant differences between optimum and control treatments, however, transcripts of the small subunit of RuBisCO were extremely down-regulated in optimum treatment condition. Reduction in RuBisCO expression at protein level was tangible under treatment condition based on the ELISA test. Therefore, it can be inferred that suppressing the expression of RuBisCO, probably resulted in higher access of expression system to free amino acids inside the cell. In this study, MeJA has been shown to be a positive factor, but the low temperature (17 °C), unlike previous studies, suppressed the expression of recombinant protein unexpectedly, probably due to the incompatibility of the viral construct with low temperature. In conclusion, the use of a suitable gene construct, which is not sensitive to temperature, is likely to result in a more favorable outcome.

Identification and validation of reference genes for gene expression studies in sweet osmanthus (Osmanthus fragrans) based on transcriptomic sequence data

Accurate normalized data is a primary requisite for quantifying gene expression using RT-qPCR technology. Despite this importance, however, suitable reference genes in Osmanthus fragrans are not available. In this study, seven potential candidate reference genes (OfL25-1, OfL25-10, OfRP2, OfTUA, OfTUB3, OfUBQ2 and Of18S) were evaluated to determine which one would be the most reliable reference genes. The expression levels of the candidate reference genes were analysed by RT-qPCR in flower, leaf, pedicel, blossom bud tissues, as well as in floral organs at different developmental stages.GeNormand NormFinderwere used to statistically analyse transcript variation.Results indicated that OfRP2 and OfL25-10 were the optimal reference genes for use in RT-qPCR when analysing different stages of floral development; while OfTUB3 and OfL25-1 were optimal across tissues. The selected reference genes were used to examineOfMYB1 expression. The results appeared to be useful for future gene expression analyses aiming to characterize developmental stages and tissues of O. fragrans.

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

Ectopic Expression of Pumpkin NAC Transcription Factor CmNAC1 Improves Multiple Abiotic Stress Tolerance in Arabidopsis

Drought, cold and salinity are the major environmental stresses that limit agricultural productivity. NAC transcription factors regulate the stress response in plants. Pumpkin (Cucurbita moschata) is an important cucurbit vegetable crop and it has strong resistance to abiotic stress; however, the biological functions of stress-related NAC genes in this crop are largely unknown. This study reports the function of CmNAC1, a stress-responsive pumpkin NAC domain protein. The CmNAC1-GFP fusion protein was transiently expressed in tobacco leaves for subcellular localization analysis, and we found that CmNAC1 is localized in the nucleus. Transactivation assay in yeast cells revealed that CmNAC1 functions as a transcription activator, and its transactivation domain is located in the C-terminus. CmNAC1 was ubiquitously expressed in different organs, and its transcript was induced by salinity, cold, dehydration, H₂O₂, and abscisic acid (ABA) treatment. Furthermore, the ectopic expression (EE) of CmNAC1 in Arabidopsis led to ABA hypersensitivity and enhanced tolerance to salinity, drought and cold stress. In addition, five ABA-responsive elements were enriched in CmNAC1 promoter. The CmNAC1-EE plants exhibited different root architecture, leaf morphology, and significantly high concentration of ABA compared with WT Arabidopsis under normal conditions. Our results indicated that CmNAC1 is a critical factor in ABA signaling pathways and it can be utilized in transgenic breeding to improve the abiotic stress tolerance of crops.

Identification of Appropriate Reference Genes for Normalization of miRNA Expression in Grafted Watermelon Plants under Different Nutrient Stresses

Watermelon (Citrullus lanatus) is a globally important crop belonging to the family Cucurbitaceae. The grafting technique is commonly used to improve its tolerance to stress, as well as to enhance its nutrient uptake and utilization. It is believed that miRNA is most likely involved in its nutrient-starvation response as a graft-transportable signal. The quantitative real-time reverse transcriptase polymerase chain reaction is the preferred method for miRNA functional analysis, in which reliable reference genes for normalization are crucial to ensure the accuracy. The purpose of this study was to select appropriate reference genes in scion (watermelon) and rootstocks (squash and bottle gourd) of grafted watermelon plants under normal growth conditions and nutrient stresses (nitrogen and phosphorus starvation). Under nutrient starvation, geNorm identified miR167c and miR167f as two most stable genes in both watermelon leaves and squash roots. miR166b was recommended by both geNorm and NormFinder as the best reference in bottle gourd roots under nutrient limitation. Expression of a new Cucurbitaceae miRNA, miR85, was used to validate the reliability of candidate reference genes under nutrient starvation. Moreover, by comparing several target genes expression in qRT-PCR analysis with those in RNA-seq data, miR166b and miR167c were proved to be the most suitable reference genes to normalize miRNA expression under normal growth condition in scion and rootstock tissues, respectively. This study represents the first comprehensive survey of the stability of miRNA reference genes in Cucurbitaceae and provides valuable information for investigating more accurate miRNA expression involving grafted watermelon plants.

Validation of Reference Genes for RT-qPCR Studies of Gene Expression in Preharvest and Postharvest Longan Fruits under Different Experimental Conditions

Reverse transcription quantitative PCR (RT-qPCR) as the accurate and sensitive method is use for gene expression analysis, but the veracity and reliability result depends on whether select appropriate reference gene or not. To date, several reliable reference gene validations have been reported in fruits trees, but none have been done on preharvest and postharvest longan fruits. In this study, 12 candidate reference genes, namely, CYP, RPL, GAPDH, TUA, TUB, Fe-SOD, Mn-SOD, Cu/Zn-SOD, 18SrRNA, Actin, Histone H3, and EF-1a, were selected. Expression stability of these genes in 150 longan samples was evaluated and analyzed using geNorm and NormFinder algorithms. Preharvest samples consisted of seven experimental sets, including different developmental stages, organs, hormone stimuli (NAA, 2,4-D, and ethephon) and abiotic stresses (bagging and girdling with defoliation). Postharvest samples consisted of different temperature treatments (4 and 22°C) and varieties. Our findings indicate that appropriate reference gene(s) should be picked for each experimental condition. Our data further showed that the commonly used reference gene Actin does not exhibit stable expression across experimental conditions in longan. Expression levels of the DlACO gene, which is a key gene involved in regulating fruit abscission under girdling with defoliation treatment, was evaluated to validate our findings. In conclusion, our data provide a useful framework for choice of suitable reference genes across different experimental conditions for RT-qPCR analysis of preharvest and postharvest longan fruits.

Expression Profile of Carotenoid Cleavage Dioxygenase Genes in Summer Squash (Cucurbita pepo L.)

Carotenoids are important dietary components that can be found in vegetable crops. The accumulation of these compounds in fruit and vegetables is altered by the activity of carotenoid cleavage dioxygenases (CCDs) enzymes that produce their degradation. The aim of this work was to study the possible implication of CCD genes in preventing carotenoid storage in the horticultural crop summer squash (Cucurbita pepo L.). The relationship between the presence of these compounds and gene expression for CCDs was studied in three varieties showing different peel and flesh colour. Expression analysis for the CCD genes CpNCED1, CpNCED2, CpNCED3, CpNCED9, CpCCD1, CpCCD4a, CpCCD4b and CpCCD8 was carried out on different organs and at several fruit developmental stages. The results showed that the CpCCD4a and CpCCD4b genes were highly expressed in the variety with lowest carotenoid content suggesting a putative role in carotenoid accumulation pattern in summer squash fruit.

Validation of endogenous reference genes in Buglossoides arvensis for normalizing RT-qPCR-based gene expression data

Selection of a stably expressed reference gene (RG) is an important step for generating reliable and reproducible quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) gene expression data. We, in this study, have sought to validate RGs for Buglossoides arvensis, a high nutraceutical value plant whose refined seed oil is entering the market under the commercial trade name Ahiflower™. This weed plant has received attention for its natural ability to significantly accumulate the poly-unsaturated fatty acid (PUFA) stearidonic acid (SDA, C18:4n-3) in its seeds, which is uncommon for most plant species. Ten candidate RGs (β-Act, 18S rRNA, EF-1a, α-Tub, UBQ, α-actin, CAC, PP2a, RUBISCO, GAPDH) were isolated from B. arvensis and TaqMan™ compliant primers/probes were designed for RT-qPCR analysis. Abundance of these gene transcripts was analyzed across different tissues and growth regimes. Two of the most widely used algorithms, geNorm and NormFinder, showed variation in expression levels of these RGs. However, combinatorial analysis of the results clearly identified CAC and α-actin as the most stable and unstable RG candidates, respectively. This study has for the first time identified and validated RGs in the non-model system B. arvensis, a weed plant projected to become an important yet sustainable source of dietary omega-3 PUFA.

Selection of reference genes for transcriptional analysis of edible tubers of potato (Solanum tuberosum L.)

Potato (Solanum tuberosum) yield has increased dramatically over the last 50 years and this has been achieved by a combination of improved agronomy and biotechnology efforts. Gene studies are taking place to improve new qualities and develop new cultivars. Reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) is a bench-marking analytical tool for gene expression analysis, but its accuracy is highly dependent on a reliable normalization strategy of an invariant reference genes. For this reason, the goal of this work was to select and validate reference genes for transcriptional analysis of edible tubers of potato. To do so, RT-qPCR primers were designed for ten genes with relatively stable expression in potato tubers as observed in RNA-Seq experiments. Primers were designed across exon boundaries to avoid genomic DNA contamination. Differences were observed in the ranking of candidate genes identified by geNorm, NormFinder and BestKeeper algorithms. The ranks determined by geNorm and NormFinder were very similar and for all samples the most stable candidates were C2, exocyst complex component sec3 (SEC3) and ATCUL3/ATCUL3A/CUL3/CUL3A (CUL3A). According to BestKeeper, the importin alpha and ubiquitin-associated/ts-n genes were the most stable. Three genes were selected as reference genes for potato edible tubers in RT-qPCR studies. The first one, called C2, was selected in common by NormFinder and geNorm, the second one is SEC3, selected by NormFinder, and the third one is CUL3A, selected by geNorm. Appropriate reference genes identified in this work will help to improve the accuracy of gene expression quantification analyses by taking into account differences that may be observed in RNA quality or reverse transcription efficiency across the samples.

Contribution of polyamines and other related metabolites to the maintenance of zucchini fruit quality during cold storage

In order to investigate the contribution of polyamines and related amino acids in the maintenance of zucchini fruit quality during cold storage, two varieties of Cucurbita pepo with different degrees of chilling tolerance were used, Natura (more tolerant) and Sinatra (moresensitive). After harvest, free putrescine levels decreased during storage at 20 °C, whereas in fruit kept at 4 °C this polyamine accumulated in both varieties, but with higher levels in the sensitive variety (Sinatra). This behavior suggests that putrescine is accumulated as a response to low temperature in zucchini fruit by stress-induced chilling injury, and not due to the postharvest storage itself. ADC activity responds quickly to chilling but sharply decreases after 14 days, whereas its expression remains high in both varieties. ODC activity takes over when the cold stress is relatively severe, as this activity was found to be much higher in Sinatra. ODCexpression also correlated with ODC activity. DAO activity increased in Natura fruit, and conversely decreased in Sinatra fruit during storage at 4 °C, whereas the proline content was higher in Natura and lower in Sinatra. Therefore, we suggest that putrescine degradation and proline accumulation contribute to the acquisition of chilling tolerance in zucchini fruit. GABA content decreased in both varieties, with a greater reduction in Natura fruit and less in Sinatra fruit. In addition, GABA transaminase showed a higher activity in Natura fruit than in Sinatra fruit during cold storage, suggesting that GABA catabolism could be involved in the tolerance to postharvest cold storage in zucchini fruit.

Selection of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in Brassica napus under various stress conditions

Data normalization is essential for reliable output of quantitative real-time reverse-transcription polymerase chain reaction (qRT-PCR) assays, as the unsuitable choice of reference gene(s), whose expression might be influenced by exogenous treatments in plant tissues, could cause misinterpretation of results. To date, no systematic studies on reference genes have been performed in stressed Brassica napus. In this study, we investigated the expression variations of nine candidate reference genes in 40 samples of B. napus leaves subjected to various exogenous treatments. Parallel analyses by geNorm and NormFinder revealed that optimal reference genes differed across the different sets of samples. The best-ranked reference genes were PP2A and TIP41 for salt stress, TIP41 and ACT7 for heavy metal (Cr(6+)) stress, PP2A and UBC21 for drought stress, F-box and SAND for cold stress, F-box and ZNF for salicylic acid stress, TIP41, ACT7, and PP2A for methyl jasmonate stress, TIP41 and ACT7 for abscisic acid stress, and TIP41, UBC21, and PP2A for Sclerotinia sclerotiorum stress. Two newly employed reference genes, TIP41 and PP2A, showed better performances, suggesting their suitability in multiple conditions. To further validate the suitability of the reference genes, the expression patterns of BnWRKY40 and BnMKS1 were studied in parallel. This study is the first systematic analysis of reference gene selection for qRT-PCR normalization in B. napus, an agriculturally important crop, under different stress conditions. The results will contribute toward more accurate and widespread use of qRT-PCR in gene analysis of the genus Brassica.

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

Watermelon is one of the major Cucurbitaceae crops and the recent availability of genome sequence greatly facilitates the fundamental researches on it. Quantitative real-time reverse transcriptase PCR (qRT–PCR) is the preferred method for gene expression analyses, and using validated reference genes for normalization is crucial to ensure the accuracy of this method. However, a systematic validation of reference genes has not been conducted on watermelon. In this study, transcripts of 15 candidate reference genes were quantified in watermelon using qRT–PCR, and the stability of these genes was compared using geNorm and NormFinder. geNorm identified ClTUA and ClACT, ClEF1α and ClACT, and ClCAC and ClTUA as the best pairs of reference genes in watermelon organs and tissues under normal growth conditions, abiotic stress, and biotic stress, respectively. NormFinder identified ClYLS8, ClUBCP, and ClCAC as the best single reference genes under the above experimental conditions, respectively. ClYLS8 and ClPP2A were identified as the best reference genes across all samples. Two to nine reference genes were required for more reliable normalization depending on the experimental conditions. The widely used watermelon reference gene 18SrRNA was less stable than the other reference genes under the experimental conditions. Catalase family genes were identified in watermelon genome, and used to validate the reliability of the identified reference genes. ClCAT1and ClCAT2 were induced and upregulated in the first 24 h, whereas ClCAT3 was downregulated in the leaves under low temperature stress. However, the expression levels of these genes were significantly overestimated and misinterpreted when 18SrRNA was used as a reference gene. These results provide a good starting point for reference gene selection in qRT–PCR analyses involving watermelon.

Screening suitable reference genes for normalization in reverse transcription quantitative real-time PCR analysis in melon

Melon (Cucumis melo. L) is not only an economically important cucurbitaceous crop but also an attractive model for studying many biological characteristics. Screening appropriate reference genes is essential to reverse transcription quantitative real-time PCR (RT-qPCR), which is key to many studies involving gene expression analysis. In this study, 14 candidate reference genes were selected, and the variations in their expression in roots and leaves of plants subjected to biotic stress, abiotic stress, and plant growth regulator treatment were assessed by RT-qPCR. The stability of the expression of the selected genes was determined and ranked using geNorm and NormFinder. geNorm identified the two most stable genes for each set of conditions: CmADP and CmUBIep across all samples, CmUBIep and CmRPL in roots, CmRAN and CmACT in leaves, CmADP and CmRPL under abiotic stress conditions, CmTUA and CmACT under biotic stress conditions, and CmRAN and CmACT under plant growth regulator treatments. NormFinder determined CmRPL to be the best reference gene in roots and under biotic stress conditions and CmADP under the other experimental conditions. CmUBC2 and CmPP2A were not found to be suitable under many experimental conditions. The catalase family genes CmCAT1, CmCAT2, and CmCAT3 were identified in melon genome and used as target genes to validate the reliability of identified reference genes. The catalase family genes showed the most upregulation 3 days after inoculation with Fusarium wilt in roots, after which they were downregulated. Their levels of expression were significantly overestimated when the unsuitable reference gene was used for normalization. These results not only provide guidelines for the selection of reference genes for gene expression analyses in melons but may also provide valuable information for studying the functions of catalase family genes in stress responses.

Selection and evaluation of novel reference genes for quantitative reverse transcription PCR (qRT-PCR) based on genome and transcriptome data in Brassica napus L

Selection of reference genes in Brassica napus, a tetraploid (4×) species, is a very difficult task without information on genome and transcriptome. By now, only several traditional reference genes which show significant expression differentiation under different conditions are used in B. napus. In the present study, based on genome and transcriptome data of the rapeseed Zhongshuang-11 cultivar, 14 candidate reference genes were screened for investigation in different tissues, cultivars, and treated conditions of B. napus. These genes were as follows: ELF5, ENTH, F-BOX7, F-BOX2, FYPP1, GDI1, GYF, MCP2d, OTP80, PPR, SPOC, Unknown1, Unknown2 and UBA. Among them, excluding GYF and FYPP1, another 12 genes, were identified to perform better than traditional reference genes ACTIN7 and GAPDH. To further validate the accuracy of the newly developed reference genes in normalization, expression levels of BnCAT1 (B. napus catalase 1) in different rapeseed tissues and seedlings under stress conditions were normalized by the three most stable reference genes PPR, GDI1, and ENTH and little difference existed in normalization results. To the best of our knowledge, this is the first time B. napus reference genes have been provided with the help of complete genome and transcriptome information. The new reference genes provided in this study are more accurate than previously reported reference genes in quantifying expression levels of B. napus genes.

Validation of suitable reference genes for gene expression analysis in the halophyte Salicornia europaea by real-time quantitative PCR

Real-time quantitative polymerase chain reaction (RT-qPCR), a reliable technique for quantifying gene expression, requires stable reference genes to normalize its data. Salicornia europaea, a stem succulent halophyte with remarkable salt resistance and high capacity for ion accumulation, has not been investigated with regards to the selection of appropriate reference genes for RT-qPCR. In this study, the expression of 11 candidate reference genes, GAPDH (Glyceraldehyde 3-phosphate dehydrogenase), Actin, α-Tub (α-tubulin), β-Tub (β-tubulin), EF1-α (Elongation factor 1-α), UBC (Ubiquitin-conjugating enzyme), UBQ (Polyubiquitin), CYP (Cyclophilin), TIP41 (TIP41-like protein), CAC (Clathrin adaptor complexes), and DNAJ (DnaJ-like protein), was analyzed in S. europaea samples, which were classified into groups according to various abiotic stresses (NaCl, nitrogen, drought, cold and heat), tissues and ages. Three commonly used software programs (geNorm, NormFinder and BestKeeper) were applied to evaluate the stability of gene expression, and comprehensive ranks of stability were generated by aggregate analysis. The results show that the relatively stable genes for each group are the following: (1) CAC and UBC for whole samples; (2) CAC and UBC for NaCl stress; (3) Actin and α-Tub for nitrogen treatment; (4) Actin and GAPDH for drought stress; (5) α-Tub and UBC for cold stress; (6) TIP41 and DNAJ for heat stress; (7) UBC and UBQ for different tissues; and (8) UBC and Actin for various developmental stages. These genes were validated by comparing transcriptome profiles. Using two stable reference genes was recommended in the normalization of RT-qPCR data. This study identifies optimal reference genes for RT-qPCR in S. europaea, which will benefit gene expression analysis under these conditions.

Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.)

BACKGROUND

We have identified a kind of parthenocarpy in zucchini squash which is associated with an incomplete andromonoecy, i.e. a partial conversion of female into bisexual flowers. Given that andromonoecy in this and other cucurbit species is caused by a reduction of ethylene production in the female flower, the associated parthenocarpic development of the fruit suggested the involvement of ethylene in fruit set and early fruit development.

RESULTS

We have compared the production of ethylene as well as the expression of 13 ethylene biosynthesis and signalling genes in pollinated and unpollinated ovaries/fruits of two cultivars, one of which is parthenocarpic (Cavili), while the other is non-parthenocarpic (Tosca). In the latter, unpollinated ovaries show an induction of ethylene biosynthesis and ethylene signal transduction pathway genes three days after anthesis, which is concomitant with the initiation of fruit abortion and senescence. Fruit set and early fruit development in pollinated flowers of both cultivars and unpollinated flowers of Cavili is coupled with low ethylene biosynthesis and signalling, which would also explain the partial andromonoecy in the parthenocarpic genotype. The reduction of ethylene production in the ovary cosegregates with parthenocarpy and partial andromonoecy in the selfing progeny of Cavili. Moreover, the induction of ethylene in anthesis (by ethephon treatments) reduced the percentage of bisexual parthenocarpic flowers in Cavili, while the inhibition of ethylene biosynthesis or response (by AVG and STS treatments) induces not only andromonoecy but also the parthenocarpic development of the fruit in both cultivars.

CONCLUSIONS

Results demonstrate that a reduction of ethylene production or signalling in the zucchini flower is able to induce fruit set and early fruit development, and therefore that ethylene is actively involved in fruit set and early fruit development. Auxin and TIBA treatments, inducing fruit set and early fruit development in this species, also inhibit ethylene production and the expression of ethylene biosynthesis and response genes. A model is presented that discusses the crosstalk between ethylene and auxin in the control of fruit set and early fruit development in zucchini squash.

Evaluation of reference genes for RT qPCR analyses of structure-specific and hormone regulated gene expression in Physcomitrella patens gametophytes

The use of the moss Physcomitrella patens as a model system to study plant development and physiology is rapidly expanding. The strategic position of P. patens within the green lineage between algae and vascular plants, the high efficiency with which transgenes are incorporated by homologous recombination, advantages associated with the haploid gametophyte representing the dominant phase of the P. patens life cycle, the simple structure of protonemata, leafy shoots and rhizoids that constitute the haploid gametophyte, as well as a readily accessible high-quality genome sequence make this moss a very attractive experimental system. The investigation of the genetic and hormonal control of P. patens development heavily depends on the analysis of gene expression patterns by real time quantitative PCR (RT qPCR). This technique requires well characterized sets of reference genes, which display minimal expression level variations under all analyzed conditions, for data normalization. Sets of suitable reference genes have been described for most widely used model systems including e.g. Arabidopsis thaliana, but not for P. patens. Here, we present a RT qPCR based comparison of transcript levels of 12 selected candidate reference genes in a range of gametophytic P. patens structures at different developmental stages, and in P. patens protonemata treated with hormones or hormone transport inhibitors. Analysis of these RT qPCR data using GeNorm and NormFinder software resulted in the identification of sets of P. patens reference genes suitable for gene expression analysis under all tested conditions, and suggested that the two best reference genes are sufficient for effective data normalization under each of these conditions.

Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses

Strawberry (Fragaria spp) is an emerging model for the development of basic genomics and recombinant DNA studies among rosaceous crops. Functional genomic and molecular studies involve relative quantification of gene expression under experimental conditions of interest. Accuracy and reliability are dependent upon the choice of an optimal reference control transcript. There is no information available on validated endogenous reference genes for use in studies testing strawberry-pathogen interactions. Thirteen potential pre-selected strawberry reference genes were tested against different tissues, strawberry cultivars, biotic stresses, ripening and senescent conditions, and SA/JA treatments. Evaluation of reference candidate's suitability was analyzed by five different methodologies, and information was merged to identify best reference transcripts. A combination of all five methods was used for selective classification of reference genes. The resulting superior reference genes, FaRIB413, FaACTIN, FaEF1α and FaGAPDH2 are strongly recommended as control genes for relative quantification of gene expression in strawberry. This report constitutes the first systematic study to identify and validate optimal reference genes for accurate normalization of gene expression in strawberry plant defense response studies.

Carotenogenic gene expression and carotenoid accumulation in three varieties of Cucurbita pepo during fruit development

The control of gene expression is a crucial regulatory mechanism in carotenoid accumulation of fruits and flowers. We investigated the role of transcriptional regulation of nine genes involved in the carotenoid biosynthesis pathway in three varieties of Cucurbita pepo with evident differences in fruit color. The transcriptional levels of the key genes involved in the carotenoid biosynthesis were higher in flower-, leaf-, and fruit skin tissues than flesh tissues. This correlated with higher concentration of carotenoid content in these tissues. The differential expression among the colored and white cultivars detected for some genes, such as LCYe, in combination with other regulatory mechanisms, could explain the large differences found in terms of carotenoid content among the three varieties. These results are a first step to elucidate carotenogenesis in C. pepo and demonstrate that, in general, regulation of the pathway genes is a critical factor that determines the accumulation of these compounds.

Evaluation of different pulverisation methods for RNA extraction in squash fruit: lyophilisation, cryogenic mill and mortar grinding

INTRODUCTION

Quality and integrity of RNA are critical for transcription studies in plant molecular biology. In squash fruit and other high water content crops, the grinding of tissue with mortar and pestle in liquid nitrogen fails to produce a homogeneous and fine powered sample desirable to ensure a good penetration of the extraction reagent.

OBJECTIVE

To develop an improved pulverisation method to facilitate the homogenisation process of squash fruit tissue prior to RNA extraction without reducing quality and yield of the extracted RNA.

METHODOLOGY

Three methods of pulverisation, each followed by the same extraction protocol, were compared. The first approach consisted of the lyophilisation of the sample in order to remove the excess of water before grinding, the second one used a cryogenic mill and the control one a mortar grinding of frozen tissue. The quality of the isolated RNA was tested by carrying out a quantitative real time downstream amplification.

RESULTS

In the three situations considered, mean values for A(260) /A(280) indicated minimal interference by proteins and RNA quality indicator (RQI) values were considered appropriate for quantitative real-time polymerase chain reaction (qRT-PCR) amplification. Successful qRT-PCR amplifications were obtained with cDNA isolated with the three protocols.

CONCLUSION

Both apparatus can improve and facilitate the grinding step in the RNA extraction process in zucchini, resulting in isolated RNA of high quality and integrity as revealed by qRT-PCR downstream application. This is apparently the first time that a cryogenic mill has been used to prepare fruit samples for RNA extraction, thereby improving the sampling strategy because the fine powder obtained represents a homogeneous mix of the organ tissue.

Selection and validation of reference genes for quantitative RT-PCR expression studies of the non-model crop Musa

Gene expression analysis by reverse transcriptase real-time or quantitative polymerase chain reaction (RT-qPCR) is becoming widely used for non-model plant species. Given the high sensitivity of this method, normalization using multiple housekeeping or reference genes is critical, and careful selection of these reference genes is one of the most important steps to obtain reliable results. In this study, reference genes commonly used for other plant species were investigated to identify genes displaying highly uniform expression patterns in different varieties, tissues, developmental stages, fungal infection, and osmotic stress conditions for the non-model crop Musa (banana and plantains). The expression stability of six candidate reference genes was tested on six different sample sets, and the results were analyzed using the publicly available algorithms geNorm and NormFinder. Our results show that variety, plant material, primer set, and gene identity can all influence the robustness and outcome of RT-qPCR analysis. In the case of Musa, a combination of three reference genes (EF1, TUB and ACT) can be used for normalization of gene expression data from greenhouse leaf samples. In the case of shoot meristem cultures, numerous combinations can be used because the investigated reference genes exhibited limited variability. In contrast, variability in expression of the reference genes was much larger among leaf samples from plants grown in vitro, for which the best combination of reference genes (L2 and ACT genes) is still suboptimal. Overall, our data confirm that the stability of candidate reference genes should be thoroughly investigated for each experimental condition under investigation. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-012-9711-1) contains supplementary material, which is available to authorized users.

Evaluation of RNA extraction methods and identification of putative reference genes for real-time quantitative polymerase chain reaction expression studies on olive (Olea europaea L.) fruits

Genome wide transcriptomic surveys together with targeted molecular studies are uncovering an ever increasing number of differentially expressed genes in relation to agriculturally relevant processes in olive (Olea europaea L). These data need to be supported by quantitative approaches enabling the precise estimation of transcript abundance. qPCR being the most widely adopted technique for mRNA quantification, preliminary work needs to be done to set up robust methods for extraction of fully functional RNA and for the identification of the best reference genes to obtain reliable quantification of transcripts. In this work, we have assessed different methods for their suitability for RNA extraction from olive fruits and leaves and we have evaluated thirteen potential candidate reference genes on 21 RNA samples belonging to fruit developmental/ripening series and to leaves subjected to wounding. By using two different algorithms, GAPDH2 and PP2A1 were identified as the best reference genes for olive fruit development and ripening, and their effectiveness for normalization of expression of two ripening marker genes was demonstrated.

Evaluation of reference genes for gene expression studies in radish (Raphanus sativus L.) using quantitative real-time PCR

Real-time quantitative reverse transcription PCR (RT-qPCR) is a rapid and reliable method for gene expression studies. Normalization based on reference genes can increase the reliability of this technique; however, recent studies have shown that almost no single reference gene is universal for all possible experimental conditions. In this study, eight frequently used reference genes were investigated, including Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Actin2/7 (ACT), Tubulin alpha-5 (TUA), Tubulin beta-1 (TUB), 18S ribosomal RNA (18SrRNA), RNA polymerase-II transcription factor (RPII), Elongation factor 1-b (EF-1b) and Translation elongation factor 2 (TEF2). Expression stability of candidate reference genes was examined across 27 radish samples, representing a range of tissue types, cultivars, photoperiodic and vernalization treatments, and developmental stages. The eight genes in these sample pools displayed a wide range of Ct values and were variably expressed. Two statistical software packages, geNorm and NormFinder showed that TEF2, RPII and ACT appeared to be relatively stable and therefore the most suitable for use as reference genes. These results facilitate selection of desirable reference genes for accurate gene expression studies in radish.

Evaluation of candidate reference genes for gene expression normalization in Brassica juncea using real time quantitative RT-PCR

The real time quantitative reverse transcription PCR (qRT-PCR) is becoming increasingly important to gain insight into function of genes. Given the increased sensitivity, ease and reproducibility of qRT-PCR, the requirement of suitable reference genes for normalization has become important and stringent. It is now known that the expression of internal control genes in living organism vary considerably during developmental stages and under different experimental conditions. For economically important Brassica crops, only a couple of reference genes are reported till date. In this study, expression stability of 12 candidate reference genes including ACT2, ELFA, GAPDH, TUA, UBQ9 (traditional housekeeping genes), ACP, CAC, SNF, TIPS-41, TMD, TSB and ZNF (new candidate reference genes), in a diverse set of 49 tissue samples representing different developmental stages, stress and hormone treated conditions and cultivars of Brassica juncea has been validated. For the normalization of vegetative stages the ELFA, ACT2, CAC and TIPS-41 combination would be appropriate whereas TIPS-41 along with CAC would be suitable for normalization of reproductive stages. A combination of GAPDH, TUA, TIPS-41 and CAC were identified as the most suitable reference genes for total developmental stages. In various stress and hormone treated samples, UBQ9 and TIPS-41 had the most stable expression. Across five cultivars of B. juncea, the expression of CAC and TIPS-41 did not vary significantly and were identified as the most stably expressed reference genes. This study provides comprehensive information that the new reference genes selected herein performed better than the traditional housekeeping genes. The selection of most suitable reference genes depends on the experimental conditions, and is tissue and cultivar-specific. Further, to attain accuracy in the results more than one reference genes are necessary for normalization.