Validation of reference genes for real-time quantitative PCR normalization in soybean developmental and germinating seeds

Differential metabolic reprogramming in developing soybean embryos in response to nutritional conditions and abscisic acid

Seed storage compound deposition is influenced by both maternal and filial tissues. Within this framework, we analyzed strategies that operate during the development and filling of soybean embryos, using in vitro culture systems combined with metabolomics and proteomics approaches. The carbon:nitrogen ratio (C:N) of the maternal supply and the hormone abscisic acid (ABA) are specific and interacting signals inducing differential metabolic reprogrammings linked to changes in the accumulation of storage macromolecules like proteins or oils. Differences in the abundance of sugars, amino acids, enzymes, transporters, transcription factors, and proteins involved in signaling were detected. Embryos adapted to the nutritional status by enhancing the metabolism of both carbon and nitrogen under lower C:N ratio condition or only carbon under higher C:N ratio condition. ABA turned off multiple pathways especially in high availability of amino acids, prioritizing the storage compounds biosynthesis. Common responses induced by ABA involved increased sucrose uptake (to increase the sink force) and oleosin (oil body structural component) accumulation. In turn, ABA differentially promoted protein degradation under lower nitrogen supply in order to sustain the metabolic demands. Further, the operation of a citrate shuttle was suggested by transcript quantification and enzymatic activity measurements. The results obtained are useful to help define biotechnological tools and technological approaches to improve oil and protein yields, with direct impact on human and animal nutrition as well as in green chemistry.

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.

Evaluation of reference genes and characterization of the MYBs in xylem radial change of Chinese fir stem

The radial change (RC) of tree stem is the process of heartwood formation involved in complex molecular mechanism. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), an evergreen species, is an important fast-growing timber tree in southern China. In this study, the top four stable genes (IDH, UBC2, RCA and H2B) were selected in RC tissues of 15 years old Chinese fir stem (RC15) and the genes (H2B, 18S, TIP41 and GAPDH) were selected in RC tissues of 30 years old Chinese fir stem (RC30). The stability of the reference genes is higher in RC30 than in RC15. Sixty-one MYB transcripts were obtained on the PacBio Sequel platform from woody tissues of one 30 years old Chinese fir stem. Based on the number of MYB DNA-binding domain and phylogenetic relationships, the ClMYB transcripts contained 21 transcripts of MYB-related proteins (1R-MYB), 39 transcripts of R2R3-MYB proteins (2R-MYB), one transcript of R1R2R3-MYB protein (3R-MYB) belonged to 18 function-annotated clades and two function-unknown clades. In RC woody tissues of 30 years old Chinese fir stem, ClMYB22 was the transcript with the greatest fold change detected by both RNA-seq and qRT-PCR. Reference genes selected in this study will be helpful for further verification of transcript abundance patterns during the heartwood formation of Chinese fir.

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.

Characterization of Two Growth Period QTLs Reveals Modification of PRR3 Genes During Soybean Domestication

Soybean yield is largely dependent on growth period. We characterized two growth period quantitative trait loci, Gp11 and Gp12, from a recombinant inbred population generated from a cross of wild (W05) and cultivated (C08) soybean. Lines carrying Gp11C08 and Gp12C08 tend to have a shorter growth period and higher expression of GmFT2a and GmFT5a. Furthermore, multiple interval mapping suggests that Gp11 and Gp12 may be genetically interacting with the E2 locus. This is consistent with the observation that GmFT2a and GmFT5a are activated by Gp11C08 and Gp12C08 at ZT4 in the recessive e2 but not the dominant E2 background. Gp11 and Gp12 are duplicated genomic regions each containing a copy of the soybean ortholog of PSEUDO RESPONSE REGULATOR 3 (GmPRR3A and GmPRR3B). GmPRR3A and GmPRR3B from C08 carry mutations that delete the CCT domain in the encoded proteins. These mutations were selected during soybean improvement and they alter the subcellular localization of GmPRR3A and GmPRR3B. Furthermore, GmPRR3A and GmPRR3B can interact with TOPLESS-related transcription factors, suggesting that they function in a transcription repressor complex. This study addresses previously unexplored components of the genetic network that probably controls the growth period of soybean and puts these loci into context with the well-characterized growth period-regulating E loci.

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.

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.

Reference gene selection for molecular studies of dormancy in wild oat (Avena fatua L.) caryopses by RT-qPCR method

Molecular studies of primary and secondary dormancy in Avena fatua L., a serious weed of cereal and other crops, are intended to reveal the species-specific details of underlying molecular mechanisms which in turn may be useable in weed management. Among others, quantitative real-time PCR (RT-qPCR) data of comparative gene expression analysis may give some insight into the involvement of particular wild oat genes in dormancy release, maintenance or induction by unfavorable conditions. To assure obtaining biologically significant results using this method, the expression stability of selected candidate reference genes in different data subsets was evaluated using four statistical algorithms i.e. geNorm, NormFinder, Best Keeper and ΔCt method. Although some discrepancies in their ranking outputs were noticed, evidently two ubiquitin-conjugating enzyme homologs, AfUBC1 and AfUBC2, as well as one homolog of glyceraldehyde 3-phosphate dehydrogenase AfGAPDH1 and TATA-binding protein AfTBP2 appeared as more stably expressed than AfEF1a (translation elongation factor 1α), AfGAPDH2 or the least stable α-tubulin homolog AfTUA1 in caryopses and seedlings of A. fatua. Gene expression analysis of a dormancy-related wild oat transcription factor VIVIPAROUS1 (AfVP1) allowed for a validation of candidate reference genes performance. Based on the obtained results it can be recommended that the normalization factor calculated as a geometric mean of Cq values of AfUBC1, AfUBC2 and AfGAPDH1 would be optimal for RT-qPCR results normalization in the experiments comprising A. fatua caryopses of different dormancy status.

Differential expression of a WRKY gene between wild and cultivated soybeans correlates to seed size

Soybean (Glycine max) probably originated from the wild soybean (Glycine soja). Glycine max has a significantly larger seed size, but the underlying genomic changes are largely unknown. Candidate regulatory genes were preliminarily proposed by data co-localizing RNA sequencing with the quantitative loci (QTLs) for seed size. The soybean gene locus SoyWRKY15a and its orthologous genes from G. max (GmWRKY15a) and G. soja (GsWRKY15a) were analyzed in detail. The coding sequences were nearly identical between the two orthologs, but GmWRKY15a was significantly more highly expressed than GsWRKY15a. Four haplotypes (H1-H4) were found and they varied in the size of a CT-core microsatellite locus in the 5'-untranslated region of this gene. H1 (with six CT-repeats) was the only allelic version found in G. max, while H3 (with five CT-repeats) was the dominant G. soja allele. Differential expression of this gene in soybean pods was correlated with CT-repeat variation, and manipulation of the CT copy number altered the reporter gene expression, suggesting a regulatory role for the simple sequence repeats. Seed weight of wild soybeans harboring H1 was significantly greater than that of soybeans having haplotypes H2, H3, or H4, and seed weight was correlated with gene expression, suggesting the influence of GsWRKY15a in controlling seed size. However, the seed size might be refractory to increased SoyWRKY15a expression in cultivated soybeans. The evolutionary significance of SoyWRKY15a variation in soybean seed domestication is discussed.

A PP2C-1 Allele Underlying a Quantitative Trait Locus Enhances Soybean 100-Seed Weight

Cultivated soybeans may lose some useful genetic loci during domestication. Introgression of genes from wild soybeans could broaden the genetic background and improve soybean agronomic traits. In this study, through whole-genome sequencing of a recombinant inbred line population derived from a cross between a wild soybean ZYD7 and a cultivated soybean HN44, and mapping of quantitative trait loci for seed weight, we discovered that a phosphatase 2C-1 (PP2C-1) allele from wild soybean ZYD7 contributes to the increase in seed weight/size. PP2C-1 may achieve this function by enhancing cell size of integument and activating a subset of seed trait-related genes. We found that PP2C-1 is associated with GmBZR1, a soybean ortholog of Arabidopsis BZR1, one of key transcription factors in brassinosteroid (BR) signaling, and facilitate accumulation of dephosphorylated GmBZR1. In contrast, the PP2C-2 allele with variations of a few amino acids at the N-terminus did not exhibit this function. Moreover, we showed that GmBZR1 could promote seed weight/size in transgenic plants. Through analysis of cultivated soybean accessions, we found that 40% of the examined accessions do not have the PP2C-1 allele, suggesting that these accessions can be improved by introduction of this allele. Taken together, our study identifies an elite allele PP2C-1, which can enhance seed weight and/or size in soybean, and pinpoints that manipulation of this allele by molecular-assisted breeding may increase production in soybean and other legumes/crops.

Differential Contribution of Malic Enzymes during Soybean and Castor Seeds Maturation

Malic enzymes (ME) catalyze the decarboxylation of malate generating pyruvate, CO2 and NADH or NADPH. In some organisms it has been established that ME is involved in lipids biosynthesis supplying carbon skeletons and reducing power. In this work we studied the MEs of soybean and castor, metabolically different oilseeds. The comparison of enzymatic activities, transcript profiles and organic acid contents suggest different metabolic strategies operating in soybean embryo and castor endosperm in order to generate precursors for lipid biosynthesis. In castor, the malate accumulation pattern agrees with a central role of this metabolite in the provision of carbon to plastids, where the biosynthesis of fatty acids occurs. In this regard, the genome of castor possesses a single gene encoding a putative plastidic NADP-ME, whose expression level is high when lipid deposition is active. On the other hand, NAD-ME showed an important contribution to the maturation of soybean embryos, perhaps driving the carbon relocation from mitochondria to plastids to support the fatty acids synthesis in the last stages of seed filling. These findings provide new insights into intermediary metabolism in oilseeds and provide new biotechnological targets to improve oil yields.

The transcriptomic signature of developing soybean seeds reveals the genetic basis of seed trait adaptation during domestication

Cultivated soybean has undergone many transformations during domestication. In this paper we report a comprehensive assessment of the evolution of gene co-expression networks based on the analysis of 40 transcriptomes from developing soybean seeds in cultivated and wild soybean accessions. We identified 2680 genes that are differentially expressed during seed maturation and established two cultivar-specific gene co-expression networks. Through analysis of the two networks and integration with quantitative trait locus data we identified two potential key drivers for seed trait formation, GA20OX and NFYA. GA20OX encodes an enzyme in a rate-limiting step of gibberellin biosynthesis, and NFYA encodes a transcription factor. Overexpression of GA20OX and NFYA enhanced seed size/weight and oil content, respectively, in seeds of transgenic plants. The two genes showed significantly higher expression in cultivated than in wild soybean, and the increases in expression were associated with genetic variations in the promoter region of each gene. Moreover, the expression of GA20OX and NFYA in seeds of soybean accessions correlated with seed weight and oil content, respectively. Our study reveals transcriptional adaptation during soybean domestication and may identify a mechanism of selection by expression for seed trait formation, providing strategies for future breeding practice.

Normalization for Relative Quantification of mRNA and microRNA in Soybean Exposed to Various Abiotic Stresses

Plant microRNAs are small non-coding, endogenic RNA molecule (containing 20-24 nucleotides) produced from miRNA precursors (pri-miRNA and pre-miRNA). Evidence suggests that up and down regulation of the miRNA targets the mRNA genes involved in resistance against biotic and abiotic stresses. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a powerful technique to analyze variations in mRNA levels. Normalizing the data using reference genes is essential for the analysis of reliable RT-qPCR data. In this study, two groups of candidate reference mRNAs and miRNAs in soybean leaves and roots treated with various abiotic stresses (PEG-simulated drought, salinity, alkalinity, salinity+alkalinity, and abscisic acid) were analyzed by RT-qPCR. We analyzed the most appropriate reference mRNA/miRNAs using the geNorm, NormFinder, and BestKeeper algorithms. According to the results, Act and EF1b were the most suitable reference mRNAs in leaf and root samples, for mRNA and miRNA precursor data normalization. The most suitable reference miRNAs found in leaf and root samples were 166a and 167a for mature miRNA data normalization. Hence the best combinations of reference mRNAs for mRNA and miRNA precursor data normalization were EF1a + Act or EF1b + Act in leaf samples, and EF1a + EF1b or 60s + EF1b in root samples. For mature miRNA data normalization, the most suitable combinations of reference miRNAs were 166a + 167d in leaf samples, and 171a + 156a or 167a + 171a in root samples. We identified potential reference mRNA/miRNAs for accurate RT-qPCR data normalization for mature miRNA, miRNA precursors, and their targeted mRNAs. Our results promote miRNA-based studies on soybean plants exposed to abiotic stress conditions.

Early Selection for Smut Resistance in Sugarcane Using Pathogen Proliferation and Changes in Physiological and Biochemical Indices

Sugarcane smut disease, caused by Sporisorium scitamineum, significantly decreases yield and use of resistant cultivars is the most cost-effective measure for disease control. Current field testing methods for identification of smut resistance are time-consuming and hindered by environmental variability. Our goal was to develop an efficient and reliable resistance identification technique that is rapid, performed in a controlled environment, and stable. Nine sugarcane cultivars with different phenotypic resistance levels were selected. TaqMan quantitative real-time polymerase chain reaction analysis was performed to measure copy number changes of smut pathogen in sugarcane buds at 0-7 days after needle puncture inoculation. There was a positive correlation between time after inoculation and the amount of smut pathogen in the sugarcane bud. This reached a peak value on 7 days, and the copy number of S. scitamineum increased in the following order: YZ03-258, FN40, YZ01-1413, GT02-467, ROC22, YT96-86, YZ03-103, FN39, LC05-136. After smut pathogen inoculation, differences in the physiological and biochemical indices of the nine cultivars were observed. Peroxidase, ascorbate peroxidase, catalase, superoxide dismutase, β-1,3-glucanase, and malondialdehyde were grouped into three main components, and the cumulative contribution rate was 80.177%, revealing that these are useful physiological and biochemical indicators of smut resistance. Subordinate function analysis indicated that the levels of smut resistance of the nine genotypes were (high to low): YZ03-258, FN40, YZ01-1413, GT02-467, ROC22, YZ03-103, YT96-86, FN39, LC05-136, which is similar to the results from copy number determination of smut pathogens. The results suggest that after artificial needle inoculation, rapid identification of physiological resistance to sugarcane smut was achieved based on copy number increases in the sugarcane smut pathogen and the physiological and biochemical changes in sugarcane bud during the early phase of infection.

Evaluation of Sorghum [Sorghum bicolor (L.)] Reference Genes in Various Tissues and under Abiotic Stress Conditions for Quantitative Real-Time PCR Data Normalization

Accurate and reliable gene expression data from qPCR depends on stable reference gene expression for potential gene functional analyses. In this study, 15 reference genes were selected and analyzed in various sample sets including abiotic stress treatments (salt, cold, water stress, heat, and abscisic acid) and tissues (leaves, roots, seedlings, panicle, and mature seeds). Statistical tools, including geNorm, NormFinder and RefFinder, were utilized to assess the suitability of reference genes based on their stability rankings for various sample groups. For abiotic stress, PP2A and CYP were identified as the most stable genes. In contrast, EIF4α was the most stable in the tissue sample set, followed by PP2A; PP2A was the most stable in all the sample set, followed by EIF4α. GAPDH, and UBC1 were the least stably expressed in the tissue and all the sample sets. These results also indicated that the use of two candidate reference genes would be sufficient for the optimization of normalization studies. To further verify the suitability of these genes for use as reference genes, SbHSF5 and SbHSF13 gene expression levels were normalized using the most and least stable sorghum reference genes in root and water stressed-leaf tissues of five sorghum varieties. This is the first systematic study of the selection of the most stable reference genes for qPCR-related assays in Sorghum bicolor that will potentially benefit future gene expression studies in sorghum and other closely related species.

Selection of reference genes for quantitative real-time PCR normalization in creeping bentgrass involved in four abiotic stresses

This study identified stable reference genes for normalization of gene expression data in qRT-PCR analysis of leaf and root tissues in creeping bentgrass under four abiotic stresses. Examination of gene expression using quantitative real-time PCR (qRT-PCR) in plant responses to abiotic stresses can provide valuable information for stress-tolerance improvement. Selecting stable reference genes for qRT-PCR analysis is critically important. The objective of this study was to determine the stability of expression for eight candidate reference genes (ACT, EF1a, TUB, UPL7, GAPDH, PP2A, PEPKR1, and CACS) in two tissues (roots and leaves) of a perennial grass species under four abiotic stresses (salt, drought, cold, and heat) using four programs (GeNorm, NormFinder, BestKeeper, and RefFinder). The results showed that (1) the combinations of CACS and UPL7 or PP2A and ACT were stably expressed in salt-treated roots or leaves; (2) the combinations of GAPDH and CACS or PP2A and PEPKR1 were stable in roots and leaves under drought stress; (3) CACS and PP2A exhibited stable expression in cold-treated roots and the combination of EF1a and UPL7 was also stable in cold-treated leaves; and (4) CACS and PP2A were the two most stable reference genes in heat-stressed roots and UPL7 combined with GAPDH and PP2A was stably expressed in heat-stressed leaves. The qRT-PCR analysis of a target gene, AsSAP expression patterns in response to salinity and drought stress, confirmed the reliability of those selected and stable reference genes. Identification of stable reference genes in creeping bentgrass will improve assay accuracy for selecting stress-tolerance genes and identifying molecular mechanisms conferring stress tolerance in this species.

Using RNA-Seq Data to Evaluate Reference Genes Suitable for Gene Expression Studies in Soybean

Differential gene expression profiles often provide important clues for gene functions. While reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is an important tool, the validity of the results depends heavily on the choice of proper reference genes. In this study, we employed new and published RNA-sequencing (RNA-Seq) datasets (26 sequencing libraries in total) to evaluate reference genes reported in previous soybean studies. In silico PCR showed that 13 out of 37 previously reported primer sets have multiple targets, and 4 of them have amplicons with different sizes. Using a probabilistic approach, we identified new and improved candidate reference genes. We further performed 2 validation tests (with 26 RNA samples) on 8 commonly used reference genes and 7 newly identified candidates, using RT-qPCR. In general, the new candidate reference genes exhibited more stable expression levels under the tested experimental conditions. The three newly identified candidate reference genes Bic-C2, F-box protein2, and VPS-like gave the best overall performance, together with the commonly used ELF1b. It is expected that the proposed probabilistic model could serve as an important tool to identify stable reference genes when more soybean RNA-Seq data from different growth stages and treatments are used.

Reference Gene Selection for qPCR Analysis in Tomato-Bipartite Begomovirus Interaction and Validation in Additional Tomato-Virus Pathosystems

Quantitative Polymerase Chain Reaction (qPCR) is currently the most sensitive technique used for absolute and relative quantification of a target gene transcript, requiring the use of appropriated reference genes for data normalization. To accurately estimate the relative expression of target tomato (Solanum lycopersicum L.) genes responsive to several virus species in reverse transcription qPCR analysis, the identification of reliable reference genes is mandatory. In the present study, ten reference genes were analyzed across a set of eight samples: two tomato contrasting genotypes ('Santa Clara', susceptible, and its near-isogenic line 'LAM 157', resistant); subjected to two treatments (inoculation with Tomato chlorotic mottle virus (ToCMoV) and its mock-inoculated control) and in two distinct times after inoculation (early and late). Reference genes stability was estimated by three statistical programs (geNorm, NormFinder and BestKeeper). To validate the results over broader experimental conditions, a set of ten samples, corresponding to additional three tomato-virus pathosystems that included tospovirus, crinivirus and tymovirus + tobamovirus, was analyzed together with the tomato-ToCMoV pathosystem dataset, using the same algorithms. Taking into account the combined analyses of the ranking order outputs from the three algorithms, TIP41 and EF1 were identified as the most stable genes for tomato-ToCMoV pathosystem, and TIP41 and EXP for the four pathosystems together, and selected to be used as reference in the forthcoming expression qPCR analysis of target genes in experimental conditions involving the aforementioned tomato-virus pathosystems.

Recommended Reference Genes for Quantitative PCR Analysis in Soybean Have Variable Stabilities during Diverse Biotic Stresses

For real-time reverse transcription-PCR (qRT-PCR) in soybean, reference genes in different tissues, developmental stages, various cultivars, and under stress conditions have been suggested but their usefulness for research on soybean under various biotic stresses occurring in North-Central U.S. is not known. Here, we investigated the expression stabilities of ten previously recommended reference genes (ABCT, CYP, EF1A, FBOX, GPDH, RPL30, TUA4, TUB4, TUA5, and UNK2) in soybean under biotic stress from Bean pod mottle virus (BPMV), powdery mildew (PMD), soybean aphid (SBA), and two-spotted spider mite (TSSM). BPMV, PMD, SBA, and TSSM are amongst the most common pest problems on soybean in North-Central U.S. and other regions. Reference gene stability was determined using three software algorithms (geNorm, NormFinder, BestKeeper) and a web-based tool (RefFinder). Reference genes showed variability in their expression as well as stability across various stressors and the best reference genes were stress-dependent. ABCT and FBOX were found to be the most stable in soybean under both BPMV and SBA stress but these genes had only minimal to moderate stability during PMD and TSSM stress. Expression of TUA4 and CYP was found to be most stable during PMD stress; TUB4 and TUA4 were stable under TSSM stress. Under various biotic stresses on soybean analyzed, GPDH expression was found to be consistently unstable. For all biotic stressors on soybean, we obtained pairwise variation (V2/3) values less than 0.15 which suggested that combined use of the two most stable reference genes would be sufficient for normalization. Further, we demonstrated the utility of normalizing the qRT-PCR data for target genes using the most stable reference genes validated in current study. Following of the recommendations from our current study will enable an accurate and reliable normalization of qRT-PCR data in soybean under biotic stress.

Comprehensive selection of reference genes for quantitative gene expression analysis during seed development in Brassica napus

MicroRNAs have higher expression stability than protein-coding genes in B. napus seeds and are therefore good reference genes for miRNA and mRNA RT-qPCR analysis. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) has become the "gold standard" to gain insight into function of genes. However, the accuracy of the technique depends on appropriate reference genes for quantification analysis in different experimental conditions. Accumulation of microRNAs (miRNAs) has also been studied by RT-qPCR, but there are no reference genes currently validated for normalization of Brassica napus miRNA expression data. In this study, we selected 43 B. napus miRNAs and 18 previously validated mRNA reference genes. The expression stability of the candidate reference genes was evaluated in different tissue samples (stages of seed development, flowers, and leaves) using geNorm, NormFinder, and RefFinder analysis. The best-ranked reference genes for expression studies during seed development (miR167-1_2, miR11-1, miR159-1 and miR168-1) were used to asses the expression of miR03-1. Since candidate miRNAs showed higher expression stability than protein-coding genes in most of the tested conditions, the expression profile of DGAT1 gene was compared when normalized by the four most stable miRNAs reference genes and by the four most stable mRNA reference genes. The expected expression pattern of DGAT1 during seed development was achieved with the use of miRNA as reference genes. In conclusion, the most stable miRNA reference genes can be employed in the normalization of RT-qPCR quantification of miRNAs and protein-coding genes. This work is the first to perform a comprehensive survey of the stability of miRNA reference genes in B. napus and provides guidelines to obtain more accurate RT-qPCR results in B. napus seeds studies.

Reference genes selection and normalization of oxidative stress responsive genes upon different temperature stress conditions in Hypericum perforatum L

Reverse transcription-quantitative real-time PCR (RT-qPCR) is a widely used technique for gene expression analysis. The reliability of this method depends largely on the suitable selection of stable reference genes for accurate data normalization. Hypericum perforatum L. (St. John's wort) is a field growing plant that is frequently exposed to a variety of adverse environmental stresses that can negatively affect its productivity. This widely known medicinal plant with broad pharmacological properties (anti-depressant, anti-tumor, anti-inflammatory, antiviral, antioxidant, anti-cancer, and antibacterial) has been overlooked with respect to the identification of reference genes suitable for RT-qPCR data normalization. In this study, 11 candidate reference genes were analyzed in H. perforatum plants subjected to cold and heat stresses. The expression stability of these genes was assessed using GeNorm, NormFinder and BestKeeper algorithms. The results revealed that the ranking of stability among the three algorithms showed only minor differences within each treatment. The best-ranked reference genes differed between cold- and heat-treated samples; nevertheless, TUB was the most stable gene in both experimental conditions. GSA and GAPDH were found to be reliable reference genes in cold-treated samples, while GAPDH showed low expression stability in heat-treated samples. 26SrRNA and H2A had the highest stabilities in the heat assay, whereas H2A was less stable in the cold assay. Finally, AOX1, AOX2, CAT1 and CHS genes, associated with plant stress responses and oxidative stress, were used as target genes to validate the reliability of identified reference genes. These target genes showed differential expression profiles over time in treated samples. This study not only is the first systematic analysis for the selection of suitable reference genes for RT-qPCR studies in H. perforatum subjected to temperature stress conditions, but may also provide valuable information about the roles of genes associated with temperature stress responses.

Identification of a soybean MOTHER OF FT AND TFL1 homolog involved in regulation of seed germination

Seed germination is an important event in the life cycle of seed plants, and is controlled by complex and coordinated genetic networks. Many genes involved in the regulation of this process have been identified in different plant species so far. Recent studies in both Arabidopsis and wheat have uncovered a new role of MOTHER OF FT AND TFL1 (MFT) in seed germination. Here, we reported a homolog of MFT in soybean (GmMFT) which strongly expressed in seeds. Detailed expression analysis showed that the mRNA level of GmMFT increased with seed development but declined during seed germination. The transcription of GmMFT also responded to exogenous application of ABA and GA3. Ectopic expression of GmMFT CDS in Arabidopsis moderately inhibited seed germination. All these evidences suggest that GmMFT may be a negative regulator of seed germination.

Evaluation of reference genes for quantitative real-time PCR in oil palm elite planting materials propagated by tissue culture

BACKGROUND

The somatic embryogenesis tissue culture process has been utilized to propagate high yielding oil palm. Due to the low callogenesis and embryogenesis rates, molecular studies were initiated to identify genes regulating the process, and their expression levels are usually quantified using reverse transcription quantitative real-time PCR (RT-qPCR). With the recent release of oil palm genome sequences, it is crucial to establish a proper strategy for gene analysis using RT-qPCR. Selection of the most suitable reference genes should be performed for accurate quantification of gene expression levels.

RESULTS

In this study, eight candidate reference genes selected from cDNA microarray study and literature review were evaluated comprehensively across 26 tissue culture samples using RT-qPCR. These samples were collected from two tissue culture lines and media treatments, which consisted of leaf explants cultures, callus and embryoids from consecutive developmental stages. Three statistical algorithms (geNorm, NormFinder and BestKeeper) confirmed that the expression stability of novel reference genes (pOP-EA01332, PD00380 and PD00569) outperformed classical housekeeping genes (GAPDH, NAD5, TUBULIN, UBIQUITIN and ACTIN). PD00380 and PD00569 were identified as the most stably expressed genes in total samples, MA2 and MA8 tissue culture lines. Their applicability to validate the expression profiles of a putative ethylene-responsive transcription factor 3-like gene demonstrated the importance of using the geometric mean of two genes for normalization.

CONCLUSIONS

Systematic selection of the most stably expressed reference genes for RT-qPCR was established in oil palm tissue culture samples. PD00380 and PD00569 were selected for accurate and reliable normalization of gene expression data from RT-qPCR. These data will be valuable to the research associated with the tissue culture process. Also, the method described here will facilitate the selection of appropriate reference genes in other oil palm tissues and in the expression profiling of genes relating to yield, biotic and abiotic stresses.

VIGS technology: an attractive tool for functional genomics studies in legumes

Legume species are among the most important crops worldwide. In recent years, six legume genomes have been completely sequenced, and there is now an urgent need for reverse-genetics tools to validate genes affecting yield and product quality. As most legumes are recalcitrant to stable genetic transformation, virus-induced gene silencing (VIGS) appears to be a powerful alternative technology for determining the function of unknown genes. VIGS technology is based on the property of plant viruses to trigger a defence mechanism related to post-transcriptional gene silencing (PTGS). Infection by a recombinant virus carrying a fragment of a plant target gene will induce homology-dependent silencing of the endogenous target gene. Several VIGS systems have been developed for legume species since 2004, including those based on Bean pod mottle virus, Pea early browning virus, and Apple latent spherical virus, and used in reverse-genetics studies of a wide variety of plant biological processes. In this work, we give an overview of the VIGS systems available for legumes, and present their successful applications in functional genomics studies. We also discuss the limitations of these VIGS systems and the future challenges to be faced in order to use VIGS to its full potential in legume species.

Evaluation and validation of reference genes for normalization of quantitative real-time PCR based gene expression studies in peanut

The quantitative real-time PCR (qPCR) based techniques have become essential for gene expression studies and high-throughput molecular characterization of transgenic events. Normalizing to reference gene in relative quantification make results from qPCR more reliable when compared to absolute quantification, but requires robust reference genes. Since, ideal reference gene should be species specific, no single internal control gene is universal for use as a reference gene across various plant developmental stages and diverse growth conditions. Here, we present validation studies of multiple stably expressed reference genes in cultivated peanut with minimal variations in temporal and spatial expression when subjected to various biotic and abiotic stresses. Stability in the expression of eight candidate reference genes including ADH3, ACT11, ATPsyn, CYP2, ELF1B, G6PD, LEC and UBC1 was compared in diverse peanut plant samples. The samples were categorized into distinct experimental sets to check the suitability of candidate genes for accurate and reliable normalization of gene expression using qPCR. Stability in expression of the references genes in eight sets of samples was determined by geNorm and NormFinder methods. While three candidate reference genes including ADH3, G6PD and ELF1B were identified to be stably expressed across experiments, LEC was observed to be the least stable, and hence must be avoided for gene expression studies in peanut. Inclusion of the former two genes gave sufficiently reliable results; nonetheless, the addition of the third reference gene ELF1B may be potentially better in a diverse set of tissue samples of peanut.

Defining reference genes for qPCR normalization to study biotic and abiotic stress responses in Vigna mungo

Expression of ACT, EF1A; H2A, EF1A, ACT and 18S, TUB showed stability under MYMIV, salinity and drought stress, respectively; these are recommended as reference genes for qPCR normalization in Vigna mungo. Accurate gene expression profiling through qPCR depends on selection of appropriate reference gene(s) for normalization. Due to lack of unanimous internal standard, suitable constitutively expressed reference genes are selected that exhibit stable expression under diverse experimental conditions. In this communication, a comparative evaluation of stability among seven V. mungo genes encoding actin (ACT), histone H2A (H2A), elongation factor 1-alpha (EF1A), 18S rRNA (18S), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), cyclophilin (CYP) and tubulin (TUB) under biotic (MYMIV) and abiotic (drought and salinity) stress conditions has been attempted. Specificity and amplification efficiency for each primer pair were verified; however, cumulative assessment of their accumulated transcripts revealed no uniformity. Therefore, individual stability and suitability of these seven candidates have been assessed in silico, by two widely used algorithms, geNorm and Normfinder. Based on the computed results, high stability was obtained for ACT and EF1A during MYMIV stress, while H2A, EFIA and ACT were found to be most suitable in salinity stress experiments and TUB and 18S during drought treatments. Combinations of ACT/TUB or ACT/EFIA were recommended for their use in the pooled analysis, while expression of 18S and CYP showed greater variations and therefore considered unsuitable as reference genes. Additionally, precise quantification of the target gene VmPRX under these stresses was shown to be a function of reference genes' stability, which tends to get affected when normalized with the least stable genes. Hence, use of these normalizers will facilitate accurate and reliable analyses of gene expression in V. mungo.