Evaluation of Appropriate Reference Genes for Gene Expression Normalization during Watermelon Fruit Development

Selection and verification of reliable internal reference genes in stem development of herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) has emerged in the cut flower market due to its beautiful appearance. The bending flower stems caused by a lack of mechanical strength is the main problem restricting the development of the cut P. lactiflora industry. So it is of great worth to reveal the basis of stem development changes in P. lactiflora to improve its cut flower quality. Quantitative research on gene expression characteristics can provide clues for understanding their biological functions, and the screening of relatively stable expression genes is a prerequisite for the quantitative study of gene expression characteristics. Thus, it is necessary to find appropriate genes during stem development so as to analyze the qRT‒PCR results. In this study, 10 genes were screened, and these genes expressed stably in stems of different stem strengths at three different developmental stages. Then, their expressions were evaluated by RefFinder, BestKeeper, NormFinder, and GeNorm programs. The results demonstrated that γ-tubulin (γ-TUB) was the most suitable gene, followed by α-tubulin (α-TUB) and β-D-glucosidase (β-GUS), whereas histone H3 (His) was the least suitable gene. Besides, the temporal and spatial expression characteristics of PlCOMT1, the key gene concerned with the synthesis of cell wall fillers in P. lactiflora, were also used to evaluate the suitability of genes. Consequently, γ-TUB and α-TUB are the two best combinations during stem development, and their combination can be used for the stem development of P. lactiflora. These findings will provide a reference for the selection of genes related to stem development and the study of molecular mechanisms related to stem development in P. lactiflora.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01325-5.

Thioredoxin h2 inhibits the MPKK5-MPK3 cascade to regulate the CBF-COR signaling pathway in Citrullus lanatus suffering chilling stress

Thioredoxins (TRXs) are ubiquitous oxidoreductases and present as a multigenic family. TRXs determine the thiol redox balance, which is crucial for plants in the response to cold stress. However, limited knowledge is available about the role of TRXs in watermelon (Citrullus lanatus), which is highly sensitive to chilling stress in agricultural practice. Here, we identified 18 genes encoding 14 typical and 4 atypical TRXs from the watermelon genome, and found that ClTRX h2 localized at the plasma membrane was largely induced by chilling. Virus-induced gene silencing of ClTRX h2 resulted in watermelon plants that were more sensitive to chilling stress. We further found that ClTRX h2 physically interacted with mitogen-activated protein kinase kinase 5 (ClMPKK5), which was confirmed to phosphorylate and activate ClMPK3 in vitro, and the activation of ClMPK3 by ClMPKK5 was blocked by a point mutation of the Cys-229 residue to Ser in ClMPKK5. Additionally, ClTRX h2 inhibited the chilling-induced activation of ClMPK3, suggesting that the ClMPKK5-ClMPK3 cascade is regulated in a redox-dependent manner. We showed that ClMPK3-silenced plants had increased tolerance to chilling, as well as enhanced transcript abundances of the C-repeat/DREB binding factor (ClCBF) and cold-responsive (ClCOR) genes. Taken together, our results indicate that redox status mediated by ClTRX h2 inhibits ClMPK3 phosphorylation through the interaction between ClTRX h2 and ClMPKK5, which subsequently regulates the CBF-COR signaling pathway when submitted to chilling stress. Hence, our results provide a link between thiol redox balance and MAPK cascade signaling, revealing a conceptual framework to understand how TRX regulates chilling stress tolerance in watermelon.

Genome-wide analysis of IQD proteins and ectopic expression of watermelon ClIQD24 in tomato suggests its important role in regulating fruit shape

The plant-specific IQ67 domain (IQD) is the largest class of calmodulin targets found in plants, and plays an important role in many biological processes, especially fruit development processes. However, the functional role of IQD proteins in the development of watermelon (Citrullus lanatus) shape remains unknown, as the IQD protein family in watermelon has not been systematically characterized. Herein, we elucidated the gene structures, chromosomal locations, evolutionary divergence, and functions of 35 IQD genes in the watermelon genome. The transcript profiles and quantitative real-time PCR analysis at different stages of fruit development showed that the ClIQD24 gene was highly expressed on 0 days after pollination. Furthermore, we found that the ectopic overexpression of ClIQD24 promoted tomato fruit elongation, thereby revealing the significance of ClIQD24 in the progression of watermelon shape. Our study will serve as a reference for further investigations on the molecular mechanisms underlying watermelon fruit shape formation.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Selection and validation of reference genes for measuring gene expression in Piper species at different life stages using RT-qPCR analysis

The secondary metabolism of Piper species is known to produce a myriad of natural products from various biosynthetic pathways which, represent a rich source of previously uncharacterized chemical compounds. The determination of gene expression profiles in multiple tissue/organ samples could provide valuable clues towards understanding the potential biological functions of chemical changes in these plants. Studies on gene expression by RT-qPCR require particularly careful selection of suitable reference genes as a control for normalization. Here, we provide a study for the identification of reliable reference genes in P. arboreum, P. gaudichaudianum, P. malacophyllum, and P. tuberculatum, at two different life stages: 2-month-old seedlings and adult plants. To do this, annotated sequences were recovered from transcriptome datasets of the above listed Piper spp. These sequences were subjected to expression analysis using RT-qPCR, followed by analysis using the geNorm and NormFinder algorithms. A set of five genes were identified showing stable expression: ACT7 (Actin-7), Cyclophilin (Peptidyl-prolyl cis-trans isomerase), EF1α (Elongation factor 1-alpha), RNABP (RNA-binding protein), and UBCE (Ubiquitin conjugating enzyme). The universality of these genes was then validated using two target genes, ADC (arginine decarboxylase) and SAMDC (S-adenosylmethionine decarboxylase), which are involved in the biosynthesis of polyamines. We showed that normalization genes varied according to Piper spp., and we provide a list of recommended pairs of the best combination for each species. This study provides the first set of suitable candidate genes for gene expression studies in the four Piper spp. assayed, and the findings will facilitate subsequent transcriptomic and functional gene research.

The branchless gene Clbl in watermelon encoding a TERMINAL FLOWER 1 protein regulates the number of lateral branches

A SNP mutation in Clbl gene encoding TERMINAL FLOWER 1 protein is responsible for watermelon branchless. Lateral branching is one of the most important traits, which directly determines plant architecture and crop productivity. Commercial watermelon has the characteristics of multiple lateral branches, and it is time-consuming and labor-costing to manually remove the lateral branches in traditional watermelon cultivation. In our present study, a lateral branchless trait was identified in watermelon material WCZ, and genetic analysis revealed that it was controlled by a single recessive gene, which named as Clbl (Citrullus lanatus branchless). A bulked segregant sequencing (BSA-seq) and linkage analysis was conducted to primarily map Clbl on watermelon chromosome 4. Next-generation sequencing-aided marker discovery and a large mapping population consisting of 1406 F₂ plants were used to further map Clbl locus into a 9011-bp candidate region, which harbored only one candidate gene Cla018392 encoding a TERMINAL FLOWER 1 protein. Sequence comparison of Cla018392 between two parental lines revealed that there was a SNP detected from C to A in the coding region in the branchless inbred line WCZ, which resulted in a mutation from alanine (GCA) to glutamate (GAA) at the fourth exon. A dCAPS marker was developed from the SNP locus, which was co-segregated with the branchless phenotype in both BC₁ and F₂ population, and it was further validated in 152 natural watermelon accessions. qRT-PCR and in situ hybridization showed that the expression level of Cla018392 was significantly reduced in the axillary bud and apical bud in branchless line WCZ. Ectopic expression of ClTFL1 in Arabidopsis showed an increased number of lateral branches. The results of this study will be helpful for better understanding the molecular mechanism of lateral branch development in watermelon and for the development of marker-assisted selection (MAS) for new branchless watermelon cultivars.

Transcriptome regulation of carotenoids in five flesh-colored watermelons (Citrullus lanatus)

BACKGROUND

Fruit flesh color in watermelon (Citrullus lanatus) is a great index for evaluating the appearance quality and a key contributor influencing consumers' preferences. But the molecular mechanism of this intricate trait remains largely unknown. Here, the carotenoids and transcriptome dynamics during the fruit development of cultivated watermelon with five different flesh colors were analyzed.

RESULTS

A total of 13 carotenoids and 16,781 differentially expressed genes (DEGs), including 1295 transcription factors (TFs), were detected in five watermelon genotypes during the fruit development. The comprehensive accumulation patterns of carotenoids were closely related to flesh color. A number of potential structural genes and transcription factors were found to be associated with the carotenoid biosynthesis pathway using comparative transcriptome analysis. The differentially expressed genes were divided into six subclusters and distributed in different GO terms and metabolic pathways. Furthermore, we performed weighted gene co-expression network analysis and predicted the hub genes in six main modules determining carotenoid contents. Cla018406 (a chaperone protein dnaJ-like protein) may be a candidate gene for β-carotene accumulation and highly expressed in orange flesh-colored fruit. Cla007686 (a zinc finger CCCH domain-containing protein) was highly expressed in the red flesh-colored watermelon, maybe a key regulator of lycopene accumulation. Cla003760 (membrane protein) and Cla021635 (photosystem I reaction center subunit II) were predicted to be the hub genes and may play an essential role in yellow flesh formation.

CONCLUSIONS

The composition and contents of carotenoids in five watermelon genotypes vary greatly. A series of candidate genes were revealed through combined analysis of metabolites and transcriptome. These results provide an important data resource for dissecting candidate genes and molecular basis governing flesh color formation in watermelon fruit.

Resistant Sources and Genetic Control of Resistance to ToLCNDV in Cucumber

Tomato leaf curl New Delhi virus (ToLCNDV) is a severe threat for cucurbit production worldwide. Resistance has been reported in several crops, but at present, there are no described accessions with resistance to ToLCNDV in cucumber (Cucumis sativus). C. sativus var. sativus accessions were mechanically inoculated with ToLCNDV and screened for resistance, by scoring symptom severity, tissue printing, and PCR (conventional and quantitative). Severe symptoms and high load of viral DNA were found in plants of a nuclear collection of Spanish landraces and in accessions of C. sativus from different geographical origins. Three Indian accessions (CGN23089, CGN23423, and CGN23633) were highly resistant to the mechanical inoculation, as well as all plants of their progenies obtained by selfing. To study the inheritance of the resistance to ToLCNDV, plants of the CGN23089 accession were crossed with the susceptible accession BGV011742, and F₁ hybrids were used to construct segregating populations (F₂ and backcrosses), which were mechanically inoculated and evaluated for symptom development and viral load by qPCR. The analysis of the genetic control fit with a recessive monogenic inheritance model, and after genotyping with SNPs distributed along the C. sativus genome, a QTL associated with ToLCNDV resistance was identified in chromosome 2 of cucumber.

Rootstock mediates transcriptional regulation of citrulline metabolism in grafted watermelon

Citrulline is a non-essential amino acid, involved in key biological functions in plants and humans. Rootstocks have a major impact on citrulline accumulation in grafted watermelon. Information regarding rootstock induced changes in citrulline metabolism is elusive. To understand the regulatory mechanism, parallel changes in the expression profiles of citrulline metabolic genes and citrulline content of watermelon were monitored during the development of self-rooted watermelon and watermelon grafted onto pumpkin, wild and bottle gourd rootstocks. Results demonstrated that rootstocks regulated the expression profiles in different ways to influence the citrulline content. GAT, NAGPR, ASS3 ASS2 and Asl2 showed the negative correlation with citrulline content in pumpkin grafted watermelon. Pumpkin rootstock promoted the citrulline content by high down-regulation and synergistic effect of ASS2, ASS3, ASL1 and ASl2 genes. In wild grafted watermelon, citrulline was accumulated as a result of down regulation of GAT, NAGS and ASL2 genes, which showed an inverse correlation with citrulline. In gourd grafted watermelon, changes in citrulline content were observed to be linked with lower expressions of GAT, NAGK, ASS2, ASS3, ASL1 and ARG which were negatively correlated with citrulline content. Our study will provide the basis to understand the molecular mechanism of citrulline accumulation in various rootstocks.

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.

Selection of Suitable Reference Genes for qPCR Gene Expression Analysis of HepG2 and L02 in Four Different Liver Cell Injured Models

Quantitative real-time PCR (qPCR) has become a widely used approach to analyze the expression level of selected genes. However, owing to variations in cell types and drug treatments, a suitable reference gene should be selected according to special experimental design. In this study, we investigated the expression level of ten candidate reference genes in hepatoma carcinoma cell (HepG2) and human hepatocyte cell line (L02) treated with ethanol (EtOH), hydrogen peroxide (H₂O₂), acetaminophen (APAP), and carbon tetrachloride (CCl₄), respectively. To analyze raw cycle threshold values (Cp values) from qPCR run, three reference gene validation programs, including Bestkeeper, geNorm, and NormFinder, were used to evaluate the stability of ten candidate reference genes. The results showed that TATA-box binding protein (TBP) and tubulin beta 2a (TUBB2a) presented the highest stability for normalization under different treatments and were regarded as the most suitable reference genes of HepG2 and L02. In addition, this study not only identified the most stable reference genes of each treatment, but also suggested that β-actin (ACTB), glyceraldehade-3-phosphate dehydrogenase (GAPDH), tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta (YWHAZ), and beta-2 microglobulin (B2M) were the least stable reference genes in HepG2 and L02. This work was the first report to systematically explore the stability of reference genes in injured models of HepG2 and L02.

A cucumber green mottle mosaic virus vector for virus-induced gene silencing in cucurbit plants

BACKGROUND

Cucurbits produce fruits or vegetables that have great dietary importance and economic significance worldwide. The published genomes of at least 11 cucurbit species are boosting gene mining and novel breeding strategies, however genetic transformation in cucurbits is impractical as a tool for gene function validation due to low transformation efficiency. Virus-induced gene silencing (VIGS) is a potential alternative tool. So far, very few ideal VIGS vectors are available for cucurbits.

RESULTS

Here, we describe a new VIGS vector derived from cucumber green mottle mosaic virus (CGMMV), a monopartite virus that infects cucurbits naturally. We show that the CGMMV vector is competent to induce efficient silencing of the phytoene desaturase (PDS) gene in the model plant Nicotiana benthamiana and in cucurbits, including watermelon, melon, cucumber and bottle gourd. Infection with the CGMMV vector harboring PDS sequences of 69-300 bp in length in the form of sense-oriented or hairpin cDNAs resulted in photobleaching phenotypes in N. benthamiana and cucurbits by PDS silencing. Additional results reflect that silencing of the PDS gene could persist for over two months and the silencing effect of CGMMV-based vectors could be passaged.

CONCLUSIONS

These results demonstrate that CGMMV vector could serve as a powerful and easy-to-use tool for characterizing gene function, controlling viral pathogens or even performing resistance breeding in cucurbits. Moreover, this study will possess considerable important reference value for developing different viral vectors.

Comparative analysis of primary metabolites and transcriptome changes between ungrafted and pumpkin-grafted watermelon during fruit development

Grafting has been reported as a factor that influences fruit quality. However, a comprehensive study of the metabolic profile related to fruit quality and the underlying molecular mechanism in grafted watermelon has not been carried out. Metabolomics and transcriptome analysis were performed on both pumpkin-grafted watermelon and ungrafted watermelon at different developmental stages. In total, 56 primary metabolites were identified with either high or low abundance between ungrafted and pumpkin-grafted watermelon. The results indicated that ornithine, arginine, lysine (amino acids), glucose, sucrose, glucosamine (sugars), malic acid, fumaric acid and succinic acid (organic acids) were among the dominant metabolites influencing fruit quality. Additionally, comparative RNA sequence analysis on grafted and ungrafted watermelon yielded 729, 174, 128 and 356 differentially expressed genes at 10, 18, 26 and 34 days after pollination (DAP), respectively. Functional annotations of these genes indicated that grafting significantly altered the biological and metabolic processes related to fruit quality. Our comparative metabolomics and transcriptome analysis revealed that FBA2, FK, SuSy, SPS, IAI, AI and sugar transporter gene (SWT3b) might play a central role in the accumulation of glucose and sucrose, whereas higher malic acid content was attributed to high down regulation of ALMT13 and ALMT8 in pumpkin-grafted watermelon. Changes in the ornithine, glutamine, alanine, tyrosine, valine, asparagine, phenylalanine, arginine and tryptophan contents were consistent with the transcript level of their metabolic genes such as NAOD, GS, AGT, TaT, aDH1, OGDH, aDC, 4CL 1, PaL, CaT and two nitrate transporter genes (NRT1) in pumpkin-grafted watermelon. This study provides the basis for understanding the graft-responsive changes in the metabolic profile and regulatory mechanism related to fruit quality.

Genetic mapping and development of molecular markers for a candidate gene locus controlling rind color in watermelon

ClCG08G017810 (ClCGMenG) encoding a 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein is associated with formation of dark green versus light green rind color in watermelon. Rind color is an important agronomic trait in watermelon [Citrullus lanatus (Thunb.) Matsum. and Nakai], but the underlying molecular mechanism for this trait is not fully known. In the present study, we identified a single locus on chromosome 8 accounting for watermelon rind color (dark green vs. light green). Genetic analysis of F₁, F₂, and BC₁ populations derived from two parental lines (9904 with dark green rind and Handel with light green rind) revealed that the watermelon rind color (dark green vs. light green) is controlled by a single locus, and dark green is dominant to light green rind. Initial mapping revealed a region of interest spanning 2.07 Mb on chromosome 8. Genetic mapping with CAPS and SNP markers narrowed down the candidate region to 31.4 kb. Gene annotation of the corresponding region in the reference genome revealed the ClCG08G017810 gene sequence encoding the 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein. The sequence alignment of the candidate gene with the two parental lines suggested a nonsynonymous SNP mutation in the coding region of ClCG08G017810, converting an arginine (R) to glycine (G). The SNP might be associated with rind color of 103 watermelon germplasm lines investigated in this study. The qRT-PCR analysis revealed higher expression of ClCG08G017810 in dark green rind than in light green rind. Therefore, ClCG08G017810 is a candidate gene associated with watermelon rind color. The present study facilitates marker-assisted selection useful for the development of cultivars with desirable rind color.

Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress

Quantitative real-time reverse-transcriptase PCR (qRT-PCR) has been frequently used for detecting gene expression. To obtain reliable results, selection of suitable reference genes is a fundamental and necessary step. Garlic (Allium sativum), a member from Alliaceae family, has been used both as a food flavoring and as a traditional medicine. In the present study, garlic plants were exposed to salt stress (200 mM NaCl) for 0, 1, 4 and 12 h, and garlic roots, bulbs, and leaves were harvested for subsequent analysis. The expression stability of eight candidate reference genes, eukaryotic translation initiation factor 4α (eIF-4α), actin (ACTIN), tubulin β-7 (TUB7), TAP42-interacting protein of 41 kDa (TIP41), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), SAND family protein (SAND), elongation factor 1 alpha (EF-1α), and protein phosphatase 2A (PP2A) were evaluated by geNorm, NormFinder, and BestKeeper. All genes tested displayed variable expression profiles under salt stress. In the leaf and root group, ACTIN was the best reference gene for normalizing gene expression. In garlic clove, ACTIN and SAND were the least variable, and were suitable for gene expression studies under salt stress; these two genes also performed well in all samples tested. Based on our results, we recommend that it is essential to use specific reference genes in different situations to obtain accurate results. Using a combination of multiple stable reference genes, such as ACTIN and SAND, to normalize gene expression is encouraged. The results from the study will be beneficial for accurate determination of gene expression in garlic and other plants.

Systematized biosynthesis and catabolism regulate citrulline accumulation in watermelon

Citrulline, a non-protein amino acid, is present in large amounts in watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai Cucurbitaceae) fruits. Amino acid profiling of various tissues of cv. Charleston Gray during plant development confirmed progressive accumulation of citrulline only in the fruit flesh and rind tissues. Citrulline content was positively correlated with precursor (ornithine) and by-product (arginine) amino acids during fruit ripening. Genetic variation in the partitioning of citrulline and related amino acids in the flesh and rind tissues was confirmed in a sub-set of watermelon cultivars. No correlation was established between morphological fruit traits (size and rind properties) and citrulline content. To understand the regulation of citrulline accumulation, we investigated the expression of genes associated with its biosynthesis and catabolism in flesh and rind tissues during fruit development. The expression of ornithine carbamoyltransferase (OTC) involved in the ultimate step of citrulline synthesis remained steady in both tissues. The expression of N-acetylornithine aminotransferase (N-AOA) involved in the production of N-acetylornithine and N-acetylornithine deacetylase (AOD-3) involved in ornithine synthesis coincided with increasing accumulation of citrulline in flesh and rind tissues during fruit development. Down-regulation N-acetylornithine-glutamate acetyltransferase (N-AOGA) suggests the subordinate role of the non-cyclic pathway in citrulline synthesis. Eccentricity between citrulline accumulation and expression of carbamoyl phosphate synthases (CPS-1, CPS-2) during fruit development suggest that the localized synthesis of carbamoyl phosphates may not be required for citrulline synthesis. Most genes involved in citrulline break-down (Argininosuccinate synthases - ASS-1, ASS-2, and ASS-3, Argininosuccinate lyases - ASL-1, Ornithine decarboxylase - ODC, Arginine decarboxylase - ADC) were consistently down-regulated during fruit development.

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 duplication improves the resistance of watermelon root to salt stress

Salinity is a major abiotic stress factor that affects crop productivity. Roots play an important role in salt stress in plants. Watermelon is a salt-sensitive crop; however, tetraploid watermelon seedlings are more tolerant to salt stress than their homogenotype diploid ancestors. To obtain insights into the reasons underlying the differences in salt tolerance with respect to the ploidy of plants, self-grafted and cross-grafted diploid and tetraploid watermelon seedlings were exposed to 300 mM NaCl for 8 days. After the treatment, the tetraploid rootstock-grafted watermelon plants showed higher salt stress tolerance than the diploid plants. There were no significant differences in the physiological effects between the rootstocks with the same ploidy. The tetraploid rootstock-grafted watermelon plants exhibited higher net photosynthetic rate, leaf stomatal conductance and transpiration rate than the diploid rootstock-grafted watermelon plants throughout the salt treatment process. The activities of antioxidant enzymes and contents of osmoregulatory compounds in the roots were higher in the tetraploid rootstock-grafted watermelon plants than in the diploid plants during the entire salt response process. Higher Na⁺/K⁺ ratio was found in all parts of diploid rootstock-grafted watermelon, especially in the roots, K⁺ and Na⁺ were preferentially accumulated in the aerial parts (leaves and stem) than in the roots, which might be driven by the Na⁺/H⁺ antiporter, as evidenced by the higher transcript levels for SOS, PMA1, HKT1 and NHX1 in the roots. Taken together, our results suggest that genome duplication improves the resistance of watermelon root to salt stress.

Comprehensive evaluation of candidate reference genes for real-time quantitative PCR (RT-qPCR) data normalization in nutri-cereal finger millet [Eleusine Coracana (L.)]

Finger millet (Eleusine coracana L.) is an annual herbaceous self-pollinating C4 cereal crop of the arid and semi-arid regions of the world. Finger millet is a food security crop proven to have resilience to changing climate and scores very high in nutrition. In the current study, we have assessed sixteen candidate reference genes for their appropriateness for the normalization studies in finger millet subjected to experimental regimes and treatments. Ten candidate reference genes (GAPDH, β-TUB, CYP, EIF4α, TIP41, UBC, G6PD, S24, MACP and MDH) were cloned and six (ACT, ELF1α, PP2A, PT, S21 and TFIID) were mined from the NCBI database as well as from the literature. Expression stability ranking of the finger millet reference genes was validated using four different statistical tools i.e., geNorm, NormFinder, BestKeeper, ΔCt and RefFinder. From the study, we endorse MACP, CYP, EIF4α to be most stable candidate reference genes in all 'tissues', whereas PT, TFIID, MACP ranked high across genotypes, β-TUB, CYP, ELF1α were found to be best under abiotic stresses and 'all samples set'. The study recommends using minimum of two reference genes for RT-qPCR data normalizations in finger millet. All in all, CYP, β-TUB, and EF1α, in combination were found to be best for robust normalizations under most experimental conditions. The best and the least stable genes were validated for confirmation by assessing their appropriateness for normalization studies using EcNAC1 gene. The report provides the first comprehensive list of suitable stable candidate reference genes for nutritional rich cereal finger millet that will be advantageous to gene expression studies in this crop.

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.

Genetic mapping reveals a marker for yellow skin in watermelon (Citrullus lanatus L.)

As a diverse species, watermelon [Citrullus lanatus (Thunb.) Matsum. &Nakai var. lanatus] has different kinds of fruit sizes, shapes, flesh colors and skin colors. Skin color is among the major objectives for breeding. Yellow skin is an important trait in watermelon, but the underlying genetic mechanism is unknown. In this study, we identified a locus for yellow skin through BSA-seq and GWAS. A segregation analysis in F₂ and BC₁ populations derived from a cross of two inbred lines ‘94E1’(yellow skin) and ‘Qingfeng’(green skin) suggested that skin color is a qualitative trait. BSA-seq mapping confirmed the locus in the F₂ population, which was detected on chromosome 4 by GWAS among 330 varieties. Several major markers, namely, 15 CAPS markers, 6 SSR markers and 2 SNP markers, were designed to delimit the region to 59.8 kb region on chromosome 4. Utilizing the two populations consisting of 10 yellow and 10 green skin watermelons, we found a tightly linked functional SNP marker for the yellow skin phenotype. The application of this marker as a selection tool in breeding programs will help to improve the breeder’s ability to make selections at early stages of growth, thus accelerating the breeding program.

Comprehensive evaluation of blood-brain barrier-forming micro-vasculatures: Reference and marker genes with cellular composition

Primary brain microvessels (BrMV) maintain the cellular characters and molecular signatures as displayed in vivo, and serve as a vital tool for biomedical research of the blood-brain barrier (BBB) and the development/optimization of brain drug delivery. The variations of relative purities or cellular composition among different BrMV samples may have significant consequences in data interpretation and research outcome, especially for experiments with high-throughput genomics and proteomics technologies. In this study, we aimed to identify suitable reference gene (RG) for accurate normalization of real-time RT-qPCR analysis, and determine the proper marker genes (MG) for relative purity assessment in BrMV samples. Out of five housekeeping genes, β-actin was selected as the most suitable RG that was validated by quantifying mRNA levels of alpha-L-iduronidase in BrMV isolated from mice with one or two expressing alleles. Four marker genes highly/selectively expressed in BBB-forming capillary endothelial cells were evaluated by RT-qPCR for purity assessment, resulting in Cldn5 and Pecam1 as most suitable MGs that were further confirmed by immunofluorescent analysis of cellular components. Plvap proved to be an indicator gene for the presence of fenestrated vessels in BrMV samples. This study may contribute to the building blocks toward overarching research needs on the blood-brain barrier.

Genetic mapping reveals a candidate gene (ClFS1) for fruit shape in watermelon (Citrullus lanatus L.)

A 159 bp deletion in ClFS1 gene encoding IQD protein is responsible for fruit shape in watermelon. Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is known for its rich diversity in fruit size and shape. Fruit shape has been one of the major objectives of watermelon breeding. However, the candidate genes and the underlying genetic mechanism for such an important trait in watermelon are unknown. In this study, we identified a locus on chromosome 3 of watermelon genome controlling fruit shape. Segregation analysis in F₂ and BC₁ populations derived from a cross between two inbred lines "Duan125" (elongate fruit) and "Zhengzhouzigua" (spherical fruit) suggests that fruit shape of watermelon is controlled by a single locus and elongate fruit (OO) is incompletely dominant to spherical fruit (oo) with the heterozygote (Oo) being oval fruit. GWAS profiles among 315 accessions identified a major locus designated on watermelon chromosome 3, which was confirmed by BSA-seq mapping in the F₂ population. The candidate gene was mapped to a region 46 kb on chromosome 3. There were only four genes present in the corresponding region in the reference genome. Four candidate genes were sequenced in this region, revealing that the CDS of Cla011257 had a 159 bp deletion which resulted in the omission of 53 amino acids in elongate watermelon. An indel marker was derived from the 159 bp deletion to test the F₂ population and 105 watermelon accessions. The results showed that Cla011257 cosegregated with watermelon fruit shape. In addition, the Cla011257 expression was the highest at ovary formation stage. The predicted protein of the Cla011257 gene fitted in IQD protein family which was reported to have association with cell arrays and Ca²⁺-CaM signaling modules. Clear understanding of the genes facilitating the fruit shape along with marker association selection will be an effective way to develop new cultivars.

Identification and validation of reference genes for quantitative real-time PCR studies in long yellow daylily, Hemerocallis citrina Borani

Gene expression analysis using reverse transcription quantitative real-time PCR (RT-qPCR) requires the use of reference gene(s) in the target species. The long yellow daylily, Hemerocallis citrina Baroni. is rich in beneficial secondary metabolites and is considered as a functional vegetable. It is widely cultivated and consumed in East Asian countries. However, reference genes for use in RT-qPCR in H. citrina are not available. In the present study, six potential reference genes, actin (ACT), AP-4 complex subunit (AP4), tubulin (TUB), ubiquitin (UBQ), 18S and 60S ribosomal RNA, were selected and their expression stability in different developmental stages, organs and accessions was evaluated using four statistical software packages (geNorm, NormFinder, BestKeeper, and RefFinder). For commercial flower buds of different landraces, the combination of 60S, TUB, and AP4 was appropriate whereas ACT and 60S was suitable for normalization of different organs. In addition, AP4 exhibited the most stable expression in flower buds among different developmental stages. UBQ was less stable than the other reference genes under the experimental conditions except under different organs was 18S. The relative expression levels of two genes, primary-amine oxidase (HcAOC3) and tyrosine aminotransferase (HcTAT) which play important roles in alkaloid biosynthesis were also examined in different organs of the 'Datong' landrace, which further confirmed the results of selected reference genes. This is the first report to evaluate the stability of reference genes in the long yellow daylily that can serve as a foundation for RT-qPCR analysis of gene expression in this species.

Identification and validation of reference genes for quantitative real-time PCR studies in long yellow daylily, Hemerocallis citrina Borani

Gene expression analysis using reverse transcription quantitative real-time PCR (RT-qPCR) requires the use of reference gene(s) in the target species. The long yellow daylily, Hemerocallis citrina Baroni. is rich in beneficial secondary metabolites and is considered as a functional vegetable. It is widely cultivated and consumed in East Asian countries. However, reference genes for use in RT-qPCR in H. citrina are not available. In the present study, six potential reference genes, actin (ACT), AP-4 complex subunit (AP4), tubulin (TUB), ubiquitin (UBQ), 18S and 60S ribosomal RNA, were selected and their expression stability in different developmental stages, organs and accessions was evaluated using four statistical software packages (geNorm, NormFinder, BestKeeper, and RefFinder). For commercial flower buds of different landraces, the combination of 60S, TUB, and AP4 was appropriate whereas ACT and 60S was suitable for normalization of different organs. In addition, AP4 exhibited the most stable expression in flower buds among different developmental stages. UBQ was less stable than the other reference genes under the experimental conditions except under different organs was 18S. The relative expression levels of two genes, primary-amine oxidase (HcAOC3) and tyrosine aminotransferase (HcTAT) which play important roles in alkaloid biosynthesis were also examined in different organs of the 'Datong' landrace, which further confirmed the results of selected reference genes. This is the first report to evaluate the stability of reference genes in the long yellow daylily that can serve as a foundation for RT-qPCR analysis of gene expression in this species.

Morphological observation, RNA-Seq quantification, and expression profiling: novel insight into grafting-responsive carotenoid biosynthesis in watermelon grafted onto pumpkin rootstock

Watermelon is an important and economical horticultural crop in China, where ~20% of the plants are grafted. The development of grafted watermelon fruit involves a diverse range of gene interactions that results in dynamic changes in fruit. However, the molecular mechanisms underlying grafting-induced fruit quality change are unclear. In the present study, we measured the lycopene content by high-performance liquid chromatography and used RNA-Seq (quantification) to perform a genome-wide transcript analysis of fruits from watermelon grafted onto pumpkin rootstock (pumpkin-grafted watermelon, PGW), self-grafted watermelon (SGW), and non-grafted watermelon (NGW). The results showed variation in the lycopene content in the flesh of PGW fruits, first increasing and then decreasing in the four stages, which was different from the trend in the flesh of NGW and SGW fruits. The transcriptome profiling data provided new information on the grafting-induced gene regulation of lycopene biosynthesis during fruit growth and development. The expression levels of 33 genes from 8 gene families (GGPS, PSY, PDS, ZDS, CRTISO, LCYb, LCYe, and CHY) related to lycopene biosynthesis, which play critical roles in fruit coloration and contribute significantly to fruit phytonutrient values, were monitored during the four periods of fruit development in watermelon. Compared with those of NGW and SGW, 14 genes were differentially expressed in PGW during fruit development, suggesting that these genes possibly help to mediate lycopene biosynthesis in grafted watermelon fruit. Our work provides some novel insights into grafting-responsive carotenoid metabolism and its potential roles during PGW fruit development and ripening.

Identification of Optimal Reference Genes for Normalization of qPCR Analysis during Pepper Fruit Development

Due to its high sensitivity and reproducibility, quantitative real-time PCR (qPCR) is practiced as a useful research tool for targeted gene expression analysis. For qPCR operations, the normalization with suitable reference genes (RGs) is a crucial step that eventually determines the reliability of the obtained results. Although pepper is considered an ideal model plant for the study of non-climacteric fruit development, at present no specific RG have been developed or validated for the qPCR analyses of pepper fruit. Therefore, this study aimed to identify stably expressed genes for their potential use as RGs in pepper fruit studies. Initially, a total of 35 putative RGs were selected by mining the pepper transcriptome data sets derived from the PGP (Pepper Genome Platform) and PGD (Pepper Genome Database). Their expression stabilities were further measured in a set of pepper (Capsicum annuum L. var. 007e) fruit samples, which represented four different fruit developmental stages (IM: Immature; MG: Mature green; B: Break; MR: Mature red) using the qPCR analysis. Then, based on the qPCR results, three different statistical algorithms, namely geNorm, Normfinder, and boxplot, were chosen to evaluate the expression stabilities of these putative RGs. It should be noted that nine genes were proven to be qualified as RGs during pepper fruit development, namely CaREV05 (CA00g79660); CaREV08 (CA06g02180); CaREV09 (CA06g05650); CaREV16 (Capana12g002666); CaREV21 (Capana10g001439); CaREV23 (Capana05g000680); CaREV26 (Capana01g002973); CaREV27 (Capana11g000123); CaREV31 (Capana04g002411); and CaREV33 (Capana08g001826). Further analysis based on geNorm suggested that the application of the two most stably expressed genes (CaREV05 and CaREV08) would provide optimal transcript normalization in the qPCR experiments. Therefore, a new and comprehensive strategy for the identification of optimal RGs was developed. This strategy allowed for the effective normalization of the qPCR analysis of the pepper fruit development at the whole pepper genome level. This study not only explored the optimal RGs specific for studying pepper fruit development, but also introduced a referable strategy of RG mining which could potentially be implicated in other plant species.

Genome-Wide Identification and Evaluation of Reference Genes for Quantitative RT-PCR Analysis during Tomato Fruit Development

Gene expression analysis in tomato fruit has drawn increasing attention nowadays. Quantitative real-time PCR (qPCR) is a routine technique for gene expression analysis. In qPCR operation, reliability of results largely depends on the choice of appropriate reference genes (RGs). Although tomato is a model for fruit biology study, few RGs for qPCR analysis in tomato fruit had yet been developed. In this study, we initially identified 38 most stably expressed genes based on tomato transcriptome data set, and their expression stabilities were further determined in a set of tomato fruit samples of four different fruit developmental stages (Immature, mature green, breaker, mature red) using qPCR analysis. Two statistical algorithms, geNorm and Normfinder, concordantly determined the superiority of these identified putative RGs. Notably, SlFRG05 (Solyc01g104170), SlFRG12 (Solyc04g009770), SlFRG16 (Solyc10g081190), SlFRG27 (Solyc06g007510), and SlFRG37 (Solyc11g005330) were proved to be suitable RGs for tomato fruit development study. Further analysis using geNorm indicate that the combined use of SlFRG03 (Solyc02g063070) and SlFRG27 would provide more reliable normalization results in qPCR experiments. The identified RGs in this study will be beneficial for future qPCR analysis of tomato fruit developmental study, as well as for the potential identification of optimal normalization controls in other plant species.

Evaluation of Reference Genes for Quantitative Real-Time PCR Analysis of the Gene Expression in Laticifers on the Basis of Latex Flow in Rubber Tree (Hevea brasiliensis Muell. Arg.)

Latex exploitation-caused latex flow is effective in enhancing latex regeneration in laticifer cells of rubber tree. It should be suitable for screening appropriate reference gene for analysis of the expression of latex regeneration-related genes by quantitative real-time PCR (qRT-PCR). In the present study, the expression stability of 23 candidate reference genes was evaluated on the basis of latex flow by using geNorm and NormFinder algorithms. Ubiquitin-protein ligase 2a (UBC2a) and ubiquitin-protein ligase 2b (UBC2b) were the two most stable genes among the selected candidate references in rubber tree clones with differential duration of latex flow. The two genes were also high-ranked in previous reference gene screening across different tissues and experimental conditions. By contrast, the transcripts of latex regeneration-related genes fluctuated significantly during latex flow. The results suggest that screening reference gene during latex flow should be an efficient and effective clue for selection of reference genes in qRT-PCR.

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

Melon (Cucumis melo L.) is an attractive model plant for investigating fruit development because of its morphological, physiological, and biochemical diversity. Quantification of gene expression by quantitative reverse transcription polymerase chain reaction (qRT-PCR) with stably expressed reference genes for normalization can effectively elucidate the biological functions of genes that regulate fruit development. However, the reference genes for data normalization in melon fruits have not yet been systematically validated. This study aims to assess the suitability of 20 genes for their potential use as reference genes in melon fruits. Expression variations of these genes were measured in 24 samples that represented different developmental stages of fertilized and parthenocarpic melon fruits by qRT-PCR analysis. GeNorm identified ribosomal protein L (CmRPL) and cytosolic ribosomal protein S15 (CmRPS15) as the best pair of reference genes, and as many as five genes including CmRPL, CmRPS15, TIP41-like family protein (CmTIP41), cyclophilin ROC7 (CmCYP7), and ADP ribosylation factor 1 (CmADP) were required for more reliable normalization. NormFinder ranked CmRPS15 as the best single reference gene, and RAN GTPase gene family (CmRAN) and TATA-box binding protein (CmTBP2) as the best combination of reference genes in melon fruits. Their effectiveness was further validated by parallel analyses on the activities of soluble acid invertase and sucrose phosphate synthase, and expression profiles of their respective encoding genes CmAIN2 and CmSPS1, as well as sucrose contents during melon fruit ripening. The validated reference genes will help to improve the accuracy of gene expression studies in melon fruits.