Selection and Validation of Reference Genes for Gene Expression Studies in Codonopsis pilosula Based on Transcriptome Sequence Data

Transcriptome-based analysis of key functional genes in the triterpenoid saponin synthesis pathway of Platycodon grandiflorum

BACKGROUND

Platycodon grandiflorum (P. grandiflorum) is a commonly used medicinal plant in China. Transcriptome sequencing studies of different tissues of P. grandiflorum have been widely conducted. However, studies on transcriptome sequencing and expression patterns of key genes in the saponin synthesis pathway of Tongcheng P. grandiflorum, a high-quality medicinal resource of different years, are relatively limited.

RESULTS

This study involved transcriptome sequencing and bioinformatics analysis of the roots from annual, biennial, and triennial P. grandiflorum in the Tongcheng area. After data filtering and assembly, we obtained 111.44 Gb of clean base data, including 742,880616 clean reads. We identified 5,156 differential expression unigenes between at least two sample groups, with differences noted among annual, biennial, and triennial P. grandiflorum plants. GO enrichment analysis annotated 3509, 1819, and 1393 DEGs in comparison TC1vsTC2, TC1vsTC3, and TC2vsTC3, respectively. Furthermore, KEGG enrichment analysis identified 16 genes encoding key enzymes in the terpene skeleton biosynthesis, sesquiterpene and triterpene biosynthesis pathways, including SE, AACT, FPPS, DXR, HMGR, HMGS, and DXS. The results of qRT-PCR experiments showed that most of the genes were most highly expressed in annual P. grandiflorum.

CONCLUSIONS

The present study provided transcriptomic data from the roots of Tongcheng P. grandiflorum of different years, which provides critical bioinformatics data on the growth and development of P. grandiflorum, laying a foundation for further research on saponins and identifying key enzymes involved in this process across different growth stages.

A protocol for identifying universal reference genes within a genus based on RNA-Seq data: a case study of poplar stem gene expression

Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) plays a crucial role in relative gene expression analysis, and accurate normalization relies on suitable reference genes (RGs). In this study, a pipeline for identifying candidate RGs from publicly available stem-related RNA-Seq data of different Populus species under various developmental and abiotic stress conditions is presented. DESeq2's median of ratios yielded the smallest coefficient of variance (CV) values in a total of 292 RNA-Seq samples and was therefore chosen as the method for sample normalization. A total of 541 stably expressed genes were retrieved based on the CV values with a cutoff of 0.3. Universal gene-specific primer pairs were designed based on the consensus sequences of the orthologous genes of each Populus RG candidate. The expression levels of 12 candidate RGs and six reported RGs in stems under different abiotic stress conditions or in different Populus species were assessed by RT-qPCR. The expression stability of selected genes was further evaluated using ΔCt, geNorm, NormFinder, and BestKeeper. All candidate RGs were stably expressed in different experiments and conditions in Populus. A test dataset containing 117 RNA-Seq samples was then used to confirm the expression stability, six candidate RGs and three reported RGs met the requirement of CV ≤ 0.3. In summary, this study was to propose a systematic and optimized protocol for the identification of constitutively and stably expressed genes based on RNA-Seq data, and Potri.001G349400 (CNOT2) was identified as the best candidate RG suitable for gene expression studies in poplar stems.

Transcriptome-based identification and validation of reference genes for corm growth stages, different tissues, and drought stress in Taro (Colocasia esculenta)

Taro is a widely utilized starch resource plant. It is essential to quantify the expression levels of functional genes associated with taro growth using real-time quantitative polymerase chain reaction (RT-qPCR). However, to obtain reliable RT-qPCR results, appropriate reference genes (RGs) are required for data normalization. In this study, we screened seven novel candidate RGs using transcriptome datasets from taro, encompassing data from growth corms and various tissues. The expression stability of these seven new RGs, along with the commonly used RGs Actin, EF1-α, and β-tubulin, was assessed using Delta Ct, BestKeeper, geNorm, and NormFinder algorithms. Furthermore, we conducted a comprehensive analysis using the RefFinder program and validated the results using the target gene, CeAGPL1. The findings revealed that ACY-1 and PIA2 were the optimal multiple RGs for normalization during corm growth, while COX10 and Armc8 were suitable for samples including various types of tissues. Furthermore, we found three RGs, Armc8, COX10 and CCX4L, were the optimal RGs for drought stress. This study assessed the suitability of RGs in taro for the first time. The identified RGs provide valuable resources for studying corm growth, diverse tissues, and drought stress. This study contributes to the advancement of our understanding of the underlying mechanisms that govern the growth of taro.

Selection and Validation of Reference Genes for Gene Expression Studies in Euonymus japonicus Based on RNA Sequencing

Euonymus japonicus is one of the most low-temperature-tolerant evergreen broad-leaved tree species in the world and is widely used in urban greening. However, there are very few molecular biology studies on its low-temperature tolerance mechanism. So far, no researcher has selected and reported on its reference genes. In this study, 21 candidate reference genes (12 traditional housekeeping genes and 9 other genes) were initially selected based on gene expression and coefficient of variation (CV) through RNA-Seq (unpublished data), and qRT-PCR was used to detect the expression levels of candidate reference genes in three different groups of samples (leaves under different temperature stresses, leaves of plants at different growth stages, and different organs). After further evaluating the expression stability of these genes using geNorm, NormFinder, Bestkeeper, and RefFind, the results show that the traditional housekeeping gene eIF5A and the new reference gene RTNLB1 have good stability in the three different groups of samples, so they are reference genes with universality. In addition, we used eIF5A and RTNLB1 as reference genes to calibrate the expression pattern of the target gene EjMAH1, which confirmed this view. This article is the first to select and report on the reference gene of E. japonicus, laying the foundation for its low-temperature tolerance mechanism and other molecular biology research.

Insights into Genus Codonopsis: From past Achievements to Future Perspectives

Codonopsis plants, as a kind of medicinal and edible herb, have a long history of application and have been widely concerned by pharmacists and biologists. In this article, the species diversity, taxonomy and distribution, ethnic medicinal records, chemical composition, pharmacological activity, and quality evaluation methods of Codonopsis species were systematically reviewed. In addition, the research progress of Codonopsis plants using biotechnology in recent years was summarized. The phytochemistry and biological activities of Codonopsis are widely evaluated. To date, more than 350 compounds have been isolated from Codonopsis. Codonopsis pilosula polysaccharides are important functional components and biomarkers. Lobetyolin, atractylenolide III, tangshenoside I, and oligosaccharide can be considered as characteristic index components to evaluate the quality of Codonopsis plants. Although recent experimental evidence has confirmed the pharmacological value of this genus, its quality control, resource development and utilization, and active ingredient synthesis mechanisms are not well studied. In particular, molecular biology research is still in its infancy, but its application prospects are broad, and it is a hot spot for future research on Codonopsis. Therefore, it is urgent to conduct a detailed study on the single level of phytochemistry, pharmacology, and molecular biology of Codonopsis to establish a scientific evaluation system and modern medication guidelines. The multi-angle, multi-level, and multi-aspect integrated association analysis is also an inevitable trend for the future in-depth study of Codonopsis plants. This research status was summarized in order to provide a broader scientific research idea and theoretical reference for the in-depth study of Codonopsis.

Screening and Validation of Internal Reference Genes for Quantitative Real-Time PCR Analysis of Leaf Color Mutants in Dendrobium officinale

Leaf color mutants (LCMs) are important resources for studying diverse metabolic processes such as chloroplast biogenesis and differentiation, pigments' biosynthesis and accumulation, and photosynthesis. However, in Dendrobium officinale, LCMs are yet to be fully studied and exploited due to the unavailability of reliable RGs (reference genes) for qRT-PCR (quantitative real-time reverse transcription PCR) normalization. Hence, this study took advantage of previously released transcriptome data to select and evaluate the suitability of ten candidate RGs, including Actin (Actin), polyubiquitin (UBQ), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), elongation factor 1-α (EF1α), β-tubulin (β-TUB), α-tubulin (α-TUB), 60S ribosomal protein L13-1 (RPL13AD), aquaporin PIP1-2 (PIP1-2), Intima protein (ALB3) and Cyclin (CYCB1-2) for normalizing leaf color-related genes' expression levels via qRT-PCR. Stability rankings analysis via common software Best-Keeper, GeNorm, and NormFinder disclosed that all ten genes met the requirements of RGs. Of them, EF1α exhibited the highest stability and was selected as the most reliable. The reliability and accuracy of EF1α were confirmed through qRT-PCR analysis of fifteen chlorophyll pathway-related genes. The expression patterns of these genes via EF1α normalization were consistent with the results by RNA-Seq. Our results offer key genetic resources for the functional characterization of leaf color-related genes and will pave the way for molecular dissection of leaf color mutations in D. officinale.

Identification and validation of reference genes in vetiver (Chrysopogon zizanioides) root transcriptome

Vetiver [Vetiveria zizanioides (L.) Roberty] is a perennial C-4 grass traditionally valued for its aromatic roots/root essential oil. Owing to its deep penetrating web-forming roots, the grass is now widely used across the globe for phytoremediation and the conservation of soil and water. This study has used the transcriptome data of vetiver roots in its two distinct geographic morphotypes (North Indian type A and South Indian type B) for reference gene(s) identification. Further, validation of reference genes using various abiotic stresses such as heat, cold, salt, and drought was carried out. The de novo assembly based on differential genes analysis gave 1,36,824 genes (PRJNA292937). Statistical tests like RefFinder, NormFinder, BestKeeper, geNorm, and Delta-Ct software were applied on 346 selected contigs. Eleven selected genes viz., GAPs, UBE2W, RP, OSCam2, MUB, RPS, Core histone 1, Core histone 2, SAMS, GRCWSP, PLDCP along with Actin were used for qRT-PCR analysis. Finally, the study identified the five best reference genes GAPs, OsCam2, MUB, Core histone 1, and SAMS along with Actin. The two optimal reference genes SAMS and Core histone 1 were identified with the help of qbase + software. The findings of the present analyses have value in the identification of suitable reference gene(s) in transcriptomic and molecular data analysis concerning various phenotypes related to abiotic stress and developmental aspects, as well as a quality control measure in gene expression experiments. Identifying reference genes in vetiver appears important as it allows for accurate normalization of gene expression data in qRT-PCR experiments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01315-7.

Identification and validation of reference genes for qRT-PCR based studies in horse gram (Macrotyloma uniflorum)

The quantitative real-time polymerase chain reaction (qRT-PCR) is the most sensitive and commonly used technique for gene expression studies in biological systems. However, the reliability of qRT-PCR results depends on the selection of reference gene(s) for data normalization. Horse gram (Macrotyloma uniflorum) is an important legume crop on which several molecular studies have been reported. However, the stability of reference genes has not been evaluated. In the present study, nine candidate reference genes were identified from horse gram RNA-seq data and evaluated in two horse gram genotypes, HPK4 and HPKM317 under six abiotic stresses viz. cold, drought, salinity, heat, abscisic acid and methyl viologen-induced oxidative stress. The results were evaluated using geNorm, Bestkeeper, Normfinder and delta-delta Ct methods and comprehensive ranking was assigned using RefFinder and RankAggreg software. The overall result showed that TCTP was one of the most stable genes in all samples and in genotype HPK4, while in HPKM317 profilin was most stably expressed. However, PSMA5 was identified as least stable in all the experimental conditions. Expression of target genes dehydrin and early response to dehydration 6 under drought stress was also validated using TCTP and profilin for data normalization, either alone or in combination, which confirmed their suitability for qRT-PCR data normalization. Thus, TCTP and profilin genes may be used for qRT-PCR data normalization for molecular and genomic studies in horse gram.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01104-0.

The transcription factor Xrp1 is required for PERK-mediated antioxidant gene induction in Drosophila

PERK is an endoplasmic reticulum (ER) transmembrane sensor that phosphorylates eIF2α to initiate the Unfolded Protein Response (UPR). eIF2α phosphorylation promotes stress-responsive gene expression most notably through the transcription factor ATF4 that contains a regulatory 5’ leader. Possible PERK effectors other than ATF4 remain poorly understood. Here, we report that the bZIP transcription factor Xrp1 is required for ATF4-independent PERK signaling. Cell-type-specific gene expression profiling in Drosophila indicated that delta-family glutathione-S-transferases (gstD) are prominently induced by the UPR-activating transgene Rh1^G69D. Perk was necessary and sufficient for such gstD induction, but ATF4 was not required. Instead, Perk and other regulators of eIF2α phosphorylation regulated Xrp1 protein levels to induce gstDs. The Xrp1 5’ leader has a conserved upstream Open Reading Frame (uORF) analogous to those that regulate ATF4 translation. The gstD-GFP reporter induction required putative Xrp1 binding sites. These results indicate that antioxidant genes are highly induced by a previously unrecognized UPR signaling axis consisting of PERK and Xrp1.

Selection and validation of reference genes for RT-qPCR analysis in Desmodium styracifolium Merr

Gene expression valuated by reverse transcription-quantitative PCR (RT-qPCR) are often applied to study the gene function. To obtain accurate and reliable results, the usage of stable reference genes is essential for RT-qPCR analysis. The traditional southern Chinese medicinal herb, Desmodium styracifolium Merr is well known for its remarkable effect on the treatment of urination disturbance, urolithiasis, edema and jaundice. However, there are no ready-made reference genes identified for D. styracifolium. In this study, 13 novel genes retrieved from transcriptome datasets of four different tissues were reported according to the coefficient of variation (CV) and maximum fold change (MFC) of gene expression. The expression stability of currently used Leguminosae ACT6 was compared to the 13 candidate reference genes in different tissues and 7-day-old seedlings under different experimental conditions, which was evaluated by five statistical algorithms (geNorm/NormFinder/BestKeeper/ΔCT/RefFinder). Our results indicated that the reference gene combinations of PP  +  UFM1, CCRP4  +  BRM and NFD6  +  NCLN1 were the most stable reference genes in leaf, stem and root tissues, respectively. The most stable reference gene combination for all tissues was CCRP4  +  CUL1. In addition, the most stable reference genes for different experimental conditions were distinct, for instance SMUP1 for MeJA treatment, ERDJ2A  +  SMUP1 for SA treatment, NCLN1  +  ERDJ2A for ABA treatment and SF3B  +  VAMP721d for salt stress, respectively. Our results lay a foundation for achieving accurate and reliable RT-qPCR results so as to correctly understand the function of genes in D. styracifolium.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02954-x.

Determination of reliable reference genes for gene expression studies in Chinese chive (Allium tuberosum) based on the transcriptome profiling

Chinese chive (Allium tuberosum) is widely cultivated around the world for its unique flavor, nutrient, and medicinal values, yet its molecular mechanism on flavor formation and other metabolic pathways remains intangible. The elucidation of these complex processes begins with investigating the expression of the genes of interest, however the appropriate reference genes (RGs) for normalizing the gene expression are still unavailable in A. tuberosum. To fill this lacuna, transcriptome-wide screening was undertaken to identify the most stable genes according to the analysis of their FPKM values. The expression stability of the RGs was further evaluated using geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The comprehensive analysis showed that GLY1 and SKP1, instead of two traditionally used RGs (eIF1α and ACT2), were the most stable genes across diverse A. tuberosum tissues, indicating the necessity to carefully validate the stability of RGs prior to their use for normalizations. As indicated by geNorm, the normalizations with at least two RGs could give more accurate results. qRT-PCR experiments were conducted with randomly selected genes, demonstrating that normalization with a combination of GLY1 and SKP1 resulted in reliable normalization results. Our finding represents the first attempt toward establishing a standardized qRT-PCR analysis in this economically important vegetable.

Transcriptome-based selection and validation of optimal house-keeping genes for skin research in goats (Capra hircus)

BACKGROUND

In quantitative real-time polymerase chain reaction (qRT-PCR) experiments, accurate and reliable target gene expression results are dependent on optimal amplification of house-keeping genes (HKGs). RNA-seq technology offers a novel approach to detect new HKGs with improved stability. Goat (Capra hircus) is an economically important livestock species and plays an indispensable role in the world animal fiber and meat industry. Unfortunately, uniform and reliable HKGs for skin research have not been identified in goat. Therefore, this study seeks to identify a set of stable HKGs for the skin tissue of C. hircus using high-throughput sequencing technology.

RESULTS

Based on the transcriptome dataset of 39 goat skin tissue samples, 8 genes (SRP68, NCBP3, RRAGA, EIF4H, CTBP2, PTPRA, CNBP, and EEF2) with relatively stable expression levels were identified and selected as new candidate HKGs. Commonly used HKGs including SDHA and YWHAZ from a previous study, and 2 conventional genes (ACTB and GAPDH) were also examined. Four different experimental variables: (1) different development stages, (2) hair follicle cycle stages, (3) breeds, and (4) sampling sites were used for determination and validation. Four algorithms (geNorm, NormFinder, BestKeeper, and ΔCt method) and a comprehensive algorithm (ComprFinder, developed in-house) were used to assess the stability of each HKG. It was shown that NCBP3 + SDHA + PTPRA were more stably expressed than previously used genes in all conditions analysis, and that this combination was effective at normalizing target gene expression. Moreover, a new algorithm for comprehensive analysis, ComprFinder, was developed and released.

CONCLUSION

This study presents the first list of candidate HKGs for C. hircus skin tissues based on an RNA-seq dataset. We propose that the NCBP3 + SDHA + PTPRA combination could be regarded as a triplet set of HKGs in skin molecular biology experiments in C. hircus and other closely related species. In addition, we also encourage researchers who perform candidate HKG evaluations and who require comprehensive analysis to adopt our new algorithm, ComprFinder.