Selection and validation of reference genes for gene expression analysis in switchgrass (Panicum virgatum) using quantitative real-time RT-PCR

Characterizing reference genes for high-fidelity gene expression analysis under different abiotic stresses and elicitor treatments in fenugreek leaves

BACKGROUND

Quantifying gene expression is a critical aspect of applied genomics research across all organisms, and real-time PCR has emerged as a powerful tool for this purpose. However, selecting appropriate internal control genes for data normalization presents specific challenges. This study aimed to identify suitable reference genes for gene expression analysis under various conditions, encompassing salinity, low and high-temperature stresses, and different elicitor treatments. These treatments included titanium dioxide, cold plasma, 24-epibrassinolide, and melatonin, resulting in a total of 13 unique treatments and 148 treatment combinations applied to fenugreek plants.

RESULTS

As per the analysis performed with the BestKeeper tool, EEF-1α, and GAPDH were recognized as the most stable reference genes under the majority of conditions. Furthermore, the GeNorm and NormFinder tools identified β-tubulin and EEF-1α as the most stable reference genes. The findings of this research demonstrated that, although the stability of three reference genes expression was acceptable in almost all evaluated treatments, fluctuations in their expression were observed under the treatments of cold stress with TiO2 NPs application, cold plasma application with salinity stress, and cold plasma application with high-temperature stress compared to others. Simultaneously, the GeNorm analysis results demonstrated that in the mentioned treatments, relying on only one reference gene is inadequate. To corroborate the results, we examined the expression profile of the SSR gene, a pivotal gene in diosgenin biosynthesis, under all investigated treatments and treatment combinations. The outcomes suggested that employing stable reference genes yielded highly consistent results.

CONCLUSIONS

The varying expression patterns of the target genes emphasize the crucial need for precise optimization of experimental conditions and selecting stable reference genes to achieve accurate results in gene expression studies utilizing real-time PCR. These findings offer valuable insights into the selection of appropriate reference genes for gene expression analysis under diverse conditions using real-time PCR.

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.

Selection of Reference Genes in Evodia rutaecarpa var. officinalis and Expression Patterns of Genes Involved in Its Limonin Biosynthesis

E. rutaecarpa var. officinalis is a traditional Chinese medicinal plant known for its therapeutic effects, which encompass the promotion of digestion, the dispelling of cold, the alleviation of pain, and the exhibition of anti-inflammatory and antibacterial properties. The principal active component of this plant, limonin, is a potent triterpene compound with notable pharmacological activities. Despite its significance, the complete biosynthesis pathway of limonin in E. rutaecarpa var. officinalis remains incompletely understood, and the underlying molecular mechanisms remain unexplored. The main purpose of this study was to screen the reference genes suitable for expression analysis in E. rutaecarpa var. officinalis, calculate the expression patterns of the genes in the limonin biosynthesis pathway, and identify the relevant enzyme genes related to limonin biosynthesis. The reference genes play a pivotal role in establishing reliable reference standards for normalizing the gene expression data, thereby ensuring precision and credibility in the biological research outcomes. In order to identify the optimal reference genes and gene expression patterns across the diverse tissues (e.g., roots, stems, leaves, and flower buds) and developmental stages (i.e., 17 July, 24 August, 1 September, and 24 October) of E. rutaecarpa var. officinalis, LC-MS was used to analyze the limonin contents in distinct tissue samples and developmental stages, and qRT-PCR technology was employed to investigate the expression patterns of the ten reference genes and eighteen genes involved in limonin biosynthesis. Utilizing a comprehensive analysis that integrated three software tools (GeNorm ver. 3.5, NormFinder ver. 0.953 and BestKeeper ver. 1.0) and Delta Ct method alongside the RefFinder website, the best reference genes were selected. Through the research, we determined that Act1 and UBQ served as the preferred reference genes for normalizing gene expression during various fruit developmental stages, while Act1 and His3 were optimal for different tissues. Using Act1 and UBQ as the reference genes, and based on the different fruit developmental stages, qRT-PCR analysis was performed on the pathway genes selected from the "full-length transcriptome + expression profile + metabolome" data in the limonin biosynthesis pathway of E. rutaecarpa var. officinalis. The findings indicated that there were consistent expression patterns of HMGCR, SQE, and CYP450 with fluctuations in the limonin contents, suggesting their potential involvement in the limonin biosynthesis of E. rutaecarpa var. officinalis. This study lays the foundation for further research on the metabolic pathway of limonin in E. rutaecarpa var. officinalis and provides reliable reference genes for other researchers to use for conducting expression analyses.

Genome-wide identification, classification, and expression analysis of heat shock transcription factor family in switchgrass (Panicum virgatum L.)

Switchgrass is one of the most promising bioenergy crops and is generally cultivated in arid climates and poor soils. Heat shock transcription factors (Hsfs) are key regulators of plant responses to abiotic and biotic stressors. However, their role and mechanism of action in switchgrass have not been elucidated. Hence, this study aimed to identify the Hsf family in switchgrass and understand its functional role in heat stress signal transduction and heat tolerance by using bioinformatics and RT-PCR analysis. Forty-eight PvHsfs were identified and divided into three main classes based on their gene structure and phylogenetic relationships: HsfA, HsfB, and HsfC. The results of the bioinformatics analysis showed a DNA-binding domain (DBD) at the N-terminal in PvHsfs, and they were not evenly distributed on all chromosomes except for chromosomes 8 N and 8 K. Many cis-elements related to plant development, stress responses, and plant hormones were identified in the promoter sequence of each PvHsf. Segmental duplication is the primary force underlying Hsf family expansion in switchgrass. The results of the expression pattern of PvHsfs in response to heat stress showed that PvHsf03 and PvHsf25 might play critical roles in the early and late stages of switchgrass response to heat stress, respectively, and HsfB mainly showed a negative response to heat stress. Ectopic expression of PvHsf03 in Arabidopsis significantly increased the heat resistance of seedlings. Overall, our research lays a notable foundation for studying the regulatory network in response to deleterious environments and for further excavating tolerance genes in switchgrass.

Genome-Wide Datasets of Chicories (Cichorium intybus L.) for Marker-Assisted Crop Breeding Applications: A Systematic Review and Meta-Analysis

Cichorium intybus L. is the most economically important species of its genus and among the most important of the Asteraceae family. In chicory, many linkage maps have been produced, several sets of mapped and unmapped markers have been developed, and dozens of genes linked to traits of agronomic interest have been investigated. This treasure trove of information, properly cataloged and organized, is of pivotal importance for the development of superior commercial products with valuable agronomic potential in terms of yield and quality, including reduced bitter taste and increased inulin production, as well as resistance or tolerance to pathogens and resilience to environmental stresses. For this reason, a systematic review was conducted based on the scientific literature published in chicory during 1980-2023. Based on the results obtained from the meta-analysis, we created two consensus maps capable of supporting marker-assisted breeding (MAB) and marker-assisted selection (MAS) programs. By taking advantage of the recently released genome of C. intybus, we built a 639 molecular marker-based consensus map collecting all the available mapped and unmapped SNP and SSR loci available for this species. In the following section, after summarizing and discussing all the genes investigated in chicory and related to traits of interest such as reproductive barriers, sesquiterpene lactone biosynthesis, inulin metabolism and stress response, we produced a second map encompassing 64 loci that could be useful for MAS purposes. With the advent of omics technologies, molecular data chaos (namely, the situation where the amount of molecular data is so complex and unmanageable that their use becomes challenging) is becoming far from a negligible issue. In this review, we have therefore tried to contribute by standardizing and organizing the molecular data produced thus far in chicory to facilitate the work of breeders.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Foxtail mosaic virus-induced gene silencing (VIGS) in switchgrass (Panicum virgatum L.)

BACKGROUND

Although the genome for the allotetraploid bioenergy crop switchgrass (Panicum virgatum) has been established, limitations in mutant resources have hampered in planta gene function studies toward crop optimization. Virus-induced gene silencing (VIGS) is a versatile technique for transient genetic studies. Here we report the implementation of foxtail mosaic virus (FoMV)-mediated gene silencing in switchgrass in above- and below-ground tissues and at different developmental stages.

RESULTS

The study demonstrated that leaf rub-inoculation is a suitable method for systemic gene silencing in switchgrass. For all three visual marker genes, Magnesium chelatase subunit D (ChlD) and I (ChlI) as well as phytoene desaturase (PDS), phenotypic changes were observed in leaves, albeit at different intensities. Gene silencing efficiency was verified by RT-PCR for all tested genes. Notably, systemic gene silencing was also observed in roots, although silencing efficiency was stronger in leaves (~ 63-94%) as compared to roots (~ 48-78%). Plants at a later developmental stage were moderately less amenable to VIGS than younger plants, but also less perturbed by the viral infection.

CONCLUSIONS

Using FoMV-mediated VIGS could be achieved in switchgrass leaves and roots, providing an alternative approach for studying gene functions and physiological traits in this important bioenergy crop.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Tibetan Medicinal Plant Saussurea Laniceps Callus under Abiotic Stresses and Hormone Treatments

Real-time quantitative PCR (RT-qPCR) is an important technique for studying gene expression analysis, but accurate and reliable results depend on the use of a stable reference gene. This study proposes to test the expression stability of candidate reference genes in the callus of Saussurea laniceps, a unique Tibetan medicinal plant. Based on the S. laniceps callus transcriptome, eleven candidate reference genes, including TUA2, TUB3, TUB8, TIF3B1, TIF3H1, ELF5A, PP2AA2, UEV1D, UBL5, UBC36, and SKIP1), were validated for RT-qPCR normalization in the callus under abiotic stress (salt, cold, and UV) and hormone treatments (abscisic acid, MeJA, and salicylic acid). The stability of the candidate genes was evaluated in all the samples of S. laniceps. Comprehensive analysis of all samples showed that the best reference genes were UBC36 and UBL5. ELF5A and TIF3B1 were ranked as the most stable genes in the sample sets under abiotic stress. For hormone stimulation, UBC36 and TIF3H1 genes had the best stability. This study provides useful guidelines and a starting point for reference gene selection for expression analysis using RT-qPCR techniques in S. laniceps.

Overexpression of PvSTK1 gene from Switchgrass (Panicum virgatum L.) affects flowering time and development of floral organ in transgenic Arabidopsis thaliana

Flowering means that the plant enters the reproductive growth stage, which is a crucial stage in the plant life cycle. Delaying flowering time to prolong vegetative growth is an important method to increase biomass yield and saccharification efficiency in switchgrass, It is of great significance to study the molecular mechanism of plant flowering and regulate the process of plant flowering in the process of biomass production. In this study, we identified 55 serine/threonine-protein kinase genes related to flower development from the switchgrass transcriptome database. Simultaneously, we cloned one of them, PvSTK1, whose expression level and differential fold were significantly higher than other members. PvSTK1 is located on chromosome 8N and its protein was in the cell membrane, cytoplasm, and nucleus. The spatio-temporal expression analysis of the PvSTK1 in switchgrass displayed that the PvSTK1 is crucial in vegetative period, however, not in the transition to reproductive period. Overexpression of PvSTK1 in Arabidopsis resulted in down-regulation of flower-promoting genes and up-regulation of flower-suppressing genes, thereby delaying flowering. In addition, PvSTK1 caused atrophy of the ovules of the florets at the base of the inflorescence, leading to sterility of the florets. The function of PvSTK1 is to inhibit the development of floral organs, and its overexpression can prolong its vegetative period. In the future, overexpression of the PvSTK1 gene in switchgrass will change the flowering time and increase yield and utilization efficiency of biomass.

Haploidy and aneuploidy in switchgrass mediated by misexpression of CENH3

Cross bred species such as switchgrass may benefit from advantageous breeding strategies requiring inbred lines. Doubled haploid production methods offer several ways that these lines can be produced that often involve uniparental genome elimination as the rate limiting step. We have used a centromere-mediated genome elimination strategy in which modified CENH3 is expressed to induce the process. Transgenic tetraploid switchgrass lines coexpressed Cas9, a poly-cistronic tRNA-gRNA tandem array containing eight guide RNAs that target two CENH3 genes, and different chimeric versions of CENH3 with alterations to the N-terminal tail region. Genotyping of CENH3 genes in transgenics identified edits including frameshift mutations and deletions in one or both copies of the two CENH3 genes. Flow cytometry of T₁ seedlings identified two T₀ lines that produced five haploid individuals representing an induction rate of 0.5% and 1.4%. Eight different T₀ lines produced aneuploids at rates ranging from 2.1 to 14.6%. A sample of aneuploid lines were sequenced at low coverage and aligned to the reference genome, revealing missing chromosomes and chromosome arms.

Reference gene selection for quantitative RT-PCR in Miscanthus sacchariflorus under abiotic stress conditions

BACKGROUND

Reference genes are necessary for quantitative real-time PCR (qRT-PCR) analysis and their stability can directly influence the accuracy of gene expression result. Miscanthus sacchariflorus, a perennial C4 grass that serves as promising biofuel plant for temperate climates, has not been explored for the identification of stable reference genes yet.

MATERIALS AND METHODS

Nine potential reference genes (ACT, EF1a, FBOX, GAPDH, PP2A, SAND, TIP41, TUB and UBC) of M. sacchariflorus under different abiotic (salinity, drought and cadmium) stresses, as well as in two tissues (roots and leaves) were evaluated. The expression stability of these genes were analyzed by four commonly used software programs (geNorm, NormFinder, BestKeeper, ΔCt method and RefFinder).

RESULTS

Our results found that FBOX and SAND are the most stable genes among all tested samples. FBOX and EF1a are suitable for gene expression normalization of cadmium-treated samples and salinity-treated leaves. FBOX and PP2A are appropriate reference genes for salt-stressed roots and PEG-treated leaves. The traditional reference gene ACT and GAPDH exhibited the most variable pattern, which would not be recommended for qRT-PCR analysis under different abiotic stresses. Furthermore, the expression levels of PIP2, NHX1 and MT2a under drought, salt and cadmium treatment were detected with above reference genes.

CONCLUSIONS

This work identified the appropriate reference genes for qRT-PCR in M. sacchariflorus and FBOX was recommended to be effective internal control for gene expression normalization in M. sacchariflorus in response to different abiotic stresses.

Establishment of reference (housekeeping) genes via quantitative real-time PCR for investigation of the genomic basis of abiotic stress resistance in Psammochloa villosa (Poaceae)

Psammochloa villosa is a desert plant growing in Northwest China with considerable resistance to abiotic stress, including drought, cold, and salt. To facilitate future studies of stress resistance in Psammochloa villosa, we sought to establish a suite of reference (or housekeeping) genes for utilization within future gene expression studies. Specifically, we selected nine candidate genes based on prior studies and new transcriptomic data for P. villosa, and we evaluated their expression stability in three different tissues of P. villosa under different treatments simulating abiotic stress conditions using four different bioinformatics assessments. Our results showed that TIP41 (TIP41-like family protein) was the most stable reference gene in drought- and salt-stressed leaves and salt-stressed stems, ELF-1α (elongation factor 1-α) was the most stable in cold-stressed leaves and drought- and salt-stressed roots, ACT (actin) was the most stable in drought-stressed stems, TUA (α-tubulin) was the most stable in cold-stressed stems, and 18S rRNA (18S ribosomal RNA) was the most stable in cold-stressed roots. Additionally, we tested the utility of these candidate reference genes to detect the expression pattern of P5CS (Δ¹-pyrroline-5-carboxylate synthetase), which is a drought-related gene. This study is the first report on selecting and validating reference genes of P. villosa under various stress conditions and will benefit future investigations of the genomic mechanisms of stress resistance in this ecologically important species.

Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase

Soil-borne microbes can establish compatible relationships with host plants, providing a large variety of nutritive and protective compounds in exchange for photosynthesized sugars. However, the molecular mechanisms mediating the establishment of these beneficial relationships remain unclear. Our previous genetic mapping and whole-genome resequencing studies identified a gene deletion event of a Populus trichocarpa lectin receptor-like kinase gene PtLecRLK1 in Populus deltoides that was associated with poor-root colonization by the ectomycorrhizal fungus Laccaria bicolor. By introducing PtLecRLK1 into a perennial grass known to be a non-host of L. bicolor, switchgrass (Panicum virgatum L.), we found that L. bicolor colonizes ZmUbipro-PtLecRLK1 transgenic switchgrass roots, which illustrates that the introduction of PtLecRLK1 has the potential to convert a non-host to a host of L. bicolor. Furthermore, transcriptomic and proteomic analyses on inoculated-transgenic switchgrass roots revealed genes/proteins overrepresented in the compatible interaction and underrepresented in the pathogenic defence pathway, consistent with the view that pathogenic defence response is down-regulated during compatible interaction. Metabolomic profiling revealed that root colonization in the transgenic switchgrass was associated with an increase in N-containing metabolites and a decrease in organic acids, sugars, and aromatic hydroxycinnamate conjugates, which are often seen in the early steps of establishing compatible interactions. These studies illustrate that PtLecRLK1 is able to render a plant susceptible to colonization by the ectomycorrhizal fungus L. bicolor and shed light on engineering mycorrhizal symbiosis into a non-host to enhance plant productivity and fitness on marginal lands.

Serendipita Fungi Modulate the Switchgrass Root Transcriptome to Circumvent Host Defenses and Establish a Symbiotic Relationship

The fungal family Serendipitaceae encompasses root-associated lineages with endophytic, ericoid, orchid, and ectomycorrhizal lifestyles. Switchgrass is an important bioenergy crop for cellulosic ethanol production owing to high biomass production on marginal soils otherwise unfit for food crop cultivation. The aim of this study was to investigate the host plant responses to Serendipita spp. colonization by characterizing the switchgrass root transcriptome during different stages of symbiosis in vitro. For this, we included a native switchgrass strain, Serendipita bescii, and a related strain, S. vermifera, isolated from Australian orchids. Serendipita colonization progresses from thin hyphae that grow between root cells to, finally, the production of large, bulbous hyphae that fill root cells during the later stages of colonization. We report that switchgrass seems to perceive both fungi prior to physical contact, leading to the activation of chemical and structural defense responses and putative host disease resistance genes. Subsequently, the host defense system appears to be quenched and carbohydrate metabolism adjusted, potentially to accommodate the fungal symbiont. In addition, prior to contact, switchgrass exhibited significant increases in root hair density and root surface area. Furthermore, genes involved in phytohormone metabolism such as gibberellin, jasmonic acid, and salicylic acid were activated during different stages of colonization. Both fungal strains induced plant gene expression in a similar manner, indicating a conserved plant response to members of this fungal order. Understanding plant responsiveness to Serendipita spp. will inform our efforts to integrate them into forages and row crops for optimal plant-microbe functioning, thus facilitating low-input, sustainable agricultural practices.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Influence of Switchgrass TDIF-like Genes on Arabidopsis Vascular Development

As a member of the CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION (CLE) family, the dodecapeptide tracheary element differentiation inhibitory factor (TDIF) has a major impact on vascular development in plants. However, the influence of polymorphisms in the TDIF peptide motif on activity remains poorly understood. The model plant, Arabidopsis provides a fast and effective tool for assaying the activity of TDIF homologs. Five TDIF homologs from a group of 93 CLE genes in switchgrass (Panicum virgatum), a perennial biomass crop, named PvTDIF-like (PvTDIFL) genes were studied. The expression levels of PvTDIFL1, PvTDIFL3 ^MR3, and PvTDIFL3 ^MR2 were relatively high and all of them were expressed at the highest levels in the rachis of switchgrass. The precursor proteins for PvTDIFL1, PvTDIFL3^MR3, and PvTDIFL3^MR2 contained one, three, and two TDIFL motifs, respectively. Treatments with exogenous PvTDIFL peptides increased the number of stele cells in the hypocotyls of Arabidopsis seedlings, with the exception of PvTDIFL_4p. Heterologous expression of PvTDIFL1 in Arabidopsis strongly inhibited plant growth, increased cell division in the vascular tissue of the hypocotyl, and disrupted the cellular organization of the hypocotyl. Although heterologous expression of PvTDIFL3 ^MR3 and PvTDIFL3 ^MR2 also affected plant growth and vascular development, PvTDIFL activity was not enhanced by the multiple TDIFL motifs encoded by PvTDIFL3 ^MR3 and PvTDIFL3 ^MR2. These data indicate that in general, PvTDIFLs are functionally similar to Arabidopsis TDIF but that the processing and activities of the PvTDIFL peptides are more complex.

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.

Expression of potential reference genes in response to macronutrient stress in rice and soybean

Given the complexity of nutrient stress responses and the availability of a few validated reference genes, we aimed to identify robust and stable reference genes for macronutrient stress in rice and soybean. Ten potential reference genes were evaluated using geNorm, NormFinder, BestKeeper, Comparative ΔCt method, and RefFinder algorithms under low and completely starved conditions of nitrogen (N), phosphorus (P), potassium (K), and sulphur (S). Results revealed distinct sets of reference gene pairs, showing stable expression under different experimental conditions. The gene pairs TIP41/UBC(9/10/18) and F-box/UBC10 were most stable in rice and soybean, respectively under N stress. Under P stress, UBC9/UBC10 in rice and F-Box/UBC10 in soybean were most stable. Similarly, TIP41/UBC10 in rice and RING FINGER/UBC9 in soybean were the best gene pairs under K stress while F-Box/TIP41 in rice and UBC9/UBC10 in soybean were the most stable gene pairs under S stress. These reference gene pairs were validated by quantifying the expression levels of high-affinity transporters like NRT2.1/NRT2.5, PT1, AKT1, and SULTR1 for N, P, K, and S stress, respectively. This study reiterates the importance of choosing reference genes based on crop species and the experimental conditions, in order to obtain concrete answers to missing links of gene regulation in response to macronutrient deficiencies.

De novo transcriptome in roots of switchgrass (Panicum virgatum L.) reveals gene expression dynamic and act network under alkaline salt stress

Background

Soil salinization is a major limiting factor for crop cultivation. Switchgrass is a perennial rhizomatous bunchgrass that is considered an ideal plant for marginal lands, including sites with saline soil. Here we investigated the physiological responses and transcriptome changes in the roots of Alamo (alkaline-tolerant genotype) and AM-314/MS-155 (alkaline-sensitive genotype) under alkaline salt stress.

Results

Alkaline salt stress significantly affected the membrane, osmotic adjustment and antioxidant systems in switchgrass roots, and the ASTTI values between Alamo and AM-314/MS-155 were divergent at different time points. A total of 108,319 unigenes were obtained after reassembly, including 73,636 unigenes in AM-314/MS-155 and 65,492 unigenes in Alamo. A total of 10,219 DEGs were identified, and the number of upregulated genes in Alamo was much greater than that in AM-314/MS-155 in both the early and late stages of alkaline salt stress. The DEGs in AM-314/MS-155 were mainly concentrated in the early stage, while Alamo showed greater advantages in the late stage. These DEGs were mainly enriched in plant-pathogen interactions, ubiquitin-mediated proteolysis and glycolysis/gluconeogenesis pathways. We characterized 1480 TF genes into 64 TF families, and the most abundant TF family was the C2H2 family, followed by the bZIP and bHLH families. A total of 1718 PKs were predicted, including CaMK, CDPK, MAPK and RLK. WGCNA revealed that the DEGs in the blue, brown, dark magenta and light steel blue 1 modules were associated with the physiological changes in roots of switchgrass under alkaline salt stress. The consistency between the qRT-PCR and RNA-Seq results confirmed the reliability of the RNA-seq sequencing data. A molecular regulatory network of the switchgrass response to alkaline salt stress was preliminarily constructed on the basis of transcriptional regulation and functional genes.

Conclusions

Alkaline salt tolerance of switchgrass may be achieved by the regulation of ion homeostasis, transport proteins, detoxification, heat shock proteins, dehydration and sugar metabolism. These findings provide a comprehensive analysis of gene expression dynamic and act network induced by alkaline salt stress in two switchgrass genotypes and contribute to the understanding of the alkaline salt tolerance mechanism of switchgrass and the improvement of switchgrass germplasm.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07368-w.

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Knowing how switchgrass (Panicum virgatum L.) responds and adapts to phosphorus (P)-limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P-stress on genome transcription and whole-plant metabolism, physiology and development in switchgrass. P-limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P-stress, and decreased root diameter with increased density and length of root hairs at severe P-stress. RNA-seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P-stress in roots and/or shoots compared to P-replete plants, with many stress-induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P-acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P-stress in shoots and roots, although this trend reversed under severe P-stress, especially in shoots.

Screening and verification of reference genes for analysis of gene expression in winter rapeseed (Brassica rapa L.) under abiotic stress

Winter rapeseed (Brassica rapa L.) is the main oilseed crop in northern China and can safely overwinter at 35 (i.e., Tianshui, China) to 48 degrees north latitude (i.e., Altai, Heilongjiang, Raohe, and Xinjiang, China). In order to identify stable reference genes to understand the molecular mechanisms of stress tolerance in winter rapeseed, internal reference genes of winter rapeseed under four abiotic stresses were analyzed using GeNorm, NormFinder, BestKeeper, and RefFinder software. The most stable combinations of internal reference genes were β-actin and SAND in cold-stressed leaves, β-actin and EF1a in cold-stressed roots, F-box and SAND in high temperature-stressed leaves, and PP2A and RPL in high temperature-stressed roots, SAND and PP2A in NaCl-stressed leaves, RPL and UBC in NaCl-stressed roots, RPL and PP2A in PEG-stressed leaves, and PP2A and RPL in PEG-stressed roots. Expression profiles of PXG3 were used to verify these results. The stable reference genes identified in this study are useful tools for identifying stress-responsive genes to understand the molecular mechanisms of stress tolerance in winter rapeseed.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR in White Clover (Trifolium repens L.) Involved in Five Abiotic Stresses

White clover (Trifolium repens L.) is a widely cultivated cool-season perennial forage legume in temperate grassland systems. Many studies have analyzed the gene expression in this grass species using quantitative real-time reverse transcription PCR (qRT-PCR). The selection of stable reference genes for qRT-PCR is crucial. However, there was no detailed study on reference genes in different tissues of white clover under various abiotic stress conditions. Herein, 14 candidate reference genes (ACT7, ACT101, TUA1109, TUB, CYP, 60SrRNA, UBQ, E3, GAPDH1, GAPDH2, PP2A, BAM3, SAMDC, and ABC) were selected and analyzed by four programs (GeNorm, NormFinder, BestKeeper, and RefFinder). Samples were taken from two tissues (leaves and roots) under five different abiotic stresses (drought, salt, heat, cold, and heavy metal stress). Our results showed that 60SrRNA and ACT101 were the two top-ranked genes for all samples. Under various experimental conditions, the most stable gene was different; however, SAMDC, UBQ, 60SrRNA, and ACT101 were always top ranked. The most suitable reference genes should be selected according to different plant tissues and growth conditions. Validation of these reference genes by expression analysis of Cyt-Cu/Zn SOD and CAT confirmed their reliability. Our study will benefit the subsequent research of gene function in this species.

Identification and Validation of Reference Genes for RT-qPCR Analysis in Switchgrass under Heavy Metal Stresses

Switchgrass (Panicum Virgatum L.) has been recognized as the new energy plant, which makes it ideal for the development of phytoremediation on heavy metal contamination in soils with great potential. This study aimed to screen the best internal reference genes for the real-time quantitative PCR (RT-qPCR) in leaves and roots of switchgrass for investigating its response to various heavy metals, such as cadmium (Cd), lead (Pb), mercury (Hg), chromium (Cr), and arsenic (As). The stability of fourteen candidate reference genes was evaluated by BestKeeper, GeNorm, NormFinder, and RefFinder software. Our results identified U2AF as the best reference gene in Cd, Hg, Cr, and As treated leaves as well as in Hg, Pb, As, and Cr stressed root tissues. In Pb treated leaf tissues, 18S rRNA was demonstrated to be the best reference gene. CYP5 was determined to be the optimal reference gene in Cd treated root tissues. The least stable reference gene was identified to be CYP2 in all tested samples except for root tissues stressed by Pb. To further validate the initial screening results, we used the different sets of combinatory internal reference genes to analyze the expression of two metal transport associated genes (PvZIP4 and PvPDB8) in young leaves and roots of switchgrass. Our results demonstrated that the relative expression of the target genes consistently changed during the treatment when CYP5/UBQ1, U2AF/ACT12, eEF1a/U2AF, or 18S rRNA/ACT12 were combined as the internal reference genes. However, the time-dependent change pattern of the target genes was significantly altered when CYP2 was used as the internal reference gene. Therefore, the selection of the internal reference genes appropriate for specific experimental conditions is critical to ensure the accuracy and reliability of RT-qPCR. Our findings established a solid foundation to further study the gene regulatory network of switchgrass in response to heavy metal stress.

Combined Proteomics and Metabolism Analysis Unravels Prominent Roles of Antioxidant System in the Prevention of Alfalfa (Medicago sativa L.) against Salt Stress

Alfalfa is the most extensively cultivated forage legume worldwide, and salinity constitutes the main environmental scourge limiting its growth and productivity. To unravel the potential molecular mechanism involved in salt tolerance in alfalfa, we accomplished a combined analysis of parallel reaction monitoring-based proteomic technique and targeted metabolism. Based on proteomic analysis, salt stress induced 226 differentially abundant proteins (DAPs). Among them, 118 DAPs related to the antioxidant system, including glutathione metabolism and oxidation-reduction pathways, were significantly up-regulated. Data are available via ProteomeXchange with identifier PXD017166. Overall, 107 determined metabolites revealed that the tricarboxylic acid (TCA) cycle, especially the malate to oxaloacetate conversion step, was strongly stimulated by salt stress. This leads to an up-regulation by about 5 times the ratio of NADPH/NADP+, as well as about 3 to 5 times in the antioxidant enzymes activities, including those of catalase and peroxidase and proline contents. However, the expression levels of DAPs related to the Calvin-Benson-Bassham (CBB) cycle and photorespiration pathway were dramatically inhibited following salt treatment. Consistently, metabolic analysis showed that the metabolite amounts related to carbon assimilation and photorespiration decreased by about 40% after exposure to 200 mM NaCl for 14 d, leading ultimately to a reduction in net photosynthesis by around 30%. Our findings highlighted also the importance of the supplied extra reducing power, thanks to the TCA cycle, in the well-functioning of glutathione to remove and scavenge the reactive oxygen species (ROS) and mitigate subsequently the oxidative deleterious effect of salt on carbon metabolism including the CBB cycle.

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

Relative gene expression analyses by RT-qPCR (reverse transcription-quantitative PCR) are highly dependent on the reference genes in normalizing the expression data of target genes. Therefore, inappropriate endogenous control genes will lead to inaccurate target gene expression profiles, and the selection and validation of suitable internal reference genes becomes essential. In this study, we retrieved the commonly used reference genes in transcriptome datasets of Codonopsis pilosula by RNA-Seq (unpublished data), and selected 15 candidate reference genes according to the coefficient of variation (CV) and fold change (FC) value of gene expression. The expression levels of candidate reference genes, which is at different growth stages, undergoing cold stress and drought stress, was determined by RT-qPCR. The expression stability of these genes was evaluated using software packages and algorithms including ΔCt, geNorm, NormFinder and Bestkeeper. Then appropriate reference genes were screened and validated by target gene-UDGPase (UDP glucose pyrophosphorylase). The optimal RGs combinations of C. pilosula, including PP2A59γ, CPY20-1, UBCE32, RPL5B and UBC18 for developmental stage, RPL5B, RPL13 and PP2A59γ for cold treatment, RPL13 and PP2A59γ for drought treatment, were found and proposed as reference genes for future work. This paper laid foundations for both the selection of reference genes and exploration in metabolic mechanism of C. pilosula.

Comparative transcriptome analysis provides key insights into seedling development in switchgrass (Panicum virgatum L.)

BACKGROUND

Switchgrass (Panicum virgatum L.), a warm-season perennial C4 plant, can be used as a forage plant, a soil and water conservation plant, a windbreak plant, and as a good source of biofuels and alternative energy with low planting costs. However, switchgrass exhibits low rates of seedling development compared to other crops, which means it is typically out-competed by weeds. There is a large variation in seedling development rates among different plantlets in switchgrass, which limits its usefulness for large-scale cultivation. Little is currently known about the molecular reasons for slow seedling growth.

RESULTS

Characterization of the seedling development process via growth indices indicated a relatively stagnant growth stage in switchgrass. A total of 678 differentially expressed genes (DEGs) were identified from the comparison of transcriptomes from slowly developed (sd) and rapidly developed (rd) switchgrass seedlings. Gene ontology and pathway enrichment analysis showed that DEGs were enriched in diterpenoid biosynthesis, thiamine metabolism, and circadian rhythm. Transcription factor enrichment and expression analyses showed MYB-related, bHLH and NAC family genes were essential for seedling growth. The transcriptome results were consistent with those of quantitative real-time polymerase chain reaction. Then, the expression profiles of maize and switchgrass were compared during seedling leaf development. A total of 128 DEGs that play key roles in seedling growth were aligned to maize genes. Transcriptional information and physiological indices suggested that several genes involved in the circadian rhythm, thiamine metabolism, energy metabolism, gibberellic acid biosynthesis, and signal transduction played important roles in seedling development.

CONCLUSIONS

The seedling development process of switchgrass was characterized, and the molecular differences between slowly developed and rapidly developed seedlings were discussed. This study provides new insights into the reasons for slow seedling development in switchgrass and will be useful for the genetic improvement of switchgrass and other crops.

Overexpression of the Lolium perenne L. delta1-pyrroline 5-carboxylate synthase (LpP5CS) gene results in morphological alterations and salinity tolerance in switchgrass (Panicum virgatum L.)

In plants, Δ1-pyrroline- 5-carboxylate synthase (P5CS) is the rate-limiting enzyme in proline biosynthesis. In this study, we introduced the LpP5CS (Lolium perenne L.) gene into switchgrass by Agrobacterium-mediated transformation. The transgenic lines (TG) were classified into two groups based on their phenotypes and proline levels. The group I lines (TG4 and TG6) had relatively high proline levels and improved biomass yield. The group II lines (TG1 and TG2) showed low proline levels, severely delayed flowering, stunted growth and reduced biomass yield. Additionally, we used RNA-seq analysis to detect the most significant molecular changes, and we analyzed differentially expressed genes, such as flowering-related and CYP450 family genes. Moreover, the biomass yield, physiological parameters, and expression levels of reactive oxygen species scavenger-related genes under salt stress all indicated that the group I plants exhibited significantly increased salt tolerance compared with that of the control plants, in contrast to the group II plants. Thus, genetic improvement of switchgrass by overexpressing LpP5CS to increase proline levels is feasible for increasing plant stress tolerance.

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.

Validation of Reference Genes for Gene Expression Normalization in RAW264.7 Cells under Different Conditions

RAW264.7 is a macrophage strain derived from mice tumour and shows a significant ability in antigen uptake. Real-time quantitative PCR (RT-qPCR) is one of the most commonly used methods in gene studies and requires suitable reference genes to normalize and quantitate the expression of gene of interest with sensitivity and specificity. However, suitable reference genes in RAW264.7 cells have not yet been identified for accurate gene expression quantification. In the current study, we evaluated expression levels of ten candidate reference genes in RAW264.7 cells under different conditions. RT-qPCR results indicated significant differences in the expression levels among the ten reference genes. Statistical analyses were carried out using geNorm, NormFinder, and BestKeeper software to further investigate the stability of the reference genes. Integrating the results from the three analytical methods, cytochrome c-1 and hydroxymethylbilane synthase were found to be the most stable and therefore more suitable reference genes, while ribosomal protein L4 and cyclophilin A were the least stable. This study emphasises the importance of identifying and selecting the most stable reference genes for normalization and provides a basis for future gene expression studies using RAW264.7 cells.

Selection of suitable reference genes for quantitive real-time PCR normalization in Miscanthus lutarioriparia

Miscanthus lutarioriparia, which is found widespread in China, has attracted great attention as a most potential bioenergy plant for years. The quantitative real time PCR (RT-qPCR) has appeared as a sensitive and powerful technique to measure gene expression in living organisms during different development stages. In this study, we evaluated ten candidate genes, including 25S ribosomal RNA gene (25S rRNA), actin1 gene (ACT1), carotenoid-binding protein 20 gene (CBP20), glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), Ubiquitin gene (UBQ), eukaryotic elongation factor 1-αgene (eEF-1α), α-tubulin gene (α-TUB), β-tubulin gene (β-TUB), eukaryotic translation initiation factor 4α-1 gene (eIF-4α) and NAC domain protein gene(NAC) in a series of 30 M. lutarioriparia samples followed by statistical algorithms geNorm and Normfinder to analyze the gene expression stability. The results indicated that eIF-4αand UBQ were the most stable expressed genes while CBP20 showed as the least stable among all the samples. Based on above research, we recommend that at least two top-ranked reference genes should be employed for expression data normalization. The best genes selected in this study will provide a starting point to select reference genes in the future in other tissues and under other experimental conditions in this energy crop candidate.

Reference gene selection for real-time quantitative PCR normalization in Hemarthria compressa and Hemarthria altissima leaf tissue

Hemarthria compressa and Hemarthria altissima are widely used as livestock forage and play important roles in tropical and subtropical grassland agricultural systems promoting healthy ecological environment and the development of animal husbandry. Leaf tissue of "Yaan" limpograss (H. compressa) and "H255" whip grass (H. altissima) were used to test the mRNA expression levels of 12 reference genes using RT-qPCR. The Delta-Ct method, BestKeeper (ver. 1.0), Genorm (ver. 3.5), Normfinder (ver. 0.953) and RefFinder were used to analyze the expression stability of the 12 reference genes under drought, salt, acid-aluminum and cold stresses to provide significant technical support for the study of gene expression under various abiotic stresses in Hemarthria. The results showed that the candidate reference genes showed divergent expression levels under various abiotic stresses. Among the genes that were selected, CL18892 showed the highest expression stability under salt stress in the leaf tissue. eEF-1α was the most stable gene under cold and acid-aluminum stresses and CL16384 was comparatively the most suitable genes under drought stress. As a whole, according to RefFinder analysis, CYP5, BMK.74327 and CL21527 were the most suitable reference genes for studying the effects of abiotic stress in Hemarthria. In general, CL16812 and CL18038 were not suitable reference genes under abiotic stress conditions that were examined in this study.

Selection and validation of reference genes for quantitative expression analysis of miRNAs and mRNAs in Poplar

BACKGROUND

Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) is a rapid and sensitive approach to identify miRNA and protein-coding gene expression in plants. However, because of the specially designated reverse transcription and shorter PCR products, very few reference genes have been identified for the quantitative analysis of miRNA expression in plants, and different internal reference genes are needed to normalize the expression of miRNAs and mRNA genes respectively. Therefore, it is particularly important to select the suitable common reference genes for normalization of quantitative PCR of miRNA and mRNA.

RESULTS

In this study, a modified reverse transcription PCR protocol was adopted for selecting and validating universal internal reference genes of mRNAs and miRNAs. Eight commonly used reference genes, four stably expressed novel genes in Populus tremula, three small noncoding RNAs and three conserved miRNAs were selected as candidate genes, and the stability of their expression was examined across a set of 38 tissue samples from four developmental stages of poplar clone 84K (Populus alba × Populus glandulosa). The expression stability of these candidate genes was evaluated systematically by four algorithms: geNorm, NormFinder, Bestkeeper and DeltaCt. The results showed that Eukaryotic initiation factor 4A III (EIF4A) and U6-2 were suitable for samples of the callus stage; U6-1 and U6-2 were best for the seedling stage; Protein phosphatase 2A-2 (PP2A-2) and U6-1 were best for the plant stage; and Protein phosphatase 2A-2 (PP2A-2) and Oligouridylate binding protein 1B (UBP) were the best reference genes in the adventitious root (AR) regeneration stage.

CONCLUSIONS

The purpose of this study was to identify the most appropriate reference genes for qRT-PCR of miRNAs and mRNAs in different tissues at several developmental stages in poplar. U6-1, EIF4A and PP2A-2 were the three most appropriate reference genes for qRT-PCR normalization of miRNAs and mRNAs during the plant regeneration process, and PP2A-2 and UBP represent the best reference genes in the AR regeneration stage of poplar. This work will benefit future studies of expression and function analysis of miRNAs and their target genes in poplar.

Evaluation and validation of housekeeping genes in two contrast species of thyme plant to drought stress using real-time PCR

To decrease errors and increase accuracy and reliability of quantitative real-time PCR (qRT-PCR) results, the use of a reference gene is inevitable. Despite the industrial importance of genus Thymus, not any validated reference gene has not been reported for T. kotschyanus and T. vulgaris which could limit such investigations. In this study, the expression stability of seven housekeeping genes including Actin, Cyclophilin-18, elongation factor-1A, glyceraldehyde-3-phosphate dehydrogenase, 18S ribosomal RNA, Cullin, and Polypyrimidine tract-binding protein were evaluated in T. kotschyanus and T. vulgaris which grown at four levels of drought stress using geNorm, NormFinder, and BestKeeper algorithms. Histone deacetylase-6 (HDA-6) gene was also used for validation of evaluated reference genes. In T. vulgaris, all of the algorithms similarly ranked elongation factor-1A and glyceraldehyde-3-phosphate dehydrogenase as the two most stably expressed genes. In T. kotschyanus, only NormFinder and BestKeeper had a similar ranking and identified Actin and glyceraldehyde-3-phosphate dehydrogenase as the two most stably expressed genes, but geNorm algorithm ranked elongation factor-1A and glyceraldehyde-3-phosphate dehydrogenase as the best two reference genes. On the other hand, all algorithms ranked 18S rRNA and Cyclophilin-18 as the least stable genes in T. kotschyanus and T. vulgaris, respectively. Validation results indicated that there was a significant change (0.53-3.19 fold change) in relative expression of HDA-6 normalized by the best stable gene compare to the least ranked gene. Our study presented the first systematic validation of reference gene(s) selection in T. kotschyanus and T. vulgaris and provided useful information to obtain more accurate qRT-PCR results in these species.

Genome-Wide Association Study in Pseudo-F2 Populations of Switchgrass Identifies Genetic Loci Affecting Heading and Anthesis Dates

Switchgrass (Panicum virgatum) is a native prairie grass and valuable bio-energy crop. The physiological change from juvenile to reproductive adult can draw important resources away from growth into producing reproductive structures, thereby limiting the growth potential of early flowering plants. Delaying the flowering of switchgrass is one approach by which to increase total biomass. The objective of this research was to identify genetic variants and candidate genes for controlling heading and anthesis in segregating switchgrass populations. Four pseudo-F₂ populations (two pairs of reciprocal crosses) were developed from lowland (late flowering) and upland (early flowering) ecotypes, and heading and anthesis dates of these populations were collected in Lafayette, IN and DeKalb, IL in 2015 and 2016. Across 2 years, there was a 34- and 73-day difference in heading and a 52- and 75-day difference in anthesis at the Lafayette and DeKalb locations, respectively. A total of 37,901 single nucleotide polymorphisms obtained by exome capture sequencing of the populations were used in a genome-wide association study (GWAS) that identified five significant signals at three loci for heading and two loci for anthesis. Among them, a homolog of FLOWERING LOCUS T on chromosome 5b associated with heading date was identified at the Lafayette location across 2 years. A homolog of ARABIDOPSIS PSEUDO-RESPONSE REGULATOR 5, a light modulator in the circadian clock associated with heading date was detected on chromosome 8a across locations and years. These results demonstrate that genetic variants related to floral development could lend themselves to a long-term goal of developing late flowering varieties of switchgrass with high biomass yield.

Functional Diversity of Diterpene Synthases in the Biofuel Crop Switchgrass

Diterpenoids constitute a diverse class of metabolites with critical functions in plant development, defense, and ecological adaptation. Major monocot crops, such as maize (Zea mays) and rice (Oryza sativa), deploy diverse blends of specialized diterpenoids as core components of biotic and abiotic stress resilience. Here, we describe the genome-wide identification and functional characterization of stress-related diterpene synthases (diTPSs) in the dedicated bioenergy crop switchgrass (Panicum virgatum). Mining of the allotetraploid switchgrass genome identified an expansive diTPS family of 31 members, and biochemical analysis of 11 diTPSs revealed a modular metabolic network producing a diverse array of diterpenoid metabolites. In addition to ent-copalyl diphosphate (CPP) and ent-kaurene synthases predictably involved in gibberellin biosynthesis, we identified syn-CPP and ent-labda-13-en-8-ol diphosphate (LPP) synthases as well as two diTPSs forming (+)-labda-8,13E-dienyl diphosphate (8,13-CPP) and ent-neo-cis-trans-clerodienyl diphosphate (CT-CLPP) scaffolds not observed previously in plants. Structure-guided mutagenesis of the (+)-8,13-CPP and ent-neo-CT-CLPP synthases revealed residue substitutions in the active sites that altered product outcome, representing potential neofunctionalization events that occurred during diversification of the switchgrass diTPS family. The conversion of ent-CPP, ent-LPP, syn-CPP, and ent-neo-CT-CLPP by promiscuous diTPSs further yielded distinct labdane-type diterpene olefins and alcohols. Of these metabolites, the formation of 9β-hydroxy-syn-pimar-15-ene and the expression of the corresponding genes were induced in roots and leaves in response to oxidative stress and ultraviolet irradiation, indicating their possible roles in abiotic stress adaptation. Together, these findings expand the known chemical space of diterpenoid metabolism in monocot crops toward systematically investigating and ultimately improving stress resilience traits in crop species.

Expression of endogenous proteins in maize hybrids in a multi-location field trial in India

Genetically modified (GM) crops undergo large scale multi-location field trials to characterize agronomics, composition, and the concentration of newly expressed protein(s) [herein referred to as transgenic protein(s)]. The concentration of transgenic proteins in different plant tissues and across the developmental stages of the plant is considered in the safety assessment of GM crops. Reference or housekeeping proteins are expected to maintain a relatively stable expression pattern in healthy plants given their role in cellular functions. Understanding the effects of genotype, growth stage and location on the concentration of endogenous housekeeping proteins may provide insight into the contribution these factors could have on transgenic protein concentrations in GM crops. The concentrations of three endogenous proteins (actin, elongation factor 1-alpha, and glyceraldehyde 3-phosphate dehydrogenase) were measured in several different maize hybrids grown across multiple field locations over 2 years. Leaf samples were collected from healthy plants at three developmental stages across the growing seasons, and protein concentrations were quantified by indirect enzyme-linked immunosorbent assay (ELISA) for each protein. In general, the concentrations of these three endogenous proteins were relatively consistent across hybrid backgrounds, when compared within one growth stage and location (2-26%CV), whereas the concentrations of proteins in the same hybrid and growth stage across different locations were more variable (12-64%CV). In general, the protein concentrations in 2013 and 2014 show similar trends in variability. Some degree of variability in protein concentrations should be expected for both transgenic and endogenous plant-expressed proteins. In the case of GM crops, the potential variation in protein concentrations due to location effects is captured in the current model of multi-location field testing.

Long non-coding RNAs of switchgrass (Panicum virgatum L.) in multiple dehydration stresses

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in plant growth and stress responses. Studies of lncRNAs in non-model plants are quite limited, especially those investigating multiple dehydration stresses. In this study, we identified novel lncRNAs and analyzed their functions in dehydration stress memory in switchgrass, an excellent biofuel feedstock and soil-conserving plant in the Gramineae family.

RESULTS

We analyzed genome-wide transcriptional profiles of leaves of 5-week-old switchgrass plantlets grown via tissue culture after primary and secondary dehydration stresses (D1 and D2) and identified 16,551 novel lncRNAs, including 4554 annotated lncRNAs (targeting 3574 genes), and 11,997 unknown lncRNAs. Gene ontology and pathway enrichment analysis of annotated genes showed that the differentially expressed lncRNAs were related to abscisic acid (ABA) and ethylene (ETH) biosynthesis and signal transduction, and to starch and sucrose metabolism. The upregulated lncRNAs and genes were related to ABA synthesis and its signal transduction, and to trehalose synthesis. Meanwhile, lncRNAs and genes related to ETH biosynthesis and signal transduction were suppressed. LncRNAs and genes involved in ABA metabolism were verified using quantitative real-time PCR, and the endogenous ABA content was determined via high performance liquid chromatography mass spectrometry (HPLC-MS). These results showed that ABA accumulated significantly during dehydration stress, especially in D2. Furthermore, we identified 307 dehydration stress memory lncRNAs, and the ratios of different memory types in switchgrass were similar to those in Arabidopsis and maize.

CONCLUSIONS

The molecular responses of switchgrass lncRNAs to multiple dehydration stresses were researched systematically, revealing novel information about their transcriptional regulatory behavior. This study provides new insights into the response mechanism to dehydration stress in plants. The lncRNAs and pathways identified in this study provide valuable information for genetic modification of switchgrass and other crops.

Selection and evaluation of housekeeping genes as endogenous controls for quantification of mRNA transcripts in Theileria parva using quantitative real-time polymerase chain reaction (qPCR)

The reliability of any quantitative real-time polymerase chain reaction (qPCR) experiment can be seriously compromised by variations between samples as well as between PCR runs. This usually result from errors in sample quantification, especially with samples that are obtained from different individuals and tissues and have been collected at various time intervals. Errors also arise from differences in qPCR efficiency between assays performed simultaneously to target multiple genes on the same plate. Consequently, the derived quantitative data for the target genes become distorted. To avoid this grievous error, an endogenous control, with relatively constant transcription levels in the target individual or tissue, is included in the qPCR assay to normalize target gene expression levels in the analysis. Several housekeeping genes (HKGs) have been used as endogenous controls in quantification studies of mRNA transcripts; however, there is no record in the literature of the evaluation of these genes for the tick-borne protozoan parasite, Theileria parva. Importantly, the expression of these genes should be invariable between different T. parva stocks, ideally under different experimental conditions, to gain extensive application in gene expression studies of this parasite. Thus, the expression of several widely used HKGs was evaluated in this study, including the genes encoding β-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 28S rRNA, cytochrome b and fructose-2.6-biphosphate aldolase (F6P) proteins. The qPCR analysis revealed that the expression of genes encoding cytochrome b, F6P and GAPDH varied considerably between the two T. parva stocks investigated, the cattle-derived T. parva Muguga and the buffalo-derived T. parva 7014. 28S rRNA and β-actin gene expression was the most stable; thus, these genes were considered suitable candidates to be used as endogenous control genes for mRNA quantification studies in T. parva.

Reference gene selection for qRT-PCR assays in Stellera chamaejasme subjected to abiotic stresses and hormone treatments based on transcriptome datasets

Background

Stellera chamaejasme Linn, an important poisonous plant of the China grassland, is toxic to humans and livestock. The rapid expansion of S. chamaejasme has greatly damaged the grassland ecology and, consequently, seriously endangered the development of animal husbandry. To draft efficient prevention and control measures, it has become more urgent to carry out research on its adaptive and expansion mechanisms in different unfavorable habitats at the genetic level. Quantitative real-time polymerase chain reaction (qRT-PCR) is a widely used technique for studying gene expression at the transcript level; however, qRT-PCR requires reference genes (RGs) as endogenous controls for data normalization and only through appropriate RG selection and qRT-PCR can we guarantee the reliability and robustness of expression studies and RNA-seq data analysis. Unfortunately, little research on the selection of RGs for gene expression data normalization in S. chamaejasme has been reported.

Method

In this study, 10 candidate RGs namely, 18S, 60S, CYP, GAPCP1, GAPDH2, EF1B, MDH, SAND, TUA1, and TUA6, were singled out from the transcriptome database of S. chamaejasme, and their expression stability under three abiotic stresses (drought, cold, and salt) and three hormone treatments (abscisic acid, ABA; gibberellin, GA; ethephon, ETH) were estimated with the programs geNorm, NormFinder, and BestKeeper.

Result

Our results showed that GAPCP1 and EF1B were the best combination for the three abiotic stresses, whereas TUA6 and SAND, TUA1 and CYP, GAPDH2 and 60S were the best choices for ABA, GA, and ETH treatment, respectively. Moreover, GAPCP1 and 60S were assessed to be the best combination for all samples, and 18S was the least stable RG for use as an internal control in all of the experimental subsets. The expression patterns of two target genes (P5CS2 and GI) further verified that the RGs that we selected were suitable for gene expression normalization.

Discussion

This work is the first attempt to comprehensively estimate the stability of RGs in S. chamaejasme. Our results provide suitable RGs for high-precision normalization in qRT-PCR analysis, thereby making it more convenient to analyze gene expression under these experimental conditions.

Reference gene selection for quantitative real-time RT-PCR normalization in Clarias magur at different larval developmental stages

Reference genes employed for normalizing quantitative PCR data are important for the accurate analysis of gene expression. To date, no reference genes have been screened for developmental gene expression studies in Clarias magur. In the present study, three commonly used and constitutively expressed genes viz. beta actin (β- actin), glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and elongation factor-alpha 1 (EFa1) were examined for their efficacy as internal control to avoid any variation during qRT-PCR expression analysis at different developmental stages of C. magur. All the selected housekeeping genes showed a variable level of mRNA expression during the developmental stages of C. magur. Using three independent statistical algorithms (delta-CT, BestKeeper and NormFinder), β-actin and GAPDH were identified as the suitable genes at different developmental stages. However, comprehensive gene stability evaluation denoted β-actin to be the most stable gene for carrying any gene expression studies. The present results, recommend β-actin as the optimal housekeeping gene for qRT-PCR analysis during different developmental stages of C. magur.

Identification and validation of reference genes for qRT-PCR analysis in mulberry (Morus alba L.)

Mulberry (Morus alba L.) is an important economic tree species in many countries. Quantitative real time PCR (qRT-PCR) has become a widely used method for gene expression studies in plants. A suitable reference gene is essential to ensure accurate and reliable results for qRT-PCR analyses. However, no reports describing the selection of reference genes have been published for mulberry. In this work, we evaluated the stability of twenty candidate reference genes in different plant tissues and under different stress conditions by qRT-PCR in mulberry using algorithms in two programs—geNorm and NormFinder. The results revealed that TUB2, UBI4, ACTIN3 and RPL4 were ranked as the most stable reference genes in the samples subsets, whereas EF1α4 and TUB3showed the least stability with both algorithms. To further validate the stability of the reference genes, the expression patterns of six genes of mulberry were analyzed by normalization with the selected reference genes. Our study will benefit future analyses of gene expression in mulberry.

Selection and validation of reference genes for quantitative real-time PCR in Artemisia sphaerocephala based on transcriptome sequence data

Artemisia sphaerocephala, a dicotyledonous perennial semi-shrub belonging to the Artemisia genus of the Compositae family, is widely distributed in northwestern China. This shrub is one of the most important pioneer plants which is capable of protecting rangelands from wind erosion. It therefore plays a vital role in maintaining desert ecosystem stability. In addition, to its use as a forage grass, it has excellent prospective applications as a source of plant oil and as a plant-based fuel. The use of internal genes is the basis for accurately assessing Real time quantitative PCR. In this study, based on transcriptome data of A. sphaerocephala, we analyzed 21 candidate internal genes to determine the optimal internal genes in this shrub. The stabilities of candidate genes were evaluated in 16 samples of A. sphaerocephala. Finally, UBC9 and TIP41-like were determined as the optimal reference genes in A. sphaerocephala by Delta Ct and three various programs. There were GeNorm, NormFinder and BestKeeper.

Validation of house-keeping genes for normalization of gene expression data during diurnal/circadian studies in rice by RT-qPCR

The circadian clock in plants is the intrinsic rhythmic expression of thousands of genes in a 24 h period, which is set by the day-night cycles in the environment. The study of the circadian clock often requires expression profiling of genes over a large number of samples for which RT-qPCR is invariably used. Reliability of the results depends largely on the house-keeping genes, which serve as control and thus should be chosen carefully to prevent erroneous results. In this study, ten house-keeping genes were chosen from rice for stability analysis with 48 tissue samples harvested from plants subjected to diurnal/circadian cycles. Although, all the genes were found to be stable, however, six of them showed cyclic expression patterns and caused major changes in the expression profiles of the target genes when used to normalize their expression data, thereby making them poor candidates for diurnal/circadian studies. In conclusion, reference genes need to be selected for diurnal/circadian studies with extra caution as more than 80% of transcriptome in plants undergoes cycling, which remains undetected by the gene stability assessment software and can severely affect the RT-qPCR based gene expression profiling. The geometric mean of two or more most stable reference genes is hence recommended for diurnal/circadian studies in plants.

Selection and evaluation of reference genes for qRT-PCR analysis in Euscaphis konishii Hayata based on transcriptome data

BACKGROUND

Quantitative real-time reverse transcription-polymerase chain reaction has been widely used in gene expression analysis, however, to have reliable and accurate results, reference genes are necessary to normalize gene expression under different experimental conditions. Several reliable reference genes have been reported in plants of Traditional Chinese Medicine, but none have been identified for Euscaphis konishii Hayata.

RESULTS

In this study, 12 candidate reference genes, including 3 common housekeeping genes and 9 novel genes based on E. konishii Hayata transcriptome data were selected and analyzed in different tissues (root, branch, leaf, capsule and seed), capsule and seed development stages. Expression stability was calculated using geNorm and NormFinder, the minimal number of reference genes required for accurate normalization was calculated by Vn/Vn + 1 using geNorm. EkEEF-5A-1 and EkADF2 were the two most stable reference genes for all samples, while EkGSTU1 and EkGAPDH were the most stable reference genes for tissue samples. For seed development stages, EkGAPDH and EkEEF-5A-1 were the most stable genes, whereas EkGSTU1 and EkGAPDH were identified as the two most stable genes in the capsule development stages. Two reference genes were sufficient to normalize gene expression across all sample sets.

CONCLUSION

Results of this study revealed that suitable reference genes should be selected for different experimental samples, and not all the common reference genes are suitable for different tissue samples and/or experimental conditions. In this study, we present the first data of reference genes selection for E. konishii Hayata based on transcriptome data, our data will facilitate further studies in molecular biology and gene function on E. konishii Hayata and other closely related species.

Transcriptional and physiological data reveal the dehydration memory behavior in switchgrass (Panicum virgatum L.)

BACKGROUND

Switchgrass (Panicum virgatum L.) is a model biofuel plant because of its high biomass, cellulose-richness, easy degradation to ethanol, and the availability of extensive genomic information. However, a little is currently known about the molecular responses of switchgrass plants to dehydration stress, especially multiple dehydration stresses.

RESULTS

Studies on the transcriptional profiles of 35-day-old tissue culture plants revealed 741 dehydration memory genes. Gene Ontology and pathway analysis showed that these genes were enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant hormone signal transduction. Further analysis of specific pathways combined with physiological data suggested that switchgrass improved its dehydration resistance by changing various aspects of its responses to secondary dehydration stress (D2), including the regulation of abscisic acid (ABA) and jasmonic acid (JA) biosynthesis and signal transduction, the biosynthesis of osmolytes (l-proline, stachyose and trehalose), energy metabolism (i.e., metabolic process relating to photosynthetic systems, glycolysis, and the TCA cycle), and lignin biosynthesis. The transcriptional data and chemical substance assays showed that ABA was significantly accumulated during both primary (D1) and secondary (D2) dehydration stresses, whereas JA accumulated during D1 but became significantly less abundant during D2. This suggests the existence of a complicated signaling network of plant hormones in response to repeated dehydration stresses. A homology analysis focusing on switchgrass, maize, and Arabidopsis revealed the conservation and species-specific distribution of dehydration memory genes.

CONCLUSIONS

The molecular responses of switchgrass plants to successive dehydration stresses have been systematically characterized, revealing a previously unknown transcriptional memory behavior. These results provide new insights into the mechanisms of dehydration stress responses in plants. The genes and pathways identified in this study will be useful for the genetic improvement of switchgrass and other crops.

Using Transcriptome Analysis to Identify Genes Involved in Switchgrass Flower Reversion

Floral reversion is a process in which differentiated floral organs revert back to vegetative organs. Although this phenomenon has been described for decades, the underlying molecular mechanisms remain unclear. In this study, we found that immature switchgrass (Panicum virgatum) inflorescences can revert to neonatal shoots when incubated on a basal medium with benzylaminopurine. We used anatomical and histological methods to verify that these shoots were formed from floret primordia through flower reversion. To further explore the gene regulation of floral reversion in switchgrass, the transcriptome of reversed, unreversed, and uncultured immature inflorescences were analyzed and 517 genes were identified as participating in flower reversion. Annotation using non-redundant databases revealed that these genes are involved in plant hormone biosynthesis and signal transduction, starch and sucrose metabolism, DNA replication and modification, and other processes crucial for switchgrass flower reversion. When four of the genes were overexpressed in Arabidopsis thaliana, vegetative growth was facilitated and reproductive growth was inhibited in transgenic plants. This study provides a basic understanding of genes regulating the floral transition in switchgrass and will promote the research of floral reversion and flower maintenance.

Selection and validation of appropriate reference genes for quantitative real-time PCR analysis in Salvia hispanica

Quantitative real-time polymerase chain reaction (qRT-PCR) has become the most popular choice for gene expression studies. For accurate expression analysis, it is pertinent to select a stable reference gene to normalize the data. It is now known that the expression of internal reference genes varies considerably during developmental stages and under different experimental conditions. For Salvia hispanica, an economically important oilseed crop, there are no reports of stable reference genes till date. In this study, we chose 13 candidate reference genes viz. Actin11 (ACT), Elongation factor 1-alpha (EF1-α), Eukaryotic translation initiation factor 3E (ETIF3E), alpha tubulin (α-TUB), beta tubulin (β-TUB), Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), Cyclophilin (CYP), Clathrin adaptor complex (CAC), Serine/threonine-protein phosphatase 2A (PP2A), FtsH protease (FtsH), 18S ribosomal RNA (18S rRNA), S-adenosyl methionine decarboxylase (SAMDC) and Rubisco activase (RCA) and the expression levels of these genes were assessed in a diverse set of tissue samples representing vegetative stages, reproductive stages and various abiotic stress treatments. Two of the widely used softwares, geNorm and Normfinder were used to evaluate the expression stabilities of these 13 candidate reference genes under different conditions. Results showed that GAPDH and CYP expression remain stable throughout in the different abiotic stress treatments, CAC and PP2A expression were relatively stable under reproductive stages and α-TUB, PP2A and ETIF3E were found to be stably expressed in vegetative stages. Further, the expression levels of Diacylglycerol acyltransferase (DGAT1), a key enzyme in triacylglycerol synthesis was analyzed to confirm the validity of reference genes identified in the study. This is the first systematic study of selection of reference genes in S. hispanica, and will benefit future expression studies in this crop.

Switchgrass PvDREB1C plays opposite roles in plant cold and salt tolerance in transgenic tobacco

BACKGROUND

The C-repeat-binding factors/DRE-binding factors (CBF/DREBs) comprise a key transcription factor family involved in plant stress tolerance. Yet, there is limited information about switchgrass DREB genes and their functional roles.

RESULTS

In this study, four cold-inducible PvDREB1s were identified from switchgrass (Panicum virgatum), among which PvDREB1C was the one responded to cold stress later than the other three PvDREB1s. Yet, ectopic overexpression of PvDREB1C led to significantly compromised, instead of improved cold tolerance in transgenic tobacco. On the other hand, PvDREB1C was transcriptionally down-regulated in response to salt stress, but overexpression of PvDREB1C improved plant salt tolerance in transgenic tobacco. The improved salt tolerance was associated with increased K+/Na+ ratio and Ca2+ content, higher cellular osmotic potential, and activation of stress-related functional genes in the leaves of transgenic plants under salt stress.

CONCLUSIONS

The current results implied that PvDREB1C played opposite roles in plant cold and salt tolerance. Although DREB1s were known as positive stress regulators, particular attentions shall be paid to their potential negative regulatory role(s).

Identification of Endogenous Reference Genes for RT-qPCR Expression Analysis in Urochloa brizantha Under Abiotic Stresses

Urochloa brizantha is one of the most important warm season forage grasses in tropical countries. Despite its importance, there are few studies on gene expression in this species under stressful conditions. Real-time (RT-qPCR) is an accurate technique for gene quantification analysis, but reference genes must be validated under the same conditions used to assess the expression of the target genes. Here, we evaluated the stability of nine reference genes: Actin 12, Eukaryotic initiation factor 4 A, Elongation factor-1 alpha, FTSH protease 4, U2 auxiliary fator, Succinol Co-enzyme A, Tubulin alfa-5, Tubulin beta-6, Ubiquitin conjugating enzyme. Total RNA was extract from leaf tissues of U. brizantha subjected to 6, 12 and 24 h of cold and heat stresses (10 and 45 °C, respectively), and drought, including moderate (−0.5 to −0.7 MPa), severe (−1.1 to −1.8 MPa) and recovery after re-watering. The RefFinder web-based tool was used to rank the most stable reference genes for each stress. Elongation factor-1 alpha, Elongation factor-1 alpha or Ubiquitin conjugating enzyme, and Eukaryotic initiation factor 4 A were the most stable genes for heat, cold and drought stress, respectively. The expression of Rubisco large subunit gene was normalized against the most stable gene selected by ReFfinder for each stress.