Selection and Verification of Appropriate Reference Genes for Expression Normalization in Cryptomeria fortunei under Abiotic Stress and Hormone Treatments

Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR

Quantitative reverse transcription PCR (qRT-PCR) can screen applicable reference genes of species, and reference genes can be used to reduce experimental errors. Sudan grass (Sorghum sudanense (Piper) Stapf) is a high-yield, abiotic-tolerant annual high-quality forage with a wide range of uses. However, no studies have reported reference genes suitable for Sudan grass. Therefore, we found eight candidate reference genes, including UBQ10, HIS3, UBQ9, Isoform0012931, PP2A, ACP2, eIF4α, and Actin, under salt stress (NaCl), drought stress (DR), acid aluminum stress (AlCl3), and methyl jasmonate treatment (MeJA). By using geNorm, NormFinder, BestKeeper, and RefFinder, we ranked eight reference genes on the basis of their expression stabilities. The results indicated that the best reference gene was PP2A under all treatments. eIF4α can be used in CK, MeJA, NaCl, and DR. HIS3 can serve as the best reference gene in AlCl3. Two target genes (Isoform0007606 and Isoform0002387) belong to drought-stress-response genes, and they are highly expressed in Sudan grass according to transcriptome data. They were used to verify eight candidate reference genes under drought stress. The expression trends of the two most stable reference genes were similar, but the trend in expression for Actin showed a significant difference. The reference genes we screened provided valuable guidance for future research on Sudan grass.

Reference gene selection for quantitative real-time PCR analysis of Hymenopellis radicata under abiotic stress

Hymenopellis radicata (H. radicata) is an edible fungus rich in protein and mineral elements, with high edible and medical value. And reference genes suitable for normalization of qRT-PCR data from this species have not been investigated. In this study, therefore, we selected 11 housekeeping genes common in biology. The expression levels of these housekeeping genes were measured in three different tissues and six different abiotic stress treatments in mycelium. They were evaluated for expression stability using online tools. The results showed that gene ACT could be stable expressed in all samples. The expressions of genes TUB and UBQ10 are the most stable under heat stress, ACT and EF are the most stable genes under salt stress, ACT and TUB are the most stable genes under oxidation stress, RPL6 and EF are the most stable genes under pH condition, ACT and RPB2 are the most stable genes under cadmium stress, and RPB2 and UBC are the most stable genes under drought condition. ACT and PP2A are the most stable genes at different tissue sites. This study is of great help to explore the gene expression pattern of H. radicata, and also provides reference for internal reference gene screening under other conditions.

Reference Genes Screening and Gene Expression Patterns Analysis Involved in Gelsenicine Biosynthesis under Different Hormone Treatments in Gelsemium elegans

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) is an accurate method for quantifying gene expression levels. Choosing appropriate reference genes to normalize the data is essential for reducing errors. Gelsemium elegans is a highly poisonous but important medicinal plant used for analgesic and anti-swelling purposes. Gelsenicine is one of the vital active ingredients, and its biosynthesis pathway remains to be determined. In this study, G. elegans leaf tissue with and without the application of one of four hormones (SA, MeJA, ETH, and ABA) known to affect gelsenicine synthesis, was analyzed using ten candidate reference genes. The gene stability was evaluated using GeNorm, NormFinder, BestKeeper, ∆CT, and RefFinder. The results showed that the optimal stable reference genes varied among the different treatments and that at least two reference genes were required for accurate quantification. The expression patterns of 15 genes related to the gelsenicine upstream biosynthesis pathway was determined by RT-qPCR using the relevant reference genes identified. Three genes 8-HGO, LAMT, and STR, were found to have a strong correlation with the amount of gelsenicine measured in the different samples. This research is the first study to examine the reference genes of G. elegans under different hormone treatments and will be useful for future molecular analyses of this medically important plant species.

Identification of Accurate Reference Genes for qRT-PCR Analysis of Gene Expression in Eremochloa ophiuroides under Multiple Stresses of Phosphorus Deficiency and/or Aluminum Toxicity

Centipedegrass (Eremochloa ophiuroides (Munro.) Hack.) is a species originating in China and is an excellent warm-season turfgrass. As a native species in southern China, it is naturally distributed in the phosphorus-deficient and aluminum-toxic acid soil areas. It is important to research the molecular mechanism of centipedegrass responses to phosphorus-deficiency and/or aluminum-toxicity stress. Quantitative Real-Time PCR (qRT-PCR) is a common method for gene expression analysis, and the accuracy of qRT-PCR results depends heavily on the stability of internal reference genes. However, there are still no reported stable and effective reference genes for qRT-PCR analysis of target genes under the acid-soil-related stresses in different organs of centipedegrass. For scientific rigor, the gene used as a reference for any plant species and/or any stress conditions should be first systematically screened and evaluated. This study is the first to provide a group of reliable reference genes to quantify the expression levels of functional genes of Eremochloa ophiuroides under multiple stresses of P deficiency and/or aluminum toxicity. In this study, centipedegrass seedlings of the acid-soil-resistant strain 'E041' and acid-soil-sensitive strain 'E089' were used for qRT-PCR analysis. A total of 11 candidate reference genes (ACT, TUB, GAPDH, TIP41, CACS, HNR, EP, EF1α, EIF4α, PP2A and actin) were detected by qRT-PCR technology, and the stability of candidate genes was evaluated with the combination of four internal stability analysis software programs. The candidate reference genes exhibited differential stability of expression in roots, stems and leaves under phosphorus-deficiency and/or aluminum-toxicity stress. On the whole, the results showed that GAPDH, TIP41 and HNR were the most stable in the total of samples. In addition, for different tissues under various stresses, the selected reference genes were also different. CACS and PP2A were identified as two stable reference genes in roots through all three stress treatments (phosphate deficiency, aluminum toxicity, and the multiple stress treatment of aluminum toxicity and phosphate deficiency). Moreover, CACS was also stable as a reference gene in roots under each treatment (phosphate deficiency, aluminum toxicity, or multiple stresses of aluminum toxicity and phosphate deficiency). In stems under all three stress treatments, GAPDH and EIF4α were the most stable reference genes; for leaves, PP2A and TIP41 showed the two highest rankings in all three stress treatments. Finally, qRT-PCR analysis of the expression patterns of the target gene ALMT1 was performed to verify the selected reference genes. The application of the reference genes identified as internal controls for qRT-PCR analysis will enable accurate analysis of the target gene expression levels and expression patterns in centipedegrass under acid-soil-related stresses.

Selection of Suitable Reference Genes for RT-qPCR Gene Expression Analysis in Centipedegrass under Different Abiotic Stress

As a C4 warm-season turfgrass, centipedegrass (Eremochloa ophiuroides (Munro) Hack.) is known for its exceptional resilience to intensive maintenance practices. In this research, the most stably expressed reference genes in the leaves of centipedegrass under different stress treatments, including salt, cold, drought, aluminum (Al), and herbicide, were screened by the quantitative real-time PCR (RT-qPCR) technique. The stability of 13 candidate reference genes was evaluated by software GeNorm V3.4, NormFinder V20, BestKeeper V1.0, and ReFinder V1.0. The results of this experiment demonstrated that the expression of the UBC (ubiquitin-conjugating enzyme) remained the most stable under cold and Al stress conditions. On the other hand, the MD (malate dehydrogenase) gene exhibited the best performance in leaf tissues subjected to salt and drought stresses. Under herbicide stress, the expression level of the RIP (60S ribosomal protein L2) gene ranked the highest. The expression levels of abiotic stress-associated genes such as PIP1, PAL, COR413, ALMT9, and BAR were assessed to validate the reliability of the selected reference genes. This study provides valuable information and reference points for gene expression under abiotic stress conditions in centipedegrass.

Two R2R3-MYB transcription factors from Chinese cedar (Cryptomeria fortunei Hooibrenk) are involved in the regulation of secondary cell wall formation

As the most abundant renewable energy source, wood comprises the secondary cell wall (SCW). SCW biosynthesis involves lignin and cellulose deposition. Increasing studies have illustrated that R2R3-MYB transcription factors (TFs) play pivotal roles in affecting lignin accumulation and SCW formation. Nevertheless, the regulatory roles of R2R3-MYBs are still unresolved in Cryptomeria fortunei Hooibrenk cambium and wood formation. To dissect the potentials of CfMYBs, we successfully cloned and intensively studied the functions of CfMYB4 and CfMYB5 in SCW formation and abiotic stress response. They both contained the conserved MYB domain capable of forming a special structure that could bind to the core motifs of downstream genes. The phylogenetic tree implied that two CfMYBs clustered into different evolutionary branches. They were predominantly expressed in the stem and were localized to the nucleus. Furthermore, CfMYB4 functioned as an activator to enhance lignin and cellulose accumulation, and increase the SCW thickness by elevating the expression levels of SCW-related genes. By contrast, CfMYB5 negatively regulated lignin and cellulose biosynthesis, and decreased SCW formation by reducing the expression of SCW biosynthetic genes. Our data not only highlight the regulatory functions of CfMYBs in lignin deposition but also provide critical insights into the development of strategies for the genetic improvement of Cryptomeria fortunei wood biomass.

Selection of reference genes for RT-qPCR analysis of rice with Rhizoctonia solani infection and biocontrol PGPR/KSi application

BACKGROUND

Rhizoctonia solani (AG1 IA) is an important pathogen of rice (Oryza sativa L.) that causes rice sheath blight (RSB). Since control of RSB by breeding and fungicides have had limited success, novel strategies like biocontrol with plant growth-promoting rhizobacteria (PGPR) can be an effective alternative.

METHOD AND RESULTS

Seven commonly used reference genes (RGs), 18SrRNA, ACT1, GAPDH2, UBC5, RPS27, eIF4a and CYP28, were evaluated for their stability in rice-R. solani-PGPR interaction for real-time quantitative PCR (RT-qPCR) analysis. Different algorithms were examined, Delta Ct, geNorm, NormFinder, BestKeeper, and comprehensive ranking by RefFinder, to evaluate RT-qPCR of rice in tissues infected with R. solani and treated with the PGPR strains, Pseudomonas saponiphilia and Pseudomonas protegens, with potassium silicate (KSi) alone or in combination with each PGPR strain. RG stability was affected for each treatment and treatment-specific RG selection was suggested. Validation analysis was done for nonexpressor of PR-1(NPR1) for each treatment.

CONCLUSION

Overall, ACT1 was the most stable RG with R. solani infection alone, GAPDH2 with R. solani infection plus KSi, UBC5 with R. solani infection plus P. saponiphilia, and eIF4a with R. solani infection plus P. protegens. Both ACT1 and RPS27 were the most stable with the combination of KSi and P. saponiphilia, while RPS27 was the most stable with the combination of KSi and P. protegens.

Physiological and gene expression changes of Cryptomeria fortunei Hooibrenk families under heat stress

Heat stress is one of the major abiotic stresses affecting plant growth and productivity. Cryptomeria fortunei (Chinese cedar) is an excellent timber and landscaping tree species in southern China thanks to its beautiful appearance, straight texture and ability to purify the air and improve the environment. In this study, we first screened 8 excellent C. fortunei families (#12, #21, #37, #38, #45, #46, #48, #54) in a second generation seed orchard. We then analyzed the electrolyte leakage (EL) and lethal temperature at 50% (LT₅₀) values under heat stress, to identify the families with the best heat resistance (#48) and the lowest heat resistance (#45) and determine the physiological and morphological response of different threshold-resistance of C. fortune to heat stress. The relative conductivity of the C. fortunei families showed an increasing trend with increasing temperature, following an "S" curve, and the half-lethal temperature ranges between 39°C and 43.2°C. The activities of SOD and POD fluctuated in the early stage of stress but decreased after 37°C. We observed the changes in the cell ultrastructure at 43°C, and the mesophyll cell structure of #48 was less damaged than that of #45. Eight heat resistance gene, including CfAPX1, CfAPX2, CfHSP11, CfHSP21, CfHSP70, CfHSFA1a, CfHSFB2a and CfHSFB4, were all up-regulated in #45 and #48, and there were significant differences between #45 and #48 under different heat stress treatments. We found a significant difference in heat tolerance between #45 and #48, such that #48 shows higher heat tolerance capability and could be exploited in breeding programs. We conclude that the strongly heat-resistant family had a more stable physiological state and a wider range of heat stress adaptations.

Selection and validation of appropriate reference genes for RT-qPCR analysis of Nitraria sibirica under various abiotic stresses

BACKGROUND

Nitraria sibirica Pall. is a halophytic shrub with strong environmental adaptability that can survive in extremely saline-alkali and drought-impacted environments. Gene expression analysis aids in the exploration of the molecular mechanisms of plant responses to abiotic stresses. RT-qPCR is the most common technique for studying gene expression. Stable reference genes are a prerequisite for obtaining accurate target gene expression results in RT-qPCR analysis.

RESULTS

In this study, a total of 10 candidate reference genes were selected from the transcriptome of N. sibirica, and their expression stability in leaves and roots under different treatment conditions (salt, alkali, drought, cold, heat and ABA) was evaluated with the geNorm, NormFinder, BestKeeper, comparative ΔCt and RefFinder programs. The results showed that the expression stability of the candidate reference genes was dependent on the tissue and experimental conditions tested. ACT7 combined with R3H, GAPDH, TUB or His were the most stable reference genes in the salt- or alkali-treated leaves, salt-treated roots and drought-treated roots, respectively; R3H and GAPDH were the most suitable combination for drought-treated leaves, heat-treated root samples and ABA-treated leaves; DIM1 and His maintained stable expression in roots under alkali stress; and TUB combined with R3H was stable in ABA-treated roots. TBCB and GAPDH exhibited stable expression in heat-treated leaves; TBCB, R3H, and ERF3A were stable in cold-treated leaves; and the three most stable reference genes for cold-treated roots were TBCB, ACT11 and DIM1. The reliability of the selected reference genes was further confirmed by evaluating the expression patterns of the NsP5CS gene under the six treatment conditions.

CONCLUSION

This study provides a theoretical reference for N. sibirica gene expression standardization and quantification under various abiotic stress conditions and will help to reveal the molecular mechanisms that confer stress tolerance to N. sibirica.

Transcriptome and metabolome changes in Chinese cedar during cold acclimation reveal the roles of flavonoids in needle discoloration and cold resistance

Cryptomeria fortunei growth and development are usually affected by low temperatures. Despite the evergreen nature of this species, most needles turn yellowish-brown in cold winters. The underlying discoloration mechanisms that cause this phenomenon in response to cold acclimation remain poorly understood. Here, we measured the pigment content and ultrastructure of normal wild-type (Wt) and evergreen mutant (GM) C. fortunei needles and performed integrated transcriptomic and metabolomic analyses to explore potential discoloration mechanisms. The results showed that the needle chlorophyll content of these two genotypes decreased in winter. Wt needles showed greater decrease in the chlorophyll content and local destruction of chloroplast ultrastructure and contained larger amounts of flavonoids than GM needles, as shown by metabolomics analysis. We subsequently identified key differentially expressed genes in the flavonoid biosynthesis pathway and observed significantly upregulated flavonol synthase expression in Wt needles compared with GM needles that significantly increased the anthoxanthin (flavones and flavonols) content, which is likely a key factor underlying the difference in needle color between these two genotypes. Therefore, flavonoid metabolism may play important roles in the cold resistance and needle discoloration of C. fortunei, and our results provide an excellent foundation for the molecular mechanism of C. fortunei in response to cold stress.

Selection of Suitable Reference Genes for Gene Expression Normalization Studies in Dendrobium huoshanense

Dendrobium huoshanense is a kind of precious herb with important medicinal and edible value in China, which is widely used in traditional Chinese medicine for various diseases. Recent studies have paid close attention to the genetic expression of the biosynthetic pathway of the main active components (polysaccharides, alkaloids, and flavonoids), and real-time polymerase chain reaction (qPCR) is one of the most widely used methods for doing so. However, so far, no reference gene selections have been reported in D. huoshanense. In this study, 15 reference gene candidates (GAPDH, eIF, EF-1α, PP2A, UBCE, RPL5, TBP, APT1, MDH, PTBP3, PEPC, CYP71, NCBP2, TIP41, and F-box) were selected and evaluated for their expression stability in D. huoshanense under various experimental conditions, including in different tissues (root, stem, and leaf), abiotic stresses (oxidative, drought, cold, and UV), and hormone treatment (methyl jasmonate) using three statistical programs (geNorm, NormFinder, and BestKeeper). Then, the RefFinder program was employed to comprehensively validate the stability of the selected reference genes. Finally, the expression profiles of the CESA and GMPP genes were further analyzed, and these results indicated that TBP, NCBP2, and CYP71 were the top three most stable reference genes after comprehensive comparison, which could be used as stable reference genes for normalizing the genes expression in D. huoshanense. This study described here provides the first data regarding on reference gene selection in D. huoshanense, which will be extremely beneficial for future research on the gene expression normalization in D. huoshanense.

Transcriptome and Metabolome Analyses Reveal Differences in Terpenoid and Flavonoid Biosynthesis in Cryptomeria fortunei Needles Across Different Seasons

Cryptomeria fortunei (Chinese cedar) has outstanding medicinal value due to its abundant flavonoid and terpenoid contents. The metabolite contents of C. fortunei needles differ across different seasons. However, the biosynthetic mechanism of these differentially synthesized metabolites (DSMs) is poorly understood. To improve our understanding of this process, we performed integrated non-targeted metabolomic liquid chromatography and gas chromatography mass spectrometry (LC-MS and GC-MS), and transcriptomic analyses of summer and winter needles. In winter, the C. fortunei needle ultrastructure was damaged, and the chlorophyll content and F _v/F _m were significantly (p < 0.05) reduced. Based on GC-MS and LC-MS, we obtained 106 and 413 DSMs, respectively; based on transcriptome analysis, we obtained a total of 41.17 Gb of clean data and assembled 33,063 unigenes, including 14,057 differentially expressed unigenes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DSMs/DEGs were significantly (p < 0.05) enriched in many biosynthesis pathways, such as terpenoids, photosynthates, and flavonoids. Integrated transcriptomic and metabonomic analyses showed that seasonal changes have the greatest impact on photosynthesis pathways, followed by terpenoid and flavonoid biosynthesis pathways. In summer Chinese cedar (SCC) needles, DXS, DXR, and ispH in the 2-methyl-pentaerythritol 4-phosphate (MEP) pathway and GGPS were highly expressed and promoted the accumulation of terpenoids, especially diterpenoids. In winter Chinese cedar (WCC) needles, 9 genes (HCT, CHS, CHI, F3H, F3'H, F3'5'H, FLS, DFR, and LAR) involved in flavonoid biosynthesis were highly expressed and promoted flavonoid accumulation. This study broadens our understanding of the metabolic and transcriptomic changes in C. fortunei needles caused by seasonal changes and provides a reference regarding the adaptive mechanisms of C. fortunei and the extraction of its metabolites.

CfAPX, a cytosolic ascorbate peroxidase gene from Cryptomeria fortunei, confers tolerance to abiotic stress in transgenic Arabidopsis

Plants subjected to biotic or abiotic stresses produce a large amount of reactive oxygen species (ROS). If ROS cannot be cleared in time, they cause a series of harmful reactions in plants. Ascorbate peroxidase (APX) is a key enzyme that removes ROS from plant cells and plays a vital role in plant stress resistance. However, to date, no studies on APX homologs in Cryptomeria fortunei have been reported. In this study, we isolated complementary DNA (cDNA) encoding APXfrom C. fortunei needles, which is referred to as CfAPX, by rapid amplification of cDNA ends (RACE). The full-length CfAPX sequence was 1226 bp in length and included a 750-bp open reading frame (ORF) encoding a protein of 249 amino acids. Phylogenetic analysis showed that APXs of different plant species have been highly evolutionarily conserved. CfAPX was shown to belong to the cytoplasmic subgroup and was more closely related to GbAPX of the gymnosperm Ginkgo biloba. CfAPX showed no transcriptional activity in yeast cells but was highly expressed in cones. To better handle abiotic stresses, compared with wild-type (WT) Arabidopsis thaliana, 35S::CfAPX transgenic Arabidopsis strongly expressed CfAPX, presented increased antioxidant enzyme activities, ascorbic acid (AsA) contents, chlorophyll levels and fluorescence parameter and reduced malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) contents. In addition, CfAPX expression in C. fortunei was mostly upregulated under stress. In summary, CfAPX confers abiotic stress responses to plants, which provides a scientific basis for subsequent breeding for increased stress resistance in C. fortunei.

Transcriptome Profiling of Maize (Zea mays L.) Leaves Reveals Key Cold-Responsive Genes, Transcription Factors, and Metabolic Pathways Regulating Cold Stress Tolerance at the Seedling Stage

Cold tolerance is a complex trait that requires a critical perspective to understand its underpinning mechanism. To unravel the molecular framework underlying maize (Zea mays L.) cold stress tolerance, we conducted a comparative transcriptome profiling of 24 cold-tolerant and 22 cold-sensitive inbred lines affected by cold stress at the seedling stage. Using the RNA-seq method, we identified 2237 differentially expressed genes (DEGs), namely 1656 and 581 annotated and unannotated DEGs, respectively. Further analysis of the 1656 annotated DEGs mined out two critical sets of cold-responsive DEGs, namely 779 and 877 DEGs, which were significantly enhanced in the tolerant and sensitive lines, respectively. Functional analysis of the 1656 DEGs highlighted the enrichment of signaling, carotenoid, lipid metabolism, transcription factors (TFs), peroxisome, and amino acid metabolism. A total of 147 TFs belonging to 32 families, including MYB, ERF, NAC, WRKY, bHLH, MIKC MADS, and C₂H₂, were strongly altered by cold stress. Moreover, the tolerant lines’ 779 enhanced DEGs were predominantly associated with carotenoid, ABC transporter, glutathione, lipid metabolism, and amino acid metabolism. In comparison, the cold-sensitive lines’ 877 enhanced DEGs were significantly enriched for MAPK signaling, peroxisome, ribosome, and carbon metabolism pathways. The biggest proportion of the unannotated DEGs was implicated in the roles of long non-coding RNAs (lncRNAs). Taken together, this study provides valuable insights that offer a deeper understanding of the molecular mechanisms underlying maize response to cold stress at the seedling stage, thus opening up possibilities for a breeding program of maize tolerance to cold stress.

Integrated Four Comparative-Omics Reveals the Mechanism of the Terpenoid Biosynthesis in Two Different Overwintering Cryptomeria fortunei Phenotypes

Chinese cedar (Cryptomeria fortunei) is a tree species with important ornamental, medicinal, and economic value. Terpenoids extracted from the essential oil of C. fortunei needles have been considered valuable ingredients in the pharmaceutical and cosmetic industries. However, the possible gene regulation mechanisms that limit terpenoid biosynthesis in this genus are poorly understood. Here, we adopted integrated metabolome analysis, transcriptome, small-RNA (sRNA), and degradome sequencing to analyze the differences in terpenoid regulatory mechanisms in two different overwintering C. fortunei phenotypes (wild-type and an evergreen mutant). A total of 1447/6219 differentially synthesized metabolites (DSMs)/unigenes (DEGs) were detected through metabolome/transcriptome analyses, and these DSMs/DEGs were significantly enriched in flavonoid and diterpenoid biosynthesis pathways. In C. fortunei needles, 587 microRNAs (miRNAs), including 67 differentially expressed miRNAs (DERs), were detected. Among them, 8346 targets of 571 miRNAs were predicted using degradome data, and a 72-miRNA-target regulatory network involved in the metabolism of terpenoids and polyketides was constructed. Forty-one targets were further confirmed to be involved in terpenoid backbone and diterpenoid biosynthesis, and target analyses revealed that two miRNAs (i.e., aly-miR168a-5p and aof-miR396a) may be related to the different phenotypes and to differential regulation of diterpenoid biosynthesis. Overall, these results reveal that C. fortunei plants with the evergreen mutation maintain high terpenoid levels in winter through miRNA-target regulation, which provides a valuable resource for essential oil-related bioengineering research.