Selection and validation of reference genes desirable for gene expression analysis by qRT-PCR in MeJA-treated ginseng hairy roots

Identification and validation of stable reference genes for RT-qPCR analyses of Kobresia littledalei seedlings

BACKGROUND

Kobreisa littledalei, belonging to the Cyperaceae family is the first Kobresia species with a reference genome and the most dominant species in Qinghai-Tibet Plateau alpine meadows. It has several resistance genes which could be used to breed improved crop varieties. Reverse Transcription Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) is a popular and accurate gene expression analysis method. Its reliability depends on the expression levels of reference genes, which vary by species, tissues and environments. However, K.littledalei lacks a stable and normalized reference gene for RT-qPCR analysis.

RESULTS

The stability of 13 potential reference genes was tested and the stable reference genes were selected for RT-qPCR normalization for the expression analysis in the different tissues of K. littledalei under two abiotic stresses (salt and drought) and two hormonal treatments (abscisic acid (ABA) and gibberellin (GA)). Five algorithms were used to assess the stability of putative reference genes. The results showed a variation amongst the methods, and the same reference genes showed tissue expression differences under the same conditions. The stability of combining two reference genes was better than a single one. The expression levels of ACTIN were stable in leaves and stems under normal conditions, in leaves under drought stress and in roots under ABA treatment. The expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression was stable in the roots under the control conditions and salt stress and in stems exposed to drought stress. Expression levels of superoxide dismutase (SOD) were stable in stems of ABA-treated plants and in the roots under drought stress. Moreover, RPL6 expression was stable in the leaves and stems under salt stress and in the stems of the GA-treated plants. EF1-alpha expression was stable in leaves under ABA and GA treatments. The expression levels of 28 S were stable in the roots under GA treatment. In general, ACTIN and GAPDH could be employed as housekeeping genes for K. littledalei under different treatments.

CONCLUSION

This study identified the best RT-qPCR reference genes for different K. littledalei tissues under five experimental conditions. ACTIN and GAPDH genes can be employed as the ideal housekeeping genes for expression analysis under different conditions. This is the first study to investigate the stable reference genes for normalized gene expression analysis of K. littledalei under different conditions. The results could aid molecular biology and gene function research on Kobresia and other related species.

Genome-wide identification and integrated analysis of TCP genes controlling ginsenoside biosynthesis in Panax ginseng

Panax ginseng is an important medicinal plant, and ginsenosides are the main bioactive molecules of ginseng. The TCP (TBI, CYC, PCF) family is a group of transcription factors (TFs) that play an important role in plant growth and development, hormone signalling and synthesis of secondary metabolites. In our study, 78 PgTCP transcripts were identified from the established ginseng transcriptome database. A phylogenetic tree analysis showed that the 67 PgTCP transcripts with complete open reading frames were classified into three subfamilies, including CIN, PCF, and CYC/TB1. Protein structure analysis showed that PgTCP genes had bHLH structures. Chromosomal localization analysis showed that 63 PgTCP genes were localized on 17 of the 24 chromosomes of the Chinese ginseng genome. Expression pattern analysis showed that PgTCP genes differed among different lineages and were spatiotemporally specific. Coexpression network analysis indicated that PgTCP genes were coexpressed and involved in plant activities or metabolic regulation in ginseng. The expression levels of PgTCP genes from class I (PCF) were significantly downregulated, while the expression levels of PgTCP genes from class II (CIN and CYC/TB1) were upregulated, suggesting that TCP genes may be involved in the regulation of secondary metabolism in ginseng. As the PgTCP26-02 gene was found to be related to ginsenoside synthesis, its predicted protein structure and expression pattern were further analysed. Our results provide new insights into the origin, differentiation, evolution and function of the PgTCP gene family in ginseng, as well as the regulation of plant secondary metabolism.

Screening of Reference Genes under Biotic Stress and Hormone Treatment of Mung Bean (Vigna radiata) by Quantitative Real-Time PCR

Mung bean (Vigna radiata) production has been greatly threatened by numerous diseases. Infection with these pathogens causes extensive changes in gene expression and the activation of hormone signal transduction. Quantitative real-time PCR (qRT-PCR) is the most common technique used for gene expression validation. Screening proper reference genes for mung bean under pathogen infection and hormone treatment is a prerequisite for ensuring the accuracy of qRT-PCR data in mung bean disease-resistance research. In this study, six candidate reference genes (Cons4, ACT, TUA, TUB, GAPDH, and EF1α) were selected to evaluate the expression stability under four soil-borne disease pathogens (Pythium myriotylum, Pythium aphanidermatum, Fusarium oxysporum, and Rhizoctonia solani) and five hormone treatments (SA, MeJA, ETH, ABA, and GA3). In the samples from different treatments, the Ct value distribution of the six candidate reference genes was different. Under the condition of hormone treatment, the Ct value ranged from a minimum of 17.87 for EF1α to a maximum of 29.63 for GAPDH. Under the condition of pathogen infection, the Ct value ranged from a minimum of 19.43 for EF1α to a maximum of 31.82 for GAPDH. After primer specificity analysis, it was found that GAPDH was not specific, so the five reference genes Cons4, ACT, TUA, TUB, and EF1α were used in subsequent experiments. The software products GeNorm, NormFinder, BestKeeper and RefFinder were used for qRT-PCR data analysis. In general, the best candidates reference genes were: TUA for SA, ABA, GA3, and Pythium myriotylum treatment; TUB for ETH treatment; ACT for MeJA and Fusarium oxysporum treatment; and EF1α for Pythium aphanidermatum and Rhizoctonia solani treatment. The most stably expressed genes in all samples were TUA, while Cons4 was the least stable reference gene. Finally, the reliability of the reference gene was further validated by analysis of the expression profiles of four mung bean genes (Vradi0146s00260, Vradi0158s00480, Vradi07g23860, and Vradi11g03350) selected from transcriptome data. Our results provide more accurate information for the normalization of qRT-PCR data in mung bean response to pathogen interaction.

Genetic and molecular dissection of ginseng (Panax ginseng Mey.) germplasm using high-density genic SNP markers, secondary metabolites, and gene expressions

Genetic and molecular knowledge of a species is crucial to its gene discovery and enhanced breeding. Here, we report the genetic and molecular dissection of ginseng, an important herb for healthy food and medicine. A mini-core collection consisting of 344 cultivars and landraces was developed for ginseng that represents the genetic variation of ginseng existing in its origin and diversity center. We sequenced the transcriptomes of all 344 cultivars and landraces; identified over 1.5 million genic SNPs, thereby revealing the genic diversity of ginseng; and analyzed them with 26,600 high-quality genic SNPs or a selection of them. Ginseng had a wide molecular diversity and was clustered into three subpopulations. Analysis of 16 ginsenosides, the major bioactive components for healthy food and medicine, showed that ginseng had a wide variation in the contents of all 16 ginsenosides and an extensive correlation of their contents, suggesting that they are synthesized through a single or multiple correlated pathways. Furthermore, we pair-wisely examined the relationships between the cultivars and landraces, revealing their relationships in gene expression, gene variation, and ginsenoside biosynthesis. These results provide new knowledge and new genetic and genic resources for advanced research and breeding of ginseng and related species.

Genome-wide identification and systematic analysis of the HD-Zip gene family and its roles in response to pH in Panax ginseng Meyer

BACKGROUND

Ginseng, Panax ginseng Meyer, is a traditional herb that is immensely valuable both for human health and medicine and for medicinal plant research. The homeodomain leucine zipper (HD-Zip) gene family is a plant-specific transcription factor gene family indispensable in the regulation of plant growth and development and plant response to environmental stresses.

RESULTS

We identified 117 HD-Zip transcripts from the transcriptome of ginseng cv. Damaya that is widely grown in Jilin, China where approximately 60% of the world's ginseng is produced. These transcripts were positioned to 64 loci in the ginseng genome and the ginseng HD-Zip genes were designated as PgHDZ genes. Identification of 82 and 83 PgHDZ genes from the ginseng acc. IR826 and cv. ChP genomes, respectively, indicated that the PgHDZ gene family consists of approximately 80 PgHDZ genes. Phylogenetic analysis showed that the gene family originated after Angiosperm split from Gymnosperm and before Dicots split from Monocots. The gene family was classified into four subfamilies and has dramatically diverged not only in gene structure and functionality but also in expression characteristics. Nevertheless, co-expression network analysis showed that the activities of the genes in the family remain significantly correlated, suggesting their functional correlation. Five hub PgHDZ genes were identified that might have central functions in ginseng biological processes and four of them were shown to be actively involved in plant response to environmental pH stress in ginseng.

CONCLUSIONS

The PgHDZ gene family was identified from ginseng and analyzed systematically. Five potential hub genes were identified and four of them were shown to be involved in ginseng response to environmental pH stress. The results provide new insights into the characteristics, diversity, evolution, and functionality of the PgHDZ gene family in ginseng and lay a foundation for comprehensive research of the gene family in plants.

Transcriptome and Phenotype Integrated Analysis Identifies Genes Controlling Ginsenoside Rb1 Biosynthesis and Reveals Their Interactions in the Process in Panax ginseng

Genes are the keys to deciphering the molecular mechanism underlying a biological trait and designing approaches desirable for plant genetic improvement. Ginseng is an important medicinal herb in which ginsenosides have been shown to be the major bioactive component; however, only a few genes involved in ginsenoside biosynthesis have been cloned through orthologue analysis. Here, we report the identification of 21 genes controlling Rb1 biosynthesis by stepwise ginseng transcriptome and Rb1 content integrated analysis. We first identified the candidate genes for Rb1 biosynthesis by integrated analysis of genes with the trait from four aspects, including gene transcript differential expression between highest- and lowest-Rb1 content cultivars, gene transcript expression-Rb1 content correlation, and biological impacts of gene mutations on Rb1 content, followed by the gene transcript co-expression network. Twenty-two candidate genes were identified, of which 21 were functionally validated for Rb1 biosynthesis by gene regulation, genetic transformation, and mutation analysis. These genes were strongly correlated in expression with the previously cloned genes encoding key enzymes for Rb1 biosynthesis. Based on the correlations, a pathway for Rb1 biosynthesis was deduced to indicate the roles of the genes in Rb1 biosynthesis. Moreover, the genes formed a strong co-expression network with the previously cloned Rb1 biosynthesis genes, and the variation in the network was associated with the variation in the Rb1 content. These results indicate that Rb1 biosynthesis is a process of correlative interactions among Rb1 biosynthesis genes. Therefore, this study provides new knowledge, 21 new genes, and 96 biomarkers for Rb1 biosynthesis useful for enhanced research and breeding in ginseng.

Evaluation and validation of reference genes in Cymbidium faberi for real-time quantitative PCR

For all organs at all Cymbidium faberi stages, ACT, UBQ3 and GAPDH can be selected as reference genes. For organs of the vegetative stage, UBQ2 and UBQ3 can be chosen for analysis of normalized gene expression. For the bud stage, ACT and UBQ3 can be used for analysis of gene expression. For the full blossom stage, ACT, UBQ3 and UBQ2 can be introduced into relative gene expression analysis. For vegetative organs, UBQ2 and ACT can be used as reference genes. For reproductive organs, ACT, UBQ3 and UBQ2 can be used as a reference for data processing. CfAG1 gene expression is more consistent when UBQ3, GAPDH and ACT are used as reference genes.

Genome-wide analysis of the C2H2 zinc finger protein gene family and its response to salt stress in ginseng, Panax ginseng Meyer

The C2H2 zinc finger protein (C2H2-ZFP) gene family plays important roles in response to environmental stresses and several other biological processes in plants. Ginseng is a precious medicinal herb cultivated in Asia and North America. However, little is known about the C2H2-ZFP gene family and its functions in ginseng. Here, we identified 115 C2H2-ZFP genes from ginseng, defined as the PgZFP gene family. It was clustered into five groups and featured with eight conserved motifs, with each gene containing one to six of them. The family genes are categorized into 17 gene ontology subcategories and have numerous regulatory elements responsive to a variety of biological process, suggesting their functional differentiation. The 115 PgZFP genes were spliced into 228 transcripts at seed setting stage and varied dramatically in expression across tissues, developmental stages, and genotypes, but they form a co-expression network, suggesting their functional correlation. Furthermore, four genes, PgZFP31, PgZFP78-01, PgZFP38, and PgZFP39-01, were identified from the gene family that were actively involved in plant response to salt stress. These results provide new knowledge on origin, differentiation, evolution, and function of the PgZFP gene family and new gene resources for C2H2-ZFP gene research and application in ginseng and other plant species.

Selection and validation of reference genes for quantitative real-time PCR of Quercus mongolica Fisch. ex Ledeb under abiotic stresses

Quercus mongolica Fisch. ex Ledeb is the main species of coniferous and broadleaved mixed forests in northeast and north China, which has high ornamental, economic, and ecological value. The appropriate reference genes must be selected for quantitative real-time PCR to reveal the molecular mechanisms of stress responses and their contribution to breeding of Q. mongolica. In the present study, we chose 11 candidate reference genes (TUA, CYP18, HIS4, RPS13, ACT97, TUB1, UBQ10, UBC5, SAND, PP2A, and SAMDC) and used four programs (GeNorm, NormFinder, BestKeeper, and RefFinder) to assess the expression stability of the above genes in roots, stems, and leaves under five abiotic stress factors (cold, salt, drought, weak light, and heavy metal). The findings revealed that under various experimental environments, the most stable genes were different; CYP18, ACT97, and RPS13 ranked the highest under most experimental environments. Moreover, two genes induced by stress, CMO and P5CS2, were chosen to demonstrate the reliability of the selected reference genes in various tissues under various stress conditions. Our research provides a significant basis for subsequent gene function studies of Q. mongolica.

Transcriptome-Wide Analysis for Ginsenoside Rb3 Synthesis-Related Genes and Study on the Expression of Methyl Jasmonate Treatment in Panax ginseng

Panax ginseng C. A. Meyer is a kind of renascent herb that belongs to the genus Panax in the family Araliaceae. It is a traditional Chinese precious herbal medicine with a long history of medicinal use. Ginsenoside Rb3 is one of the important active ingredients in ginseng and has important physiological activity in the treatment of many diseases. In this study, we screened and systematically analyzed the candidate genes related to ginsenoside Rb3 synthesis through bioinformatics methods; discussed the functions, expression patterns, and interactions of the genes related to ginsenoside Rb3 synthesis; and finally, selected seven genes, mainly PgRb3, that directly contribute to the synthesis of ginsenoside Rb3. This study provides a reference for revealing the expression rules of ginsenoside Rb3 synthesis-related genes and elucidating the regulatory mechanism of methyl jasmonate, lays a theoretical foundation for the research of ginsenoside Rb3 synthesis, and provides theoretical and technical support for the factory production of ginsenoside monomer saponins.

Identification of Specific Glycosyltransferases Involved in Flavonol Glucoside Biosynthesis in Ginseng Using Integrative Metabolite Profiles, DIA Proteomics, and Phylogenetic Analysis

Ginseng contains a variety of flavonol glycosides that possess diverse biological activities; however, scant information of flavonoid glycosylation was reported in ginseng. We found that panasenoside and kaempferol 3-O-glucoside were commonly accumulated along with cultivation years in leaves. In order to explore the procedure of flavonol glycosylation in ginseng, 50 UDP-glycosyltransferases (UGTs) were screened out using differentiated data-independent acquisition (DIA) proteomics and phylogenetic analysis. UGT92A10 and UGT94Q4 were found contributing to the formation of kaempferol 3-O-glucoside. UGT73A18, UGT74T4, and UGT75W1 could catalyze galactosylation of kaempferol 3-O-glucoside. Ser278, Trp335, Gln338, and Val339 were found forming hydrogen bonds with UDP-galactose in UGT75W1 by docking. MeJA induced transcripts of UGT73A18 and UGT74T4 by over fourfold, consistent with the decrease of kaempferol 3-O-glucoside, which indicated that these genes may be related to resisting adversity stress in ginseng. These results highlight the significance of integrative metabolite profiles, proteomics, and phylogenetic analysis for exploring flavonol glycosylation in ginseng.

Selection of suitable reference genes for qRT-PCR expression analysis of Codonopsis pilosula under different experimental conditions

Codonopsis pilosula is a well-known medicinal plant. Although its transcriptome sequence has been published, suitable reference genes have not been systematically identified for conducting expression analyses via quantitative real-time polymerase chain reaction (qRT-PCR). To screen appropriate genes for use with this species, we applied four different methods-GeNorm, NormFinder, BestKeeper, and RefFinder-to evaluate the stability of 13 candidates: CpiEF1Bb, CpiCACS, CpiF-Box, Cpiβ-Tubulin, CpiGAPDH, CpiActin2, CpiAPT1, CpiActin7, CpiActin8, CpiRPL6, CpiHAF1, CpiTubulin6, and CpiUBQ12. Expression was examined by qRT-PCR for various tissue types, chemical treatments, and developmental stages. For all tested samples, CpiGAPDH proved to be the most stable. Comprehensive analysis indicated that the most stable internal reference genes were CpiF-Box and CpiCACS in different tissues and at different developmental stages, respectively. Under NaCl stress, CpiAPT1 was the best internal reference gene. For methyl jasmonate and abscisic acid treatments, CpiGAPDH and CpiF-Box, respectively, presented the highest degree of expression stability. Based on these findings, we chose CpiSPL9 as the target gene for validating the suitability of these selected reference genes. All of these results provide a foundation for accurate quantification of expression levels by genes of interest in C. pilosula.

Transcriptome-Wide Identification, Evolutionary Analysis, and GA Stress Response of the GRAS Gene Family in Panax ginseng C. A. Meyer

GRAS transcription factors are a kind of plant-specific transcription factor that have been found in a variety of plants. According to previous studies, GRAS proteins are widely involved in the physiological processes of plant signal transduction, stress, growth and development. The Jilin ginseng (Panax ginseng C.A. Meyer) is a heterogeneous tetraploid perennial herb of the Araliaceae family, ginseng genus. Important information regarding the GRAS transcription factors has not been reported in ginseng. In this study, 59 Panax ginseng GRAS (PgGRAS) genes were obtained from the Jilin ginseng transcriptome data and divided into 13 sub-families according to the classification of Arabidopsis thaliana. Through systematic evolution, structural variation, function and gene expression analysis, we further reveal GRAS’s potential function in plant growth processes and its stress response. The expression of PgGRAS genes responding to gibberellin acids (GAs) suggests that these genes could be activated after application concentration of GA. The qPCR analysis result shows that four PgGRAS genes belonging to the DELLA sub-family potentially have important roles in the GA stress response of ginseng hairy roots. This study provides not only a preliminary exploration of the potential functions of the GRAS genes in ginseng, but also valuable data for further exploration of the candidate PgGRAS genes of GA signaling in Jilin ginseng, especially their roles in ginseng hairy root development and GA stress response.