Selection of reliable reference genes for RT-qPCR analysis during developmental stages and abiotic stress in Setaria viridis

Elevated CO2 increases biomass of Sorghum bicolor green prop roots under drought conditions via soluble sugar accumulation and photosynthetic activity

Elevated [CO2 ] (E[CO2 ]) mitigates agricultural losses of C4 plants under drought. Although several studies have described the molecular responses of the C4 plant species Sorghum bicolor during drought exposure, few have reported the combined effects of drought and E[CO2 ] (E[CO2 ]/D) on the roots. A previous study showed that, among plant organs, green prop roots (GPRs) under E[CO2 ]/D presented the second highest increase in biomass after leaves compared with ambient [CO2 ]/D. GPRs are photosynthetically active and sensitive to drought. To understand which mechanisms are involved in the increase in biomass of GPRs, we performed transcriptome analyses of GPRs under E[CO2 ]/D. Whole-transcriptome analysis revealed several pathways altered under E[CO2 ]/D, among which photosynthesis was strongly affected. We also used previous metabolome data to support our transcriptome data. Activities associated with photosynthesis and central metabolism increased, as seen by the upregulation of photosynthesis-related genes, a rise in glucose and polyol contents, and increased contents of chlorophyll a and carotenoids. Protein-protein interaction networks revealed that proliferation, biogenesis, and homeostasis categories were enriched and contained mainly upregulated genes. The findings suggest that the previously reported increase in GPR biomass of plants grown under E[CO2 ]/D is mainly attributed to glucose and polyol accumulation, as well as photosynthesis activity and carbon provided by respiratory CO2 refixation. Our findings reveal that an intriguing and complex metabolic process occurs in GPRs under E[CO2 ]/D, showing the crucial role of these organs in plant drought /tolerance.

How ambient temperature affects the heading date of foxtail millet (Setaria italica)

Foxtail millet (Setaria italica), a short-day plant, is one of the important crops for food security encountering climate change, particularly in regions where it is a staple food. Under the short-day condition in Taiwan, the heading dates (HDs) of foxtail millet accessions varied by genotypes and ambient temperature (AT). The allelic polymorphisms in flowering time (FT)-related genes were associated with HD variations. AT, in the range of 13°C-30°C that was based on field studies at three different latitudes in Taiwan and observations in the phytotron at four different AT regimes, was positively correlated with growth rate, and high AT promoted HD. To elucidate the molecular mechanism of foxtail millet HD, the expression of 14 key FT-related genes in four accessions at different ATs was assessed. We found that the expression levels of SiPRR95, SiPRR1, SiPRR59, SiGhd7-2, SiPHYB, and SiGhd7 were negatively correlated with AT, whereas the expression levels of SiEhd1, SiFT11, and SiCO4 were positively correlated with AT. Furthermore, the expression levels of SiGhd7-2, SiEhd1, SiFT, and SiFT11 were significantly associated with HD. A coexpression regulatory network was identified that shown genes involved in the circadian clock, light and temperature signaling, and regulation of flowering, but not those involved in photoperiod pathway, interacted and were influenced by AT. The results reveal how gene × temperature and gene × gene interactions affect the HD in foxtail millet and could serve as a foundation for breeding foxtail millet cultivars for shift production to increase yield in response to global warming.

Clinical efficacy of azacytidine and venetoclax and prognostic impact of Tim-3 and galectin-9 in acute myeloid leukemia and high-risk myelodysplastic syndromes: A single-center real-life experience

Treatment of acute myeloid leukemia (AML) and high-risk myelodysplastic syndromes (MDS) is difficult in older patients with comorbidities and high-risk disease factors. Venetoclax, the first-in-class Bcl-2 inhibitor, has proven efficacy and safety in combination with azacytidine for treatment of high-risk myeloid diseases. In this single-center real-life retrospective study, a total of 27 consecutive patients treated with azacytidine plus venetoclax were included, and clinical outcomes, hematological improvements, and biomarkers of responsiveness to therapy were compared to those observed in an historical cohort of 95 consecutive patients treated with azacytidine as single agent. Azacytidine plus venetoclax was effective and safe in older and frail AML and high-risk MDS patients, with median overall survival of 22.3 months, higher than that reported in phase III trial (14.7 months), and higher than that of historical cohort (5.94 months). Progression-free survival was higher in patients treated with the drug combination compared to those treated with azacytidine as single agent (p = 0.0065). Clinical benefits might increase when azacytidine and venetoclax are administered as upfront therapy (p = 0.0500). We showed that Tim-3 expression could be a promising therapeutic target in refractory/relapsed patients, and galectin-9 a biomarker of responsiveness to therapy. Moreover, patients treated with azacytidine and venetoclax displayed a higher overall survival regardless the presence of negative prognostic markers at diagnosis (e.g., increased WT1 copies and/or normalized blast count). These encouraging results in a real-world setting supported efficacy and safety of azacytidine plus venetoclax as upfront therapy in AML and high-risk MDS, with clinical outcomes comparable to those of clinical trials when an appropriate venetoclax management with bone marrow assessment at every first, second, fourth, and eighth cycle, and dose adjustments for toxicities are performed.

Multi-omics intervention in Setaria to dissect climate-resilient traits: Progress and prospects

Millets constitute a significant proportion of underutilized grasses and are well known for their climate resilience as well as excellent nutritional profiles. Among millets, foxtail millet (Setaria italica) and its wild relative green foxtail (S. viridis) are collectively regarded as models for studying broad-spectrum traits, including abiotic stress tolerance, C₄ photosynthesis, biofuel, and nutritional traits. Since the genome sequence release, the crop has seen an exponential increase in omics studies to dissect agronomic, nutritional, biofuel, and climate-resilience traits. These studies have provided first-hand information on the structure, organization, evolution, and expression of several genes; however, knowledge of the precise roles of such genes and their products remains elusive. Several open-access databases have also been instituted to enable advanced scientific research on these important crops. In this context, the current review enumerates the contemporary trend of research on understanding the climate resilience and other essential traits in Setaria, the knowledge gap, and how the information could be translated for the crop improvement of related millets, biofuel crops, and cereals. Also, the review provides a roadmap for studying other underutilized crop species using Setaria as a model.

MicroRNAs expression profiles in early responses to different levels of water deficit in Setaria viridis

Water deficit is a major constraint for crops of economic importance in almost all agricultural regions. However, plants have an active defense system to adapt to these adverse conditions, acting in the reprogramming of gene expression responsible for encoding microRNAs (miRNAs). These miRNAs promote the regulation to the target gene expression by the post-transcriptional (PTGS) and transcriptional gene silencing (TGS), modulating several pathways including defense response to water deficit. The broader knowledge of the miRNA expression profile and its regulatory networks in response to water deficit can provide evidence for the development of new biotechnological tools for genetic improvement of several important crops. In this study, we used Setaria viridis accession A10.1 as a C4 model plant to widely investigate the miRNA expression profile in early responses to different levels of water deficit. Ecophysiological studies in Setaria viridis under water deficit and after rewatering demonstrated a drought tolerant accession, capable of a rapid recovery from the stress. Deep small RNA sequencing and degradome studies were performed in plants submitted to drought to identify differentially expressed miRNA genes and their predicted targets, using in silico analysis. Our findings showed that several miRNAs were differentially modulated in response to distinctive levels of water deficit and after rewatering. The predicted mRNA targets mainly corresponded to genes related to cell wall remodeling, antioxidant system and drought-related transcription factors, indicating that these genes are rapidly regulated in early responses to drought stress. The implications of these modulations are extensively discussed, and higher-effect miRNAs are suggested as major players for potential use in genetic engineering to improve drought tolerance in economically important crops, such as sugarcane, maize, and sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01226-z.

Characterizing serotonin biosynthesis in Setaria viridis leaves and its effect on aphids

A combined transcriptomic and metabolic analysis of Setaria viridis leaves responding to aphid infestation was used to identify genes related to serotonin biosynthesis. Setaria viridis (green foxtail), a short life-cycle C4 plant in the Poaceae family, is the wild ancestor of Setaria italica (foxtail millet), a resilient crop that provides good yields in dry and marginal land. Although S. viridis has been studied extensively in the last decade, the molecular mechanisms of insect resistance in this species remain under-investigated. To address this issue, we performed a metabolic analysis of S. viridis and discovered that these plants accumulate the tryptophan-derived compounds tryptamine and serotonin. To elucidate the defensive functions of serotonin, Rhophalosiphum padi (bird cherry-oat aphids) were exposed to this compound, either by exogenous application to the plant medium or with artificial diet bioassays. In both cases, exposure to serotonin increased aphid mortality. To identify genes that are involved in serotonin biosynthesis, we conducted a transcriptome analysis and identified several predicted S. viridis tryptophan decarboxylase (TDC) and tryptamine 5-hydroxylase (T5H) genes. Two candidate genes were ectopically expressed in Nicotiana tabacum, where SvTDC1 (Sevir.6G066200) had tryptophan decarboxylase activity, and SvT5H1 (Sevir.8G219600) had tryptamine hydroxylase activity. Moreover, the function of the SvTDC1 gene was validated using virus-induced gene silencing in S. italica, which caused a reduction in serotonin levels. This study provides the first evidence of serotonin biosynthesis in Setaria leaves. The biosynthesis of serotonin may play an important role in defense responses and could prove to be useful for developing more pest-tolerant Setaria italica cultivars.

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 Oil Palm Early Nodulin 93 Protein Gene (EgENOD93) Enhances In Vitro Shoot Regeneration in Arabidopsis thaliana

EgENOD93 was first identified in a cDNA microarray study of oil palm tissue culture where it was highly expressed in leaf explants with embryogenic potential. Functional characterization via an RNA interference study of its orthologue in Medicago truncatula demonstrated a significant role of this gene in somatic embryo formation. In this study, EgENOD93 was overexpressed in the important model plant Arabidopsis thaliana to investigate the embryogenic potential of EgENOD93 transgenic Arabidopsis explants compared to explants from control plants (pMDC140 and WT). Experiments using leaf explants revealed higher numbers of regenerated shoots at day 27 in all the homozygous transgenic Arabidopsis cultures (Tg01, Tg02 and Tg03) compared to controls. The expression level of EgENOD93 in Arabidopsis cultures was quantified using reverse transcription quantitative real-time PCR (RT-qPCR). The results supported the overexpression of this gene in transgenic Arabidopsis cultures, with 6 and 10 times higher expression of EgENOD93 in callus at Day 9 and Day 20, respectively. Overall, the results support the role of EgENOD93 in the enhancement of shoot regeneration in transgenic Arabidopsis. This together with the previous results observed in oil palm and Medicago truncatula suggests that ENOD93 plays a key role in the induction of somatic embryogenesis. A similarity to early nodulation-like ontogeny is possible.

Validation of reference genes for quantitative real-time PCR in chemical exposed and at different age's brackish water flea Diaphanosoma celebensis

Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), a primary approach for evaluating gene expression, requires an appropriate normalization strategy to confirm relative gene expression levels by comparison, and rule out variations that might occur in analytical procedures. The best option is to use a reference gene whose expression level is stable across various experimental conditions to compare the mRNA levels of a target gene. However, there is limited information on how the reference gene is differentially expressed at different ages (growth) in small invertebrates with notable changes such as molting. In this study, expression profiles of nine candidate reference genes from the brackish water flea, Diaphanosoma celebensis, were evaluated under diverse exposure to toxicants and according to growth. As a result, four different algorithms showed similar stabilities of genes for chemical exposures in the case of limited conditions using the same developmental stage (H2A was stable, whereas Act was fairly unstable in adults), while the results according to age showed a significantly different pattern in suite of candidate reference genes. This affected the results of genes EcRA and GST, which are involved in development and detoxification mechanisms, respectively. Our finding is the first step towards establishing a standardized real-time qRT-PCR analysis of this environmentally important invertebrate that has potential for aquatic ecotoxicology, particularly in estuarine environments.

The Characteristics of Mercury Flux at the Interfaces between Two Typical Plants and the Air in Leymus chinensis Grasslands

Mercury is a global pollutant. The mercury exchanges between vegetation and the atmosphere are important for the global mercury cycle. Grassland ecosystems occupy more than 25% of the global land area and have different succession processes and ecological functions. The current research regarding mercury exchanges between forests and the atmosphere have attracted much attention, but the research regarding grasslands tends to be rare. To reveal the characteristics of mercury exchanges in grasslands, this study conducted field in-situ monitoring experiments in a Leymus meadow grassland regions of the Songnen Plains in northeastern China. The exchange flux values of the GEM (gaseous element mercury) between the plants and the atmosphere were measured using a dynamic flux bag method (DFB). The experiments were conducted for the purpose of assessing the mercury flux levels between the vegetation and the atmosphere in a typical Leymus chinensis meadow. The goal was to further the understanding of the change characteristics and influential factors and to describe the source and sink actions and dynamics between the grassland vegetation and the atmosphere. The diurnal variation characteristics were as follows: High during the day and low at night, with peaks generally appearing at noon. The growing period was characterized by absorption peaks of atmospheric mercury by the plants. The breeding period was characterized by the peak release of atmospheric mercury by the plants. The change characteristics were as follows: During the growing period, the duration of the plants in a mercury absorption state exceeded 96.5%, which was represented as the net sink of the atmospheric mercury. During the breeding period, the time of mercury release ranged between 46.4% and 66.8%, making the breeding period the net source of atmospheric mercury. The results of this study's analysis indicated that each environmental factor was correlated with the mercury flux, and the environmental factors had different effects on the mercury flux during the different stages of plant growth. The atmospheric mercury concentration levels were the main factor during the growing period. Atmospheric humidity was the main factor during the breeding period. Solar radiation was the decisive factor during the entire experimental period.

Efficient genome editing and gene knockout in Setaria viridis with CRISPR/Cas9 directed gene editing by the non-homologous end-joining pathway

The CRISPR/Cas9 system has been used for genome editing in several organisms, including higher plants. This system induces site-specific mutations in the genome based on the nucleotide sequence of engineered guide RNAs. The complex genomes of C4 grasses makes genome editing a challenge in key grass crops like maize (Zea mays), sorghum (Sorghum bicolor), Brachiaria spp., switchgrass (Panicum virgatum), and sugarcane (Saccharum spp.). Setaria viridis is a diploid C4 grass widely used as a model for these C4 crop plants. Here, an optimized CRISPR/Cas9 binary vector that exploits the non-homologous end joining (NHEJ) system was used to knockout a green fluorescent protein (gfp) transgene in S. viridis accession A10.1. Transformation of embryogenic callus by A. tumefaciens generated ten glufosinate-ammonium resistant transgenic events. In the T0 generation, 60% of the events were biallelic mutants in the gfp transgene with no detectable accumulation of GFP protein and without insertions or deletions in predicted off-target sites. The gfp mutations generated by CRISPR/Cas9 were stable and displayed Mendelian segregation in the T1 generation. Altogether, the system described here is a highly efficient genome editing system for S. viridis, an important model plant for functional genomics studies in C4 grasses. Also, this system is a potential tool for improvement of agronomic traits in C4 crop plants with complex genomes.

Comparative Transcriptome Profiling of Kappaphycus alvarezii (Rhodophyta, Solieriaceae) in Response to Light of Different Wavelengths and Carbon Dioxide Enrichment

Rhodophyta (red algae) comprises over 6000 species, however, there have only been a few comparative transcriptomic studies due to their under-representation in genomic databases. Kappaphycus alvarezii, a Gigartinales algae, is a valuable source of carrageenan and is extensively cultivated in many countries. The majority of seaweed farming in Southeast Asia is done in intertidal zones under varying light (i.e., spectra and irradiance) and carbon dioxide (CO₂) conditions, which affects the rate of photosynthesis. This study conducted transcriptome profiling to investigate the photosynthetic mechanisms in K. alvarezii exposed to different wavelengths of light (i.e., blue, green, and red light, in comparison to white light) and CO₂ availability. We analyzed the responses of photosynthetic protein complexes to light and observed that light of different wavelengths regulates a similar set of photosynthetic apparatuses. Under CO₂ enrichment, genes encoding C₃ and C₄ enzymes were found to be actively transcribed, suggesting the likely shift in the carbon metabolism pathway or the involvement of these genes in adaptive physiological processes. This study contributes to the understanding of the regulatory mechanisms of photosynthetic carbon metabolism in red algae and has implications for the culture and commercial production of these economically valuable macroalgae.

Suppression of a BAHD acyltransferase decreases p-coumaroyl on arabinoxylan and improves biomass digestibility in the model grass Setaria viridis

Grass cell walls have hydroxycinnamic acids attached to arabinosyl residues of arabinoxylan (AX), and certain BAHD acyltransferases are involved in their addition. In this study, we characterized one of these BAHD genes in the cell wall of the model grass Setaria viridis. RNAi silenced lines of S. viridis (SvBAHD05) presented a decrease of up to 42% of ester-linked p-coumarate (pCA) and 50% of pCA-arabinofuranosyl, across three generations. Biomass from SvBAHD05 silenced plants exhibited up to 32% increase in biomass saccharification after acid pre-treatment, with no change in total lignin. Molecular dynamics simulations suggested that SvBAHD05 is a p-coumaroyl coenzyme A transferase (PAT) mainly involved in the addition of pCA to the arabinofuranosyl residues of AX in Setaria. Thus, our results provide evidence of p-coumaroylation of AX promoted by SvBAHD05 acyltransferase in the cell wall of the model grass S. viridis. Furthermore, SvBAHD05 is a promising biotechnological target to engineer crops for improved biomass digestibility for biofuels, biorefineries and animal feeding.

Effect of Nitrogen Sources on Omega-3 Polyunsaturated Fatty Acid Biosynthesis and Gene Expression in Thraustochytriidae sp

The molecular mechanism that contributes to nitrogen source dependent omega-3 polyunsaturated fatty acid (n-3 PUFA) synthesis in marine oleaginous protists Thraustochytriidae sp., was explored in this study. The fatty acid (FA) synthesis was significantly influenced by the supplement of various levels of sodium nitrate (SN) (1–50 mM) or urea (1–50 mM). Compared with SN (50 mM) cultivation, cells from urea (50 mM) cultivation accumulated 1.16-fold more n-3 PUFAs (49.49% docosahexaenoic acid (DHA) (w/w, of total FAs) and 5.28% docosapentaenoic acid (DPA) (w/w, of total FAs)). Strikingly higher quantities of short chain FAs (<18 carbons) (52.22-fold of that in urea cultivation) were produced from SN cultivation. Ten candidate reference genes (RGs) were screened by using four statistical methods (geNorm, NormFinder, Bestkeeper and RefFinder). MFT (Mitochondrial folate transporter) and NUC (Nucleolin) were determined as the stable RGs to normalize the RT-qPCR (real-time quantitative polymerase chain reaction) data of essential genes related to n-3 PUFAs-synthesis. Our results elucidated that the gene transcripts of delta(3,5)-delta(2,4)-dienoyl-CoA isomerase, enoyl-CoA hydratase, fatty acid elongase 3, long-chain fatty acid acyl-CoA ligase, and acetyl-CoA carboxylase were up-regulated under urea cultivation, contributing to the extension and unsaturated bond formation. These findings indicated that regulation of the specific genes through nitrogen source could greatly stimulate n-3 PUFA production in Thraustochytriidae sp.

Identification of Ginger (Zingiber officinale Roscoe) Reference Genes for Gene Expression Analysis

Quantitative real-time PCR (qRT-PCR) is widely used in the detection of gene expression level. However, there is no suitable ginger reference gene for qPCR analysis. Therefore, it is the primary task to select and validate the appropriate ginger reference gene to normalize the expression of target genes. In this study, 14 candidate reference genes were selected and analyzed in different tissues (leaf, and rhizome), different development stages, different varieties, and abiotic stress (ABA and salt stress). Expression stability was calculated using geNorm and NormFinder, Bestkeeper, and RefFinder. For abiotic stress and total conditions, 28S and COX were identified as the most stable genes. In addition, RPII was the most stable in the different development stages and different varieties. TEF2 and RPL2 were the least stably expressed in the tissue and all the conditions. In order to verify the feasibility of these genes as reference genes, we used the most stable and least stable reference genes to normalize the expression levels of ZoSPS genes under different conditions. This work can provide theoretical support for future research on ginger gene expression.

Analysis of reference genes stability and histidine kinase expression under cold stress in Cordyceps militaris

The development of fungal fruiting bodies from a hyphal thallus is inducible under low temperature (cold stress). The molecular mechanism has been subject to surprisingly few studies. Analysis of gene expression level has become an important means to study gene function and its regulation mechanism. But identification of reference genes (RGs) stability under cold stress have not been reported in famous medicinal mushroom-forming fungi Cordyceps militaris. Herein, 12 candidate RGs had been systematically validated under cold stress in C. militaris. Three different algorithms, geNorm, NormFinder and BestKeeper were applied to evaluate the expression stability of the RGs. Our results showed that UBC and UBQ were the most stable RGs for cold treatments in short and long periods, respectively. 2 RGs (UBC and PP2A) and 3 RGs (UBQ, TUB and CYP) were the suitable RGs for cold treatments in short and long periods, respectively. Moreover, target genes, two-component-system histidine kinase genes, were selected to validate the most and least stable RGs under cold treatment, which indicated that use of unstable expressed genes as RGs leads to biased results. Our results provide a good starting point for accurate reverse transcriptase quantitative polymerase chain reaction normalization by using UBC and UBQ in C. militaris under cold stress and better support for understanding the mechanism of response to cold stress and fruiting body formation in C. militaris and other mushroom-forming fungi in future research.

Plant science's next top models

BACKGROUND

Model organisms are at the core of life science research. Notable examples include the mouse as a model for humans, baker's yeast for eukaryotic unicellular life and simple genetics, or the enterobacteria phage λ in virology. Plant research was an exception to this rule, with researchers relying on a variety of non-model plants until the eventual adoption of Arabidopsis thaliana as primary plant model in the 1980s. This proved to be an unprecedented success, and several secondary plant models have since been established. Currently, we are experiencing another wave of expansion in the set of plant models.

SCOPE

Since the 2000s, new model plants have been established to study numerous aspects of plant biology, such as the evolution of land plants, grasses, invasive and parasitic plant life, adaptation to environmental challenges, and the development of morphological diversity. Concurrent with the establishment of new plant models, the advent of the 'omics' era in biology has led to a resurgence of the more complex non-model plants. With this review, we introduce some of the new and fascinating plant models, outline why they are interesting subjects to study, the questions they will help to answer, and the molecular tools that have been established and are available to researchers.

CONCLUSIONS

Understanding the molecular mechanisms underlying all aspects of plant biology can only be achieved with the adoption of a comprehensive set of models, each of which allows the assessment of at least one aspect of plant life. The model plants described here represent a step forward towards our goal to explore and comprehend the diversity of plant form and function. Still, several questions remain unanswered, but the constant development of novel technologies in molecular biology and bioinformatics is already paving the way for the next generation of plant models.

An improved protocol for efficient transformation and regeneration of Setaria italica

An efficient and improved transformation method for functional genetics studies in S. italica, being a boon for the Setaria scientific community and for crop improvement. Foxtail millet (Setaria italica) is a short life cycle C4 plant, with sequenced genome, and a potential model plant for C4 species. S. italica is also important on a global food security and healthiness context due to its importance in arid and semi-arid areas. However, despite its importance, there are just few transformation protocols directed to this species. The current protocols reached about 5.5-9% of efficiency, which do not make it a valuable model organism. Different types of explants were used in the above mentioned methods, such as immature and mature inflorescence and shoot apex. However, these techniques have many limitations, such as unavailability of explants throughout the year and a crucial, laborious and considerable time-consuming selection. Aiming a simplified and efficient methodology, we adopted dry mature seeds as explants, which are available in abundance, are constant along the year and well responsive to tissue culture, in addition to a differentiated approach that reaches on an average 19.2% transformation efficiency of S. italica. Thus, we propose a protocol that optimizes the transformation efficiency of this cereal crop allowing a high increase on transformation and regeneration rates. Our transformation protocol provides an interesting tool for Setaria community research as well as enables new strategies for breeding enhanced productivity in the species.

Systematic Identification of Housekeeping Genes Possibly Used as References in Caenorhabditis elegans by Large-Scale Data Integration

For accurate gene expression quantification, normalization of gene expression data against reliable reference genes is required. It is known that the expression levels of commonly used reference genes vary considerably under different experimental conditions, and therefore, their use for data normalization is limited. In this study, an unbiased identification of reference genes in Caenorhabditis elegans was performed based on 145 microarray datasets (2296 gene array samples) covering different developmental stages, different tissues, drug treatments, lifestyle, and various stresses. As a result, thirteen housekeeping genes (rps-23, rps-26, rps-27, rps-16, rps-2, rps-4, rps-17, rpl-24.1, rpl-27, rpl-33, rpl-36, rpl-35, and rpl-15) with enhanced stability were comprehensively identified by using six popular normalization algorithms and RankAggreg method. Functional enrichment analysis revealed that these genes were significantly overrepresented in GO terms or KEGG pathways related to ribosomes. Validation analysis using recently published datasets revealed that the expressions of newly identified candidate reference genes were more stable than the commonly used reference genes. Based on the results, we recommended using rpl-33 and rps-26 as the optimal reference genes for microarray and rps-2 and rps-4 for RNA-sequencing data validation. More importantly, the most stable rps-23 should be a promising reference gene for both data types. This study, for the first time, successfully displays a large-scale microarray data driven genome-wide identification of stable reference genes for normalizing gene expression data and provides a potential guideline on the selection of universal internal reference genes in C. elegans, for quantitative gene expression analysis.

Comprehensive transcriptome analysis of reference genes for fruit development of Euscaphis konishii

Background

Quantitativereal-time reverse transcriptase polymerase chain reaction is the common method to quantify relative gene expression. Normalizating using reliable genes is critical in correctly interpreting expression data from qRT-PCR. Euscaphis konishii is a medicinal plant with a long history in China, which has various chemical compounds in fruit. However, there is no report describing the selection of reference genes in fruit development of Euscaphis konishii.

Methods

We selected eight candidate reference genes based on RNA-seq database analysis, and ranked expression stability using statistical algorithms GeNorm, NormFinder, BestKeeper and ReFinder. Finally, The nine genes related to the anthocyanin synthesis pathway of Euscaphis konishii were used to verify the suitability of reference gene.

Results

The results showed that the stability of EkUBC23, EkCYP38 and EkGAPDH2 was better, and the low expression reference genes (EkUBC23 and EkCYP38) were favourable for quantifying low expression target genes, while the high expression reference gene (EkGAPDH2) was beneficial for quantifying high expression genes. In this study, we present the suitable reference genes for fruit development of Euscaphis konishii based on transcriptome data, our study will contribute to further studies in molecular biology and gene function on Euscaphis konishii and other closely related species.

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.

Selection and validation of reference genes for gene expression studies in Pseudomonas brassicacearum GS20 using real-time quantitative reverse transcription PCR

Pseudomonas brassicacearum GS20 is an antagonistic strain of bacteria recently isolated from the rhizosphere of Codonopsis pilosula. No validated reference gene has been indentified from P. brassicacearum to use in the normalization of real-time quantitative reverse transcription-PCR data. Therefore, in this study, nine candidate reference genes (recA, gyrA, rpoD, proC, gmk, rho, 16S, ftsz, and secA) were assessed at different growth phases and under various growth conditions. The expression stability of these candidate genes was evaluated using BestKeeper, NormFinder and GeNorm. In general, the results showed rho, rpoD and gyrA were the most suitable reference genes for P. brassicacearum GS20. The relative expression of iron-regulated gene (fhu) was normalized to verify the reliability of the proposed reference genes under iron-replete and iron-limited conditions. The trend in relative expression was consistent with the change in siderophore production under different iron conditions. This study presents reliable reference genes for transcriptional studies in P. brassicacearum GS20 under the chosen experimental conditions.

Genome-wide characterization of the laccase gene family in Setaria viridis reveals members potentially involved in lignification

Five laccase genes are potentially involved in developmental lignification in the model C4 grass Setaria viridis and their different tissue specificities suggest subfunctionalization events. Plant laccases are copper-containing glycoproteins involved in monolignol oxidation and, therefore, their activity is essential for lignin polymerization. Although these enzymes belong to large multigene families with highly redundant members, not all of them are thought to be involved in lignin metabolism. Here, we report on the genome-wide characterization of the laccase gene family in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol oxidation. A total of 52 genes encoding laccases (SvLAC1 to SvLAC52) were found in the genome of S. viridis, and phylogenetic analyses showed that these genes were heterogeneously distributed among the characteristic six subclades of the family and are under relaxed selective constraints. The observed expansion in the total number of genes in this species was mainly caused by tandem duplications within subclade V, which accounts for 68% of the whole family. Comparative phylogenetic analyses showed that the expansion of subclade V is specifically observed for the Paniceae tribe within the Panicoideae subfamily in grasses. Five SvLAC genes (SvLAC9, SvLAC13, SvLAC15, SvLAC50, and SvLAC52) fulfilled the criteria established to identify lignin-related candidates: (1) phylogenetic proximity to previously characterized lignin-related laccases from other species, (2) similar expression pattern to that observed for lignin biosynthetic genes in the S. viridis elongating internode, and (3) high expression in S. viridis tissues undergoing active lignification. In addition, in situ hybridization experiments not only confirmed that these selected SvLAC genes were expressed in lignifying cells, but also that their expression showed different tissue specificities, suggesting subfunctionalization events within the family. These five laccase genes are strong candidates to be involved in lignin polymerization in S. viridis and might be good targets for lignin bioengineering strategies.

Identification of Appropriate Reference Genes for Normalizing miRNA Expression in Citrus Infected by Xanthomonas citri subsp. citri

MicroRNAs (miRNAs) are short noncoding RNA molecules that regulate gene expression at the posttranscriptional level. Reverse transcription-quantitative PCR (RT-qPCR) is one of the most common methods used for quantification of miRNA expression, and the levels of expression are normalized by comparing with reference genes. Thus, the selection of reference genes is critically important for accurate quantification. The present study was intended to identify appropriate miRNA reference genes for normalizing the level of miRNA expression in Citrus sinensis L. Osbeck and Citrus reticulata Blanco infected by Xanthomonas citri subsp. citri, which caused citrus canker disease. Five algorithms (Delta Ct, geNorm, NormFinder, BestKeeper and RefFinder) were used for screening reference genes, and two quantification approaches, poly(A) extension RT-qPCR and stem-loop RT-qPCR, were used to determine the most appropriate method for detecting expression patterns of miRNA. An overall comprehensive ranking output derived from the multi-algorithms showed that poly(A)-tailed miR162-3p/miR472 were the best reference gene combination for miRNA RT-qPCR normalization in citrus canker research. Candidate reference gene expression profiles determined by poly(A) RT-qPCR were more consistent in the two citrus species. To the best of our knowledge, this is the first systematic comparison of two miRNA quantification methods for evaluating reference genes. These results highlight the importance of rigorously assessing candidate reference genes and clarify some contradictory results in miRNA research on citrus.

Stable reference genes for RT-qPCR analysis of gene expression in the Musa acuminata-Pseudocercospora musae interaction

Leaf pathogens are limiting factors in banana (Musa spp.) production, with Pseudocercospora spp. responsible for the important Sigatoka disease complex. In order to investigate cellular processes and genes involved in host defence responses, quantitative real-time PCR (RT-qPCR) is an analytical technique for gene expression quantification. Reliable RT-qPCR data, however, requires that reference genes for normalization of mRNA levels in samples are validated under the conditions employed for expression analysis of target genes. We evaluated the stability of potential reference genes ACT1, α-TUB, UBQ1, UBQ2, GAPDH, EF1α, APT and RAN. Total RNA was extracted from leaf tissues of Musa acuminata genotypes Calcutta 4 (resistant) and Cavendish Grande Naine (susceptible), both subjected to P. musae infection. Expression stability was determined with NormFinder, BestKeeper, geNorm and RefFinder algorithms. UBQ2 and RAN were the most stable across all M. acuminata samples, whereas when considering inoculated and non-inoculated leaf samples, APT and UBQ2 were appropriate for normalization in Calcutta 4, with RAN and α-TUB most stable in Cavendish Grande Naine. This first study of reference genes for relative quantification of target gene expression in the M. acuminata-P. musae interaction will enable reliable analysis of gene expression in this pathosystem, benefiting elucidation of disease resistance mechanisms.

The lignin toolbox of the model grass Setaria viridis

The core set of biosynthetic genes potentially involved in developmental lignification was identified in the model C4 grass Setaria viridis. Lignin has been recognized as a major recalcitrant factor negatively affecting the processing of plant biomass into bioproducts. However, the efficient manipulation of lignin deposition in order to generate optimized crops for the biorefinery requires a fundamental knowledge of several aspects of lignin metabolism, including regulation, biosynthesis and polymerization. The current availability of an annotated genome for the model grass Setaria viridis allows the genome-wide characterization of genes involved in the metabolic pathway leading to the production of monolignols, the main building blocks of lignin. Here we performed a comprehensive study of monolignol biosynthetic genes as an initial step into the characterization of lignin metabolism in S. viridis. A total of 56 genes encoding bona fide enzymes catalyzing the consecutive ten steps of the monolignol biosynthetic pathway were identified in the S. viridis genome. A combination of comparative phylogenetic studies, high-throughput expression analysis and quantitative RT-PCR analysis was further employed to identify the family members potentially involved in developmental lignification. Accordingly, 14 genes clustered with genes from closely related species with a known function in lignification and showed an expression pattern that correlates with lignin deposition. These genes were considered the "core lignin toolbox" responsible for the constitutive, developmental lignification in S. viridis. These results provide the basis for further understanding lignin deposition in C4 grasses and will ultimately allow the validation of biotechnological strategies to produce crops with enhanced processing properties.

Selection and Validation of Reference Genes for qRT-PCR in Lentinula edodes under Different Experimental Conditions

Lentinula edodes is the most consumed mushroom in Asia due to its nutritional and medicinal values, and the optimal reference gene is crucial for normalization of its gene expression analysis. Here, the expression stability of 18 candidate reference genes (CRGs) in L. edodes was analyzed by three statistical algorithms (geNorm, NormFinder and BestKeeper) under different stresses (heat, cadmium excess and Trichoderma atroviride infection), different substrates (straw, sawdust and corn stalk) and different development stages (mycelia, primordia and fruit bodies). Among the 18 CRGs, 28S, Actin and α-tub exhibited the highest expression stability in L. edodes under all conditions, while GPD, SPRYP and MSF showed the least stable expression. The best reference gene in different conditions was different. The pairwise variation values showed that two genes would be sufficient for accurate normalization under different conditions of L. edodes. This study will contribute to more accurate estimation of the gene relative expression levels under different conditions using the optimal reference gene in qRT-PCR (quantitative reverse transcription polymerase chain reaction) analysis.

Reference genes selection for Calotropis procera under different salt stress conditions

Calotropis procera is a perennial Asian shrub with significant adaptation to adverse climate conditions and poor soils. Given its increased salt and drought stress tolerance, C. procera stands out as a powerful candidate to provide alternative genetic resources for biotechnological approaches. The qPCR (real-time quantitative polymerase chain reaction), widely recognized among the most accurate methods for quantifying gene expression, demands suitable reference genes (RGs) to avoid over- or underestimations of the relative expression and incorrect interpretation. This study aimed at evaluating the stability of ten RGs for normalization of gene expression of root and leaf of C. procera under different salt stress conditions and different collection times. The selected RGs were used on expression analysis of three target genes. Three independent experiments were carried out in greenhouse with young plants: i) Leaf₁₀₀ = leaf samples collected 30 min, 2 h, 8 h and 45 days after NaCl-stress (100 mM NaCl); ii) Root₅₀ and iii) Root₂₀₀ = root samples collected 30 min, 2 h, 8 h and 1day after NaCl-stress (50 and 200 mM NaCl, respectively). Stability rank among the three algorithms used showed high agreement for the four most stable RGs. The four most stable RGs showed high congruence among all combination of collection time, for each software studied, with minor disagreements. CYP23 was the best RG (rank of top four) for all experimental conditions (Leaf₁₀₀, Root_50, and Root₂₀₀). Using appropriated RGs, we validated the relative expression level of three differentially expressed target genes (NAC78, CNBL4, and ND1) in Leaf₁₀₀ and Root₂₀₀ samples. This study provides the first selection of stable reference genes for C. procera under salinity. Our results emphasize the need for caution when evaluating the stability RGs under different amplitude of variable factors.

Identification and characterization of core abscisic acid (ABA) signaling components and their gene expression profile in response to abiotic stresses in Setaria viridis

Abscisic acid (ABA) is an essential phytohormone that regulates growth, development and adaptation of plants to environmental stresses. In Arabidopsis and other higher plants, ABA signal transduction involves three core components namely PYR/PYL/RCAR ABA receptors (PYLs), type 2C protein phosphatases (PP2Cs) and class III SNF-1-related protein kinase 2 (SnRK2s). In the present study, we reported the identification and characterization of the core ABA signaling components in Setaria viridis, an emerging model plant for cereals and feedstock crops presenting C4 metabolism, leading to the identification of eight PYL (SvPYL1 to 8), twelve PP2C (SvPP2C1 to 12) and eleven SnRK2 (SvSnRK2.1 through SvSnRK2.11) genes. In order to study the expression profiles of these genes, two different S. viridis accessions (A10.1 and Ast-1) were submitted to drought, salinity and cold stresses, in addition to application of exogenous ABA. Differential gene expression profiles were observed in each treatment and plant genotype, demonstrating variations of ABA stress responses within the same species. These differential responses to stresses were also assessed by physiological measurements such as photosynthesis, stomatal conductance and transpiration rate. This study allows a detailed analysis of gene expression of the core ABA signaling components in Setaria viridis submitted to different treatments and provides suitable targets for genetic engineering of C4 plants aiming tolerance to abiotic stresses.

Selection and Validation of Reference Genes for Gene Expression Studies by RT-PCR in Dalbergia odorifera

Perennial tree Dalbergia odorifera T. Chen could form the precious heartwood used to produce chinese traditional medicine, rosewood furniture and fragrances. However the formation of heartwood is time-consuming and low efficient, leading to the severe destruction of its wild resources. Thus, it is urgent to study the molecular mechanism of heartwood formation in D. odorifera. But till now, there is no report about the reference gene selection in this species. In this study, the expression stability of nine candidate reference genes were evaluated across different tissues and stems treated by wound and chemical stimulators. Four algorithms were applied to obtain the robust genes. The results support HIS2, GAPDH, and CYP to be the most stable reference genes in samples under different wound treatments while DNAj was the least stable. In different tissues, HIS2, UBQ, and RPL were the most stable reference genes while DNAj was the least stable. The selected reference genes were validated through the normalization of the qRT-PCR data of six heartwood related genes in terpene biosynthesis pathway and ethylene signal pathway. The results showed that their expression levels were accurate when they were normalized by the most stable reference gene HIS2, or by the combination of the two or three most stable reference genes. These results demonstrated that these selected reference genes are reliable.

SiSTL1, encoding a large subunit of ribonucleotide reductase, is crucial for plant growth, chloroplast biogenesis, and cell cycle progression in Setaria italica

The activity of ribonucleotide reductase (RNR), which catalyses the transformation of four ribonucleoside diphosphates (NDPs) to their corresponding deoxyribonucleoside diphosphates (dNDPs), is the main determiner of the cellular concentration of dNTP pools and should be tightly coordinated with DNA synthesis and cell-cycle progression. Constitutively increased or decreased RNR activity interferes with DNA replication and leads to arrested cell cycle progression; however, the mechanisms underlying these disruptive effects in higher plants remain to be uncovered. In this study, we identified a RNR large subunit mutant, sistl1, in Setaria italica (foxtail millet), which exhibited growth retardation as well as striped leaf phenotype, i.e. irregularly reduced leaf vein distances and decreased chloroplast biogenesis. We determined that a Gly737 to Glu substitution occurring in the C-terminus of the SiSTL1 protein slightly affected its optimal function, leading in turn to the reduced expression of genes variously involved in the assembly and activation of the DNA pre-replicative complex, elongation of replication forks and S phase entry. Our study provides new insights into how SiSTL1 regulates plant growth, chloroplast biogenesis, and cell cycle progression in Poaceae crops.

SiSTL1, encoding a large subunit of ribonucleotide reductase, is crucial for plant growth, chloroplast biogenesis, and cell cycle progression in Setaria italica

The activity of ribonucleotide reductase (RNR), which catalyses the transformation of four ribonucleoside diphosphates (NDPs) to their corresponding deoxyribonucleoside diphosphates (dNDPs), is the main determiner of the cellular concentration of dNTP pools and should be tightly coordinated with DNA synthesis and cell-cycle progression. Constitutively increased or decreased RNR activity interferes with DNA replication and leads to arrested cell cycle progression; however, the mechanisms underlying these disruptive effects in higher plants remain to be uncovered. In this study, we identified a RNR large subunit mutant, sistl1, in Setaria italica (foxtail millet), which exhibited growth retardation as well as striped leaf phenotype, i.e. irregularly reduced leaf vein distances and decreased chloroplast biogenesis. We determined that a Gly737 to Glu substitution occurring in the C-terminus of the SiSTL1 protein slightly affected its optimal function, leading in turn to the reduced expression of genes variously involved in the assembly and activation of the DNA pre-replicative complex, elongation of replication forks and S phase entry. Our study provides new insights into how SiSTL1 regulates plant growth, chloroplast biogenesis, and cell cycle progression in Poaceae crops.

Genome-Wide Identification and Evaluation of New Reference Genes for Gene Expression Analysis Under Temperature and Salinity Stresses in Ciona savignyi

Rapid adaptation/accommodation to changing environments largely contributes to maximal survival of invaders during biological invasions, usually leading to success in crossing multiple barriers and finally in varied environments in recipient habitats. Gene expression is one of the most important and rapid ways during responses to environmental stresses. Selection of proper reference genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). Here we identified eight candidate novel reference genes from the RNA-Seq data in an invasive model ascidian Ciona savignyi under temperature and salinity stresses. Subsequently, the expression stability of these eight novel reference genes, as well as other six traditionally used reference genes, was evaluated using RT-qPCR and comprehensive tool RefFinder. Under the temperature stress, two traditional reference genes, ribosomal proteins S15 and L17 (RPS15, RPL17), and one novel gene Ras homolog A (RhoA), were recommended as the top three stable genes, which can be used to normalize target genes with a high and moderate expression level, respectively. Under the salinity stress, transmembrane 9 superfamily member (TMN), MOB kinase activator 1A-like gene (MOB) and ubiquitin-conjugating enzyme (UBQ2) were suggested as the top three stable genes. On the other hand, several commonly used reference genes such as α-tubulin (TubA), β-tubulin (TubB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) showed unstable expressions, thus these genes should not be used as internal controls for gene expression analysis. We also tested the expression level of an important stress response gene, large proline-rich protein bag6-like gene (BAG) using different reference genes. As expected, we observed different results and conclusions when using different normalization methods, thus suggesting the importance of selection of proper reference genes and associated normalization methods. Our results provide a valuable reference gene resource for the normalization of gene expression in the study of environmental adaptation/accommodation during biological invasions using C. savignyi as a model.

Selection and validation of reference genes for target gene analysis with quantitative real-time PCR in the leaves and roots of Carex rigescens under abiotic stress

Carex rigescens is an ornamental turfgrass in northern China which has a relatively low maintenance cost and robust tolerance to many adverse environmental conditions, so it could be considered a new material for researching into plant stress resistance. However, suitable reference genes are vacant for obtaining reliable results in quantitative real-time PCR (qRT-PCR) analysis of C. rigescens in adversity research. In this study, we tested the expression stability of nine potential reference genes in leaves and roots under five different abiotic stress conditions, including cold, salt, heat, osmotic and cadmium (Cd). We then selected the best reference genes according to the analysis results calculated by three algorithmic programs (geNorm, NormFinder and BestKeeper) and used the RankAggreg package to merge the outputted data. The results showed that combinations of at least two reference genes should be used for reliable normalization except in heat-treated root samples, which require three reference genes. eIF-4α, GADPH, SAND and PEPKR1 and their combination were found to be the most stably expressed reference genes, while SAM, TUA-α and UPL7 were the three least stable reference genes among most of experimental samples. In addition, five stress-induced genes (Cu-Zn SOD, P5CS, LEA, GST, and APX) were chosen to verify the stability of the selected reference genes in various tissues and under various stress conditions. The results of this study will provide an important fundamental basis both for gene expression verification for transcriptomic and proteomic analyses and for gene expression analysis for future gene function research in C. rigescens.

Selection and Validation of Novel RT-qPCR Reference Genes under Hormonal Stimuli and in Different Tissues of Santalum album

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used technique to investigate gene expression levels due to its high throughput, specificity, and sensitivity. An appropriate reference gene is essential for RT-qPCR analysis to obtain accurate and reliable results. To date, no reliable reference gene has been validated for the economically tropical tree, sandalwood (Santalum album L.). In this study, 13 candidate reference genes, including 12 novel putative reference genes selected from a large set of S. album transcriptome data, as well as the currently used β-actin gene (ACT), were validated in different tissues (stem, leaf, root and callus), as well as callus tissue under salicylic acid (SA), jasmonic acid methyl ester (MeJA), and gibberellin (GA) treatments using geNorm, NormFinder, BestKeeper, Delta Ct and comprehensive RefFinder algorithms. Several novel candidate reference genes were much more stable than the currently used traditional gene ACT. ODD paired with Fbp1 for SA treatment, CSA and Fbp3 for MeJA treatment, PP2C and Fbp2 for GA treatment, as well as Fbp1 combined with Fbp2 for the total of three hormone treatments were the most accurate reference genes, respectively. FAB1A, when combined with PP2C, was identified as the most suitable reference gene combination for the four tissues tested, while the combination of HLMt, PPR and FAB1A were the most optimal reference genes for all of the experimental samples. In addition, to verify our results, the relative expression level of the SaSSy gene was evaluated by the validated reference genes and their combinations in the three S. album tissues and under MeJA treatment. The evaluated reference genes in this study will improve the accuracy of RT-qPCR analysis and will benefit S. album functional genomics studies in different tissues and under hormone stimuli in the future.

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.

Comprehensive evaluation of candidate reference genes for real-time quantitative PCR (RT-qPCR) data normalization in nutri-cereal finger millet [Eleusine Coracana (L.)]

Finger millet (Eleusine coracana L.) is an annual herbaceous self-pollinating C4 cereal crop of the arid and semi-arid regions of the world. Finger millet is a food security crop proven to have resilience to changing climate and scores very high in nutrition. In the current study, we have assessed sixteen candidate reference genes for their appropriateness for the normalization studies in finger millet subjected to experimental regimes and treatments. Ten candidate reference genes (GAPDH, β-TUB, CYP, EIF4α, TIP41, UBC, G6PD, S24, MACP and MDH) were cloned and six (ACT, ELF1α, PP2A, PT, S21 and TFIID) were mined from the NCBI database as well as from the literature. Expression stability ranking of the finger millet reference genes was validated using four different statistical tools i.e., geNorm, NormFinder, BestKeeper, ΔCt and RefFinder. From the study, we endorse MACP, CYP, EIF4α to be most stable candidate reference genes in all 'tissues', whereas PT, TFIID, MACP ranked high across genotypes, β-TUB, CYP, ELF1α were found to be best under abiotic stresses and 'all samples set'. The study recommends using minimum of two reference genes for RT-qPCR data normalizations in finger millet. All in all, CYP, β-TUB, and EF1α, in combination were found to be best for robust normalizations under most experimental conditions. The best and the least stable genes were validated for confirmation by assessing their appropriateness for normalization studies using EcNAC1 gene. The report provides the first comprehensive list of suitable stable candidate reference genes for nutritional rich cereal finger millet that will be advantageous to gene expression studies in this crop.

Validation of Internal Control Genes for Quantitative Real-Time PCR Gene Expression Analysis in Morchella

The reliability of qRT-PCR results depend on the stability of reference genes used for normalization, suggesting the necessity of identification of reference genes before gene expression analysis. Morels are edible mushrooms well-known across the world and highly prized by many culinary kitchens. Here, several candidate genes were selected and designed according to the Morchella importuna transcriptome data. The stability of the candidate genes was evaluated with geNorm and NormFinder under three different experimental conditions, and several genes with excellent stability were selected. The extensive adaptability of the selected genes was tested in ten Morchella species. Results from the three experimental conditions revealed that ACT1 and INTF7 were the most prominent genes in Morchella, CYC3 was the most stable gene in different development stages, INTF4/AEF3 were the top-ranked genes across carbon sources, while INTF3/CYC3 pair showed the robust stability for temperature stress treatment. We suggest using ACT1, AEF3, CYC3, INTF3, INTF4 and INTF7 as reference genes for gene expression analysis studies for any of the 10 Morchella strains tested in this study. The stability and practicality of the gene, vacuolar protein sorting (INTF3), vacuolar ATP synthase (INTF4) and14-3-3 protein (INTF7) involving the basic biological processes were validated for the first time as the candidate reference genes for quantitative PCR. Furthermore, the stability of the reference genes was found to vary under the three different experimental conditions, indicating the importance of identifying specific reference genes for particular conditions.

Selection of Reliable Reference Genes for RT-qPCR Analysis of Bursaphelenchus mucronatus Gene Expression From Different Habitats and Developmental Stages

Quantitative reverse transcription polymerase chain reaction (RT-qPCR), a sensitive technique for gene expression analysis, depends on the stability of the reference genes used for data normalization under different experimental conditions. Bursaphelenchus mucronatus, a pine-parasitic nematode varying in virulence, is widely distributed in natural pine forests throughout the northern hemisphere, but has not been investigated with respect to the identification of reference genes suitable for the normalization of RT-qPCR data. In the present study, eight candidate reference genes were analyzed in B. mucronatus under different habitat conditions and at different developmental stages. The expression stability of these genes was assessed by geNorm, NormFinder, BestKeeper, delta Cq, and RefFinder algorithms. In general, our results identified encoding beta-tubulin as the most stable gene. Moreover, pairwise analysis showed that three reference genes were sufficient to normalize the gene expression data under each set of conditions, with genes encoding beta-tubulin, 18S ribosomal RNA and ubiquitin-conjugating enzyme being the most suitable reference genes for different habitat conditions, whereas genes encoding beta-tubulin, histone, and 18S ribosomal RNA exhibited the most stable expression at different developmental stages. Validation of the selected reference genes was performed by profiling the expression of the fatty acid- and retinol-binding protein gene in different habitats, and by profiling the expression of the arginine kinase gene at different developmental stages. This first systematic analysis for the selection of suitable reference genes for RT-qPCR in B. mucronatus will facilitate future functional analyses and deep mining of genetic resources in this nematode.

Evaluation of suitable reference genes for qRT-PCR normalization in strawberry (Fragaria × ananassa) under different experimental conditions

Background

Strawberry has received much attention due to its nutritional value, unique flavor, and attractive appearance. The availability of the whole genome sequence and multiple transcriptome databases allows the great possibility to explore gene functions, comprehensively. Gene expression profiles of a target gene can provide clues towards the understanding of its biological function. Quantitative real-time PCR (qRT-PCR) is a preferred method for rapid quantification of gene expression. The accuracy of the results obtained by this method requires the reference genes with consistently stable expression to normalize its data.

Results

In present study, the expression stability of seven candidate reference genes in diverse sample subsets of different tissues and fruit developmental stages, and plant subjected to light quality and low temperature treatments was evaluated using three statistical algorithms, geNorm, NormFinder, and BestKeeper. Our data indicated that the expression stability of reference genes varied under different experimental conditions. Overall, DBP, HISTH4, ACTIN1 and GAPDH expressed much more stably. PIRUV, ACTIN2 and 18S were not recommended for normalization in given experimental conditions due to low stability. In addition, the relative expression pattern of HY5 (ELONGATED HYPOCOTYL5) was conducted to further confirm the reliability of the reference genes, which demonstrated the correct adoption of reference genes was of great importance in qRT-PCR analysis.

Conclusions

Expression stability of reference genes from strawberry varied across selected experimental conditions. Systematic validation of reference genes prior to calculation of target gene expression level should be done to improve the accuracy and consistency of qRT-PCR analysis.

Electronic supplementary material

The online version of this article (10.1186/s12867-018-0109-4) contains supplementary material, which is available to authorized users.

Comprehensive evaluation of blood-brain barrier-forming micro-vasculatures: Reference and marker genes with cellular composition

Primary brain microvessels (BrMV) maintain the cellular characters and molecular signatures as displayed in vivo, and serve as a vital tool for biomedical research of the blood-brain barrier (BBB) and the development/optimization of brain drug delivery. The variations of relative purities or cellular composition among different BrMV samples may have significant consequences in data interpretation and research outcome, especially for experiments with high-throughput genomics and proteomics technologies. In this study, we aimed to identify suitable reference gene (RG) for accurate normalization of real-time RT-qPCR analysis, and determine the proper marker genes (MG) for relative purity assessment in BrMV samples. Out of five housekeeping genes, β-actin was selected as the most suitable RG that was validated by quantifying mRNA levels of alpha-L-iduronidase in BrMV isolated from mice with one or two expressing alleles. Four marker genes highly/selectively expressed in BBB-forming capillary endothelial cells were evaluated by RT-qPCR for purity assessment, resulting in Cldn5 and Pecam1 as most suitable MGs that were further confirmed by immunofluorescent analysis of cellular components. Plvap proved to be an indicator gene for the presence of fenestrated vessels in BrMV samples. This study may contribute to the building blocks toward overarching research needs on the blood-brain barrier.

Suppression of a single BAHD gene in Setaria viridis causes large, stable decreases in cell wall feruloylation and increases biomass digestibility

Feruloylation of arabinoxylan (AX) in grass cell walls is a key determinant of recalcitrance to enzyme attack, making it a target for improvement of grass crops, and of interest in grass evolution. Definitive evidence on the genes responsible is lacking so we studied a candidate gene that we identified within the BAHD acyl-CoA transferase family. We used RNA interference (RNAi) silencing of orthologs in the model grasses Setaria viridis (SvBAHD01) and Brachypodium distachyon (BdBAHD01) and determined effects on AX feruloylation. Silencing of SvBAHD01 in Setaria resulted in a c. 60% decrease in AX feruloylation in stems consistently across four generations. Silencing of BdBAHD01 in Brachypodium stems decreased feruloylation much less, possibly due to higher expression of functionally redundant genes. Setaria SvBAHD01 RNAi plants showed: no decrease in total lignin, approximately doubled arabinose acylated by p-coumarate, changes in two-dimensional NMR spectra of unfractionated cell walls consistent with biochemical estimates, no effect on total biomass production and an increase in biomass saccharification efficiency of 40-60%. We provide the first strong evidence for a key role of the BAHD01 gene in AX feruloylation and demonstrate that it is a promising target for improvement of grass crops for biofuel, biorefining and animal nutrition applications.

Copper toxicity and date palm (Phoenix dactylifera) seedling tolerance: Monitoring of related biomarkers

Date palm (Phoenix dactylifera) seeds were exposed to different copper (Cu) solutions to examine plant stress responses. Low Cu concentrations (0.02 and 0.2 mM) caused an increase of seed germination, whereas higher Cu amounts (2 mM) significantly inhibited seed germination, delayed hypocotyl elongation, increased seedling mortality, and reduced the germination index by more than 90%. Metal-related toxicity symptoms appeared after 15 d of 2 mM of Cu exposure. Biochemical activities such as amylase activity and redox balance elements were examined to study the relationship between external Cu amount and internal plant response. The present study showed that amylolytic activity was dose- and time-dependent. Likewise, H₂ O₂ production increased after exposure to Cu, which was correlated with thiobarbituric acid reactive substance (TBARS) accumulation. Furthermore at low Cu concentrations, superoxide dismutase (SOD) and catalase (CAT) activities increased, suggesting that date palm seed stimulated its metal homeostasis networks. However, the highest cupric ion amounts increased cell oxidant accumulation and reduced enzyme production. Gene expression level measures of P. dactylifera phytochelatin synthase (Pdpcs) and P. dactylifera metallothionein (Pdmt) encoding genes have been carried out to investigate the implication of PdPCS and PdMT proteins in Cu homeostasis and/or its sequestration. Phoenix dactylifera metallothionein induction reached a peak after 30 d of exposure to 0.2 mM of Cu. However, it was down-regulated in plants exposed to higher Cu concentrations. In the same conditions, Pdpcs was overexpressed during 1 mo of exposure before it decreased thereafter. These observations provide a new insight into date palm cell response to Cu, a metal that can be toxic but that is also an essential element. Environ Toxicol Chem 2018;37:797-806. © 2017 SETAC.

SiYGL2 Is Involved in the Regulation of Leaf Senescence and Photosystem II Efficiency in Setaria italica (L.) P. Beauv

A yellow-green leaf mutant was isolated from EMS-mutagenized lines of Setaria italica variety Yugu1. Map-based cloning revealed the mutant gene is a homolog of Arabidopsis thaliana AtEGY1. EGY1 (ethylene-dependent gravitropism-deficient and yellow-green 1) is an ATP-independent metalloprotease (MP) that is required for chloroplast development, photosystem protein accumulation, hypocotyl gravitropism, leaf senescence, and ABA signal response in A. thaliana. However, the function of EGY1 in monocotyledonous C₄ plants has not yet been described. The siygl2 mutant is phenotypically characterized by chlorotic organs, premature senescence, and damaged PS II function. Sequence comparisons of the AtEGY1 and SiYGL2 proteins reveals the potential for SiYGL2 to encode a partially functional protein. Phenotypic characterization and gene expression analysis suggested that SiYGL2 participates in the regulation of chlorophyll content, leaf senescence progression, and PS II function. Additionally, our research will contribute to further characterization of the mechanisms regulating leaf senescence and photosynthesis in S. italica, and in C₄ plants in general.

SiSTL2 Is Required for Cell Cycle, Leaf Organ Development, Chloroplast Biogenesis, and Has Effects on C4 Photosynthesis in Setaria italica (L.) P. Beauv

Deoxycytidine monophosphate deaminase (DCD) is a key enzyme in the de novo dTTP biosynthesis pathway. Previous studies have indicated that DCD plays key roles in the maintenance of the balance of dNTP pools, cell cycle progression, and plant development. However, few studies have elucidated the functions of the DCD gene in Panicoideae plants. Setaria has been proposed as an ideal model of Panicoideae grasses, especially for C₄ photosynthesis research. Here, a Setaria italica stripe leaf mutant (sistl2) was isolated from EMS-induced lines of "Yugu1," the wild-type parent. The sistl2 mutant exhibited semi-dwarf, striped leaves, abnormal chloroplast ultrastructure, and delayed cell cycle progression compared with Yugu1. High-throughput sequencing and map-based cloning identified the causal gene SiSTL2, which encodes a DCD protein. The occurrence of a single-base G to A substitution in the fifth intron introduced alternative splicing, which led to the early termination of translation. Further physiological and transcriptomic investigation indicated that SiSTL2 plays an essential role in the regulation of chloroplast biogenesis, cell cycle, and DNA replication, which suggested that the gene has conserved functions in both foxtail millet and rice. Remarkably, in contrast to DCD mutants in C₃ rice, sistl2 showed a significant reduction in leaf cell size and affected C₄ photosynthetic capacity in foxtail millet. qPCR showed that SiSTL2 had a similar expression pattern to typical C₄ genes in response to a low CO₂ environment. Moreover, the loss of function of SiSTL2 resulted in a reduction of leaf ¹³C content and the enrichment of DEGs in photosynthetic carbon fixation. Our research provides in-depth knowledge of the role of DCD in the C₄ photosynthesis model S. italica and proposed new directions for further study of the function of DCD.

Stability evaluation of reference genes for gene expression analysis by RT-qPCR in soybean under different conditions

Real-time quantitative reverse transcription PCR is a sensitive and widely used technique to quantify gene expression. To achieve a reliable result, appropriate reference genes are highly required for normalization of transcripts in different samples. In this study, 9 previously published reference genes (60S, Fbox, ELF1A, ELF1B, ACT11, TUA5, UBC4, G6PD, CYP2) of soybean [Glycine max (L.) Merr.] were selected. The expression stability of the 9 genes was evaluated under conditions of biotic stress caused by infection with soybean mosaic virus, nitrogen stress, across different cultivars and developmental stages. ΔCt and geNorm algorithms were used to evaluate and rank the expression stability of the 9 reference genes. Results obtained from two algorithms showed high consistency. Moreover, results of pairwise variation showed that two reference genes were sufficient to normalize the expression levels of target genes under each experimental setting. For virus infection, ELF1A and ELF1B were the most stable reference genes for accurate normalization. For different developmental stages, Fbox and G6PD had the highest expression stability between two soybean cultivars (Tanlong No. 1 and Tanlong No. 2). ELF1B and ACT11 were identified as the most stably expressed reference genes both under nitrogen stress and among different cultivars. The results showed that none of the candidate reference genes were uniformly expressed at different conditions, and selecting appropriate reference genes was pivotal for gene expression studies with particular condition and tissue. The most stable combination of genes identified in this study will help to achieve more accurate and reliable results in a wide variety of samples in soybean.

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.

Selection and Validation of Appropriate Reference Genes for qRT-PCR Analysis in Isatis indigotica Fort

Due to its sensitivity and specificity, real-time quantitative PCR (qRT-PCR) is a popular technique for investigating gene expression levels in plants. Based on the Minimum Information for Publication of Real-Time Quantitative PCR Experiments (MIQE) guidelines, it is necessary to select and validate putative appropriate reference genes for qRT-PCR normalization. In the current study, three algorithms, geNorm, NormFinder, and BestKeeper, were applied to assess the expression stability of 10 candidate reference genes across five different tissues and three different abiotic stresses in Isatis indigotica Fort. Additionally, the IiYUC6 gene associated with IAA biosynthesis was applied to validate the candidate reference genes. The analysis results of the geNorm, NormFinder, and BestKeeper algorithms indicated certain differences for the different sample sets and different experiment conditions. Considering all of the algorithms, PP2A-4 and TUB4 were recommended as the most stable reference genes for total and different tissue samples, respectively. Moreover, RPL15 and PP2A-4 were considered to be the most suitable reference genes for abiotic stress treatments. The obtained experimental results might contribute to improved accuracy and credibility for the expression levels of target genes by qRT-PCR normalization in I. indigotica.

Validation of optimal reference genes for quantitative real time PCR in muscle and adipose tissue for obesity and diabetes research

The global incidence of obesity has led to an increasing need for understanding the molecular mechanisms that drive this epidemic and its comorbidities. Quantitative real-time RT-PCR (RT-qPCR) is the most reliable and widely used method for gene expression analysis. The selection of suitable reference genes (RGs) is critical for obtaining accurate gene expression information. The current study aimed to identify optimal RGs to perform quantitative transcriptomic analysis based on RT-qPCR for obesity and diabetes research, employing in vitro and mouse models, and human tissue samples. Using the ReFinder program we evaluated the stability of a total of 15 RGs. The impact of choosing the most suitable RGs versus less suitable RGs on RT-qPCR results was assessed. Optimal RGs differed between tissue and cell type, species, and experimental conditions. By employing different sets of RGs to normalize the mRNA expression of peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC1α), we show that sub-optimal RGs can markedly alter the PGC1α gene expression profile. Our study demonstrates the importance of validating RGs prior to normalizing transcriptional expression levels of target genes and identifies optimal RG pairs for reliable RT-qPCR normalization in cells and in human and murine muscle and adipose tissue for obesity/diabetes research.

Overexpression of BdMATE Gene Improves Aluminum Tolerance in Setaria viridis

Acidic soils are distributed worldwide, predominantly in tropical and subtropical areas, reaching around 50% of the arable soil. This type of soil strongly reduces crop production, mainly because of the presence of aluminum, which has its solubility increased at low pH levels. A well-known physiological mechanism used by plants to cope with Al stress involves activation of membrane transporters responsible for organic acid anions secretion from the root apex to the rhizosphere, which chelate Al, preventing its absorption by roots. In sorghum, a membrane transporter gene belonging to multidrug and toxic compound extrusion (MATE) family was identified and characterized as an aluminum-activated citrate transporter gene responsible for Al tolerance in this crop. Setaria viridis is an emerging model for C4 species and it is an important model to validate some genes for further C4 crops transformation, such as sugarcane, maize, and wheat. In the present work, Setaria viridis was used as a model plant to overexpress a newly identified MATE gene from Brachypodium distachyon (BdMATE), closely related to SbMATE, for aluminum tolerance assays. Transgenic S. viridis plants overexpressing a BdMATE presented an improved Al tolerance phenotype, characterized by sustained root growth and exclusion of aluminum from the root apex in transgenic plants, as confirmed by hematoxylin assay. In addition, transgenic plants showed higher root citrate exudation into the rhizosphere, suggesting that Al tolerance improvement in these plants could be related to the chelation of the metal by the organic acid anion. These results suggest that BdMATE gene can be used to transform C4 crops of economic importance with improved aluminum tolerance.