Selection and evaluation of reference genes for qRT-PCR of Scutellaria baicalensis Georgi under different experimental conditions

Metabolomic and transcriptomic reveal flavonoid biosynthesis and regulation mechanism in Phlomoides rotata from different habitats

Phlomoides rotata is a traditional medical plant at 3100-5200 m altitude in the Tibet Plateau. In this study, flavonoid metabolites were investigated in P. rotata from Henan County (HN), Guoluo County (GL), Yushu County (YS), and Chengduo County (CD) habitats in Qinghai. The level of kaempferol 3-neohesperidoside, sakuranetin, and biochanin A was high in HN. The content of limocitrin and isoquercetin was high in YS. The levels of ikarisoside A and chrysosplenol D in GL were high. Schaftoside, miquelianin, malvidin chloride, and glabrene in CD exhibited high levels. The results showed a significant correlation between 59 flavonoids and 29 DEGs. Eleven flavonoids increased with altitude. PAL2, UFGT6, COMT1, HCT2, 4CL4, and HCT3 genes were crucial in regulating flavonoid biosynthesis. Three enzymes CHS, 4CL, and UFGT, were crucial in regulating flavonoid biosynthesis. This study provided biological and chemical evidence for the different uses of various regional plants of P. rotata.

Polyethylene Glycol Priming Enhances the Seed Germination and Seedling Growth of Scutellaria baicalensis Georgi under Salt Stress

Seed priming has become a practical pre-sowing strategy to deal with abiotic stresses. This study aims to explore the effects of polyethylene glycol (PEG) priming on seed germination and seedling growth of Scutellaria baicalensis Georgi under salt stress. Regardless of seed priming, salt stress significantly inhibited the seed germination and seedling growth of S. baicalensis. PEG priming significantly alleviates the inhibitory effects of salt stress on seed germination and seedling growth when compared to non-priming and water priming. Among all treatments, PEG priming exhibited the highest germination rate, germination potential, seed vigor index, fresh weight, dry weight, and plant length; the highest contents of proline, soluble sugar, and soluble protein; the highest K+/Na+ ratio and relative water content; the highest antioxidant activities and contents; but the lowest H2O2, malondialdehyde (MDA) content, and relative electrical conductivity in response to salt stress. In addition, PEG priming had the highest transcript levels of antioxidant-related genes among all treatments under NaCl stress. Taken together, the results demonstrated that seed priming with PEG could be recommended as an effective practice to enhance the germination and early seedling growth of S. baicalensis under saline conditions.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Systematic identification of reference genes for qRT-PCR of Ardisia kteniophylla A. DC under different experimental conditions and for anthocyanin-related genes studies

Ardisia kteniophylla A. DC, widely known as folk medicinal herb and ornamental plant, has been extensively investigated due to its unique leaf color, anti-cancer and other pharmacological activities. The quantitative real-time PCR (qRT-PCR) was an excellent tool for the analysis of gene expression with its high sensitivity and quantitative properties. Normalizing gene expression with stable reference genes was essential for qRT-PCR accuracy. In addition, no studies have yet been performed on the selection, verification and stability of internal reference genes suitable for A. kteniophylla, which has greatly hindered the biomolecular researches of this species. In this study, 29 candidate genes were successfully screened according to stable expression patterns of large-scale RNA seq data that from a variety of tissues and the roots of different growth stages in A. kteniophylla. The candidates were then further determined via qRT-PCR in various experimental samples, including MeJA, ABA, SA, NaCl, CuSO4, AgNO3, MnSO4, CoCl2, drought, low temperature, heat, waterlogging, wounding and oxidative stress. To assess the stability of the candidates, five commonly used strategies were employed: delta-CT, geNorm, BestKeeper, NormFinder, and the comprehensive tool RefFinder. In summary, UBC2 and UBA1 were found to be effective in accurately normalizing target gene expression in A. kteniophella regardless of experimental conditions, while PP2A-2 had the lowest stability. Additionally, to verify the reliability of the recommended reference genes under different colored leaf samples, we examined the expression patterns of six genes associated with anthocyanin synthesis and regulation. Our findings suggested that PAP1 and ANS3 may be involved in leaf color change in A. kteniphella. This study successfully identified the ideal reference gene for qRT-PCR analysis in A. kteniphella, providing a foundation for future research on gene function, particularly in the biosynthesis of anthocyanins.

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.

Evaluation of Candidate Reference Genes for Gene Expression Analysis in Wild Lamiophlomis rotata

Lamiophlomis rotata (Benth.) Kudo is a perennial and unique medicinal plant of the Qinghai–Tibet Plateau. It has the effects of diminishing inflammation, activating blood circulation, removing blood stasis, reducing swelling, and relieving pain. However, thus far, reliable reference gene identifications have not been reported in wild L. rotata. In this study, we identified suitable reference genes for the analysis of gene expression related to the medicinal compound synthesis in wild L. rotata subjected to five different-altitude habitats. Based on the RNA-Seq data of wild L. rotata from five different regions, the stability of 15 candidate internal reference genes was analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. TFIIS, EF-1α, and CYP22 were the most suitable internal reference genes in the leaves of L. rotata from different regions, while OBP, TFIIS, and CYP22 were the optimal reference genes in the roots of L. rotata. The reference genes identified here would be very useful for gene expression studies with different tissues in L. rotata from different habitats.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

Selection and validation of reference genes for RT-qPCR-based analyses of Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae)

RT-qPCR remains a vital approach for molecular biology studies aimed at quantifying gene expression in a range of physiological or pathological settings. However, the use of appropriate reference genes is essential to attain meaningful RT-qPCR results. Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae) is an important egg parasitoid wasp and natural enemy of fruit bugs and forest caterpillars. While recent transcriptomic studies have analyzed gene expression profiles in A. japonicus specimens, offering a robust foundation for functional research focused on this parasitoid, no validated A. japonicus reference genes have yet been established, hampering further research efforts. Accordingly, this study sought to address this issue by screening for the most stable internal reference genes in A. japonicus samples to permit reliable RT-qPCR analyses. The utility of eight candidate reference genes (ACTIN, TATA, GAPDH, TUB, RPL13, RPS6, EF1α, RPS3a) was assessed under four different conditions by comparing developmental stages (larvae, pupae, adults), tissues (abdomen, chest, head), sex (male or female adults), or diapause states (diapause induction for 25, 35, 45, or 55 days, or diapause termination). RefFinder was used to calculate gene stability based on the integration of four algorithms (BestKeeper, Normfinder, geNorm, and ΔCt method) to determine the optimal RT-qPCR reference gene. Based on this approach, RPS6 and RPL13 were found to be the most reliable reference genes when assessing different stages of development, while ACTIN and EF1α were optimal when comparing adults of different sexes, RPL13 and EF1α were optimal when analyzing different tissues, and TATA and ACTIN were optimal for different diapause states. These results provide a valuable foundation for future RT-qPCR analyses of A. japonicus gene expression and function under a range of experimental conditions.

Validation of Appropriate Reference Genes for qRT–PCR Normalization in Oat (Avena sativa L.) under UV-B and High-Light Stresses

Oat is a food and forage crop species widely cultivated worldwide, and it is also an important forage grass in plateau regions of China, where there is a high level of ultraviolet radiation and sunlight. Screening suitable reference genes for oat under UV-B and high-light stresses is a prerequisite for ensuring the accuracy of real-time quantitative PCR (qRT–PCR) data used in plant adaptation research. In this study, eight candidate reference genes (sulfite oxidase, SUOX; victorin binding protein, VBP; actin-encoding, Actin1; protein PSK SIMULATOR 1-like, PSKS1; TATA-binding protein 2-like, TBP2; ubiquitin-conjugating enzyme E2, UBC2; elongation factor 1-alpha, EF1-α; glyceraldehyde-3-phosphate dehydrogenase 1, GAPDH1;) were selected based on previous studies and our oat transcriptome data. The expression stability of these reference genes in oat roots, stems, and leaves under UV-B and high-light stresses was first calculated using three frequently used statistical software (geNorm, NormFinder, and BestKeeper), and then the comprehensive stability of these genes was evaluated using RefFinder. The results showed that the most stably expressed reference genes in the roots, stems, and leaves of oat under UV-B stress were EF1-α, TBP2, and PSKS1, respectively; the most stably expressed reference genes in the roots, stems, and leaves under high-light stress were PSKS1, UBC2, and PSKS1, respectively. PSKS1 was the most stably expressed reference gene in all the samples. The reliability of the selected reference genes was further validated by analysis of the expression of the phenylalanine ammonia-lyase (PAL) gene. This study highlights reference genes for accurate quantitative analysis of gene expression in different tissues of oat under UV-B and high-light stresses.

Systematic Analysis and Functional Characterization of R2R3-MYB Genes in Scutellaria baicalensis Georgi

R2R3-MYB transcription factors participate in multiple critical biological processes, particularly as relates to the regulation of secondary metabolites. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine and possesses various bioactive attributes including anti-inflammation, anti-HIV, and anti-COVID-19 properties due to its flavonoids. In the current study, a total of 95 R2R3-MYB genes were identified in S. baicalensis and classified into 34 subgroups, as supported by similar exon–intron structures and conserved motifs. Among them, 93 R2R3-SbMYBs were mapped onto nine chromosomes. Collinear analysis revealed that segmental duplications were primarily responsible for driving the evolution and expansion of the R2R3-SbMYB gene family. Synteny analyses showed that the ortholog numbers of the R2R3-MYB genes between S. baicalensis and other dicotyledons had a higher proportion compared to that which is found from the monocotyledons. RNA-seq data indicated that the expression patterns of R2R3-SbMYBs in different tissues were different. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that 36 R2R3-SbMYBs from different subgroups exhibited specific expression profiles under various conditions, including hormone stimuli treatments (methyl jasmonate and abscisic acid) and abiotic stresses (drought and cold shock treatments). Further investigation revealed that SbMYB18/32/46/60/70/74 localized in the nucleus, and SbMYB18/32/60/70 possessed transcriptional activation activity, implying their potential roles in the regulatory mechanisms of various biological processes. This study provides a comprehensive understanding of the R2R3-SbMYBs gene family and lays the foundation for further investigation of their biological function.

Selection of suitable candidate genes for mRNA expression normalization in bulbil development of Pinellia ternata

Pinellia ternata (Thunb.) Breit. (Abbreviated as P. ternata). It is a commonly prescribed Chinese traditional medicinal herb for the treatment of phlegm, cough, and morning sick. Bulbil reproduction is one of the main reproductive methods of P. ternata. The accurate quantification of gene expression patterns associated with bulbil development might be helpful to explore the molecular mechanism involved in P. ternata reproduction. Quantitative real-time PCR was the most preferred method for expression profile and function analysis of mRNA. However, the reference genes in different tissues of P. ternata in different periods of bulbil development have not been studied in detail. In present study, the expression stability of eight candidate reference genes were determined with programs: geNorm, NormFinder, BestKeeper, and refFinder. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was identified as the top- rated reference gene in all samples of P. ternata, while different combinations of reference gene proved to be the most stable depending on development stage and tissue type. Furthermore, the reliability of GAPDH expression was verified by six P. ternata related genes in hormone and nutrient biosynthesis pathways, and the expression profiles of these genes were agreed with the results of RNA-seq digital gene expression analysis. These results can contribute to studies of gene expression patterns and functional analysis of P. ternata involved in bulbil development.

Baicalin ameliorates chronic unpredictable mild stress-induced depression through the BDNF/ERK/CREB signaling pathway

Brain-derived neurotrophic factor (BDNF) can activate the extracellular regulated protein kinase (ERK)/cAMP response element binding protein (CREB) cascade revealing an important role in antidepressant effects. Here, we studied the neuroprotective effect of baicalin (BA) in mice with chronic unpredictable mild stress (CUMS)-induced via a BDNF/ERK/CREB signaling pathway. Depression was induced via six weeks of CUMS in male ICR mice, and drug therapy was given simultaneously for the last three weeks. Cognitive dysfunctions were then evaluated via sucrose preference test (SPT), open field test (OFT), Morris water maze test (MWM), tail suspension test (TST), and novelty suppressed feeding test (NSF). Western blot and real-time PCR were then used to detect the relative expression of ERK, CREB, p-ERK, and p-CREB. Integrated optical density (IOD) tests of p-ERK and p-CREB were then evaluated via immunofluorescence. The behavior results showed that the cognitive dysfunctions increased in the CUMS group versus the control (CON) group (p < 0.01). There were decreases in fluoxetine (FLU) and BA groups (p < 0.05, p < 0.01). The protein ratios of p-ERK/ERK, p-CREB/CREB and ERK mRNA, and CREB mRNA expression decreased in the CUMS group (p < 0.01) and markedly increased in the FLU and BA groups (p < 0.05, p < 0.01). The IOD value of the p-ERK and p-CREB in the CUMS group was decreased versus the CON group (p < 0.01), and these changes were improved via BA and FLU treatment (p < 0.05, p < 0.01). This study indicated that BA can improve cognitive functions and has antidepressant effects in mice, which may be associated with activation of the BDNF/ERK/CREB signaling pathway in the hippocampus.