Validation of reference genes for gene expression analysis in chicory (Cichorium intybus) using quantitative real-time PCR

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

Expression assay of the COLQ in a family with congenital myasthenic syndrome and symptomatic carriers

Congenital myasthenic syndromes-5 (CMS5) is a rare autosomal recessive heterogeneous disorder, caused by pathogenic variants in the COLQ that lead to skeletal muscle weakness and abnormal fatigability. The onset is usually from birth to childhood. Disease-causing variants in the collagen-like tail subunit are the most explained etiology in synaptic CMS, causing defected acetylcholinesterase. In this study whole-exome sequencing (WES) was performed in an affected boy with muscle weakness, ophthalmoplegia, and bilateral ptosis and gene expression assay by qRT-PCR was performed in entire family. A homozygous nonsense variant in the COLQ [NM_005677.4:c.679C>T], (p.Arg227Ter) was identified in the proband. Segregation analysis by Sanger sequencing confirmed the homozygous state in the proband and heterozygous state in his parents and four of the siblings. The mRNA expression level in the proband was 0.02 of a healthy person, and in the carriers were 0.42 of a healthy person. This study presents an Iranian family with two affected children and eight symptomatic carriers with attenuated mRNA expression. This study provides evidence that carriers of the COLQ disease-causing variants could become symptomatic with some yet unknown pathogenesis mechanism and underscore the importance of further investigations to elucidate this mechanism.

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

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

CiXTH29 and CiLEA4 Role in Water Stress Tolerance in Cichorium intybus Varieties

Drought causes massive crop quality and yield losses. Limiting the adverse effects of water deficits on crop yield is an urgent goal for a more sustainable agriculture. With this aim, six chicory varieties were subjected to drought conditions during seed germination and at the six week-old plant growth stage, in order to identify some morphological and/or molecular markers of drought resistance. Selvatica, Zuccherina di Trieste and Galatina varieties, with a high vegetative development, showed a major germination index, greater seedling development (6 days of growth) and a greater dehydration resistance (6 weeks of growth plus 10 days without water) than the other ones (Brindisina, Esportazione and Rossa Italiana). Due to the reported involvement, in the abiotic stress response, of xyloglucan endotransglucosylase/hydrolases (XTHs) and late embryogenesis abundant (LEA) multigene families, XTH29 and LEA4 expression profiles were investigated under stress conditions for all analyzed chicory varieties. We showed evidence that chicory varieties with high CiXTH29 and CiLEA4 basal expression and vegetative development levels better tolerate drought stress conditions than varieties that show overexpression of the two genes only in response to drought. Other specific morphological traits characterized almost all chicory varieties during dehydration, i.e., the appearance of lysigen cavities and a general increase of the amount of xyloglucans in the cell walls of bundle xylem vessels. Our results highlighted that high CiXTH29 and CiLEA4 basal expression, associated with a high level of vegetative growth, is a potential marker for drought stress tolerance.

Reference genes selection of Paeonia ostii 'Fengdan' under osmotic stresses and hormone treatments by RT-qPCR

BACKGROUND

Tree peony possess significant ornamental, medicinal and oil values. Osmotic stresses including dehydratiuon and salinity limit the expansion of cultivation area of tree peony. Information on reference genes selection under osmotic stress and hormone stimulation of tree peony still limited. This study aimed to determine the stable reference genes suitable for tree peony under osmotic stresses and hormone treatments, and provide a theoretical basis for the molecular biology research.

METHODS AND RESULTS

Twelve candidate reference genes were evaluated in Paeonia ostii 'Fengdan' under osmotic stress and hormone treatments by RT-qPCR. Delta Ct method, geNorm, and NormFinder were used for the comprehensive expression stability ranking comparison. The results revealed that tubulin-α was the preferred internal reference genes for drought and ABA treatment, tubulin-β was identified as the most suitable reference gene under drought and OPDA induction, 18s-rRNA was regarded as the most stable gene for salinity and JA treatment, eIF-5 A was listed as the most stable gene for JA and MeJA treatments. The experiments also displayed that EF1-α were comparatively unstable under ABA and BR hormone treatments.

CONCLUSION

These preferred reference genes could be useful in qPCR studies involving osmotic or hormonal stresses in Paeonia ostii 'Fengdan'. It is anticipated that the results will benefit tree peony functional genomics studies and molecular breeding research in the future.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of European eel, Anguilla anguilla

Eels are important aquaculture species for which an increasing number of reference genes are being identified and applied. In this study, five housekeeping genes [RPL7 (ribosomal protein L7), 18 S (18 S ribosomal RNA), EF1A (elongation factor 1α), ACTB (β-actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase)] were chosen to evaluate their reliability as reference genes for quantitative real-time PCR (qPCR) for the study of Anguilla anguilla. The expression of the selected genes in different eel tissues was determined using qPCR at different growth stages or upon challenge by Anguillid herpesvirus (AngHV), and the expression levels of these genes were then compared and evaluated using the geNorm and NormFinder algorithms. Then, RefFinder was used to comprehensively rank the examined housekeeping genes. Interestingly, the expression of the evaluated housekeeping genes exhibited tissue-dependent and treatment-dependent variations. In different growth periods A. anguilla tissues, the most stable genes were the following: ACTB in mucus; 18 S in skin and kidney; RPL7 in muscle, gill, intestine and brain; EF1A in heart and liver; and GAPDH in spleen. In contrast, in AngHV-challenged A. anguilla tissues, the most stable genes were the following: 18 S in mucus; RPL7 in skin, gill, heart, spleen, kidney and intestine; EF1A in muscle and liver; and ACTB in brain. Further comparison analysis indicated that the expression of RPL7 and EF1A was stable in multiple A. anguilla tissues in different growth periods and in eels challenged by AngHV. Nonetheless, the expression level of GAPDH in eel tissues was lower, and it was unstable in several tissues. These results indicated that the selection of reference genes for qPCR analysis in A. anguilla should be made in accordance with experimental parameters, and both RPL7 and EF1A could be used as reference genes for qPCR study of A. anguilla at different growth stages or upon challenge by AngHV. The reference genes identified in this study could improve the accuracy of qPCR data and facilitate further studies aimed at understanding the biology of eels.

HPRT1 Most Suitable Reference Gene for Accurate Normalization of mRNA Expression in Canine Dermal Tissues with Radiation Therapy

Reference genes are crucial in molecular biological studies as an internal control for gene re-search as they exhibit consistent expression patterns across many tissue types. In canines, radiation therapy is the most important therapeutic tool to cure various diseases like cancer. However, when using radiation for therapeutic strategy, radiation exposure to healthy tissues leads to some possible side effects such as acute radiation-induced skin injury and alters gene expression. Therefore, the analysis of a change in reference gene expression during the skin recovery process after radiation therapy is essential in healthy canine tissue. In the present study, we analyzed eight reference genes (ACTB, GAPDH, YWHAZ, GUSB, HPRT1, RPL4, RPS5, and TBP) in canine dermal tissues at 0, 1, 2, 3, 4, 5, 7, and 9 weeks of radiation exposure that affected the skin condition of canines. The stability of reference genes is determined by evaluating radiation therapy’s effect on healthy canine dermal tissue. Epidermal marker, Keratin 10 expression varies each week after irradiation, and HPRT1 is found to be the most suitable for normalization of mRNA expression in radiation-exposed canine dermal tissues. Changes in the gene expression level were evaluated by using a reliable tool such as quantitative real-time polymerase chain reaction (qRT-PCR). In order to achieve a valid qRT-PCR result, the most stable reference genes used for normalization after the radiation exposure process are important. Therefore, the current study was designed to evaluate the most stable reference gene for the post-irradiation canine tissues. After radiation exposure, the alternation of reference gene expression was estimated by three algorithms (geNorm, Normfinder, and Bestkeeper). The RG validation programs (GeNorm and NormFinder) suggested that HPRT1, RPL4, and TBP were suitable for normalization in qRT-PCR. Furthermore, three algorithms suggested that HPRT1 was the most stable reference gene for normalization with qRT-PCR results, regardless of before and after radiation exposure. Whereas GAPDH was found to be the most unstable reference gene. In addition, the use of stable or unstable reference genes for the normalization of Keratin 10 expression showed statistical differences. Therefore, we observed that, to obtain accurate and suitable PCR results of the canine tissues with and without radiation exposure, the HPRT1 reference gene is recommended for normalization with its high stability. Additionally, the use of RGs such as HPRT1, RPL4, and TBP for normalization in qRT-PCR experiments is recommended for post-radiation canine tissues to generate more accurate and reliable data. These results will provide fundamental information regarding internal controls for gene expression studies and can be used for the analysis of gene patterns in regenerative medicine.

Validation of Reference Genes in a Population of Blueberry (Vaccinium corymbosum) Plants Regenerated in Colchicine

For the first time we report the validation of reference genes in plants from a population of blueberry (Vaccinium corymbosum) clones cultured in vitro on a colchicine-supplemented medium. Nodal segment explants of the cultivar Duke were regenerated by organogenesis under different periods of colchicine 1 mg/L exposure (1, 2, 3, 5, 30 days). The clones selected for the study showed variability for phenotypic traits after 2 years of adaptation to field conditions, compared to plants of the donor genotype that were regenerated on a medium without colchicine. Vaccinium myrtillus (GAPDH) and Vaccinium macrocarpon (ATP1, NADH, RPOB and COX2) were used as reference genomes for primer design. The results show that colchicine treatments can cause genomic changes in blueberry plants. At the molecular level, exposure of plants to colchicine in early periods could promote an increase in gene expression of specific genes such as ATP1, COX2, GAPDH, MATK, NADH and RPOB. However, prolonged exposure (30 days) could decrease gene expression of the genes studied. For qPCR assays, the primers designed for ATP1, COX2, GAPDH and MATK genes showed high efficiency. In addition, the GAPDH, ATP1, NADH and COX2 genes showed high stability and could be recommended as potential reference genes for gene expression assays in Vaccinium.

Validation of Appropriate Reference Genes for qRT–PCR Normalization in Oat (Avena sativa L.) under UV-B and High-Light Stresses

Oat is a food and forage crop species widely cultivated worldwide, and it is also an important forage grass in plateau regions of China, where there is a high level of ultraviolet radiation and sunlight. Screening suitable reference genes for oat under UV-B and high-light stresses is a prerequisite for ensuring the accuracy of real-time quantitative PCR (qRT–PCR) data used in plant adaptation research. In this study, eight candidate reference genes (sulfite oxidase, SUOX; victorin binding protein, VBP; actin-encoding, Actin1; protein PSK SIMULATOR 1-like, PSKS1; TATA-binding protein 2-like, TBP2; ubiquitin-conjugating enzyme E2, UBC2; elongation factor 1-alpha, EF1-α; glyceraldehyde-3-phosphate dehydrogenase 1, GAPDH1;) were selected based on previous studies and our oat transcriptome data. The expression stability of these reference genes in oat roots, stems, and leaves under UV-B and high-light stresses was first calculated using three frequently used statistical software (geNorm, NormFinder, and BestKeeper), and then the comprehensive stability of these genes was evaluated using RefFinder. The results showed that the most stably expressed reference genes in the roots, stems, and leaves of oat under UV-B stress were EF1-α, TBP2, and PSKS1, respectively; the most stably expressed reference genes in the roots, stems, and leaves under high-light stress were PSKS1, UBC2, and PSKS1, respectively. PSKS1 was the most stably expressed reference gene in all the samples. The reliability of the selected reference genes was further validated by analysis of the expression of the phenylalanine ammonia-lyase (PAL) gene. This study highlights reference genes for accurate quantitative analysis of gene expression in different tissues of oat under UV-B and high-light stresses.

Gapdh Shows Altered Gene Expression in Alcohol Models

Background and objectives. Though glyceraldehyde-3-phosphate dehydrogenase (Gapdh) is one of the most commonly used housekeeping genes for comparison of gene expression data, studies have revealed that Gapdh expression is not constant at different developmental stages, and is modulated by many factors including ethanol. In view of this, in the present study, we investigated the effect of ethanol on Gapdh stability and expression levels in different model systems.

Methods. Evaluation of Gapdh stability was determined by comparison with three commonly used housekeeping genes in alcohol-exposed and control mouse blastocysts, embryos, and placentas (ActB, Hsp90ab1, and Atp5b), and in alcohol-dependent and healthy humans (ACTB, ATP5B, and HSPCB), and was ranked by the software program RefFinder. To detect the Gapdh expression patterns as a target gene, qRT-PCR analysis was applied.

Results. The order of expression stabilities obtained by the RefFinder ranked Gapdh as the most unstable reference gene in studied groups. The study of the expression dynamics of Gapdh showed significant upregulation in ethanol-exposed mouse blastocysts and embryos (FC = 1.56, p = 0.05 and FC = 6.0, p = 0.01) and downregulation in placentas (FC = 0.60, p = 0.01).

Conclusions. Our results revealed that the expression of Gapdh can vary in different model systems under ethanol exposure and imply the need for caution when using it as an internal control. This is the first report for statistically significant change in Gapdh expression after ethanol exposure during preimplantation mouse development.

Selection of suitable candidate genes for mRNA expression normalization in bulbil development of Pinellia ternata

Pinellia ternata (Thunb.) Breit. (Abbreviated as P. ternata). It is a commonly prescribed Chinese traditional medicinal herb for the treatment of phlegm, cough, and morning sick. Bulbil reproduction is one of the main reproductive methods of P. ternata. The accurate quantification of gene expression patterns associated with bulbil development might be helpful to explore the molecular mechanism involved in P. ternata reproduction. Quantitative real-time PCR was the most preferred method for expression profile and function analysis of mRNA. However, the reference genes in different tissues of P. ternata in different periods of bulbil development have not been studied in detail. In present study, the expression stability of eight candidate reference genes were determined with programs: geNorm, NormFinder, BestKeeper, and refFinder. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was identified as the top- rated reference gene in all samples of P. ternata, while different combinations of reference gene proved to be the most stable depending on development stage and tissue type. Furthermore, the reliability of GAPDH expression was verified by six P. ternata related genes in hormone and nutrient biosynthesis pathways, and the expression profiles of these genes were agreed with the results of RNA-seq digital gene expression analysis. These results can contribute to studies of gene expression patterns and functional analysis of P. ternata involved in bulbil development.

Stable reference gene selection for quantitative real-time PCR normalization in passion fruit (Passiflora edulis Sims.)

BACKGROUND

Passiflora edulis is a tropical fruit with high nutrient and medicinal values that is widely planted in southern China. However, the molecular biology of P. edulis has not been well studied. There are few reports regarding the choice of reference genes for gene expression studies of passion fruit.

METHODS AND RESULTS

By using three algorithms, implemented in geNorm, NormFinder and BestKeeper, we have selected ten candidate reference genes to explore their transcriptional expression stability in various tissues and under cold stress conditions. EF1 and HIS were stably expressed in five tissues. Ts and OTU were stably in vegetative organs. 50 S and Liom were stably in reproductive organs. The transcriptional abundance of EF1 and UBQ was stable in cold-treated and recovery treated leaf samples of P. edulis. In all samples, EF1 and Ts exhibited the highest expression stability. Evaluation of selected genes using simple statistical methods (ANOVA and post hoc analysis). Overall, EF1 emerged as the optimum reference gene for qRT-PCR normalize in P. edulis. In addition, the qRT-PCR analysis revealed that expression of ICE1 increases with the duration of cold treatment.

CONCLUSIONS

In this study, we successfully screened stable reference genes from 10 candidates in P. edulis and verified the results by analyzing the expression level of ICE1. The results provide reliable and effective reference genes for future research on gene expression analysis in P. edulis, and lay a foundation for follow-up research on functional genes in P. edulis.

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

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

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

Validation of Reference Genes for Quantitative Real-Time PCR Normalization in Ananas comosus var. bracteatus During Chimeric Leaf Development and Response to Hormone Stimuli

Reverse transcription quantitative real-time PCR (RT-qPCR) is a common way to study gene regulation at the transcriptional level due to its sensibility and specificity, but it needs appropriate reference genes to normalize data. Ananas comosus var. bracteatus, with white-green chimeric leaves, is an important pantropical ornamental plant. Up to date, no reference genes have been evaluated in Ananas comosus var. bracteatus. In this work, we used five common statistics tools (geNorm, NormFinder, BestKeeper, ΔCt method, RefFinder) to evaluate 10 candidate reference genes. The results showed that Unigene.16454 and Unigene.16459 were the optimal reference genes for different tissues, Unigene.16454 and zinc finger ran-binding domain-containing protein 2 (ZRANB2) for chimeric leaf at different developmental stages, isocitrate dehydrogenase [NADP] (IDH) and triacylglycerol lipase SDP1-like (SDP) for seedlings under different hormone treatments. The comprehensive results showed IDH, pentatricopeptide repeat-containing protein (PPRC), Unigene.16454, and caffeoyl-CoA O methyltransferase 5-like (CCOAOMT) are the top-ranked stable genes across all the samples. The stability of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the least during all experiments. Furthermore, the reliability of recommended reference gene was validated by the detection of porphobilinogen deaminase (HEMC) expression levels in chimeric leaves. Overall, this study provides appropriate reference genes under three specific experimental conditions and will be useful for future research on spatial and temporal regulation of gene expression and multiple hormone regulation pathways in Ananas comosus var. bracteatus.

Selection of the reference genes for quantitative gene expression by RT-qPCR in the desert plant Stipagrostis pennata

The desert pioneer plant Stipagrostis pennata plays an important role in sand fixation, wind prevention, and desert ecosystem recovery. An absence of reference genes greatly limits investigations into the regulatory mechanism by which S. pennata adapts to adverse desert environments at the molecular and genetic levels. In this study, eight candidate reference genes were identified from rhizosheath development transcriptome data from S. pennata, and their expression stability in the rhizosheaths at different development stages, in a variety of plant tissues, and under drought stress was evaluated using four procedures, including geNorm, NormFinder, BestKeeper, and RefFinder. The results showed that GAPDH and elF were the most stable reference genes under drought stress and in rhizosheath development, and ARP6 and ALDH were relatively stable in all plant tissues. In addition, elF was the most suitable reference gene for all treatments. Analysis of the consistency between the reverse transcription-quantitative PCR (RT-qPCR) and RNA sequencing data showed that the identified elF and GAPDH reference genes were stable during rhizosheath development. These results provide reliable reference genes for assuring the accuracy of RT-qPCR and offer a foundation for further investigations into the genetic responses of S. pennata to abiotic stress.

Comparative transcriptome analysis provides insights into steviol glycoside synthesis in stevia (Stevia rebaudiana Bertoni) leaves under nitrogen deficiency

Transcriptome analysis revealed the potential mechanism of nitrogen regulating steviol glycosides synthesis via shifting of leaf carbon metabolic flux or inducing certain transcription factors. Nitrogen (N) plays key regulatory roles in both stevia (Stevia rebaudiana) growth and the synthesis of its functional metabolite steviol glycosides (SGs), but the mechanism by which this nutrient regulates SGs synthesis remains to be elucidated. To address this question, a pot experiment was performed in a greenhouse where stevia plants fertilized with N (the control as CK plants) and compared with plants without the supply of N. Physiological and biochemical analyses were conducted to test the growth and metabolic responses of plants to N regimes. Our results showed that N deficiency significantly inhibited plant growth and leaf photosynthesis, while increased leaf SGs contents in stevia (49.97, 46.64 and 84.80% respectively for rebaudioside A, stevioside, and rebaudioside C), which may be partly due to "concentration effect". Then, transcriptome analysis was conducted to understand the underlying mechanisms. A total of 535 differentially expressed genes were identified, and carbon metabolism-related events were highlighted by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Many of these genes were significantly upregulated by N-deficiency, including those involved in "phenylpropanoid biosynthesis", "flavonoid biosynthesis" and "starch and sucrose metabolism". Our study also analyzed the expression patterns of SGs synthesis-related genes under two N regimes and the potential transcription factors linking N nutrition and SG metabolism. N-deficiency may promote SGs synthesis by changing the carbon metabolism flux or inducing certain transcription factors. Our results provide deeper insight into the relationship between N nutrition and SGs synthesis in stevia plants.

Reference genes for quantitative real-time PCR normalization of Cenostigma pyramidale roots under salt stress and mycorrhizal association

Cenostigma pyramidale is a native legume of the Brazilian semiarid region which performs symbiotic association with arbuscular mycorrhizal fungi (AMF), being an excellent model for studying genes associated with tolerance against abiotic and biotic stresses. In RT-qPCR approach, the use of reference genes is mandatory to avoid incorrect interpretation of the relative expression. This study evaluated the stability of ten candidate reference genes (CRGs) from C. pyramidale root tissues under salt stress (three collection times) and associated with AMF (three different times of salinity). The de novo transcriptome was obtained via RNA-Seq sequencing. Three algorithms were used to calculate the stability of CRGs under different conditions: (i) global (Salt, Salt+AMF, AMF and Control, and collection times), (ii) only non-inoculated plants, and (iii) AMF (only inoculated plants). HAG2, SAC1, aRP3 were the most stable CRGs for global and AMF assays, whereas HAG2, SAC1, RHS1 were the best for salt stress assay. This CRGs were used to validate the relative expression of two up-regulated transcripts in Salt2h (RAP2-3 and PIN8). Our study provides the first set of reference genes for C. pyramidale under salinity and AMF, supporting future researches on gene expression with this species.

Identification and selection of reference genes for gene expression analysis by quantitative real-time PCR in Suaeda glauca's response to salinity

Quantitative real-time polymerase chain reaction (qPCR) using a stable reference gene is widely used for gene expression research. Suaeda glauca L. is a succulent halophyte and medicinal plant that is extensively used for phytoremediation and extraction of medicinal compounds. It thrives under high-salt conditions, which promote the accumulation of high-value secondary metabolites. However, a suitable reference gene has not been identified for gene expression standardization in S. glauca under saline conditions. Here, 10 candidate reference genes, ACT7, ACT11, CCD1, TUA5, UPL1, PP2A, DREB1D, V-H⁺-ATPase, MPK6, and PHT4;5, were selected from S. glauca transcriptome data. Five statistical algorithms (ΔCq, geNorm, NormFinder, BestKeeper, and RefFinder) were applied to determine the expression stabilities of these genes in 72 samples at different salt concentrations in different tissues. PP2A and TUA5 were the most stable reference genes in different tissues and salt treatments, whereas DREB1D was the least stable. The two reference genes were sufficient to normalize gene expression across all sample sets. The suitability of identified reference genes was validated with MYB and AP2 in germinating seeds of S. glauca exposed to different NaCl concentrations. Our study provides a foundational framework for standardizing qPCR analyses, enabling accurate gene expression profiling in S. glauca.

Candidate reference genes for quantitative gene expression analysis in Lagerstroemia indica

Quantitative gene expression analysis by qPCR requires reference genes for normalization. Lagerstroemia indica (crape myrtle) is a popular ornamental plant in the world, but suitable endogenous reference genes are lacking. To find suitable reference genes, we evaluated the stabilities of nine candidate genes in six experimental data sets: six different tissues, three leaf colors, nine flower colors, and under three abiotic stresses (salt, drought, cold) using four statistical algorithms. A target gene LiMYB56 (homolog of Arabidopsis MYB56) was used to verify the authenticity and accuracy of the candidate reference genes. The results showed that the combination of two stably expressed reference genes, rather than a single reference gene, improved the accuracy of the qPCR. LiEF1α-2 + LiEF1α-3 was best for the tissue, salt treatment, and drought treatment sets; LiEF1α-2 + LiEF1α-1 was optimal for leaf color; LiEF1α-2 + LiACT7 was optimal for cold treatment; and LiUBC + LiEF1α-1 was best for the flower color set. Notably, LiEF1α-2 had high expression stability in all six experimental sets, implying it may be a good reference gene for expression studies in L. indica. Our results will facilitate future gene expression studies in L. indica.

Identification and evaluation of antioxidant and reference genes for quantitative real-time PCR in blood of Caiman latirostris

The quantitative real-time polymerase chain reaction (qPCR) has been one of the most promising approaches to perform rapid and accurate quantification of DNA in various biological systems. The aim of this study was to standardized the qPCR technique for the analysis of important genes involved in the main routes of antioxidant defense against reactive oxygen species (catalase: cat and superoxide dismutase: sod) and evaluate the stability of different reference genes in blood of Caiman latirostris hatchlings. The stability of the reference genes, β-actin, glyceraldehyde 3-phosphate dehydrogenase (gapdh) and ribosomal protein L8 (rpl8) was determined using the comparative ΔCt, NormFinder, geNorm, BestKeeper and RefFinder. Then, cat and sod genes were normalized with each reference gene and their mRNA abundances were determined through the qPCR. Stability of genes was ranked through the different methods in the following order: β-actin, rpl8 and gapdh , under normal physiological conditions. The results reveal that cat and sod genes present a similar relative mRNA abundance with β-actin and rpl8. This is the first report of the analysis of antioxidant mRNA as potential biomarkers of oxidative stress in blood for all crocodilians species. Besides, we determined the stability of different reference genes that can be used for normalization of mRNA abundance patterns in blood of C. latirostris, without the need to sacrifice the animals.

Asparagine accumulation in chicory storage roots is controlled by translocation and feedback regulation of asparagine biosynthesis in leaves

The presence of acrylamide (AA), a potentially carcinogenic and neurotoxic compound, in food has become a major concern for public health. AA in plant-derived food mainly arises from the reaction of the amino acid asparagine (Asn) and reducing sugars during processing of foodstuffs at high temperature. Using a selection of genotypes from the chicory (Cichorium intybus L.) germplasm, we performed Asn measurements in storage roots and leaves to identify genotypes contrasting for Asn accumulation. We combined molecular analysis and grafting experiments to show that leaf to root translocation controls Asn biosynthesis and accumulation in chicory storage roots. We could demonstrate that Asn accumulation in storage roots depends on Asn biosynthesis and transport from the leaf, and that a negative feedback loop by Asn on CiASN1 expression impacts Asn biosynthesis in leaves. Our results provide a new model for Asn biosynthesis in root crop species and highlight the importance of characterizing and manipulating Asn transport to reduce AA content in processed plant-based foodstuffs.

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

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

Expression profiling of the genes encoding ABA route components and the ACC oxidase isozymes in the senescing leaves of Populus tremula

Abscisic acid (ABA) triggers and regulates, while ethylene modulates autumn leaf senescence. The expression profiles of genes encoding ABA route components and the ACC oxidase isozymes were investigated in Populus tremula during the early and moderate stages of autumn leaf senescence. The targets of interest were Ptre-HAB1-like genes (Ptre-HAB1, Ptre-HAB3a and Ptre-HAB3b), the subclass 3 of Ptre-SnRK2s genes (Ptre-SnRK2.6a, Ptre-SnRK2.6b and Ptre-SnRK2.6b) and Ptre-RbohD1, Ptre-RbohF1, and Ptre-RbohF2 genes encoding the poplar components, which are counterparts of the ABA route key regulators or the counterparts of its secondary messengers, such as Homology to ABA-insensitive 1 (HAB1), Sucrose non-fermenting 1-related protein kinases 2 (SnRK2s) or Respiratory burst oxidase D and Respiratory burst oxidase F (RbohD and RbohF, respectively) in Arabidopsis, and Ptre-ACO3, Ptre-ACO5, and Ptre-ACO6 genes encoding ACC oxidase isozymes involved in ethylene biosynthesis. The fold change in their expression levels enabled to distinguish the distinct expression patterns for the following pairs of genes: Ptre-HAB3a and Ptre-SnRK2.6a, Ptre-HAB3b and Ptre-SnRK2.2, and Ptre-HAB1 and Ptre-SnRK2.6b, where each pair involves the genes encoding the negative and positive regulators of ABA route, respectively. Among the investigated genes, the fold change of expression was the highest for Ptre-ACO3, Ptre-ACO6, and Ptre-SnRK2.6b genes during both the studied stages, and additionally for Ptre-HAB1 and Ptre-RbohD1 genes during the moderate stage. In contrast, the Ptre-RbohF1 and Ptre-RbohF2 genes exhibited only the transient upregulation at the early stage of senescence. In an in vitro study, the ability of protein kinases Ptre-SnRK2.6a and Ptre-SnRK2.6b to phosphorylate the N-terminal regions of Ptre-RbohD1 and Ptre-RbohF2 was studied; the activity of Ptre-SnRK2.6b against the studied Ptre-Rbohs was noticeably lower than that exhibited by Ptre-SnRK2.6a. It seems that despite the high similarity of their polypeptides, Ptre-SnRK2.6a and Ptre-SnRK2.6b may play different biological roles; nonetheless, it requires in vivo confirmation. Surprisingly, the highest protein kinase activity against the Ptre-Rbohs was detected in the heterologous reaction with AT-SnRK2.6/OST1 which suggests that the discussed interactions are evolutionary conserved.

Reference gene selection for quantitative real-time PCR (qRT-PCR) expression analysis in Galium aparine L

To accurately evaluate expression levels of target genes, stable internal reference genes is required for normalization of quantitative real-time PCR (qRT-PCR) data. However, there have been no systematical investigation on the stability of reference genes used in the bedstraw weed, Galium aparine L. (BGA). In this study, the expression profiles of seven traditionally used reference genes, namely 18S, 28S, ACT, GAPDH, EF1α, RPL7 and TBP in BGA were assessed under both biotic (developmental time and tissue), and abiotic (temperature, regions and herbicide) conditions. Four analytical algorithms (geNorm, Normfinder, BestKeeper and the ΔCt method) were used to analyze the suitability of these genes as internal reference genes. RefFinder, a comprehensive analytical software, was used to rank the overall stability of the candidate genes. The optimal normalization internal control genes were ranked as: 28S and RPL7 were best for all the different experimental conditions (developmental stages, tissues, temperature, regions and herbicide treatment); 28S and RPL7 for developmental stages; TBP and GAPDH for different tissues; 28S and GAPDH were relatively stable for different temperature; 28S and TBP were suitable for herbicide treatment. A specific set of reference genes were recommended for each experimental condition in BGA.

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

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

Identification and Selection of Reference Genes for Quantitative Transcript Analysis in Corydalis yanhusuo

Corydalis yanhusuo is a medicinal plant frequently used in traditional Chinese medicine, which has effective medical effects in many aspects. Real-time polymerase chain reaction (RT-PCR) has been one of the most widely used methods in biosynthesis research due to its high sensitivity and quantitative properties in gene expression analysis. To obtain accurate normalization, reference genes are often selected in advance; however, no reference genes are available in C. yanhusuo. Herein, 12 reference gene candidates, named cyclophilin 2 (CYP2), elongation factor 1-α (EF1-α), protein phosphatase 2 (PP2A), SAND protein family (SAND), polypyrimidine tract-binding protein (PTBP), TIP41-like protein (TIP41), lyceraldehyde-3-phosphate hydrogenase (GAPDH), ubiquitin-conjugating enzyme 9 (UBC9), cyclophilin 1 (CYP1), tubulin beta (TUBA), thioredoxin (YLS8), and polyubiquitin 10 (UBQ10), were selected for stability analysis. After being treated with hormone, UV, salt, metal, oxidative, drought, cold (4 °C), and hot stresses (40 °C), the qRT-PCR data of the selected genes was analyzed with NormFinder, geNorm, and BestKeeper. The result indicated that GAPDH, SNAD, and PP2A were the top three most stable reference genes under most treatments. This study selected and validated reliable reference genes in C. yanhusuo under various environmental conditions, which can provide great help for future research on gene expression normalization in C. yanhusuo.

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

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

Transcription Factor Networks in Leaves of Cichorium endivia: New Insights into the Relationship Between Photosynthesis and Leaf Development

Cichorium endivia is a leafy crop closely related to Lactuca sativa that comprises two major botanical varieties characterized by a high degree of intraspecific morphological variation: var. latifolium with broad leaves (escarole) and var. crispum with narrow crisp curly leaves (endive). To investigate the relationship between leaf morphology and photosynthetic activity, escaroles and endives were used as a crop model due to the striking morphological diversity of their leaves. We constructed a leaf database for transcription factors (TFs) and photosynthesis-related genes from a refined C. endivia transcriptome and used RNA-seq transcriptomic data from leaves of four commercial endive and escarole cultivars to explore transcription factor regulatory networks. Cluster and gene co-expression network (GCN) analyses identified two main anticorrelated modules that control photosynthesis. Analysis of the GCN network topological properties identified known and novel hub genes controlling photosynthesis, and candidate developmental genes at the boundaries between shape and function. Differential expression analysis between broad and curly leaves suggested three novel TFs putatively involved in leaf shape diversity. Physiological analysis of the photosynthesis properties and gene expression studies on broad and curly leaves provided new insights into the relationship between leaf shape and function.

Selection and validation of appropriate reference genes for real-time quantitative PCR analysis in Momordica charantia

Real time quantitative reverse transcription PCR (RT-qPCR) has been attracting more attention for its high sensitivity in gene expression analysis. Given the widely use of RT-qPCR in normalization, it is playing a pivotal role for seeking suitable reference genes in different species. In current work, 12 candidate reference genes including Actin 2 (ACT2), Cyclophilin 2 (CYP2), Glyceraldehyde-3-phosphate dehydrogenase C2 (GAPC2), Elongation factor 1-α (EF1-α), Nuclear cap binding protein 20 (NCBP20), Serine/threonine-protein phosphatase PP2A (PP2A), Polypyrimidine tract-binding protein 1 (PTBP1), SAND family protein (SNAD), TIP41-like protein (TIP41), Tubulin beta-6 (TUB6), Ubiquitin-conjugating enzyme 9 (UBC9) and Glyceraldehyde-3-phosphatedehydrogenase (GAPDH) were screened from the transcriptome datasets of M. charantia. Afterwards, GeNorm, NormFinder and BestKeeper algorithms were applied to assess the expression stability of these 12 genes under different abiotic stresses including drought, cold, high-salt, hormone, UV, oxidative and metal stress. The results indicated that 12 selected genes exhibited various stability across the samples under different external stress conditions, but TIP41, PTBP1 and PP2A presented high stability among all the reference genes. To validate the suitability of the identified reference genes, the results of hormone subset were compared with RNA sequencing (RNA-seq) data, and the relative abundance of Ascorbate peroxidase 1（APX1）was used to confirm the reliability of the results. This work assesses the stability of reference genes in M. charantia under different abiotic stress conditions, which will be beneficent for accurate normalization of target genes in M. charantia.

Evaluation of reference genes for quantitative real-time PCR analysis of messenger RNAs and microRNAs in rainbow trout Oncorhynchus mykiss under heat stress

We assessed the expression stability of several messenger (m)RNAs and micro (mi)RNAs from liver and head kidney of rainbow trout Oncorhynchus mykiss using high-throughput RNA sequencing (RNA-seq) and miRNA-seq data. Additionally, four commonly used reference genes and one small non-coding RNA (u6) were also selected to identify ideal reference mRNAs and miRNAs for quantitative real-time (qrt)-PCR analysis of heat stress responses. GeNorm, NormFinder, BestKeeper and comparative ΔC_t were employed for analysis of qrt-PCR data to systematically assess the expression stability of candidate mRNAs and miRNAs and stability was ranked using geometric means. β-actin and ef1-α were the most stably expressed reference mRNAs in liver and head kidney, respectively and ssa-mir-26a-5p and ssa-mir-462b-5p were the most stably expressed miRNAs in these tissues. This is the first identification of appropriate reference mRNAs and miRNAs for qrt-PCR analysis of O. mykiss under heat stress.

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

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

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

Selection of reference genes for reverse transcription-qPCR analysis in the biomonitor macrophyte Bidens laevis L

The RT-qPCR has been the method used to analyze gene expression in plants but its benefits have not been completely exploited in the field of plants ecotoxicology when used as molecular biomarkers. The correct use of RT-qPCR demands to establish a certain number of reference genes (RG) which are expected to be invariable in their expression although it does not always happen. The main goals of this work were to: (1) analyze the stability of six potential RG, (2) establish the optimum number of RG, (3) select the most suitable RG to be applied in Bidens laevis under different test conditions and tissues and (4) confirm its convenience by normalizing the expression of one gene of interest under three different challenges. When all data were pooled together, the geNorm algorithm pointed out beta-actin and beta-tubulin (TUB) as the optimal RG pair while NormFinder algorithm selected nicotinamide adenine dinucleotide dehydrogenase (NADHD) and histone 3 (H3) as possessing the most invariable levels of expression. On the other hand, when data were grouped by tissues, ANOVA test selected H3 and TUB, while data grouped by conditions indicated that H3 and NADHD were the most stable RG under this analysis. Therefore, for a general-purpose set of RG, the overall analysis showed that a set of three RG would be optimum, and H3, TUB and NADHD were the selected ones. On the other hand, as RG can vary depending on the tissues or conditions, results achieved with ANOVA would be more reliable. Thus, appropriate normalization process would clearly need more than one RG.

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.

Evaluation of reference genes for RT-qPCR analysis in wild and cultivated Cannabis

RT-qPCR has been widely used for gene expression analysis in recent years. The accuracy of this technique largely depends on the selection of suitable reference genes. In order to facilitate gene expression analysis in wild and cultivated Cannabis, the expression stability of seven candidate reference genes (ACT2, 18S rRNA, GAPDH, UBQ, TUB, PP2A and EF1α) were assessed in leaves samples of different development stages and different organs of both wild and cultivated Cannabis in the present study. Their expression stabilities were evaluated through three software packages (GeNorm, Normfinder and Bestkeeper). Results showed that UBQ and EF1α were the highly ranked genes in different leaves samples, and PP2A was the most stable reference gene in different organs, while GAPDH was the least stable one. And the validation of the reference genes selected was further confirmed by the expression patterns of MDS and OLS.

Genome-wide identification of differentially expressed miRNAs induced by ethephon treatment in abscission layer cells of cotton (Gossypium hirsutum)

Cotton is one of the most important economic crops and its production is influenced by various adverse factors. Boll shedding at an inappropriate period causes the severe loss of cotton yield. Ethephon can promote the formation of abscission layer cells, resulting in boll shedding. To genome-widely investigate the differentially expressed miRNAs involved in the formation of abscission layer cells, two sRNA libraries were constructed using abscission tissues of cotton pedicels treated with ethephon or water. A total of 55,264,755 and 59,069,866 clean reads were generated in two libraries. 460 mature miRNAs were identified, including 301 known miRNAs and 159 novel miRNAs. Among them, one known gra-MIR530b and seven novel miRNAs (gar-novel-miR88-3p/gra-novel-miR8-3p, gar-novel-miR16-5p, gar-novel-miR100-3p, gar-novel-miR75-3p/gra-novel-miR92-3p,gar-novel-miR32-5p/gra-novel-miR30-5p, gar-novel-miR9-5p and gra-novel-miR61-3p) were expressed differentially in abscission tissues. Seven-two genes were predicted as potential targets of the eight miRNAs. Gene ontology (GO) analyses revealed that most of these target genes are auxin response factors (ARFs), GTPase activators, and the regulators of oxidative phosphorylation and photosynthesis. Hence, we speculate that these miRNAs might regulate cell separation and aging to affect the formation of bscission layer cells. Our findings will provide a new insight into the regulatory mechanism of miRNAs involved in boll shedding.

Inulin chain length modification using a transgenic approach opening new perspectives for chicory

Chicory capable of synthesizing long-chain inulin is of great interest. During the growing season, the sucrose-sucrose 1-fructosyltransferase (1-SST) activity is vital for production of long-chain inulin in chicory. With the purpose to increase inulin chain length, we employed Agrobacterium-mediated transformation method. Transgenic chicory plants (Cichorium intybus L. var. sativum) cv. 'Melci' has been developed to overexpress sucrose-sucrose 1-fructosyltransferase (1-SST) under the control of the CaMV 35S promoter. The integration of the T-DNA into the plant genome was confirmed by PCR on genomic DNA using gene-specific primers. Quantification of the 1-SST transcript expression level revealed that transgenic plants showed higher 1-SST expression than those in non-transgenic plants. Further analyses proved that the fructan content of the roots significantly increased in the transgenic plants. These results revealed that overexpression of the 1-SST, the key gene in inulin biosynthesis in chicory, might serve as a novel approach to develop plants with the long-chain inulin content.

Selection and Validation of Reference Genes for mRNA Expression by Quantitative Real-Time PCR Analysis in Neolamarckia cadamba

Neolamarckia cadamba is an economically-important fast-growing tree species in South China and Southeast Asia. As a prerequisite first step for future gene expression studies, we have identified and characterized a series of stable reference genes that can be used as controls for quantitative real time PCR (qRT-PCR) expression analysis in this study. The expression stability of 15 candidate reference genes in various tissues and mature leaves under different conditions was evaluated using four different algorithms, i.e., geNorm, NormFinder, BestKeeper and RefFinder. Our results showed that SAMDC was the most stable of the selected reference genes across the set of all samples, mature leaves at different photosynthetic cycles and under drought stress, whereas RPL10A had the most stable expression in various tissues. PGK and RPS25 were considered the most suitable reference for mature leaves at different developmental stages and under cold treatment, respectively. Additionally, the gene expression profiles of sucrose transporter 4 (NcSUT4), and 9-cis-epoxycarotenoid dioxygenase 3 (NcNCED3) were used to confirm the validity of candidate reference genes. Collectively, our study is the first report to validate the optimal reference genes for normalization under various conditions in N. cadamba and will benefit the future discovery of gene function in this species.

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.

Strategies to enhance the production of pinoresinol and its glucosides by endophytic fungus (Phomopsis sp. XP-8) isolated from Tu-chung bark

To improve the production yield of (+)-pinoresinol (Pin), (+)-pinoresinol monoglucoside (PMG), and (+)-pinoresinol diglucoside (PDG), different methods were conducted, including co-culture with resveratrol-producing Alternaria sp. MG1 spores and addition of Tu-chung in a medium at the start of cultivation, ultrasound treatment (40 kHZ, 10 min) on 5-day culture, and addition of ethanol and sodium butyrate on Day 3, followed by cultivation for an additional period of 2 days. At the end of the cultivation period (5 days), the liquid phase was collected for product analysis. Cells were collected for the determination of gene expression levels and then used in bioconversion using resting cells for another period of 2 days. The liquid phase was measured to determine the output of the target products and the expression levels of the key genes related to the biosynthesis of these compounds. Consequently, co-culture with Alternaria MG1 and addition of Tu-chung bark in the medium efficiently increased Pin, PMG, and PDG production yield in the biosynthesis systems using potato dextrose broth medium and resting cells of Phomopsis sp. XP-8. The key genes related to the biosynthesis of these compounds were significantly upregulated. However, in the majority of cases, the addition of ethanol and sodium butyrate, and ultrasound treatment decreased the production yield of Pin, PMG, and PDG. The change in production yield was not consistently accompanied by a change in gene expression.

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

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

Xylella fastidiosa induces differential expression of lignification related-genes and lignin accumulation in tolerant olive trees cv. Leccino

Recently, Xylella fastidiosa was reported in Italy, associated with the "Olive Quick Decline Syndrome". The cv. Leccino exhibits an evident tolerance with a slow disease progression compared with the other cultivars. Between the mechanisms proposed to explain the putative tolerance of some hosts to X. fastidiosa diseases, lignin deposition plays an important role. Analysis of phenolic compounds in healthy and infected Leccino and Cellina di Nardò leaves showed, in the two cultivars, a reduction of hydroxytyrosol glucoside (usually associated with drought and cold stress) and, only in Leccino, an increase of quinic acid, precursor of lignin. To determine if lignin biosynthesis is involved in defence response, we investigated the expression of genes coding for entry-point enzymes in different branches of the phenylpropanoid pathway. In stems of Cellina di Nardò infected plants, Cinnamate-4-Hydroxylase (C4H) and 4-Coumarate:CoA Ligase (4CL) resulted strongly down-regulated, indicating a plant disease response since the inhibition of C4H is reported to promote the accumulation of benzoic acid and salicylic acid as defence signals. Instead, in the cv. Leccino, Cinnamoyl-CoA Reductase (CCR, reported to be strongly induced during the formation of lignin defence response associated) was up-regulated in the stem of infected plants; moreover, Polyphenol oxidase (PPO), coding for an enzyme involved in the hydroxytyrosol biosynthesis, was down-regulated. The quantification of lignin in healthy and infected branches of both cultivars, showed a significant increase of total lignin in infected Leccino compared with the sensitive cultivar; moreover, histochemical observations of stem sections exhibited a different lignin distribution in the sclerenchyma and in the xylem tissue of infected Leccino plants compared to sections of healthy ones. Results suggest a critical role for lignin in X. fastidiosa tolerance of cv. Leccino.

Selection of reference genes for quantitative real-time PCR analysis in chicken ovary following silver nanoparticle treatment

Reliable results of quantitative real time PCR (qPCR) analysis require normalization of target gene expression level using reference genes (RGs). However housekeeping genes expression may vary under experimental conditions, so selection of the proper RGs is a crucial step in a qPCR analysis. Several algorithms have been developed to address this problem: geNorm, NormFinder and BestKeeper. In this study, we have used these three tools to evaluate the stability of RGs in the ovarian tissues of hens treated with silver nanoparticles. Eight genes were selected for the validation: HPRT, HMBS, VIM, SDHA, TBP, RPL13, GAPDH and 18S rRNA. According to geNorm the best combination of reference genes is SDHA and TPP. NormFinder also selected SDHA as the most suitable gene, but in combination with RPL13. Analysis in BestKeeper showed that SDHA, RPL13 might be the best choice in gene expression studies using the chicken ovary. In conclusion, the results obtained depend on the algorithm used and it arises from the diverse calculation strategies used in these programs. The outcome from the NormFinder is considered to be the most trustworthy and used in further qPCR analysis.

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

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

Identification of TMEM208 and PQLC2 as reference genes for normalizing mRNA expression in colorectal cancer treated with aspirin

Numerous evidences indicate that aspirin usage causes a significant reduction in colorectal cancer. However, the molecular mechanisms about aspirin preventing colon cancer are largely unknown. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a most frequently used method to identify the target molecules regulated by certain compound. However, this method needs stable internal reference genes to analyze the expression change of the targets. In this study, the transcriptional stabilities of several traditional reference genes were evaluated in colon cancer cells treated with aspirin, and also, the suitable internal reference genes were screened by using a microarray and were further identified by using the geNorm and NormFinder softwares, and then were validated in more cell lines and xenografts. We have showed that three traditional internal reference genes, β-actin, GAPDH and α-tubulin, are not suitable for studying gene transcription in colon cancer cells treated with aspirin, and we have identified and validated TMEM208 and PQLC2 as the ideal internal reference genes for detecting the molecular targets of aspirin in colon cancer in vitro and in vivo. This study reveals stable internal reference genes for studying the target genes of aspirin in colon cancer, which will contribute to identify the molecular mechanism behind aspirin preventing colon cancer.

Genome-Wide Analysis of the Biosynthesis and Deactivation of Gibberellin-Dioxygenases Gene Family in Camellia sinensis (L.) O. Kuntze

Gibberellins (GAs), a class of diterpenoid phytohormones, play a key role in regulating diverse processes throughout the life cycle of plants. Bioactive GA levels are rapidly regulated by Gibberellin-dioxygenases (GAox), which are involved in the biosynthesis and deactivation of gibberellin. In this manuscript, a comprehensive genome-wide analysis was carried out to find all GAox in Camellia sinensis. For the first time in a tea plant, 14 CsGAox genes, containing two domains, DIOX_N (PF14226) and 2OG-FeII_Oxy, were identified (PF03171). These genes all belong to 2-oxoglutarate-dependent dioxygenases (2-ODD), including four CsGA20ox (EC: 1.14.11.12), three CsGA3ox (EC: 1.14.11.15), and seven CsGA2ox (EC: 1.14.11.13). According to the phylogenetic classification as in Arabidopsis, the CsGAox genes spanned five subgroups. Each CsGAox shows tissue-specific expression patterns, although these vary greatly. Some candidate genes, which may play an important role in response to external abiotic stresses, have been identified with regards to patterns, such as CsGA20ox2, CsGA3ox2, CsGA3ox3, CsGA2ox1, CsGA2ox2, and CsGA2ox4. The bioactive GA levels may be closely related to the GA20ox, GA3ox and GA2ox genes. In addition, the candidate genes could be used as marker genes for abiotic stress resistance breeding in tea plants.