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

Selection and Validation of Reference Genes in Virus-Infected Sweet Potato Plants

Quantitative real-time PCR (qRT-PCR) in sweet potatoes requires accurate data normalization; however, there are insufficient studies on appropriate reference genes for gene expression analysis. We examined variations in the expression of eight candidate reference genes in the leaf and root tissues of sweet potatoes (eight nonvirus-infected or eight virus-infected samples). Parallel analyses with geNorm, NormFinder, and Best-Keeper show that different viral infections and origin tissues affect the expression levels of these genes. Based on the results of the evaluation of the three software, the adenosine diphosphate-ribosylation factor is suitable for nonvirus or virus-infected sweet potato leaves. Cyclophilin and ubiquitin extension proteins are suitable for nonvirus-infected sweet potato leaves. Phospholipase D1 alpha is suitable for virus-infected sweet potato leaves. Actin is suitable for roots of nonvirus-infected sweet potatoes. Glyceraldehyde-3-phosphate dehydrogenase is suitable for virus-infected sweet potato roots. The research provides appropriate reference genes for further analysis in leaf and root samples of viruses in sweet potatoes.

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.

Genome-wide characterization of ascorbate peroxidase gene family in pepper (Capsicum annuum L.) in response to multiple abiotic stresses

Pepper is widely grown all over the world, so it faces many abiotic stresses, such as drought, high temperature, low temperature, salt damage, and so on. Stresses causing the accumulation of reactive oxidative species (ROS) in plants are removed by antioxidant defense systems, and ascorbate peroxidase (APX) is an important antioxidant enzyme. Therefore, the present study performed genome-wide identification of the APX gene family in pepper. We identified nine members of the APX gene family in the pepper genome according to the APX proteins' conserved domain in Arabidopsis thaliana. The physicochemical property analysis showed that CaAPX3 had the longest protein sequence and the largest molecular weight of all genes, while CaAPX9 had the shortest protein sequence and the smallest MW. The gene structure analysis showed that CaAPXs were composed of seven to 10 introns. The CaAPX genes were divided into four groups. The APX genes of groups I and IV were localized in the peroxisomes and chloroplasts, respectively; the group II genes were localized in the chloroplasts and mitochondria; and the group III genes were located in the cytoplasm and extracell. The conservative motif analysis showed that all APX genes in the pepper had motif 2, motif 3, and motif 5. The APX gene family members were distributed on five chromosomes (Chr. 2, 4, 6, 8, and 9). The cis-acting element analysis showed that most CaAPX genes contain a variety of cis-elements related to plant hormones and abiotic stress. RNA-seq expression analysis showed that the expression patterns of nine APXs were different in vegetative and reproductive organs at different growth and development stages. In addition, the qRT-PCR analysis of the CaAPX genes revealed significant differential expression in response to high temperature, low temperature, and salinity stresses in leaf tissue. In conclusion, our study identified the APX gene family members in the pepper and predicted the functions of this gene family, which would provide resources for further functional characterization of CaAPX genes.

Microscopic and metabolic investigations disclose the factors that lead to skin cracking in chili-type pepper fruit varieties

The hydrophobic cuticle encasing the fruit skin surface plays critical roles during fruit development and post-harvest. Skin failure often results in the fruit surface cracking and forming a wound-periderm tissue made of suberin and lignin. The factors that make the fruit skin susceptible to cracking have yet to be fully understood. Herein, we investigated two varieties of chili peppers (Capsicum annuum L.), Numex Garnet, whose fruit has intact skin, and Vezena Slatka, whose fruit has cracked skin. Microscopical observations, gas chromatography-mass spectrometry, biochemical and gene expression assays revealed that Vezena Slatka fruit form a thicker cuticle with greater levels of cutin monomers and hydroxycinnamic acids, and highly express key cutin-related genes. The skin of these fruit also had a lower epidermal cell density due to cells with very large perimeters, and highly express genes involved in epidermal cell differentiation. We demonstrate that skin cracking in the Vezena Slatka fruit is accompanied by a spatial accumulation of lignin-like polyphenolic compounds, without the formation of a typical wound-periderm tissues made of suberized cells. Lastly, we establish that skin cracking in chili-type pepper significantly affects fruit quality during post-harvest storage in a temperature-dependent manner. In conclusion, our data highlight cuticle thickness and epidermal cell density as two critical factors determining fruit skin susceptibility to cracking in chili-type pepper fruit.

Biocontrol of Diseases Caused by Phytophthora capsici and P. parasitica in Pepper Plants

The main objective of this study was to evaluate the ability of Trichoderma aggressivum f. europaeum, T. longibrachiatum, Paecilomyces variotii, and T. saturnisporum as biological control agents (BCAs) against diseases caused by P. capsici and P. parasitica in pepper. For this purpose, their antagonistic activities were evaluated both in vitro and in vivo. We analysed the expression patterns of five defence related genes, CaBGLU, CaRGA1, CaBPR1, CaPTI1, and CaSAR8.2, in leaves. All BCAs showed a high in vitro antagonistic activity, significantly reducing the mycelial growth of P. capsici and P. parasitica. The treatments with T. aggressivum f. europaeum, T. longibrachiatum, and P. variotii substantially reduced the severity of the disease caused by P. capsici by 54, 76, and 70%, respectively, and of the disease caused by P. parasitica by 66, 55, and 64%, respectively. T. saturnisporum had the lowest values of disease reduction. Reinoculation with the four BCAs increased the control of both plant pathogens. Markedly different expression patterns were observed in the genes CaBGLU, CaRGA1, and CaSAR8.2. Based on the results, all four BCAs under study could be used as a biological alternative to chemicals for the control of P. capsici and P. parasitica in pepper with a high success rate.

Red LED Light Improves Pepper (Capsicum annuum L.) Seed Radicle Emergence and Growth through the Modulation of Aquaporins, Hormone Homeostasis, and Metabolite Remobilization

Red LED light (R LED) is an efficient tool to improve seed germination and plant growth under controlled environments since it is more readily absorbed by photoreceptors' phytochromes compared to other wavelengths of the spectrum. In this work, the effect of R LED on the radicle emergence and growth (Phase III of germination) of pepper seeds was evaluated. Thus, the impact of R LED on water transport through different intrinsic membrane proteins, via aquaporin (AQP) isoforms, was determined. In addition, the remobilization of distinct metabolites such as amino acids, sugars, organic acids, and hormones was analysed. R LED induced a higher germination speed index, regulated by an increased water uptake. PIP2;3 and PIP2;5 aquaporin isoforms were highly expressed and could contribute to a faster and more effective hydration of embryo tissues, leading to a reduction of the germination time. By contrast, TIP1;7, TIP1;8, TIP3;1 and TIP3;2 gene expressions were reduced in R LED-treated seeds, pointing to a lower need for protein remobilization. NIP4;5 and XIP1;1 were also involved in radicle growth but their role needs to be elucidated. In addition, R LED induced changes in amino acids and organic acids as well as sugars. Therefore, an advanced metabolome oriented to a higher energetic metabolism was observed, conditioning better seed germination performance together with a rapid water flux.

Evaluation of reference genes for transcript normalization in Fragaria chiloensis fruit and vegetative tissues

Quantitative real-time PCR (RT-qPCR) is used extensively in gene expression studies. For adequate comparisons, the identification and use of reliable reference genes are crucial. Nevertheless, the availability of such genes in strawberry species is limited and has yet to be described for the Chilean strawberry, Fragaria chiloensis. In this study, the expression dynamics of a set of 10 candidate reference genes were analyzed in various F. chiloensis vegetative tissues (root, runners, stem, leaf, and flower), and fruits at different ripening stages or subjected to different hormonal treatments (ABA, auxin). The expression stability of candidate genes was examined by a series of algorithms, such as geNorm, NormFinder, BestKeeper, and ΔCt, for comparisons and rankings. Finally, by using RefFinder, a comprehensive and comparative ranking of the four methods was achieved. The results highlight that the expression stability of candidate reference genes fluctuates depending on tissue type, fruit stage, and hormonal treatment. As reference genes, the use of FcCHP2 and FcACTIN1 is recommended for F. chiloensis vegetative tissues; FcDBP and FcCHP1 for fruit ripening stages; FcGAPDH and FcDBP for fruit subjected to ABA and NDGA treatments; FcCHP1 and FcCHP2 for fruit under AUXIN and TIBA treatments; and FcDBP and FcCHP2 when all fruit stages and hormonal treatments are compared. If just one reference gene is employed as a normalizer, FcDBP should be chosen as it is the most stable internal control in most conditions. Therefore, the present study delivers a set of reliable reference genes for RT-qPCR expression analysis in F. chiloensis tissues and fruits subjected to several hormonal treatments.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01227-y.

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

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

Gene expression analysis for selection and validation of suitable housekeeping gene(s) in cadmium exposed pigeonpea plants inoculated with arbuscular mycorrhizae

The expression stability of six commonly used housekeeping genes (18S rRNA-18S ribosomal RNA, EF1α-elongation factor 1α, ACT1-Actin 1, GAPDH-Glyceraldehyde-3-phosphate dehydrogenase, TUB6-Tubulin/FtsZ family and UBC-Ubiquitin-conjugating enzyme) were scrutinized in leaves and roots of Cd stressed pigeonpea plants inoculated with arbuscular mycorrhizal (AM) species- Rhizoglomus intraradices (Ri), Funneliformis mosseae (Fm), Claroideoglomus etunicatum (Ce), C. claroideum (Cc). The stability profile of each gene was assessed using ΔC_t, BestKeeper, NormFinder, RefFinder and geNorm algorithmic programs, which ranked different genes as most and least stable according to the tissues analysed. All the statistical algorithms ranked TUB6 as most stable and EF1α least stable housekeeping (HK) genes in both the plant tissues. The selected HK genes were verified using metallothionein (CcMT1) i.e. a stress responsive gene, whose expression altered under conditions of metal stress and AM inoculation. The expression pattern of CcMT1 varied highly when least stable reference gene was used for normalization as compared to most stable gene, under different treatments. Thus, there is a need of selecting suitable reference gene to achieve reliable results in gene expression studies using quantitative real time PCR (qRT-PCR). The study conducted will help future gene expression analysis in pigeonpea under specific stress.

Temporal expression patterns of fruit-specific α- EXPANSINS during cell expansion in bell pepper (Capsicum annuum L.)

BACKGROUND

Expansins (EXPs) facilitate non-enzymatic cell wall loosening during several phases of plant growth and development including fruit growth, internode expansion, pollen tube growth, leaf and root development, and during abiotic stress responses. In this study, the spatial and temporal expression patterns of C. annuum α- EXPANSIN (CaEXPA) genes were characterized. Additionally, fruit-specific CaEXPA expression was correlated with the rate of cell expansion during bell pepper fruit development.

RESULTS

Spatial expression patterns revealed that CaEXPA13 was up-regulated in vegetative tissues and flowers, with the most abundant expression in mature leaves. Expression of CaEXPA4 was associated with stems and roots. CaEXPA3 was expressed abundantly in flower at anthesis suggesting a role for CaEXPA3 in flower development. Temporal expression analysis revealed that 9 out of the 21 genes were highly expressed during fruit development. Of these, expression of six genes, CaEXPA5, CaEXPA7, CaEXPA12, CaEXPA14 CaEXPA17 and CaEXPA19 were abundant 7 to 21 days after anthesis (DAA), whereas CaEXPA6 was strongly expressed between 14 and 28 DAA. Further, this study revealed that fruit growth and cell expansion occur throughout bell pepper development until ripening, with highest rates of fruit growth and cell expansion occurring between 7 and 14 DAA. The expression of CaEXPA14 and CaEXPA19 positively correlated with the rate of cell expansion, suggesting their role in post-mitotic cell expansion-mediated growth of the bell pepper fruit. In this study, a ripening specific EXP transcript, CaEXPA9 was identified, suggesting its role in cell wall disassembly during ripening.

CONCLUSIONS

This is the first genome-wide study of CaEXPA expression during fruit growth and development. Identification of fruit-specific EXPAs suggest their importance in facilitating cell expansion during growth and cell wall loosening during ripening in bell pepper. These EXPA genes could be important targets for future manipulation of fruit size and ripening characteristics.

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.

Genome-Wide Identification and Characterization of the Mitochondrial Transcription Termination Factors (mTERFs) in Capsicum annuum L

Mitochondrial transcription termination factors (mTERFs) regulate the expression of mitochondrial genes and are closely related to the function of the mitochondrion and chloroplast. In this study, the mTERF gene family in capsicum (Capsicum annuum L.) was identified and characterized through genomic and bioinformatic analyses. Capsicum was found to possess at least 35 mTERF genes (CamTERFs), which were divided into eight major groups following phylogenetic analysis. Analysis of CamTERF promoters revealed the presence of many cis-elements related to the regulation of cellular respiration and photosynthesis. In addition, CamTERF promoters contained cis-elements related to phytohormone regulation and stress responses. Differentially expressed genes in different tissues and developmental phases were identified using RNA-seq data, which revealed that CamTERFs exhibit various expression and co-expression patterns. Gene ontology (GO) annotations associated CamTERFs primarily with mitochondrion and chloroplast function and composition. These results contribute towards understanding the role of mTERFs in capsicum growth, development, and stress responses. Moreover, our data assist in the identification of CamTERFs with important functions, which opens avenues for future studies.

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.

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.

Candidate Reference Genes Selection and Application for RT-qPCR Analysis in Kenaf with Cytoplasmic Male Sterility Background

Cytoplasmic male sterility (CMS) is a maternally inherited trait that results in the production of dysfunctional pollen. Based on reliable reference gene-normalized real-time quantitative PCR (RT-qPCR) data, examining gene expression profile can provide valuable information on the molecular mechanism of kenaf CMS. However, studies have not been conducted regarding selection of reference genes for normalizing RT-qPCR data in the CMS and maintainer lines of kenaf crop. Therefore, we studied 10 candidate reference genes (ACT3, ELF1A, G6PD, PEPKR1, TUB, TUA, CYP, GAPDH, H3, and 18S) to assess their expression stability at three stages of pollen development in CMS line 722A and maintainer line 722B of kenaf. Five computational statistical approaches (GeNorm, NormFinder, ΔCt, BestKeeper, and RefFinder) were used to evaluate the expression stability levels of these genes. According to RefFinder and GeNorm, the combination of TUB, CYP, and PEPKR1 was identified as an internal control for the accurate normalization across all sample set, which was further confirmed by validating the expression of HcPDIL5-2a. Furthermore, the combination of TUB, CYP, and PEPKR1 was used to differentiate the expression pattern of five mitochondria F₁F₀-ATPase subunit genes (atp1, atp4, atp6, atp8, and atp9) by RT-qPCR during pollen development in CMS line 722A and maintainer line 722B. We found that atp1, atp6, and atp9 exhibited significantly different expression patterns during pollen development in line 722A compared with line 722B. This is the first systematic study of reference genes selection for CMS and will provide useful information for future research on the gene expressions and molecular mechanisms underlying CMS in kenaf.