Evaluation of candidate reference genes for normalization of quantitative RT-PCR in soybean tissues under various abiotic stress conditions

Meta-analysis of plant growth-promoting rhizobacteria interaction with host plants: implications for drought stress response gene expression

Introduction

The molecular and physiological mechanisms activated in plants during drought stress tolerance are regulated by several key genes with both metabolic and regulatory roles. Studies focusing on crop gene expression following plant growth-promoting rhizobacteria (PGPR) inoculation may help understand which bioinoculant is closely related to the induction of abiotic stress responses.

Methods

Here, we performed a meta-analysis following Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines to summarise information regarding plant-PGPR interactions, focusing on the regulation of nine genes involved in plant drought stress response. The literature research yielded 3,338 reports, of which only 41 were included in the meta-analysis based on the chosen inclusion criteria. The meta-analysis was performed on four genes (ACO, APX, ACS and DREB2); the other five genes (ERD15, MYB, MYC, acdS, WRKY) had an insufficient number of eligible articles.

Results

Forest plots obtained through each meta-analysis showed that the overexpression of ACO, APX, ACS and DREB2 genes was not statistically significant. Unlike the other genes, DREB2 showed statistically significant results in both the presence and absence of PGPR. Considering I2>75 %, the results showed a high heterogeneity among the studies included, and the cause for this was examined using subgroup analysis. Moreover, the funnel plot and Egger's test showed that the analyses were affected by strong publication bias.

Discussion

This study argues that the presence of PGPR may not significantly influence the expression of drought stress response-related crop genes. This finding may be due to high heterogeneity, lack of data on the genes examined, and significant publication bias.

H2O2-dependent oxidation of the transcription factor GmNTL1 promotes salt tolerance in soybean

Reactive oxygen species (ROS) play an essential role in plant growth and responses to environmental stresses. Plant cells sense and transduce ROS signaling directly via hydrogen peroxide (H2O2)-mediated posttranslational modifications (PTMs) on protein cysteine residues. Here, we show that the H2O2-mediated cysteine oxidation of NAC WITH TRANS-MEMBRANE MOTIF1-LIKE 1 (GmNTL1) in soybean (Glycine max) during salt stress promotes its release from the endoplasmic reticulum (ER) membrane and translocation to the nucleus. We further show that an oxidative posttranslational modification on GmNTL1 residue Cys-247 steers downstream amplification of ROS production by binding to and activating the promoters of RESPIRATORY BURST OXIDASE HOMOLOG B (GmRbohB) genes, thereby creating a feed-forward loop to fine-tune GmNTL1 activity. In addition, oxidation of GmNTL1 Cys-247 directly promotes the expression of CATION H+ EXCHANGER 1 (GmCHX1)/SALT TOLERANCE-ASSOCIATED GENE ON CHROMOSOME 3 (GmSALT3) and Na+/H+ Antiporter 1 (GmNHX1). Accordingly, transgenic overexpression of GmNTL1 in soybean increases the H2O2 levels and K+/Na+ ratio in the cell, promotes salt tolerance, and increases yield under salt stress, while an RNA interference-mediated knockdown of GmNTL1 elicits the opposite effects. Our results reveal that the salt-induced oxidation of GmNTL1 promotes its relocation and transcriptional activity through an H2O2-mediated posttranslational modification on cysteine that improves resilience of soybean against salt stress.

Transcriptomic analysis of humic acid in relieving the inhibitory effect of high nitrogen on soybean nodulation

Introduction

Nitrogen fertilizer intake promotes soybean growth before the formation of nodules, but excess nitrogen has an inhibitory effect on soybean nodulation. It is important to balance nitrogen levels to meet both growth and nodulation needs.

Methods

the nitrogen level suitable for soybean growth and nodulation was studied, the role of humic acid (HA) in alleviating the inhibition of high nitrogen on soybean nodulation was analyzed, and transcriptomic analysis was performed to understand its mechanism.

Results

The results showed that a lower level of nitrogen with 36.4 mg urea per pot could increase the number of nodules of soybean, and a higher level of nitrogen with 145.9 mg urea per pot (U4 group) had the best growth indicators but inhibited nodulation significantly. HA relieved the inhibitory effect at high nitrogen level, and the number of nodules increased by 122.1% when 1.29 g HA was added (H2 group) compared with the U4 group. The transcriptome analysis was subsequently performed on the H2 and U4 groups, showing that there were 2995 differentially expressed genes (DEGs) on the 25th day, accounting for 6.678% of the total annotated genes (44,848) under the test conditions. These DEGs were enriched in mitogen-activated protein kinase signaling pathway-plant, flavonoid biosynthesis, and plant hormone signal transduction based on the -log10 (P adjusted) value in the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG).

Discussion

HA balanced the nitrogen level through the above pathways in soybean planting to control the number of nodules.

KARRIKIN UPREGULATED F-BOX 1 negatively regulates drought tolerance in Arabidopsis

The karrikin (KAR) receptor and several related signaling components have been identified by forward genetic screening, but only a few studies have reported on upstream and downstream KAR signaling components and their roles in drought tolerance. Here, we characterized the functions of KAR UPREGULATED F-BOX 1 (KUF1) in drought tolerance using a reverse genetics approach in Arabidopsis (Arabidopsis thaliana). We observed that kuf1 mutant plants were more tolerant to drought stress than wild-type (WT) plants. To clarify the mechanisms by which KUF1 negatively regulates drought tolerance, we performed physiological, transcriptome, and morphological analyses. We found that kuf1 plants limited leaf water loss by reducing stomatal aperture and cuticular permeability. In addition, kuf1 plants showed increased sensitivity of stomatal closure, seed germination, primary root growth, and leaf senescence to abscisic acid (ABA). Genome-wide transcriptome comparisons of kuf1 and WT rosette leaves before and after dehydration showed that the differences in various drought tolerance-related traits were accompanied by differences in the expression of genes associated with stomatal closure (e.g. OPEN STOMATA 1), lipid and fatty acid metabolism (e.g. WAX ESTER SYNTHASE), and ABA responsiveness (e.g. ABA-RESPONSIVE ELEMENT 3). The kuf1 mutant plants had higher root/shoot ratios and root hair densities than WT plants, suggesting that they could absorb more water than WT plants. Together, these results demonstrate that KUF1 negatively regulates drought tolerance by modulating various physiological traits, morphological adjustments, and ABA responses and that the genetic manipulation of KUF1 in crops is a potential means of enhancing their drought tolerance.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

Reference Genes and Expression Analysis of Seed Desaturases Genes in Soybean Mutant Accessions

Soybean oil is the second most-produced vegetable oil worldwide. To enhance the nutritional quality and oxidative stability of soybean oil, many soybean breeding programs are trying to increase oleic acid content and reduce linoleic and linolenic acid contents. The fatty acid profile of soybean oil is controlled by many genes, including those which code for omega-3 and omega-6 desaturases. Mutations in GmFAD2-1 and GmFAD3 genes are widely studied and their combinations can produce soybean oil with high oleic and low linoleic and linolenic content. However, few studies evaluate the effect of these mutations on gene expression. Therefore, the present study sought to identify reference genes, evaluate the expression of GmFAD2-1 and GmFAD3 seed desaturase genes in thirteen wild-type and mutated soybean accessions, and associate the expression patterns with fatty acid composition and with the GmFAD2-1 and GmFAD3 genotypes. GmCONS7 and GmUKN2 were identified as the best reference genes for combined use to normalize data. The GmFAD2-1A mutation of PI603452 accession was associated with a decrease in gene expression of GmFAD2-1A; however, downregulation may not be due to the truncated enzyme structure alone. These results suggested that there are factors other than GmFAD2-1A and GmFAD2-1B that have a considerable effect on oleic content, at least in soybeans with mutations in these two genes.

Genome-Wide Identification of Wild Soybean Mitochondrial Calcium Uniporter Family Genes and Their Responses to Cold and Carbonate Alkaline Stresses

The mitochondrial calcium uniporter (MCU), as an important component of the Ca²⁺ channel uniporter complex, plays a regulatory role in intracellular Ca²⁺ signal transduction. However, only a few studies to date have investigated plant MCU genes. In this study, we identified the MCU family genes in wild soybean and investigated their expression under cold and carbonate alkaline stresses. Eleven Glycine soja MCU genes (GsMCUs) were identified and clustered into two subgroups (subgroups I and II), and subgroup II could be further divided into two branches (MCU5 and MCU6). A total of 21 pairs of GsMCUs were characterized as duplicated genes, and displayed a similar exon-intron architecture. All GsMCU proteins contained one conserved MCU domain, within which two transmembrane domains were found. An analysis of the conserved motifs further supported that the GsMCUs showed high conservation in protein sequence and structure. Moreover, we found that all GsMCUs were expressed ubiquitously in different tissues and organs, and GsMCUs from the same subgroup displayed varied tissue expression profiles. In addition, based on RNA-seq and qRT-PCR assays, six and nine GsMCUs were differentially expressed under cold and carbonate alkaline stress, respectively. Promoter analysis also uncovered the existence of two canonical cold-related cis-acting elements, LTR and DRE/CRT, as well as stress-related phytohormone-responsive elements. Our results provide valuable information about the MCU family in soybean responses to cold and carbonate alkaline stress, which will be helpful in further characterizing their biological roles in response to abiotic stress.

Comprehensive Profiling of Tubby-Like Proteins in Soybean and Roles of the GmTLP8 Gene in Abiotic Stress Responses

Tubby-like proteins (TLPs) are transcription factors that are widely present in eukaryotes and generally participate in growth and developmental processes. Using genome databases, a total of 22 putative TLP genes were identified in the soybean genome, and unevenly distributed across 13 chromosomes. Phylogenetic analysis demonstrated that the predicted GmTLP proteins were divided into five groups (I-V). Gene structure, protein motifs, and conserved domains were analyzed to identify differences and common features among the GmTLPs. A three-dimensional protein model was built to show the typical structure of TLPs. Analysis of publicly available gene expression data showed that GmTLP genes were differentially expressed in response to abiotic stresses. Based on those data, GmTLP8 was selected to further explore the role of TLPs in soybean drought and salt stress responses. GmTLP8 overexpressors had improved tolerance to drought and salt stresses, whereas the opposite was true of GmTLP8-RNAi lines. 3,3-diaminobenzidine and nitro blue tetrazolium staining and physiological indexes also showed that overexpression of GmTLP8 enhanced the tolerance of soybean to drought and salt stresses; in addition, downstream stress-responsive genes were upregulated in response to drought and salt stresses. This study provides new insights into the function of GmTLPs in response to abiotic stresses.

Metabolic adjustment and regulation of gene expression are essential for increased resistance to severe water deficit and resilience post-stress in soybean

Background

Soybean is the main oilseed crop grown in the world; however, drought stress affects its growth and physiology, reducing its yield. The objective of this study was to characterize the physiological, metabolic, and genetic aspects that determine differential resistance to water deficit in soybean genotypes.

Methods

Three soybean genotypes were used in this study, two lineages (L11644 and L13241), and one cultivar (EMBRAPA 48-C48). Plants were grown in pots containing 8 kg of a mixture of soil and sand (2:1) in a greenhouse under sunlight. Soil moisture in the pots was maintained at field capacity until the plants reached the stage of development V4 (third fully expanded leaf). At this time, plants were subjected to three water treatments: Well-Watered (WW) (plants kept under daily irrigation); Water Deficit (WD) (withholding irrigation until plants reached the leaf water potential at predawn of -1.5 ± 0.2 MPa); Rewatered (RW) (plants rehydrated for three days after reached the water deficit). The WW and WD water treatments were evaluated on the eighth day for genotypes L11644 and C48, and on the tenth day for L13241, after interruption of irrigation. For the three genotypes, the treatment RW was evaluated after three days of resumption of irrigation. Physiological, metabolic and gene expression analyses were performed.

Results

Water deficit inhibited growth and gas exchange in all genotypes. The accumulation of osmolytes and the concentrations of chlorophylls and abscisic acid (ABA) were higher in L13241 under stress. The metabolic adjustment of lineages in response to WD occurred in order to accumulate amino acids, carbohydrates, and polyamines in leaves. The expression of genes involved in drought resistance responses was more strongly induced in L13241. In general, rehydration provided recovery of plants to similar conditions of control treatment. Although the C48 and L11644 genotypes have shown some tolerance and resilience responses to severe water deficit, greater efficiency was observed in the L13241 genotype through adjustments in morphological, physiological, genetic and metabolic characteristics that are combined in the same plant. This study contributes to the advancement in the knowledge about the resistance to drought in cultivated plants and provides bases for the genetic improvement of the soybean culture.

Proteomic and lipidomics analyses of high fatty acid AhDGAT3 transgenic soybean reveals the key lipase gene associated with the lipid internal mechanism

Vegetable oil is one of the most important components of human nutrition. Soybean (Glycine max) is an important oil crop worldwide and contains rich unsaturated fatty acids. Diacylglycerol acyltransferase (DGAT) is a key rate-limiting enzyme in the Kennedy pathway from diacylglycerol (DAG) to triacylglycerol (TAG). In this study, we conducted further research using T3 AhDGAT3 transgenic soybean. A high-performance gas chromatography flame ionization detector showed that oleic acid (18:1) content and total fatty acid content of transgenic soybean were significantly higher than those of the wild type (WT). However, linoleic acid (18:2) was much lower than that in the WT. For further mechanistic studies, 20 differentially expressed proteins (DEPs) and 119 differentially expressed metabolites (DEMs) were identified between WT (JACK) and AhDGAT3 transgenic soybean mature seeds using proteomic and lipidomics analyses. Combined proteomic and lipidomics analyses showed that the upregulation of the key DEP (lipase GDSL domain-containing protein) in lipid transport and metabolic process induced an increase in the total fatty acid and 18:1 composition, but a decrease in the 18:2 composition of fatty acids. Our study provides new insights into the deep study of molecular mechanism underlying the enhancement of fatty acids in transgenic soybeans, especially oleic acid and total fatty acid, which are enhanced by over-expression of AhDGAT3.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Analysis of Development and Tissue-Dependent Flower Color Formation in Cymbidium lowianum

The development and tissue-dependent color formation of the horticultural plant results in various color pattern flowers. Anthocyanins and carotenoids contribute to the red and yellow colors, respectively. In this study, quantitative real-time polymerase chain reaction (qRT-PCR) is used to analyze the expression profiles of anthocyanin and carotenoids biosynthesis genes in Cymbidium lowianum (Rchb.f.) Rchb.f. Appropriate reference gene selection and validation are required before normalization of gene expression in qRT-PCR analysis. Thus, we firstly selected 12 candidate reference genes from transcriptome data, and used geNorm and Normfinder to evaluate their expression stability in lip (divided into abaxial and adaxial), petal, and sepal of the bud and flower of C. lowianum. Our results show that the two most stable reference genes in different tissues of C. lowianum bud and flower are EF1δ and 60S, the most unstable reference gene is 26S. The expression profiles of the CHS and BCH genes were similar to FPKM value profiles after normalization to the two most stable reference genes, EF1δ and 60S, with the upregulated CHS and BCH expression in flower stage, indicating that the ABP and CBP were activated across the stages of flower development. However, when the most unstable reference gene, 26S, was used to normalize the qRT-PCR data, the expression profiles of CHS and BCH differed from FPKM value profiles, indicating the necessity of selecting stable reference genes. Moreover, CHS and BCH expression was highest in the abaxial lip and adaxial lip, respectively, indicating that the ABP and CBP were activated in abaxial and adaxial lip, respectively, resulting in a presence of red or yellow segments in abaxial and adaxial lip. This study is the first to provide reference genes in C. lowianum, and also provide useful information for studies that aim to understand the molecular mechanisms of flower color formation in C. lowianum.

Identification and validation of reference genes for qRT-PCR based studies in horse gram (Macrotyloma uniflorum)

The quantitative real-time polymerase chain reaction (qRT-PCR) is the most sensitive and commonly used technique for gene expression studies in biological systems. However, the reliability of qRT-PCR results depends on the selection of reference gene(s) for data normalization. Horse gram (Macrotyloma uniflorum) is an important legume crop on which several molecular studies have been reported. However, the stability of reference genes has not been evaluated. In the present study, nine candidate reference genes were identified from horse gram RNA-seq data and evaluated in two horse gram genotypes, HPK4 and HPKM317 under six abiotic stresses viz. cold, drought, salinity, heat, abscisic acid and methyl viologen-induced oxidative stress. The results were evaluated using geNorm, Bestkeeper, Normfinder and delta-delta Ct methods and comprehensive ranking was assigned using RefFinder and RankAggreg software. The overall result showed that TCTP was one of the most stable genes in all samples and in genotype HPK4, while in HPKM317 profilin was most stably expressed. However, PSMA5 was identified as least stable in all the experimental conditions. Expression of target genes dehydrin and early response to dehydration 6 under drought stress was also validated using TCTP and profilin for data normalization, either alone or in combination, which confirmed their suitability for qRT-PCR data normalization. Thus, TCTP and profilin genes may be used for qRT-PCR data normalization for molecular and genomic studies in horse gram.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01104-0.

Quantitative proteomic and lipidomics analyses of high oil content GmDGAT1-2 transgenic soybean illustrate the regulatory mechanism of lipoxygenase and oleosin

KEY MESSAGE

Proteomic and lipidomics analyses of WT and GmDGAT1-2 transgenic soybeans showed that GmDGAT1-2 over-expression induced lipoxygenase down-regulatation and oleoin up-regulatation, which significantly changed the compositions and total fatty acid. The main goal of soybean breeding is to increase the oil content. Diacylglycerol acyltransferase (DGAT) is a key rate-limiting enzyme in fatty acid metabolism and may regulate oil content. Herein, 10 GmDGAT genes were isolated from soybean and transferred into wild-type (WT) Arabidopsis. The total fatty acid was 1.2 times higher in T3 GmDGAT1-2 transgenic Arabidopsis seeds than in WT. Therefore, GmDGAT1-2 was transferred into WT soybean (JACK), and four T3 transgenic soybean lines were obtained. The results of high-performance gas chromatography and Soxhlet extractor showed that, compared with those of JACK, oleic acid (18:1), and total fatty acid levels in transgenic soybean plants were much higher, but linoleic acid (18:2) was lower than WT. Palmitic acid (16:0), stearic acid (18:0), and linolenic acid (18:3) were not significantly different. For mechanistic studies, 436 differentially expressed proteins (DEPs) and 180 differentially expressed metabolites (DEMs) were identified between WT (JACK) and transgenic soybean pods using proteomic and lipidomics analyses. Four lipoxygenase proteins were down-regulated in linoleic acid metabolism while four oleosin proteins were up-regulated in the final oil formation. The results showed an increase in the total fatty acid and 18:1 composition, and a decrease in the 18:2 composition of fatty acid. Our study brings new insights into soybean genetic transformation and the deep study of molecular mechanism that changes the total fatty acid, 18:1, and 18:2 compositions in GmDGAT1-2 transgenic soybean.

Screening of stable internal reference gene of Quinoa under hormone treatment and abiotic stress

Real-time quantitative polymerase chain reaction is the most commonly used method to accurately detect gene expression patterns. The method requires stable internal reference genes to standardize the data. However, studies have shown that there is no stable expression of internal reference genes in different tissues and under different treatments. Therefore, in order to study the optimal reference genes of quinoa under different hormones and abiotic stress, leaves and stems from quinoa seedlings treated with low temperature (4 °C), salt (200 mmol/L) and abscisic acid (200 mmol/L) were used as experimental materials. Using ACT-1, eIF, EF1α, GAPDH, TUA, TUB-9, TUB-1, H2A and L8-1 as candidate reference genes, the expression stability of these 9 quinoa candidate reference genes under different hormone treatment and abiotic stress was evaluated by using geNorm, NormFinder and BestKeeper software. The results showed that TUB-1 gene under salt stress, L8-1 gene under low temperature stress, EF-1α gene induced by ABA. PLIM2c WLIM1and WLIM2b were selected to verify the candidate internal reference genes, and finally the expression of GAPDH was most unstable under the three treatments, which was not suitable to be the internal reference gene of quinoa under specific conditions, while EF1α showed good stability under the three different treatments and was suitable to be used as the internal reference gene. In conclusion, the results of this study could provide an important reference for quantifying the expression level of reference genes in quinoa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01094-z.

The Drought-Mediated Soybean GmNAC085 Functions as a Positive Regulator of Plant Response to Salinity

Abiotic stress factors, such as drought and salinity, are known to negatively affect plant growth and development. To cope with these adverse conditions, plants have utilized certain defense mechanisms involved in various aspects, including morphological, biochemical and molecular alterations. Particularly, a great deal of evidence for the biological importance of the plant-specific NAM, ATAF1/2, CUC2 (NAC) transcription factors (TFs) in plant adaptation to abiotic stress conditions has been reported. A previous in planta study conducted by our research group demonstrated that soybean (Glycine max) GmNAC085 mediated drought resistance in transgenic Arabidopsis plants. In this study, further characterization of GmNAC085 function in association with salt stress was performed. The findings revealed that under this condition, transgenic soybean plants overexpressing GmNAC085 displayed better germination rates than wild-type plants. In addition, biochemical and transcriptional analyses showed that the transgenic plants acquired a better defense system against salinity-induced oxidative stress, with higher activities of antioxidant enzymes responsible for scavenging hydrogen peroxide or superoxide radicals. Higher transcript levels of several key stress-responsive genes involved in the proline biosynthetic pathway, sodium ion transporter and accumulation of dehydrins were also observed, indicating better osmoprotection and more efficient ion regulation capacity in the transgenic lines. Taken together, these findings and our previous report indicate that GmNAC085 may play a role as a positive regulator in plant adaptation to drought and salinity conditions.

Reference gene identification for gene expression analysis in rice under different metal stress

Real-time quantitative polymerase chain reaction (RT-qPCR) is the most common approach to quantify changes in gene expression. Appropriate internal reference genes are essential for normalization of data of RT-qPCR. In the present study, we identified suitable reference genes for analysis of gene expression in rice seedlings subjected to different heavy metal stresses such as deficiencies of iron and zinc and toxicities of cobalt, cadmium and nickel. First, from publically available RNA-Seq data we identified 10 candidate genes having stable expression. We also included commonly used house-keeping gene OsUBQ5 (Ubiquitin 5) in our analysis. Expression stability of all the 11 genes was determined by two independent tools, NormFinder and geNorm. Our results show that selected candidate reference genes have higher stability in their expression compared to that of OsUBQ5. Genes with locus ID LOC_Os03g16690, encoding an oxysterol-binding protein (OsOBP) and LOC_Os01g56580, encoding Casein Kinase_1a.3 (OsCK1a.3) were identified to be the most stably expressed reference genes under most of the conditions tested. Finally, the study reveals that it is better to use a specific reference gene for a specific heavy metal stress condition rather than using a common reference gene for multiple heavy metal stress conditions. The reference genes identified here would be very useful for gene expression studies under heavy metal stresses in rice.

Identification and validation of new reference genes for accurate quantitative reverse transcriptase-PCR normalization in the Antarctic plant Colobanthus quitensis under abiotic stress conditions

The Antarctic ecotype of Colobanthus quitensis is a vascular plant highly adapted to the harsh environmental conditions of Maritime Antarctica which is now facing with the rapid local warming experienced in the Antarctic Peninsula during the last decades. Thus, the identification of the molecular mechanisms leading to the adaptation to this warming trend is a new target for modern cell physiology. The selection of suitable reference genes for quantification of key stress-responsive genes through quantitative Reverse Transcriptase-Polymerase Chain Reaction (qRT-PCR) is important to ensure accurate and reliable results. In this study, we evaluated the expression stability of eleven candidate genes in C. quitensis under different abiotic stress conditions using geNorm and RefFinder tools. The statistical analysis showed that the appropriate reference genes varied depending on the experimental conditions, even if EF1α and PP2Acs ranked as the most stable reference genes when all stress conditions were considered. To further validate the stability of the selected reference genes, the expression patterns of C. quitensis catalase gene (CqCAT) was analyzed. The reference genes validated in this study will be useful for improving the accuracy of qRT-PCR analysis for gene expression studies of the Antarctic ecotype of C. quitensis and could be extended to other ecotypes adapted to low temperatures.

Phosphate or nitrate imbalance induces stronger molecular responses than combined nutrient deprivation in roots and leaves of chickpea plants

The negative effects of phosphate (Pi) and/or nitrate (NO₃^- ) fertilizers on the environment have raised an urgent need to develop crop varieties with higher Pi and/or nitrogen use efficiencies for cultivation in low-fertility soils. Achieving this goal depends upon research that focuses on the identification of genes involved in plant responses to Pi and/or NO₃^- starvation. Although plant responses to individual deficiency in either Pi (-Pi/+NO₃^- ) or NO₃^- (+Pi/-NO₃^- ) have been separately studied, our understanding of plant responses to combined Pi and NO₃^- deficiency (-Pi/-NO₃^- ) is still very limited. Using RNA-sequencing approach, transcriptome changes in the roots and leaves of chickpea cultivated under -Pi/+NO₃^- , +Pi/-NO₃^- or -Pi/-NO₃^- conditions were investigated in a comparative manner. -Pi/-NO₃^- treatment displayed lesser effect on expression changes of genes related to Pi or NO₃^- transport, signalling networks, lipid remodelling, nitrogen and Pi scavenging/remobilization/recycling, carbon metabolism and hormone metabolism than -Pi/+NO₃^- or +Pi/-NO₃^- treatments. Therefore, the plant response to -Pi/-NO₃^- is not simply an additive result of plant responses to -Pi/+NO₃^- and +Pi/-NO₃^- treatments. Our results indicate that nutrient imbalance is a stronger stimulus for molecular reprogramming than an overall deficiency.

Medicago sativa and Medicago truncatula Show Contrasting Root Metabolic Responses to Drought

Drought is an environmental stressor that affects crop yield worldwide. Understanding plant physiological responses to stress conditions is needed to secure food in future climate conditions. In this study, we applied a combination of plant physiology and metabolomic techniques to understand plant responses to progressive water deficit focusing on the root system. We chose two legume plants with contrasting tolerance to drought, the widely cultivated alfalfa Medicago sativa (Ms) and the model legume Medicago truncatula (Mt) for comparative analysis. Ms taproot (tapR) and Mt fibrous root (fibR) biomass increased during drought, while a progressive decline in water content was observed in both species. Metabolomic analysis allowed the identification of key metabolites in the different tissues tested. Under drought, carbohydrates, abscisic acid, and proline predominantly accumulated in leaves and tapRs, whereas flavonoids increased in fibRs in both species. Raffinose-family related metabolites accumulated during drought. Along with an accumulation of root sucrose in plants subjected to drought, both species showed a decrease in sucrose synthase (SUS) activity related to a reduction in the transcript level of SUS1, the main SUS gene. This study highlights the relevance of root carbon metabolism during drought conditions and provides evidence on the specific accumulation of metabolites throughout the root system.

Integrative omic and transgenic analyses reveal the positive effect of ultraviolet-B irradiation on salvianolic acid biosynthesis through upregulation of SmNAC1

Salvianolic acids (SalAs), a group of secondary metabolites in Salvia miltiorrhiza, are widely used for treating cerebrovascular diseases. Their biosynthesis is modulated by a variety of abiotic factors, including ultraviolet-B (UV-B) irradiation; however, the underlying mechanisms remain largely unknown. Here, an integrated metabolomic, proteomic, and transcriptomic approach coupled with transgenic analyses was employed to dissect the mechanisms underlying UV-B irradiation-induced SalA biosynthesis. Results of metabolomics showed that 28 metabolites, including 12 SalAs, were elevated in leaves of UV-B-treated S. miltiorrhiza. Meanwhile, the contents of several phytohormones, including jasmonic acid and salicylic acid, which positively modulate the biosynthesis of SalAs, also increased in UV-B-treated S. miltiorrhiza. Consistently, 20 core biosynthetic enzymes and numerous transcription factors that are involved in SalA biosynthesis were elevated in treated samples as indicated by a comprehensive proteomic analysis. Correlation and gene expression analyses demonstrated that the NAC1 gene, encoding a NAC transcriptional factor, was positively involved in UV-B-induced SalA biosynthesis. Accordingly, overexpression and RNA interference of NAC1 increased and decreased SalA contents, respectively, through regulation of key biosynthetic enzymes. Furthermore, ChIP-qPCR and Dual-LUC assays showed that NAC1 could directly bind to the CATGTG and CATGTC motifs present in the promoters of the SalA biosynthesis-related genes PAL3 and TAT3, respectively, and activate their expression. Our results collectively demonstrate that NAC1 plays a crucial role in UV-B irradiation-induced SalA biosynthesis. Taken together, our findings provide mechanistic insights into the UV-B-induced SalA biosynthesis in S. miltiorrhiza, and shed light on a great potential for the development of SalA-abundant varieties through genetic engineering.

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.

ABAS1 from soybean is a 1R-subtype MYB transcriptional repressor that enhances ABA sensitivity

Transcription factors (TFs) help plants respond to environmental stresses by regulating gene expression. Up till now, studies on the MYB family of TFs have mainly focused on the highly abundant R2R3-subtype. While the less well-known 1R-subtype has been generally shown to enhance abscisic acid (ABA) sensitivity by acting as transcriptional activators, the mechanisms of their functions are unclear. Here we identified an ABA sensitivity-associated gene from soybean, ABA-Sensitive 1 (GmABAS1), of the 1R-subtype of MYB. Using the GFP-GmABAS1 fusion protein, we demonstrated that GmABAS1 is localized in the nucleus, and with yeast reporter systems, we showed that it is a transcriptional repressor. We then identified the target gene of GmABAS1 to be Glyma.01G060300, an annotated ABI five-binding protein 3 and showed that GmABAS1 binds to the promoter of Glyma.01G060300 both in vitro and in vivo. Furthermore, Glyma.01G060300 and GmABAS1 exhibited reciprocal expression patterns under osmotic stress, inferring that GmABAS1 is a transcriptional repressor of Glyma.01G060300. As a further confirmation, AtAFP2, an orthologue of Glyma.01G060300, was down-regulated in GmABAS1-transgenic Arabidopsis thaliana, enhancing the plant's sensitivity to ABA. This is the first time a 1R-subtype of MYB from soybean has been reported to enhance ABA sensitivity by acting as a transcriptional repressor.

Identification of soybean trans-factors associated with plastid RNA editing sites

RNA editing is a posttranscriptional process that changes nucleotide sequences, among which cytosine-to-uracil by a deamination reaction can revert non-neutral codon mutations. Pentatricopeptide repeat (PPR) proteins comprise a family of RNA-binding proteins, with members acting as editing trans-factors that recognize specific RNA cis-elements and perform the deamination reaction. PPR proteins are classified into P and PLS subfamilies. In this work, we have designed RNA biotinylated probes based in soybean plastid RNA editing sites to perform trans-factor specific protein isolation. Soybean cis-elements from these three different RNA probes show differences in respect to other species. Pulldown samples were submitted to mass spectrometry for protein identification. Among detected proteins, five corresponded to PPR proteins. More than one PPR protein, with distinct functional domains, was pulled down with each one of the RNA probes. Comparison of the soybean PPR proteins to Arabidopsis allowed identification of the closest homologous. Differential gene expression analysis demonstrated that the PPR locus Glyma.02G174500 doubled its expression under salt stress, which correlates with the increase of its potential rps14 editing. The present study represents the first identification of RNA editing trans-factors in soybean. Data also indicated that potential multiple trans-factors should interact with RNA cis-elements to perform the RNA editing.

Identification of candidate reference genes for quantitative RT-PCR in Miscanthus sinensis subjected to various abiotic stresses

Quantitative real-time PCR (qRT-PCR) has been widely used for studying gene expression at the transcript level. Its accuracy usually relies on the reference genes that are utilized for data normalization. Miscanthus sinensis, a perennial C4 grass with high biomass and strong resistance to adversities, is often utilized as a high value energy crop. However, no reliable reference genes have been investigated for normalizing gene expression for this species. In this study, 12 candidate reference genes were selected to identify their stability under five different abiotic stress treatments (drought, salt, cadmium, chromium and arsenic) by using geNorm, NormFinder, BestKeeper and RefFinder softwares. The results showed that 18S rRNA and Unigene33312 were the best reference genes under drought treatments. Unigene33312 and Unigene33024 were found to be the most stably expressed genes under salt stress and Cd stress. Moreover, Unigene33024 and PP2A were the most suitable reference genes under Cr stress and Unigene33024 and Sb09g019750 were deemed more suitable reference genes under As stress. In total, considering all the samples, Unigene33024 and PP2A were the most stable genes while ACTIN and Unigene26576 were the least stable reference genes for internal control. The expression patterns of two target genes (Cu/Zn SOD and CAT) were used to further verify those selected reference genes under different conditions. The results showed that the most and the least stable reference genes had clearly different expression patterns. This work comprehensively estimated the stability of reference genes in M. sinensis which may give insight to the reference genes selection in other tissues as well as other related varieties. These suggested reference genes would assist in further putative gene expression validation in M. sinensis.

Selection of Reference Genes for qPCR Analyses of Gene Expression in Ramie Leaves and Roots across Eleven Abiotic/Biotic Treatments

Quantitative real-time PCR (qPCR) is commonly used for deciphering gene functions. For effective qPCR analyses, suitable reference genes are needed for normalization. The objective of this study is to identify the appropriate reference gene(s) for qPCR analyses of the leaves and roots of ramie (Boehmeria nivea L.), an important natural fiber crop. To accomplish this goal, we investigated the expression patterns of eight common plant qPCR reference genes in ramie leaves and roots under five abiotic stresses, five hormonal treatments, and one biotic stress. The relative expression stabilities of the eight genes were evaluated using four common but different approaches: geNorm, NormFinder, BestKeeper, and RefFinder. Across the 11 tested conditions, ACT1 was the most stably expressed among the eight genes while GAPDH displayed the biggest variation. Overall, while variations in the suggested reference genes were found for different tissue x treatment combinations, our analyses revealed that together, genes ACT1, CYP2, and UBQ can provide robust references for gene expression studies of ramie leaves under most conditions, while genes EF-1α, TUB, and ACT1 can be used for similar studies of ramie roots. Our results should help future functional studies of the genes in ramie genome across tissues and environmental conditions.

A GmSIN1/GmNCED3s/GmRbohBs Feed-Forward Loop Acts as a Signal Amplifier That Regulates Root Growth in Soybean Exposed to Salt Stress

Abscisic acid (ABA) and reactive oxygen species (ROS) act as key signaling molecules in the plant response to salt stress; however, how these signals are transduced and amplified remains unclear. Here, a soybean (Glycine max) salinity-induced NAM/ATAF1/2/CUC2 (NAC) transcription factor encoded by SALT INDUCED NAC1 (GmSIN1) was shown to be a key component of this process. Overexpression of GmSIN1 in soybean promoted root growth and salt tolerance and increased yield under salt stress; RNA interference-mediated knockdown of GmSIN1 had the opposite effect. The rapid induction of GmSIN1 in response to salinity required ABA and ROS, and the effect of GmSIN1 on root elongation and salt tolerance was achieved by boosting cellular ABA and ROS contents. GmSIN1 upregulated 9-cis-epoxycarotenoid dioxygenase coding genes in soybean (GmNCED3s, associated with ABA synthesis) and Respiratory burst oxidase homolog B genes in soybean (GmRbohBs, associated with ROS generation) by binding to their promoters at a site that has not been described to date. Together, GmSIN1, GmNCED3s, and GmRbohBs constitute a positive feed-forward system that enables the rapid accumulation of ABA and ROS, effectively amplifying the initial salt stress signal. These findings suggest that the combined modulation of ABA and ROS contents enhances soybean salt tolerance.

Genome-Wide Analysis of the WOX Gene Family and Function Exploration of GmWOX18 in Soybean

WUSCHEL-related homeobox (WOX) is a family of transcription factors that are unique to plants and is characterized by the presence of a homeodomain. The WOX transcription factor plays an important role in regulating plant growth and development and the response to abiotic stress. Soybean is one of the most important oil crops worldwide. In this study, based on the available genome data of soybean, the WOX gene family was identified by bioinformatics analysis. The chromosome distribution, gene and protein structures, phylogenetic relationship and gene expression patterns of this family were comprehensively compared. The results showed that a total of 33 putative WOX genes in the soybean genome were found and then designated as GmWOX1- GmWOX33, which were distributed across 19 chromosomes except chromosome 16. Multiple sequence analysis of the GmWOX gene family revealed a highly conserved homeodomain. Phylogenetic tree analysis showed that 33 WOX genes could be divided into three major clades (modern/WUS, intermediate and ancient) in soybean. Of these 33 WOX genes, some showed differential expression patterns in the tested tissues (leaves, pods, unopen and open flowers, nodules, seed, roots, root hairs, stems, shoot apical meristems and shoot tips). In addition, the expression profile and qRT-PCR analysis showed that most of the GmWOX genes responded to different abiotic stress treatments (cold and drought). According to the expression pattern of GmWOX genes in the high regeneration capacity soybean material P3, overexpression of GmWOX18 was selected for function analysis. The overexpression of GmWOX18 increased the regeneration ability of clustered buds. The results will provide valuable information for further studies on the roles of WOX genes in regulating soybean growth, development and responses to abiotic stress, as well as a basis for the functional identification and analysis of WOX genes in soybean.

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.

Water Deficit Elicits a Transcriptional Response of Genes Governing d-pinitol Biosynthesis in Soybean (Glycine max)

d-pinitol is the most commonly accumulated sugar alcohol in the Leguminosae family and has been observed to increase significantly in response to abiotic stress. While previous studies have identified genes involved in d-pinitol synthesis, no study has investigated transcript expression in planta. The present study quantified the expression of several genes involved in d-pinitol synthesis in different plant tissues and investigated the accumulation of d-pinitol, myo-inositol and other metabolites in response to a progressive soil drought in soybean (Glycine max). Expression of myo-inositol 1-phosphate synthase (INPS), the gene responsible for the conversion of glucose-6-phosphate to myo-inositol-1-phosphate, was significantly up regulated in response to a water deficit for the first two sampling weeks. Expression of myo-inositol O-methyl transferase (IMT1), the gene responsible for the conversion of myo-inositol into d-ononitol was only up regulated in stems at sampling week 3. Assessment of metabolites showed significant changes in their concentration in leaves and stems. d-Pinitol concentration was significantly higher in all organs sampled from water deficit plants for all three sampling weeks. In contrast, myo-inositol, had significantly lower concentrations in leaf samples despite up regulation of INPS suggesting the transcriptionally regulated flux of carbon through the myo-inositol pool is important during water deficit.

Validation of reference genes for qRT-PCR data normalisation in lentil (Lens culinaris) under leaf developmental stages and abiotic stresses

Lentil (Lens culinaris) is one of the most important staple food crops of developing countries. Transcriptome based global gene expression profiling followed by validation of expression of important genes through quantitative real time-PCR (qRT-PCR) has achieved significance in recent years. However, there is a severe scarcity of information regarding stable reference genes in lentil, which is mandatory for qRT-PCR data normalisation. Hence, the present study was under-taken to identify the most stable reference gene(s) in lentil. Expression stability of eight candidate genes viz. ribulose 1,5-bisphosphate carboxylase large subunit (Rbcl), ribosomal protein L2 (RPL2), 18S rRNA, tubulin (Tub), elongation factor 1α (EF1α), glyceraldehydes-3-phosphate dehydrogenase (GAPDH), heat shock protein (HSP70), and Maturase (mat K) was evaluated in five varieties of lentil at three different stages of leaf development and abiotic stress conditions using qRT-PCR. The results were analysed using four types of statistical software viz., geNorm, BestKeeper, NormFinder and RefFinder; all softwares identified RPL2 as most stable under abiotic stress conditions and developmental stages followed by Tub and Rbcl; while, HSP70 was identified as least stable. Relative expression of the target genes, defensin and PR4, was evaluated under abiotic stress conditions and data normalisation was done using two stable reference genes, RPL2 and Tub, either alone or in combination and with two least stable genes, HSP70 and 18S. The present work provides a list of potential reference genes in lentil, which will help in selection of appropriate reference gene for qRT-PCR data normalization depending upon the experiment.

Identification of suitable internal control genes for transcriptional studies in Eleusine coracana under different abiotic stress conditions

Finger millet [Eleusine coracana (L.) Gaertn] is an excellent food and forage crop of arid and semiarid areas in Africa and Asia. It is well adapted to drought, heat, high salinity, poor soil fertility and low pH with an efficient C₄ carbon fixation mechanism for high yield potential. To normalize the target gene expression data, the identification of suitable reference genes is essential. Ten candidate reference genes were selected and their expression stability was analyzed in various samples treated with different abiotic stress conditions. Five different statistical algorithms: geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to determine the stability of these genes. Our results revealed GAPDH, EEF1a, ACT and CYC as highly stable reference genes and PP2A and eIF4A as least stable reference genes across all the samples and suggesting that these genes could be used for accurate transcript normalization under abiotic stress. To the best of our knowledge, this is the first report on identification of suitable reference genes for accurate transcript normalization using qRT-PCR in finger millet under abiotic stress.

Selection of reference genes for quantitative real-time PCR normalization in Narcissus pseudonarcissu in different cultivars and different organs

Quantitative real-time PCR (qRT-PCR) has been a widely used accurate technique for gene expression analysis in various species. However, its results require data normalization by reliable reference genes. Despite the horticultural importance of Narcissus pseudonarcissus, and genome sequence has become available for the species, no gene expression study based on the stability of reference genes in qRT-PCR has been conducted. To boost the use of qRT-PCR in N. pseudonarcissus, we uncovered eight commonly used candidate reference genes for their stability. The expression levels of the eight genes were detected for the normalization in five different organs (bulbs, scapes, leaves, perianths and coronas) of three N. pseudonarcissus cultivars (‘Marieke’, ‘Pinza’ and ‘Slim Whitman’) by qRT-PCR. Subsequently, three commonly used computational programs were applied for evaluating the stability of the candidate reference gene's expressions. It turned out that for all the samples and most subgroups, ACT and GAPDH were the most suitable reference genes for normalization. However, the best reference genes were found not always the same one across diverse samples by different computational programs. Our study was the first reference gene evaluation in N. pseudonarcissus and will promote future studies on gene expression levels of N. pseudonarcissus.

Validation of reference genes for accurate normalization of gene expression with quantitative real-time PCR in Haloxylon ammodendron under different abiotic stresses

Haloxylon ammodendron plays an important role in maintaining the structure and function of the entire ecosystem where it grows. No suitable reference genes have been reported in H. ammodendron plants to date. In this study, a total of 8 reference genes (18S, ACT1, ACT7, UBC18, TUA5, GAPDH, EF-1α and UBQ10) were selected from the available trancriptome database, and the expression stability of these 8 candidate genes was validated under different abiotic stress with three different statistical algorithms (geNorm, NormFinder and BestKeeper). The results produced from different models were in agreement with each other essentially: 18S and TUA5 were the most stable genes under drought stress, 18S, the most stable gene under heat stress and mechanical damage, ACT7 and UBC18, stable under salt stress while TUA5 and GAPDH expressed constantly under mechanical damage, and ACT1 expressed steadily under cold conditions. Expression profiles of several stress response genes, including FT-5, FT-9, DREB2A and DREB2C, were further confirmed with various candidate reference genes. None of the candidate genes showed a constant expression among all tested samples. Hence, it's essential to use more than one reference gene in order to guarantee the accuracy of quantitative real-time PCR. The results of this study will contribute to the accuracy and reliability in transcripts quantification, which is of significance to transcription-based studies and applications in this important shrub H. ammodendron.

Characterization of a new GmFAD3A allele in Brazilian CS303TNKCA soybean cultivar

KEY MESSAGE

We molecularly characterized a new mutation in the GmFAD3A gene associated with low linolenic content in the Brazilian soybean cultivar CS303TNKCA and developed a molecular marker to select this mutation. Soybean is one of the most important crops cultivated worldwide. Soybean oil has 13% palmitic acid, 4% stearic acid, 20% oleic acid, 55% linoleic acid and 8% linolenic acid. Breeding programs are developing varieties with high oleic and low polyunsaturated fatty acids (linoleic and linolenic) to improve the oil oxidative stability and make the varieties more attractive for the soy industry. The main goal of this study was to characterize the low linoleic acid trait in CS303TNKCA cultivar. We sequenced CS303TNKCA GmFAD3A, GmFAD3B and GmFAD3C genes and identified an adenine point deletion in the GmFAD3A exon 5 (delA). This alteration creates a premature stop codon, leading to a truncated protein with just 207 residues that result in a non-functional enzyme. Analysis of enzymatic activity by heterologous expression in yeast support delA as the cause of low linolenic acid content in CS303TNKCA. Thus, we developed a TaqMan genotyping assay to associate delA with low linolenic acid content in segregating populations. Lines homozygous for delA had a linolenic acid content of 3.3 to 4.4%, and the variation at this locus accounted for 50.83 to 73.70% of the phenotypic variation. This molecular marker is a new tool to introgress the low linolenic acid trait into elite soybean cultivars and can be used to combine with high oleic trait markers to produce soybean with enhanced economic value. The advantage of using CS303TNKCA compared to other lines available in the literature is that this cultivar has good agronomic characteristics and is adapted to Brazilian conditions.

Selection and validation of reference genes for quantitative gene expression analyses in black locust (Robinia pseudoacacia L.) using real-time quantitative PCR

Black locust (Robinia pseudoacacia L.) is an easy to raise, fast growing, medium-sized deciduous tree species highly tolerant to harsh eco-conditions, i.e., drought and harsh winters, and it is widely adaptable to sandy, loamy, and marshy soils. The basis for this adaptability remains to be investigated at the transcriptomic level using real-time quantitative PCR (qPCR). Selection of a reliable gene for the normalization of qPCR data is important for obtaining accurate results in gene expression. The goal of this study was to identify an appropriate reference gene from 12 candidate genes for gene expression analysis in black locust exposed to various stressors such as abscisic acid (ABA), NaCl, polyethylene glycol (PEG) and varying temperatures. In GeNorm and NormFinder analyses, ACT (actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase) gene expression were the most stable in all conditions except heat stress, but in BestKeeper analysis, GAPDH and helicase gene expression were the most stable under NaCl and heat stress. In contrast, ACT and GAPDH were highest under abscisic acid (ABA), GAPDH and βTUB (beta tubulin) under cold stress, and helicase and EF1α (elongation factor 1 alpha) under PEG stress. We found that the most stable reference gene combination for all conditions was ACT and GAPDH. Additionally, the expression pattern of NAC2 (a transcription factor) and BGL2 in different tissues and under different stress conditions was analyzed relative to ACT and GAPDH and UBQ (ubiquitin) the least stably expressed gene. NAC2 and BGL2 both had highest expression in flowers and pods under ABA stress at 48h. This study provides useful reference genes for future gene expression studies in black locust.

Plant reference genes for development and stress response studies

Many reference genes are used by different laboratories for gene expression analyses to indicate the relative amount of input RNA/DNA in the experiment. These reference genes are supposed to show least variation among the treatments and with the control sets in a given experiment. However, expression of reference genes varies significantly from one set of experiment to the other. Thus, selection of reference genes depends on the experimental conditions. Sometimes the average expression of two or three reference genes is taken as standard. This review consolidated the details of about 120 genes attempted for normalization during comparative expression analysis in 16 different plants. Plant species included in this review are Arabidopsis thaliana, cotton (Gossypium hirsutum), tobacco (Nicotiana benthamiana and N. tabacum), soybean (Glycine max), rice (Oryza sativa), blueberry (Vaccinium corymbosum), tomato (Solanum lycopersicum), wheat (Triticum aestivum), potato (Solanum tuberosum), sugar cane (Saccharum sp.), carrot (Daucus carota), coffee (Coffea arabica), cucumber (Cucumis sativus), kiwi (Actinidia deliciosa) and grape (Vitis vinifera). The list includes model and cultivated crop plants from both monocot and dicot classes. We have categorized plant-wise the reference genes that have been used for expression analyses in any or all of the four different conditions such as biotic stress, abiotic stress, developmental stages and various organs and tissues, reported till date. This review serves as a guide during the reference gene hunt for gene expression analysis studies.

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.

Comparative transcriptome analysis of nodules of two Mesorhizobium-chickpea associations with differential symbiotic efficiency under phosphate deficiency

Phosphate (Pi) deficiency is known to be a major limitation for symbiotic nitrogen fixation (SNF), and hence legume crop productivity globally. However, very little information is available on the adaptive mechanisms, particularly in the important legume crop chickpea (Cicer arietinum L.), which enable nodules to respond to low-Pi availability. Thus, to elucidate these mechanisms in chickpea nodules at molecular level, we used an RNA sequencing approach to investigate transcriptomes of the nodules in Mesorhizobium mediterraneum SWRI9-(MmSWRI9)-chickpea and M. ciceri CP-31-(McCP-31)-chickpea associations under Pi-sufficient and Pi-deficient conditions, of which the McCP-31-chickpea association has a better SNF capacity than the MmSWRI9-chickpea association during Pi starvation. Our investigation revealed that more genes showed altered expression patterns in MmSWRI9-induced nodules than in McCP-31-induced nodules (540 vs. 225) under Pi deficiency, suggesting that the Pi-starvation-more-sensitive MmSWRI9-induced nodules required expression change in a larger number of genes to cope with low-Pi stress than the Pi-starvation-less-sensitive McCP-31-induced nodules. The functional classification of differentially expressed genes (DEGs) was examined to gain an understanding of how chickpea nodules respond to Pi starvation, caused by soil Pi deficiency. As a result, more DEGs involved in nodulation, detoxification, nutrient/ion transport, transcriptional factors, key metabolic pathways, Pi remobilization and signalling were found in Pi-starved MmSWRI9-induced nodules than in Pi-starved McCP-31-induced nodules. Our findings have enabled the identification of molecular processes that play important roles in the acclimation of nodules to Pi deficiency, ultimately leading to the development of Pi-efficient chickpea symbiotic associations suitable for Pi-deficient soils.

Identification and comprehensive evaluation of reference genes for RT-qPCR analysis of host gene-expression in Brassica juncea-aphid interaction using microarray data

Brassica juncea is a chief oil yielding crop in many parts of the world including India. With advancement of molecular techniques, RT-qPCR based study of gene-expression has become an integral part of experimentations in crop breeding. In RT-qPCR, use of appropriate reference gene(s) is pivotal. The virtue of the reference genes, being constant in expression throughout the experimental treatments, needs to be validated case by case. Appropriate reference gene(s) for normalization of gene-expression data in B. juncea during the biotic stress of aphid infestation is not known. In the present investigation, 11 reference genes identified from microarray database of Arabidopsis-aphid interaction at a cut off FDR ≤0.1, along with two known reference genes of B. juncea, were analyzed for their expression stability upon aphid infestation. These included 6 frequently used and 5 newly identified reference genes. Ranking orders of the reference genes in terms of expression stability were calculated using advanced statistical approaches such as geNorm, NormFinder, delta Ct and BestKeeper. The analysis suggested CAC, TUA and DUF179 as the most suitable reference genes. Further, normalization of the gene-expression data of STP4 and PR1 by the most and the least stable reference gene, respectively has demonstrated importance and applicability of the recommended reference genes in aphid infested samples of B. juncea.

Identification of suitable reference genes for real-time quantitative PCR analysis of hydrogen peroxide-treated human umbilical vein endothelial cells

Background

Oxidative stress can induce cell injury in vascular endothelial cells, which is the initial event in the development of atherosclerosis. Although quantitative real-time polymerase chain reaction (qRT-PCR) has been widely used in gene expression studies in oxidative stress injuries, using carefully validated reference genes has not received sufficient attention in related studies. The objective of this study, therefore, was to select a set of stably expressed reference genes for use in qRT-PCR normalization in oxidative stress injuries in human umbilical vein endothelial cells (HUVECs) induced by hydrogen peroxide (H₂O₂).

Results

Using geNorm analysis, we found that five stably expressed reference genes were sufficient for normalization in qRT-PCR analysis in HUVECs treated with H₂O₂. Genes with the most stable expression according to geNorm were U6, TFRC, RPLP0, GAPDH, and ACTB, and according to NormFinder were ALAS1, TFRC, U6, GAPDH, and ACTB.

Conclusion

Taken together, our study demonstrated that the expression stability of reference genes may differ according to the statistical program used. U6, TFRC, RPLP0, GAPDH, and ACTB was the optimal set of reference genes for studies on gene expression performed by qRT-PCR assays in HUVECs under oxidative stress study.

Selection of suitable reference genes for quantitative real-time PCR gene expression analysis in Salix matsudana under different abiotic stresses

Salix matsudana is a deciduous, rapidly growing willow species commonly cultivated in China, which can tolerate drought, salt, and heavy metal stress conditions. Selection of suitable reference genes for quantitative real-time PCR is important for normalizing the expression of the key genes associated with various stresses. To validate suitable reference genes, we selected 11 candidate reference genes (five traditional housekeeping genes and six novel genes) and analyzed their expression stability in various samples, including different tissues and under different abiotic stress treatments. The expression of these genes was determined using five programs-geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. The results showed that α-TUB2 (alpha-tubulin 2) and DnaJ (chaperone protein DnaJ 49) were the most stable reference genes across all the tested samples. We measured the expression profiles of the defense response gene SmCAT (catalase) using the two most stable and one least stable reference genes in all samples of S. matsudana. The relative quantification of SmCAT varied greatly according to the different reference genes. We propose that α-TUB2 and DnaJ should be the preferred reference genes for normalization and quantification of transcript levels in future gene expression studies in willow species under various abiotic stress conditions.

Evaluation of Reference Genes for Normalization of Gene Expression Using Quantitative RT-PCR under Aluminum, Cadmium, and Heat Stresses in Soybean

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is widely used to analyze the relative gene expression level, however, the accuracy of qRT-PCR is greatly affected by the stability of reference genes, which is tissue- and environment- dependent. Therefore, choosing the most stable reference gene in a specific tissue and environment is critical to interpret gene expression patterns. Aluminum (Al), cadmium (Cd), and heat stresses are three important abiotic factors limiting soybean (Glycine max) production in southern China. To identify the suitable reference genes for normalizing the expression levels of target genes by qRT-PCR in soybean response to Al, Cd and heat stresses, we studied the expression stability of ten commonly used housekeeping genes in soybean roots and leaves under these three abiotic stresses, using five approaches, BestKeeper, Delta Ct, geNorm, NormFinder and RefFinder. We found TUA4 is the most stable reference gene in soybean root tips under Al stress. Under Cd stress, Fbox and UKN2 are the most stable reference genes in roots and leaves, respectively, while 60S is the most suitable reference gene when analyzing both roots and leaves together. For heat stress, TUA4 and UKN2 are the most stable housekeeping genes in roots and leaves, respectively, and UKN2 is the best reference gene for analysis of roots and leaves together. To validate the reference genes, we quantified the relative expression levels of six target genes that were involved in soybean response to Al, Cd or heat stresses, respectively. The expression patterns of these target genes differed between using the most and least stable reference genes, suggesting the selection of a suitable reference gene is critical for gene expression studies.

Selection of Candidate Reference Genes for Gene Expression Analysis in Kentucky Bluegrass (Poa pratensis L.) under Abiotic Stress

Kentucky bluegrass (Poa pratensis L.) belong to Gramineae and is widely used in lawns, golf courses, landscapes, and sport fields as a prominent cool-season grass. Gene expression patterns during different stages of plant development can provide clues toward the understanding of its biological functions. The selection and validation of reference genes are the first steps in any real-time quantitative PCR gene expression study. Therefore, suitable reference genes are necessary for obtaining reliable results in real-time quantitative PCR analyses of Kentucky bluegrass. In the present study, 9 candidate reference genes were chosen, and their expression stability in the leaves and roots of Kentucky bluegrass under different stresses (drought, salt, heat, and cold) were evaluated using the GeNorm, NormFinder, BestKeeper, and RefFinder programs. The results showed that the expression stability of the candidate reference genes was dependent on the experimental conditions. The combination of SAM with GAPDH was the most stable in leaves under salt stress and cold stress, while TUB combined with ACT or GAPDH was stable in roots under salt or cold stress, respectively. ACT and SAM maintained stable expression in drought-treated leaves, and GAPDH combined with ACT was stable in drought-treated roots. SAM and TUB exhibited stable expression in heat-treated leaves. ACT and RPL were stable in heat-treated roots. In addition, the expression patterns of PpFEH in response to drought and cold stress were used to confirm the reliability of the selected reference genes, indicating that the use of an inappropriate reference gene as the internal control will cause erroneous results. This work is the first study on the expression stability of reference genes in Kentucky bluegrass and will be particularly useful in the selection of stress-tolerance genes and the identification of the molecular mechanisms conferring stress tolerance in this species.

Identification and evaluation of reference genes for expression studies by RT-qPCR during embryonic development of the emerging model organism, Macrobrachium olfersii

RT-qPCR is a sensitive and highly efficient technique that is widely used in gene expression analysis and to provide insight into the molecular mechanisms underlying embryonic development. The freshwater prawn, Macrobrachium olfersii is an emerging model organism, but, the stable reference genes of this species need to be identified and validated for RT-qPCR analysis. Thus, the aim of this study was to evaluate the expression stability of six genes (β-act, GAPDH, EF-1α, RpL8, RpS6, AK) in embryos and in adult tissues (cerebral ganglia, muscle and hepatopancreas) of M. olfersii. The expression stabilities of these genes were evaluated using geNorm, NormFinder, BestKeeper, ΔCt method and integrated tool RefFinder. In the general ranking, RpL8 and RpS6 were the most stable genes in embryos, while RpS6 and RpL8 were the most stable in a combined adult tissue analysis. Analysis of the adult tissues revealed that β-act and AK were the most stable genes in cerebral ganglia, RpL8 and AK in muscle, and RpS6 and β-act in hepatopancreas. EF-1α and GAPDH were the least stable genes and as normalizer genes in RT-qPCR affected expression of the Distal-less gene during M. olfersii development. This study provides suitable reference genes for RT-qPCR analysis and allows future studies of the gene expression in M. olfersii for understanding the molecular mechanisms of their development. To our knowledge, this is the first published study that identifies and evaluates reference genes for RT-qPCR analysis in M. olfersii and could be useful as basis for evaluations of reference genes in other prawns.

Identification of reference genes for quantitative RT-PCR analysis of microRNAs and mRNAs in castor bean (Ricinus communis L.) under drought stress

Quantitative real-time PCR (RT-qPCR) is one of the most powerful and sensitive techniques to the study of gene expression. Several factors influence RT-qPCR performance though, including the stability of the reference genes used for data normalization. While the selection of appropriate reference genes is crucial for accurate and reliable gene expression analysis, no suitable reference genes have been previously identified in castor bean under drought stress. In this study, the expression stability of eleven mRNAs, thirteen microRNAs (miRNAs) and one small nuclear RNA were analyzed in roots and leaves across different levels of water deficit. Three different algorithms were employed to analyze the RT-qPCR data, and the resulting outputs were merged using a non-weighted unsupervised rank aggregation method. Our analysis indicated that the Elongation factor 1-beta (EF1B), Protein phosphatase 2A (PP2A) and ADP-ribosylation factor (ADP) ranked as the best candidates across diverse samples submitted to different levels of drought conditions. EF1B and Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and EF1B and SKP1/ASK-interacting protein 16 (SKIP16) were found as the most suitable reference genes for expression analysis in roots and leaves, respectively. In addition, miRNAs miR168, miR160 and miR397 were selected as optimal reference genes across all tissues and treatments. miR168 and miR156 were recommended as reference for roots, while miR168 and miR160 were recommended for leaves. Together, our results constitute the first attempt to identify and validate the most suitable reference genes for accurate normalization of gene expression in castor bean under drought stress.

Selection of Suitable Reference Genes for RT-qPCR Normalization under Abiotic Stresses and Hormone Stimulation in Persimmon (Diospyros kaki Thunb)

The success of quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) to quantify gene expression depends on the stability of the reference genes used for data normalization. To date, systematic screening for reference genes in persimmon (Diospyros kaki Thunb) has never been reported. In this study, 13 candidate reference genes were cloned from 'Nantongxiaofangshi' using information available in the transcriptome database. Their expression stability was assessed by geNorm and NormFinder algorithms under abiotic stress and hormone stimulation. Our results showed that the most suitable reference genes across all samples were UBC and GAPDH, and not the commonly used persimmon reference gene ACT. In addition, UBC combined with RPII or TUA were found to be appropriate for the "abiotic stress" group and α-TUB combined with PP2A were found to be appropriate for the "hormone stimuli" group. For further validation, the transcript level of the DkDREB2C homologue under heat stress was studied with the selected genes (CYP, GAPDH, TUA, UBC, α-TUB, and EF1-α). The results suggested that it is necessary to choose appropriate reference genes according to the test materials or experimental conditions. Our study will be useful for future studies on gene expression in persimmon.