Selection of reference genes for quantitative real-time PCR expression studies in the apomictic and sexual grass Brachiaria brizantha

Brown-top millet: an overview of breeding, genetic, and genomic resources development for crop improvement

MAIN CONCLUSION

Brown-top millet is a lesser-known millet with a high grain nutrient value, early maturation, and drought tolerance that needs basic research to understand and conserve food security. Brown-top millet [Urochloa ramosa (L.)] is currently cultivated in some developing countries (especially in India) for food and fodder, although it is less known among the small millets. Like other millets, it contains macro- and micronutrients, vitamins, minerals, proteins, and fiber, all of which have rich health benefits. The nutritional importance and health benefits of brown-top millet are still unknown to many people due to a lack of awareness, wide cultivation, and research. Hence, this millet is currently overshadowed by other major cereals. This review article aims to present the nutritional, breeding, genetic, and genomic resources of brown-top millet to inform millet and other plant researchers. It is important to note that genetic and genomic resources have not yet been created for this millet. To date, there are no genomic and transcriptomic resources for brown-top millet to develop single nucleotide polymorphisms (SNP) and insertion/Deletions (InDels) for breeding studies. Furthermore, studies regarding nutritional significance and health benefits are required to investigate the exact nutritional contents and health benefits of the brown-top millet. The present review delves into the nutritional value and health advantages of brown-top millet, as supported by the available literature. The limitations of producing brown-top millet have been enumerated. We also cover the status of marker-assisted breeding and functional genomics research on closely related species. Lastly, we draw insights for further research such as developing omics resources and applying genome editing to study and improve brown-top millet. This review will help to start breeding and other molecular studies to increase the growth and development of this cereal.

Identification of Accurate Reference Genes for qRT-PCR Analysis of Gene Expression in Eremochloa ophiuroides under Multiple Stresses of Phosphorus Deficiency and/or Aluminum Toxicity

Centipedegrass (Eremochloa ophiuroides (Munro.) Hack.) is a species originating in China and is an excellent warm-season turfgrass. As a native species in southern China, it is naturally distributed in the phosphorus-deficient and aluminum-toxic acid soil areas. It is important to research the molecular mechanism of centipedegrass responses to phosphorus-deficiency and/or aluminum-toxicity stress. Quantitative Real-Time PCR (qRT-PCR) is a common method for gene expression analysis, and the accuracy of qRT-PCR results depends heavily on the stability of internal reference genes. However, there are still no reported stable and effective reference genes for qRT-PCR analysis of target genes under the acid-soil-related stresses in different organs of centipedegrass. For scientific rigor, the gene used as a reference for any plant species and/or any stress conditions should be first systematically screened and evaluated. This study is the first to provide a group of reliable reference genes to quantify the expression levels of functional genes of Eremochloa ophiuroides under multiple stresses of P deficiency and/or aluminum toxicity. In this study, centipedegrass seedlings of the acid-soil-resistant strain 'E041' and acid-soil-sensitive strain 'E089' were used for qRT-PCR analysis. A total of 11 candidate reference genes (ACT, TUB, GAPDH, TIP41, CACS, HNR, EP, EF1α, EIF4α, PP2A and actin) were detected by qRT-PCR technology, and the stability of candidate genes was evaluated with the combination of four internal stability analysis software programs. The candidate reference genes exhibited differential stability of expression in roots, stems and leaves under phosphorus-deficiency and/or aluminum-toxicity stress. On the whole, the results showed that GAPDH, TIP41 and HNR were the most stable in the total of samples. In addition, for different tissues under various stresses, the selected reference genes were also different. CACS and PP2A were identified as two stable reference genes in roots through all three stress treatments (phosphate deficiency, aluminum toxicity, and the multiple stress treatment of aluminum toxicity and phosphate deficiency). Moreover, CACS was also stable as a reference gene in roots under each treatment (phosphate deficiency, aluminum toxicity, or multiple stresses of aluminum toxicity and phosphate deficiency). In stems under all three stress treatments, GAPDH and EIF4α were the most stable reference genes; for leaves, PP2A and TIP41 showed the two highest rankings in all three stress treatments. Finally, qRT-PCR analysis of the expression patterns of the target gene ALMT1 was performed to verify the selected reference genes. The application of the reference genes identified as internal controls for qRT-PCR analysis will enable accurate analysis of the target gene expression levels and expression patterns in centipedegrass under acid-soil-related stresses.

An exonuclease V homologue is expressed predominantly during early megasporogenesis in apomictic Brachiaria brizantha

MAIN CONCLUSION

An exonuclease V homologue from apomictic Brachiaria brizantha is expressed and localized in nucellar cells at key moments when these cells differentiate to give rise to unreduced gametophytes. Brachiaria is a genus of forage grasses with economical and agricultural importance to Brazil. Brachiaria reproduces by aposporic apomixis, in which unreduced embryo sacs, derived from nucellar cells, other than the megaspore mother cell (MMC), are formed. The unreduced embryo sacs produce an embryo without fertilization resulting in clones of the mother plant. Comparative gene expression analysis in ovaries of sexual and apomictic Brachiaria spp. revealed a sequence from B. brizantha that showed a distinct pattern of expression in ovaries of sexual and apomictic plants. In this work, we describe a gene named BbrizExoV with strong identity to exonuclease V (Exo V) genes from other grasses. Sequence analysis in signal prediction tools showed that BbrizExoV might have dual localization, depending on the translation point. A longer form to the nucleus and a shorter form which would be directed to the chloroplast. This is also the case for monocot sequences analyzed from other species. The long form of BbrizExoV protein localizes to the nucleus of onion epidermal cells. Analysis of ExoV proteins from dicot species, with exception of Arabidopsis thaliana ExoVL protein, showed only one localization. Using a template-based AlphaFold 2 modelling approach the structure of BbrizExoV in complex with metal and ssDNA was predicted based on the holo structure of the human counterpart. Features predicted to define ssDNA binding but a lack of sequence specificity are shared between the human enzyme and BbrizExoV. Expression analyses indicated the precise site and timing of transcript accumulation during ovule development, which coincides with the differentiation of nucelar cells to form the typical aposporic four-celled unreduced gametophyte. A putative function for this protein is proposed based on its homology and expression pattern.

Identification of IPT9 in Brachiaria brizantha (syn. Urochloa brizantha) and expression analyses during ovule development in sexual and apomictic plants

BACKGROUND

In Brachiaria sexual reproduction, during ovule development, a nucellar cell differentiates into a megaspore mother cell (MMC) that, through meiosis and mitosis, gives rise to a reduced embryo sac. In aposporic apomictic Brachiaria, next to the MMC, other nucellar cells differentiate into aposporic initials that enter mitosis directly forming an unreduced embryo sac. The IPT (isopentenyltransferase) family comprises key genes in the cytokinin (CK) pathway which are expressed in Arabidopsis during ovule development. BbrizIPT9, a B. brizantha (syn. Urochloa brizantha) IPT9 gene, highly similar to genes of other Poaceae plants, also shows similarity with Arabidopsis IPT9, AtIPT9. In this work, we aimed to investigate association of BbrizIPT9 with ovule development in sexual and apomictic plants.

METHODS AND RESULTS

RT-qPCR showed higher BbrizIPT9 expression in the ovaries of sexual than in the apomictic B. brizantha. Results of in-situ hybridization showed strong signal of BbrizIPT9 in the MMC of both plants, at the onset of megasporogenesis. By analyzing AtIPT9 knockdown mutants, we verified enlarged nucellar cell, next to the MMC, in a percentage significantly higher than in the wild type, suggesting that knockout of AtIPT9 gene triggered the differentiation of extra MMC-like cells.

CONCLUSIONS

Our results indicate that AtIPT9 might be involved in the proper differentiation of a single MMC during ovule development. The expression of a BbrizIPT9, localized in male and female sporocytes, and lower in apomicts than in sexuals, and effect of IPT9 knockout in Arabidopsis, suggest involvement of IPT9 in early ovule development.

Ovule Transcriptome Analysis Discloses Deregulation of Genes and Pathways in Sexual and Apomictic Limonium Species (Plumbaginaceae)

The genus Limonium Mill. (sea lavenders) includes species with sexual and apomixis reproductive strategies, although the genes involved in these processes are unknown. To explore the mechanisms beyond these reproduction modes, transcriptome profiling of sexual, male sterile, and facultative apomictic species was carried out using ovules from different developmental stages. In total, 15,166 unigenes were found to be differentially expressed with apomictic vs. sexual reproduction, of which 4275 were uniquely annotated using an Arabidopsis thaliana database, with different regulations according to each stage and/or species compared. Gene ontology (GO) enrichment analysis indicated that genes related to tubulin, actin, the ubiquitin degradation process, reactive oxygen species scavenging, hormone signaling such as the ethylene signaling pathway and gibberellic acid-dependent signal, and transcription factors were found among differentially expressed genes (DEGs) between apomictic and sexual plants. We found that 24% of uniquely annotated DEGs were likely to be implicated in flower development, male sterility, pollen formation, pollen-stigma interactions, and pollen tube formation. The present study identifies candidate genes that are highly associated with distinct reproductive modes and sheds light on the molecular mechanisms of apomixis expression in Limonium sp.

Screening and validating of endogenous reference genes in Chlorella sp. TLD 6B under abiotic stress

Chlorella sp. TLD 6B, a microalgae growing in the Taklamakan Desert, Xinjiang of China, is a good model material for studying the physiological and environmental adaptation mechanisms of plants in their arid habitats, as its adaptation to the harsh desert environment has led to its strong resistance. However, when using real-time quantitative polymerase chain reaction (RT-qPCR) to analyze the gene expression of this algae under abiotic stress, it is essential to find the suitable endogenous reference genes so to obtain reliable results. This study assessed the expression stability of 9 endogenous reference genes of Chlorella sp. TLD 6B under four abiotic stresses (drought, salt, cold and heat). These genes were selected based on the analysis results calculated by the three algorithmic procedures of geNorm, NormFinder, and BestKeeper, which were ranked by refinder. Our research showed that 18S and GTP under drought stress, 18S and IDH under salt stress, CYP and 18S under cold stress, GTP and IDH under heat stress were the most stable endogenous reference genes. Moreover, UBC and 18S were the most suitable endogenous reference gene combinations for all samples. In contrast, GAPDH and α-TUB were the two least stable endogenous reference genes in all experimental samples. Additionally, the selected genes have been verified to be durable and reliable by detecting POD and PXG3 genes using above endogenous reference genes. The identification of reliable endogenous reference genes guarantees more accurate RT-qPCR quantification for Chlorella sp. TLD 6B, facilitating functional genomics studies of deserts Chlorella as well as the mining of resistance genes.

Selection of suitable internal control gene for assaying gene expression in rice through qRT-PCR during sheath blight infection

qRT-PCR is a globally accepted technique for assaying gene expression in relative terms which compares the difference between critical threshold (C_T) values of a gene calculated form two independently isolated RNA samples. Independent RNA isolations, however, include error due to batch effect which must be normalized for error-free calculation of relative gene expression. Hence, C_T values of internal control (IC) genes are used for normalization during the calculation of expression fold-change in gene expression analysis. The expression of ICs genes expected to be stable in all the experimental conditions. However, it is almost impossible to find such a gene which do not depict expression fluctuation in response to the changes in experimental conditions. Hence, it is necessary to identify suitable IC gene(s) for any given experimental condition before conducting any particular gene expression study. Here, we examined the suitability of eight candidate IC genes, namely glyceraldehyde 3-phosphate dehydrogenase (GAPDH), eukaryotic elongation factor-1 (eEF-1α), 25 S rRNA (25 S), 18 S rRNA (18 S), ubiquitin C E2 ligase (UBC), Actin (Act), ubiquitin 5 (UBQ5) and ubiquitin 10 (UBQ10), for assaying gene expression in rice during sheath blight infection. Our analysis suggest that GAPDH might be the IC of choice when expression studies include contrasting genotypes differing in their tolerance to sheath blight pathogen as well as progressive infection time. While if expression analysis have to be performed only in one genotype but under progressive sheath blight infection, UBQ5 might be chosen as IC because of its high expression stability under the proposed experimental setup.

Diverged subpopulations in tropical Urochloa (Brachiaria) forage species indicate a role for facultative apomixis and varying ploidy in their population structure and evolution

BACKGROUND

Urochloa (syn. Brachiaria) is a genus of tropical grasses sown as forage feedstock, particularly in marginal soils. Here we aimed to clarify the genetic diversity and population structure in Urochloa species to understand better how population evolution relates to ploidy level and occurrence of apomictic reproduction.

METHODS

We explored the genetic diversity of 111 accessions from the five Urochloa species used to develop commercial cultivars. These accessions were conserved from wild materials collected at their centre of origin in Africa, and they tentatively represent the complete Urochloa gene pool used in breeding programmes. We used RNA-sequencing to generate 1.1 million single nucleotide polymorphism loci. We employed genetic admixture, principal component and phylogenetic analyses to define subpopulations.

RESULTS

We observed three highly differentiated subpopulations in U. brizantha, which were unrelated to ploidy: one intermixed with U. decumbens, and two diverged from the former and the other species in the complex. We also observed two subpopulations in U. humidicola, unrelated to ploidy; one subpopulation had fewer accessions but included the only characterized sexual accession in the species. Our results also supported a division of U. decumbens between diploids and polyploids, and no subpopulations within U. ruziziensis and U. maxima.

CONCLUSIONS

Polyploid U. decumbens are more closely related to polyploid U. brizantha than to diploid U. decumbens, which supports the divergence of both polyploid groups from a common tetraploid ancestor and provides evidence for the hybridization barrier of ploidy. The three differentiated subpopulations of apomictic polyploid U. brizantha accessions constitute diverged ecotypes, which can probably be utilized in hybrid breeding. Subpopulations were not observed in non-apomictic U. ruziziensis. Sexual Urochloa polyploids were not found (U. brizantha, U. decumbens) or were limited to small subpopulations (U. humidicola). The subpopulation structure observed in the Urochloa sexual-apomictic multiploidy complexes supports geographical parthenogenesis, where the polyploid genotypes exploit the evolutionary advantage of apomixis, i.e. uniparental reproduction and clonality, to occupy extensive geographical areas.

Trends in Apomixis Research: The 10 Most Cited Research Articles Published in the Pregenomic and Genomic Eras

Apomixis, or asexual reproduction by seed, represents an easy shortcut for life cycle renewal based on maternal embryo production without ploidy reduction (meiosis) and ploidy restitution (syngamy). Although the first studies officially published on this topic in scientific journals date back to the early 1930s, the identification and introduction of genes involved in asexual reproduction in species of agronomic interest still represent a major challenge. Through a bibliometric analysis of the research programs implemented in apomixis over the last 40 years, the present study was aimed to discuss not only the main findings achieved but also the investigational methods and model species used. We split the critical survey of the most cited original articles into pregenomic and genomic eras to identify potential trends and depict scenarios that have emerged in the scientific community working on apomixis, as well as to determine any correlation between the exponential increase in acquired basic knowledge and the development of advanced analytical technologies. This review found a substantial stagnation in the use of the same model species, with few exceptions, for at least 40 years. In contrast, the development of new molecular techniques, genomic platforms, and repositories has directly affected the approaches used in research, which has been directed toward an increasingly focused study of the genetic and epigenetic determinants of apomixis.

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.

An Overview of the Genetics and Genomics of the Urochloa Species Most Commonly Used in Pastures

Pastures based on perennial monocotyledonous plants are the principal source of nutrition for ruminant livestock in tropical and subtropical areas across the globe. The Urochloa genus comprises important species used in pastures, and these mainly include Urochloa brizantha, Urochloa decumbens, Urochloa humidicola, and Urochloa ruziziensis. Despite their economic relevance, there is an absence of genomic-level information for these species, and this lack is mainly due to genomic complexity, including polyploidy, high heterozygosity, and genomes with a high repeat content, which hinders advances in molecular approaches to genetic improvement. Next-generation sequencing techniques have enabled the recent release of reference genomes, genetic linkage maps, and transcriptome sequences, and this information helps improve our understanding of the genetic architecture and molecular mechanisms involved in relevant traits, such as the apomictic reproductive mode. However, more concerted research efforts are still needed to characterize germplasm resources and identify molecular markers and genes associated with target traits. In addition, the implementation of genomic selection and gene editing is needed to reduce the breeding time and expenditure. In this review, we highlight the importance and characteristics of the four main species of Urochloa used in pastures and discuss the current findings from genetic and genomic studies and research gaps that should be addressed in future research.

Expression stability of internal reference gene in response to Trichoderma polysporum infection in Avena fatua L

Trichoderma polysporum was a pathogenic fungi which showed strong pathogenicity to Avena fatua L. in recent study. The stress response of A. fatua to T. polysporum is mediated by the regulation of gene expression. Quantification of the expression of genes requires normalizing RT-qPCR data using reference genes with stable expression in the system studied as internal standards. To construct a RT-qPCR system suitable for response of A. fatua to T. polysporum, and screen stable internal reference genes, GeNorm, NormFinder, BestKeeper and RefFinde were used to perform SYBR Green-based RT-qPCR analysis on eight candidate internal reference genes (18S, 28S, TUA, UBC, ACT, GAPDH, TBP and EF-1α) in A. fatua samples after inoculation of T. polysporum Strain HZ-31. The results showed that TBP, 18S and UBC were the most stable internal reference genes, TBP and TUA, TBP and GAPDH, 18S and TBP, UBC and 18S were the most suitable combination of the two internal reference genes, which could be used as internal reference genes for functional gene expression analysis during the interaction between T. polysporum and A. fatua. This is the first study describing a set of reference genes with a stable expression under fungi stress in A. fatua. These genes are also candidate reference genes of choice for studies seeking to identify stress-responsive genes in A. fatua.

Selection and validation of appropriate reference genes for RT-qPCR analysis of flowering stages and different genotypes of Iris germanica L

Iris germanica L. is a perennial herbaceous plant that has been widely cultivated worldwide and is popular for its elegant and vibrantly colorful flowers. Selection of appropriate reference genes is the prerequisite for accurate normalization of target gene expression by quantitative real-time PCR. However, to date, the most suitable reference genes for flowering stages have not been elucidated in I. germanica. In this study, eight candidate reference genes were examined for the normalization of RT-qPCR in three I. germanica cultivars, and their stability were evaluated by four different algorithms (GeNorm, NormFinder, BestKeeper, and Ref-finder). The results revealed that IgUBC and IgGAPDH were the most stable reference genes in '00246' and 'Elizabeth', and IgTUB and IgUBC showed stable expression in '2010200'. IgUBC and IgGAPDH were the most stable in all samples, while IgUBQ showed the least stability. Finally, to validate the reliability of the selected reference genes, the expression patterns of IgFT (Flowering Locus T gene) was analyzed and emphasized the importance of appropriate reference gene selection. This work presented the first systematic study of reference genes selection during flower bud development and provided guidance to research of the molecular mechanisms of flowering stages in I. germanica.

Evaluation of suitable reference genes for normalization of quantitative real-time PCR analysis in rice plants under Xanthomonas oryzae pv. oryzae--infection and melatonin supplementation

Exogenous melatonin (MT) was found to be an interesting tool for enhancing the resistance of rice to Xanthomonasoryzaepv. oryzae (Xoo)-caused bacterial blight (BB). However, the accurate comparison of the expression levels across samples was a challenging task. In this work, the stability of 10 common used housekeeping genes under Xoo-infection and MT supplementation in rice was analyzed using quantitative real-time PCR (qRT-PCR), and algorithms geNorm, NormFinder and BestKeeper. Our results indicated that most reference genes remained stable in Xoo-infected rice plants, while a number of reference genes were affected by MT supplementation. Among all studied genes, the transcript levels of 18S(18S ribosomal RNA) and UBC (Ubiquitin-conjugating enzyme E2) remained unaltered by Xoo infection, while UBC and UBQ5(Ubiquitin 5) were the most stable genes when examining simultaneous Xoo-infection and MT supplementation, demonstrating that UBC is a suitable reference gene for qRT-PCR data normalization in rice under Xoo-infection and MT supplementation.

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

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

SERK genes identification and expression analysis during somatic embryogenesis and sporogenesis of sexual and apomictic Brachiaria brizantha (Syn. Urochloa brizantha)

In Brachiaria brizantha BbrizSERK1, BbrizSERK2 and BbrizSERK3 were identified. SERK expression marks somatic embryogenesis, sexual MMC, and sexual and apomictic PMC. BbrizSERK3 might have a regulatory role in reproductive development. Somatic embryogenesis receptor-like kinase (SERK) consists of plasma membrane receptor genes that have been characterized in various species, associated with several aspects of plant development, including reproduction. SERK genes are involved in anther development and in early embryo development in sexual and asexual seed formation. To comprehend the complexity of the SERK genes and their function in Brachiaria reproduction, we performed a homology-based search in a genomic database of a sexual B. brizantha and identified sequences of three SERK genes, BbrizSERK1, BbrizSERK2, and BbrizSERK3. RNASeq data showed equivalent abundance of BbrizSERK1 and BbrizSERK2 transcripts in ovaries at early megasporogenesis of sexuals and apomicts, while BbrizSERK3 transcripts were more abundant in ovaries of sexuals than in apomicts. BbrizSERK3 results in three coding sequences due to alternative splicing, among them Variant 1 results in a protein with all the predicted domains of a SERK. BbrizSERK transcripts were detected in male reproductive tissues of both sexual and apomictic plants, suggesting a role in controlling anther development. BbrizSERK transcripts were detected early in ovule development, in the integuments, and in the megaspore mother cell of the sexual plant, but not in the cells that give rise to apomictic embryo sacs, suggesting a role in female reproductive development of sexuals. This paper provides evidences that SERK genes plays a role in the onset and establishment of somatic embryogenesis and in the reproductive development of B. brizantha and suggests a distinct role of BbrizSERK in apomixis initiation.

Public-transcriptome-database-assisted selection and validation of reliable reference genes for qRT-PCR in rice

Accurate quantitative reverse transcription PCR (qRT-PCR) requires reliable reference genes whose expression does not vary in different tissues and developmental stages. However, few reliable reference genes are available for qRT-PCR in rice (Oryza sativa). Here, we established an effective strategy for identifying novel reference genes (NRGs) for reliable normalization of qRT-PCR data in various rice organs and developmental stages. We selected candidate NRGs using the Information Commons for Rice Database and confirmed their expression in Rice Expression Profile Database (RiceXPro) data. Genes with low variation (<2.5 cycle quantification) across tissues and developmental stages, and little fluctuation in expression in heatmaps from RiceXPro data were considered stable NRGs. To validate this strategy, we selected 11 candidate NRGs and calculated their expression stability in different spatio-temporal conditions using five programs, and compared these genes with five established reference genes (ERGs). Only one of the ERGs (UBQ5) was reliable and 10 of the candidate NRGs were more stable than the four remaining ERGs. Therefore, public transcriptomic databases are useful for identifying NRGs. We selected two NRGs, UFC1 (Homolog of UFM1-Conjugating Enzyme 1) and FhaB (Homolog of Adhesin FhaB) for qRT-PCR analysis in rice; their homologs might be suitable for other monocot plants.

Selection of Suitable Reference Genes for RT-qPCR Gene Expression Analysis in Siberian Wild Rye (Elymus sibiricus) under Different Experimental Conditions

Elymus sibiricus, which is a perennial and self-pollinated grass, is the typical species of the genus Elymus, which plays an important role in forage production and ecological restoration. No reports have, so far, systematically described the selection of optimal reference genes for reverse transcriptase quantitative real-time polymerase chain reaction (RT-qPCR) analysis in E. sibiricus. The goals of this study were to evaluate the expression stability of 13 candidate reference genes in different experimental conditions, and to determine the appropriate reference genes for gene expression analysis in E. sibiricus. Five methods including Delta Ct (ΔCt), BestKeeper, NormFinder, geNorm, and RefFinder were used to assess the expression stability of 13 potential reference genes. The results of the RefFinder analysis showed that TBP2 and HIS3 were the most stable reference genes in different genotypes. TUA2 and PP2A had the most stable expression in different developmental stages. TBP2 and PP2A were suitable reference genes in different tissues. Under salt stress, ACT2 and TBP2 were identified as the most stable reference genes. ACT2 and TUA2 showed the most stability under heat stress. For cold stress, PP2A and ACT2 presented the highest degree of expression stability. DNAJ and U2AF were considered as the most stable reference genes under osmotic stress. The optimal reference genes were selected to investigate the expression pattern of target gene CSLE6 in different conditions. This study provides suitable reference genes for further gene expression analysis using RT-qPCR in E. sibiricus.

Comparison of Reliable Reference Genes Following Different Hormone Treatments by Various Algorithms for qRT-PCR Analysis of Metasequoia

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is the most sensitive technique for evaluating gene expression levels. Choosing appropriate reference genes for normalizing target gene expression is important for verifying expression changes. Metasequoia is a high-quality and economically important wood species. However, few systematic studies have examined reference genes in Metasequoia. Here, the expression stability of 14 candidate reference genes in different tissues and following different hormone treatments were analyzed using six algorithms. Candidate reference genes were used to normalize the expression pattern of FLOWERING LOCUS T and pyrabactin resistance-like 8. Analysis using the GrayNorm algorithm showed that ACT2 (Actin 2), HIS (histone superfamily protein H3) and TATA (TATA binding protein) were stably expressed in different tissues. ACT2, EF1α (elongation factor-1 alpha) and HIS were optimal for leaves treated with the flowering induction hormone solution, while Cpn60β (60-kDa chaperonin β-subunit), GAPDH (glyceraldehyde-3-phosphate dehydrogenase) and HIS were the best reference genes for treated buds. EF1α, HIS and TATA were useful reference genes for accurate normalization in abscisic acid-response signaling. Our results emphasize the importance of validating reference genes for qRT-PCR analysis in Metasequoia. To avoid errors, suitable reference genes should be used for different tissues and hormone treatments to increase normalization accuracy. Our study provides a foundation for reference gene normalization when analyzing gene expression in Metasequoia.

Characterization of Laguncularia racemosa transcriptome and molecular response to oil pollution

Mangroves are ecosystems of economic and ecological importance. Laguncularia racemosa (Combretaceae), popularly known as white mangrove, is a species that greatly contributes to the community structure of neotropical and West African mangrove forests. Despite the significance of these ecosystems, they have been destroyed by oil spills that can cause yellowing of leaves, increased sensitivity to other stresses and death of trees. However, the molecular response of plants to oil stress is poorly known. In this work, Illumina reads were de novo assembled into 46,944 transcripts of L. racemosa roots and leaves, including putative isoform variants. In addition to improving the genomic information available for mangroves, the L. racemosa assembled transcriptome allowed us to identify reference genes to normalize quantitative real-time PCR (qPCR) expression data from oil-stressed mangrove plants, which were used in RNASeq validation. The analysis of expression changes induced by the oil exposure revealed 310 and 286 responsive transcripts of leaves and roots, respectively, mainly up-regulated. Enriched GO categories related to chloroplasts and photosynthesis were found among both leaf and root oil-responsive transcripts, while "response to heat" and "response to hypoxia" were exclusively enriched in leaves and roots, respectively. The comparison of L. racemosa 12-h-oil-stressed leaf expression profile to previous Arabidopsis heat-stress studies and co-expression evidence also pointed to similarities between the heat and oil responses, in which the HSP-coding genes seem to play a key role. A subset of the L. racemosa oil-responsive root genes exhibited similar up-regulation profiles to their Arabidopsis homologs involved in hypoxia responses, including the HRA1 and LBD41 TF-coding genes. Genes linked to the ethylene pathway such as those coding for ERF TFs were also modulated during the L. racemosa root response to oil stress. Taken together, these results show that oil contamination affects photosynthesis, protein metabolism, hypoxia response and the ethylene pathway in L. racemosa 12-h-oil-exposed leaves and roots.

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.

A banana NAC transcription factor (MusaSNAC1) impart drought tolerance by modulating stomatal closure and H2O2 content

MusaSNAC1 function in H₂O₂ mediated stomatal closure and promote drought tolerance by directly binding to CGT[A/G] motif in regulatory region of multiple stress-related genes. Drought is a abiotic stress-condition, causing reduced plant growth and diminished crop yield. Guard cells of the stomata control photosynthesis and transpiration by regulating CO₂ exchange and water loss, thus affecting growth and crop yield. Roles of NAC (NAM, ATAF1/2 and CUC2) protein in regulation of stress-conditions has been well documented however, their control over stomatal aperture is largely unknown. In this study we report a banana NAC protein, MusaSNAC1 which induced stomatal closure by elevating H₂O₂ content in guard cells during drought stress. Overexpression of MusaSNAC1 in banana resulted in higher number of stomata closure causing reduced water loss and thus elevated drought-tolerance. During drought, expression of GUS (β-glucuronidase) under P _MusaSNAC1 was remarkably elevated in guard cells of stomata which correlated with its function as a transcription factor regulating stomatal aperture closing. MusaSNAC1 is a transcriptional activator belonging to SNAC subgroup and its 5'-upstream region contain multiple Dof1 elements as well as stress-associated cis-elements. Moreover, MusaSNAC1 also regulate multiple stress-related genes by binding to core site of NAC-proteins CGT[A/G] in their 5'-upstream region. Results indicated an interesting mechanism of drought tolerance through stomatal closure by H₂O₂ generation in guard cells, regulated by a NAC-protein in banana.

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

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

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.

GID1 expression is associated with ovule development of sexual and apomictic plants

BbrizGID1 is expressed in the nucellus of apomictic Brachiaria brizantha, previous to aposporous initial differentiation. AtGID1a overexpression triggers differentiation of Arabidopsis thaliana MMC-like cells, suggesting its involvement in ovule development. GIBBERELLIN-INSENSITIVE DWARF1 (GID1) is a gibberellin receptor previously identified in plants and associated with reproductive development, including ovule formation. In this work, we characterized the Brachiaria brizantha GID1 gene (BbrizGID1). BbrizGID1 showed up to 92% similarity to GID1-like gibberellin receptors of other plants of the Poaceae family and around 58% to GID1-like gibberellin receptors of Arabidopsis thaliana. BbrizGID1 was more expressed in ovaries at megasporogenesis than in ovaries at megagametogenesis of both sexual and apomictic plants. In ovules, BbrizGID1 transcripts were detected in the megaspore mother cell (MMC) of sexual and apomictic B. brizantha. Only in the apomictic plants, expression was also observed in the surrounding nucellar cells, a region in which aposporous initial cells differentiate to form the aposporic embryo sac. AtGID1a ectopic expression in Arabidopsis determines the formation of MMC-like cells in the nucellus, close to the MMC, that did not own MMC identity. Our results suggest that GID1 might be involved in the proper differentiation of a single MMC during ovule development and provide valuable information on the role of GID1 in sexual and apomictic reproduction.

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

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

Evaluation of reference genes for reverse transcription quantitative real-time PCR (RT-qPCR) studies in Silene vulgaris considering the method of cDNA preparation

Accurate gene expression measurements are essential in studies of both crop and wild plants. Reverse transcription quantitative real-time PCR (RT-qPCR) has become a preferred tool for gene expression estimation. A selection of suitable reference genes for the normalization of transcript levels is an essential prerequisite of accurate RT-qPCR results. We evaluated the expression stability of eight candidate reference genes across roots, leaves, flower buds and pollen of Silene vulgaris (bladder campion), a model plant for the study of gynodioecy. As random priming of cDNA is recommended for the study of organellar transcripts and poly(A) selection is indicated for nuclear transcripts, we estimated gene expression with both random-primed and oligo(dT)-primed cDNA. Accordingly, we determined reference genes that perform well with oligo(dT)- and random-primed cDNA, making it possible to estimate levels of nucleus-derived transcripts in the same cDNA samples as used for organellar transcripts, a key benefit in studies of cyto-nuclear interactions. Gene expression variance was estimated by RefFinder, which integrates four different analytical tools. The SvACT and SvGAPDH genes were the most stable candidates across various organs of S. vulgaris, regardless of whether pollen was included or not.

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

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

Banana NAC transcription factor MusaNAC042 is positively associated with drought and salinity tolerance

Banana is an important fruit crop and its yield is hampered by multiple abiotic stress conditions encountered during its growth. The NAC (NAM, ATAF, and CUC) transcription factors are involved in plant response to biotic and abiotic stresses. In the present study, we studied the induction of banana NAC042 transcription factor in drought and high salinity conditions and its overexpression in transgenic banana to improve drought and salinity tolerance. MusaNAC042 expression was positively associated with stress conditions like salinity and drought and it encoded a nuclear localized protein. Transgenic lines of banana cultivar Rasthali overexpressing MusaNAC042 were generated by Agrobacterium-mediated transformation of banana embryogenic cells and T-DNA insertion was confirmed by PCR and Southern blot analysis. Our results using leaf disc assay indicated that transgenic banana lines were able to tolerate drought and high salinity stress better than the control plants and retained higher level of total chlorophyll and lower level of MDA content (malondialdehyde). Transgenic lines analyzed for salinity (250 mM NaCl) and drought (Soil gravimetric water content 0.15) tolerance showed higher proline content, better Fv/Fm ratio, and lower levels of MDA content than control suggesting that MusaNAC042 may be involved in responses to higher salinity and drought stresses in banana. Expression of several abiotic stress-related genes like those coding for CBF/DREB, LEA, and WRKY factors was altered in transgenic lines indicating that MusaNAC042 is an efficient modulator of abiotic stress response in banana.

Validation of Reference Genes for Quantitative Real-Time PCR during Bicolor Tepal Development in Asiatic Hybrid Lilies (Lilium spp.)

Quantitative real-time PCR (qRT-PCR) is a reliable and high-throughput technique for gene expression studies, but its accuracy depends on the expression stability of reference genes. To date, several reliable reference gene identifications have been reported in Lilium spp., but none has been obtained for lily tepals at different developmental stages. In this study, ten candidate reference genes were selected and evaluated for their expression stability in Lilium 'Tiny Padhye' during the process of bicolor tepal development. The expression stability of these candidates was evaluated by three software programs (geNorm, NormFinder, and BestKeeper) and the comparative ΔCt method, and comprehensive stability rankings were generated by RefFinder. As a result, TIP41-like family gene (TIP41) and actin (ACT) were the best combination of reference genes for tepals at different developmental stages; TIP41 and F-box family gene (F-box) for tepals under shading treatment; ACT, actin11 (ACT11), and elongation factor 1-α (EF1-α) for different tissues; and ACT, TIP41, and ACT11 for all samples. The selected optimal reference genes were further verified by analyzing the expression levels of flavonoid 3'-hydroxylase (LhF3'H) and anthocyanidin 3-O-glucosyltransfersae (LhUFGT) in tepals at different developmental stages. This study provides useful information for gene expression characterization in lilies under different experimental conditions, and can serve as a basis for similar research in other closely related species.

Validation of Suitable Reference Genes for RT-qPCR Data in Achyranthes bidentata Blume under Different Experimental Conditions

Real-time quantitative polymerase chain reaction (RT-qPCR) is a sensitive technique for gene expression studies. However, choosing the appropriate reference gene is essential to obtain reliable results for RT-qPCR assays. In the present work, the expression of eight candidate reference genes, EF1-α (elongation factor 1-α), GAPDH (glyceraldehyde 3-phosphate dehydrogenase), UBC (ubiquitin-conjugating enzyme), UBQ (polyubiquitin), ACT (actin), β-TUB (β-tubulin), APT1 (adenine phosphoribosyltransferase 1), and 18S rRNA (18S ribosomal RNA), was evaluated in Achyranthes bidentata samples using two algorithms, geNorm and NormFinder. The samples were classified into groups according to developmental stages, various tissues, stresses (cold, heat, drought, NaCl), and hormone treatments (MeJA, IBA, SA). Suitable combination of reference genes for RT-qPCR normalization should be applied according to different experimental conditions. In this study, EF1-α, UBC, and ACT genes were verified as the suitable reference genes across all tested samples. To validate the suitability of the reference genes, we evaluated the relative expression of CAS, which is a gene that may be involved in phytosterol synthesis. Our results provide the foundation for gene expression analysis in A. bidentata and other species of Amaranthaceae.

Candidate Reference Genes Selection and Application for RT-qPCR Analysis in Kenaf with Cytoplasmic Male Sterility Background

Cytoplasmic male sterility (CMS) is a maternally inherited trait that results in the production of dysfunctional pollen. Based on reliable reference gene-normalized real-time quantitative PCR (RT-qPCR) data, examining gene expression profile can provide valuable information on the molecular mechanism of kenaf CMS. However, studies have not been conducted regarding selection of reference genes for normalizing RT-qPCR data in the CMS and maintainer lines of kenaf crop. Therefore, we studied 10 candidate reference genes (ACT3, ELF1A, G6PD, PEPKR1, TUB, TUA, CYP, GAPDH, H3, and 18S) to assess their expression stability at three stages of pollen development in CMS line 722A and maintainer line 722B of kenaf. Five computational statistical approaches (GeNorm, NormFinder, ΔCt, BestKeeper, and RefFinder) were used to evaluate the expression stability levels of these genes. According to RefFinder and GeNorm, the combination of TUB, CYP, and PEPKR1 was identified as an internal control for the accurate normalization across all sample set, which was further confirmed by validating the expression of HcPDIL5-2a. Furthermore, the combination of TUB, CYP, and PEPKR1 was used to differentiate the expression pattern of five mitochondria F₁F₀-ATPase subunit genes (atp1, atp4, atp6, atp8, and atp9) by RT-qPCR during pollen development in CMS line 722A and maintainer line 722B. We found that atp1, atp6, and atp9 exhibited significantly different expression patterns during pollen development in line 722A compared with line 722B. This is the first systematic study of reference genes selection for CMS and will provide useful information for future research on the gene expressions and molecular mechanisms underlying CMS in kenaf.

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

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

Leaf transcriptome of two highly divergent genotypes of Urochloa humidicola (Poaceae), a tropical polyploid forage grass adapted to acidic soils and temporary flooding areas

BACKGROUND

Urochloa humidicola (Koronivia grass) is a polyploid (6x to 9x) species that is used as forage in the tropics. Facultative apospory apomixis is present in most of the genotypes of this species, although one individual has been described as sexual. Molecular studies have been restricted to molecular marker approaches for genetic diversity estimations and linkage map construction. The objectives of the present study were to describe and compare the leaf transcriptome of two important genotypes that are highly divergent in terms of their phenotypes and reproduction modes: the sexual BH031 and the aposporous apomictic cultivar BRS Tupi.

RESULTS

We sequenced the leaf transcriptome of Koronivia grass using an Illumina GAIIx system, which produced 13.09 Gb of data that consisted of 163,575,526 paired-end reads between the two libraries. We de novo-assembled 76,196 transcripts with an average length of 1,152 bp and filtered 35,093 non-redundant unigenes. A similarity search against the non-redundant National Center of Biotechnology Information (NCBI) protein database returned 65 % hits. We annotated 24,133 unigenes in the Phytozome database and 14,082 unigenes in the UniProtKB/Swiss-Prot database, assigned 108,334 gene ontology terms to 17,255 unigenes and identified 5,324 unigenes in 327 known metabolic pathways. Comparisons with other grasses via a reciprocal BLAST search revealed a larger number of orthologous genes for the Panicum species. The unigenes were involved in C4 photosynthesis, lignocellulose biosynthesis and flooding stress responses. A search for functional molecular markers revealed 4,489 microsatellites and 560,298 single nucleotide polymorphisms (SNPs). A quantitative real-time PCR analysis validated the RNA-seq expression analysis and allowed for the identification of transcriptomic differences between the two evaluated genotypes. Moreover, 192 unannotated sequences were classified as containing complete open reading frames, suggesting that the new, potentially exclusive genes should be further investigated.

CONCLUSION

The present study represents the first whole-transcriptome sequencing of U. humidicola leaves, providing an important public information source of transcripts and functional molecular markers. The qPCR analysis indicated that the expression of certain transcripts confirmed the differential expression observed in silico, which demonstrated that RNA-seq is useful for identifying differentially expressed and unique genes. These results corroborate the findings from previous studies and suggest a hybrid origin for BH031.

Screening and Validation of Housekeeping Genes of the Root and Cotyledon of Cunninghamia lanceolata under Abiotic Stresses by Using Quantitative Real-Time PCR

Cunninghamia lanceolata (Chinese fir) is a fast-growing and commercially important conifer of the Cupressaceae family. Due to the unavailability of complete genome sequences and relatively poor genetic background information of the Chinese fir, it is necessary to identify and analyze the expression levels of suitable housekeeping genes (HKGs) as internal reference for precise analysis. Based on the results of database analysis and transcriptome sequencing, we have chosen five candidate HKGs (Actin, GAPDH, EF1a, 18S rRNA, and UBQ) with conservative sequences in the Chinese fir and related species for quantitative analysis. The expression levels of these HKGs in roots and cotyledons under five different abiotic stresses in different time intervals were measured by qRT-PCR. The data were statistically analyzed using the following algorithms: NormFinder, BestKeeper, and geNorm. Finally, RankAggreg was applied to merge the sequences generated from three programs and rank these according to consensus sequences. The expression levels of these HKGs showed variable stabilities under different abiotic stresses. Among these, Actin was the most stable internal control in root, and GAPDH was the most stable housekeeping gene in cotyledon. We have also described an experimental procedure for selecting HKGs based on the de novo sequencing database of other non-model plants.

Validation of Reference Genes for Gene Expression by Quantitative Real-Time RT-PCR in Stem Segments Spanning Primary to Secondary Growth in Populus tomentosa

The vertical segments of Populus stems are an ideal experimental system for analyzing the gene expression patterns involved in primary and secondary growth during wood formation. Suitable internal control genes are indispensable to quantitative real time PCR (qRT-PCR) assays of gene expression. In this study, the expression stability of eight candidate reference genes was evaluated in a series of vertical stem segments of Populus tomentosa. Analysis through software packages geNorm, NormFinder and BestKeeper showed that genes ribosomal protein (RP) and tubulin beta (TUBB) were the most unstable across the developmental stages of P. tomentosa stems, and the combination of the three reference genes, eukaryotic translation initiation factor 5A (eIF5A), Actin (ACT6) and elongation factor 1-beta (EF1-beta) can provide accurate and reliable normalization of qRT-PCR analysis for target gene expression in stem segments undergoing primary and secondary growth in P. tomentosa. These results provide crucial information for transcriptional analysis in the P. tomentosa stem, which may help to improve the quality of gene expression data in these vertical stem segments, which constitute an excellent plant system for the study of wood formation.

Selection and Verification of Candidate Reference Genes for Mature MicroRNA Expression by Quantitative RT-PCR in the Tea Plant (Camellia sinensis)

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is a rapid and sensitive method for analyzing microRNA (miRNA) expression. However, accurate qRT-PCR results depend on the selection of reliable reference genes as internal positive controls. To date, few studies have identified reliable reference genes for differential expression analysis of miRNAs among tissues, and among experimental conditions in plants. In this study, three miRNAs and four non-coding small RNAs (ncRNA) were selected as reference candidates, and the stability of their expression was evaluated among different tissues and under different experimental conditions in the tea plant (Camellia sinensis) using the geNorm and NormFinder programs. It was shown that miR159a was the best single reference gene in the bud to the fifth leaf, 5S rRNA was the most suitable gene in different organs, miR6149 was the most stable gene when the leaves were attacked by Ectropis oblique and U4, miR5368n and miR159a were the best genes when the leaves were treated by methyl jasmonate (MeJA), salicylic acid (SA) and abscisic acid (ABA), respectively. Our results provide suitable reference genes for future investigations on miRNA functions in tea plants.

Selection and Validation of Reference Genes for qRT-PCR in Cycas elongata

Quantitative reverse transcription PCR (qRT-PCR) is a sensitive technique used in gene expression studies. To achieve a reliable quantification of transcripts, appropriate reference genes are required for comparison of transcripts in different samples. However, few reference genes are available for non-model taxa, and to date, reliable reference genes in Cycas elongata have not been well characterized. In this study, 13 reference genes (ACT7, TUB, UBQ, EIF4, EF1, CLATHRIN1, PP2A, RPB2, GAPC2, TIP41, MAPK, SAMDC and CYP) were chosen from the transcriptome database of C. elongata, and these genes were evaluated in 8 different organ samples. Three software programs, NormFinder, GeNorm and BestKeeper, were used to validate the stability of the potential reference genes. Results obtained from these three programs suggested that CeGAPC2 and CeRPB2 are the most stable reference genes, while CeACT7 is the least stable one among the 13 tested genes. Further confirmation of the identified reference genes was established by the relative expression of AGAMOUSE gene of C. elongata (CeAG). While our stable reference genes generated consistent expression patterns in eight tissues, we note that our results indicate that an inappropriate reference gene might cause erroneous results. Our systematic analysis for stable reference genes of C. elongata facilitates further gene expression studies and functional analyses of this species.

Selection of Reference Genes for Real-Time Quantitative PCR in Pinus massoniana Post Nematode Inoculation

Pinus massoniaia Lamb has gained more and more attention as the most important tree species for timber and forestation in South China. Gene expression studies are of great importance to identify new and elite cultivars. Real-time quantitative PCR, a highly sensitive and specific method, is commonly used in the analysis of gene expression. The appropriate reference genes must be employed to normalize the calculation program for ascertaining repeatable and significant results. Herein, eleven housekeeping genes were evaluated during different stages of P. massoniana post nematode inoculation in this study. Three statistical approaches such as geNorm, NormFinder and BestKeeper were selected to analyze the stability of candidate genes. The results indicated that U2af and β-TUB were the most stable reference genes. These two genes could be used for the normalization in most of the experiments of P. massoniana, while Histone and AK were the least stable ones. In addition, EF expressed at the lowest average Ct value was the most abundant candidate gene. As an important gene associated with defense mechanisms, ABC transporter was analyzed by qRT-PCR, and the results were used to confirm the reliability of two genes. The selected reference genes in the present study will be conducive to future gene expression normalized by qRT-PCR in P. massoniana.

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.

Identification of a Novel Nematotoxic Protein by Challenging the Model Mushroom Coprinopsis cinerea with a Fungivorous Nematode

The dung of herbivores, the natural habitat of the model mushroom Coprinopsis cinerea, is a nutrient-rich but also very competitive environment for a saprophytic fungus. We showed previously that C. cinerea expresses constitutive, tissue-specific armories against antagonists such as animal predators and bacterial competitors. In order to dissect the inducible armories against such antagonists, we sequenced the poly(A)-positive transcriptome of C. cinerea vegetative mycelium upon challenge with fungivorous and bacterivorous nematodes, Gram-negative and Gram-positive bacteria and mechanical damage. As a response to the fungivorous nematode Aphelenchus avenae, C. cinerea was found to specifically induce the transcription of several genes encoding previously characterized nematotoxic lectins. In addition, a previously not characterized gene encoding a cytoplasmic protein with several predicted Ricin B-fold domains, was found to be strongly upregulated under this condition. Functional analysis of the recombinant protein revealed a high toxicity toward the bacterivorous nematode Caenorhabditis elegans. Challenge of the mycelium with A. avenae also lead to the induction of several genes encoding putative antibacterial proteins. Some of these genes were also induced upon challenge of the mycelium with the bacteria Escherichia coli and Bacillus subtilis. These results suggest that fungi have the ability to induce specific innate defense responses similar to plants and animals.

Identification and Evaluation of Reference Genes for Accurate Transcription Normalization in Safflower under Different Experimental Conditions

Safflower (Carthamus tinctorius L.) has received a significant amount of attention as a medicinal plant and oilseed crop. Gene expression studies provide a theoretical molecular biology foundation for improving new traits and developing new cultivars. Real-time quantitative PCR (RT-qPCR) has become a crucial approach for gene expression analysis. In addition, appropriate reference genes (RGs) are essential for accurate and rapid relative quantification analysis of gene expression. In this study, fifteen candidate RGs involved in multiple metabolic pathways of plants were finally selected and validated under different experimental treatments, at different seed development stages and in different cultivars and tissues for real-time PCR experiments. These genes were ABCS, 60SRPL10, RANBP1, UBCL, MFC, UBCE2, EIF5A, COA, EF1-β, EF1, GAPDH, ATPS, MBF1, GTPB and GST. The suitability evaluation was executed by the geNorm and NormFinder programs. Overall, EF1, UBCE2, EIF5A, ATPS and 60SRPL10 were the most stable genes, and MBF1, as well as MFC, were the most unstable genes by geNorm and NormFinder software in all experimental samples. To verify the validation of RGs selected by the two programs, the expression analysis of 7 CtFAD2 genes in safflower seeds at different developmental stages under cold stress was executed using different RGs in RT-qPCR experiments for normalization. The results showed similar expression patterns when the most stable RGs selected by geNorm or NormFinder software were used. However, the differences were detected using the most unstable reference genes. The most stable combination of genes selected in this study will help to achieve more accurate and reliable results in a wide variety of samples in safflower.

Selection and validation of reference genes for target gene analysis with quantitative RT-PCR in leaves and roots of bermudagrass under four different abiotic stresses

Quantitative real-time reverse transcriptase polymerase chain reaction (qRT-PCR) is an effective method for quantifying expression levels of target genes. The accuracy of qRT-PCR results is largely dependent on the selection of stable reference genes. The stability of reference gene expression may vary with plant species and environmental conditions. The objective of this study was to select stable reference genes for qRT-PCR analysis of target genes in different organs under different abiotic stresses for a perennial grass species, bermudagrass (Cynodon dactylon). The stability of eight potential reference genes (TUB, ACT, GAPDH, EF1α, TIP41, PP2A, CACS and UPL7) was evaluated under four different abiotic stresses (salt, drought, cold and heat) and in leaves and roots of bermudagrass. Four programs (geNorm, NormFinder, BestKeeper and RefFinder) were employed to evaluate the stability of reference gene expression and to identify the most stable reference genes for bermudagrass. Eight potential reference genes exhibited differential expression stability in leaves and roots under salt, drought, cold and heat stress. The expression levels of PP2A and CACS were stable in roots and leaves under salt stress, in leaves under drought stress and in roots exposed to cold and heat stress. EF1α and TIP41 expression was stable in roots of drought-stressed plants. UPL7, TUB and GAPDH were stably expressed in leaves under cold stress. Expression levels of PP2A and TIP41 were stable in leaves under heat stress. The use of the reference genes identified as internal controls for examination of gene expression patterns and quantification of expression levels of target genes will enable accurate qRT-PCR analysis in bermudagrass.

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.

Reference gene selection for quantitative gene expression studies during biological invasions: A test on multiple genes and tissues in a model ascidian Ciona savignyi

As invasive species have successfully colonized a wide range of dramatically different local environments, they offer a good opportunity to study interactions between species and rapidly changing environments. Gene expression represents one of the primary and crucial mechanisms for rapid adaptation to local environments. Here, we aim to select reference genes for quantitative gene expression analysis based on quantitative Real-Time PCR (qRT-PCR) for a model invasive ascidian, Ciona savignyi. We analyzed the stability of ten candidate reference genes in three tissues (siphon, pharynx and intestine) under two key environmental stresses (temperature and salinity) in the marine realm based on three programs (geNorm, NormFinder and delta Ct method). Our results demonstrated only minor difference for stability rankings among the three methods. The use of different single reference gene might influence the data interpretation, while multiple reference genes could minimize possible errors. Therefore, reference gene combinations were recommended for different tissues - the optimal reference gene combination for siphon was RPS15 and RPL17 under temperature stress, and RPL17, UBQ and TubA under salinity treatment; for pharynx, TubB, TubA and RPL17 were the most stable genes under temperature stress, while TubB, TubA and UBQ were the best under salinity stress; for intestine, UBQ, RPS15 and RPL17 were the most reliable reference genes under both treatments. Our results suggest that the necessity of selection and test of reference genes for different tissues under varying environmental stresses. The results obtained here are expected to reveal mechanisms of gene expression-mediated invasion success using C. savignyi as a model species.

Evaluation of Appropriate Reference Genes for Reverse Transcription-Quantitative PCR Studies in Different Tissues of a Desert Poplar via Comparision of Different Algorithms

Despite the unshakable status of reverse transcription-quantitative PCR in gene expression analysis, it has certain disadvantages, including that the results are highly dependent on the reference genes selected for data normalization. Since inappropriate endogenous control genes will lead to inaccurate target gene expression profiles, the validation of suitable internal reference genes is essential. Given the increasing interest in functional genes and genomics of Populus euphratica, a desert poplar showing extraordinary adaptation to salt stress, we evaluated the expression stability of ten candidate reference genes in P. euphratica roots, stems, and leaves under salt stress conditions. We used five algorithms, namely, ΔC_t, NormFinder, geNorm, GrayNorm, and a rank aggregation method (RankAggreg) to identify suitable normalizers. To support the suitability of the identified reference genes and to compare the relative merits of these different algorithms, we analyzed and compared the relative expression levels of nine P. euphratica functional genes in different tissues. Our results indicate that a combination of multiple reference genes recommended by GrayNorm algorithm (e.g., a combination of Actin, EF1α, GAPDH, RP, UBQ in root) should be used instead of a single reference gene. These results are valuable for research of gene identification in different P. euphratica tissues.

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.

Validating Internal Control Genes for the Accurate Normalization of qPCR Expression Analysis of the Novel Model Plant Setaria viridis

Employing reference genes to normalize the data generated with quantitative PCR (qPCR) can increase the accuracy and reliability of this method. Previous results have shown that no single housekeeping gene can be universally applied to all experiments. Thus, the identification of a suitable reference gene represents a critical step of any qPCR analysis. Setaria viridis has recently been proposed as a model system for the study of Panicoid grasses, a crop family of major agronomic importance. Therefore, this paper aims to identify suitable S. viridis reference genes that can enhance the analysis of gene expression in this novel model plant. The first aim of this study was the identification of a suitable RNA extraction method that could retrieve a high quality and yield of RNA. After this, two distinct algorithms were used to assess the gene expression of fifteen different candidate genes in eighteen different samples, which were divided into two major datasets, the developmental and the leaf gradient. The best-ranked pair of reference genes from the developmental dataset included genes that encoded a phosphoglucomutase and a folylpolyglutamate synthase; genes that encoded a cullin and the same phosphoglucomutase as above were the most stable genes in the leaf gradient dataset. Additionally, the expression pattern of two target genes, a SvAP3/PI MADS-box transcription factor and the carbon-fixation enzyme PEPC, were assessed to illustrate the reliability of the chosen reference genes. This study has shown that novel reference genes may perform better than traditional housekeeping genes, a phenomenon which has been previously reported. These results illustrate the importance of carefully validating reference gene candidates for each experimental set before employing them as universal standards. Additionally, the robustness of the expression of the target genes may increase the utility of S. viridis as a model for Panicoid grasses.