Validation of Appropriate Reference Genes for Real-Time Quantitative PCR Gene Expression Analysis in Rice Plants Exposed to Metal Stresses

Environmental pollution by toxic heavy metals may lead to the possible contamination of the rice plant (Oryza sativa L.). Although gene expression analysis through real-time quantitative PCR (RT-qPCR) has increased our knowledge about biological responses to heavy metals, gene network that mediates rice plant responses to heavy metal stress remains elusive. In such scenario, validation of reference gene is a major requirement for successful analyzes involving RT-qPCR. In this study, we analyzed the expression stability of eight commonly used housekeeping genes (GAPDH, Actin, eIF-4α, UBQ 5, UBQ 10, UBC, EF-1α and β-TUB) in rice leaves exposed to four kinds of heavy metals (Zn, Cu, Cd and Pb). The expression stability of these genes was determined using geNorm, NormFinder, BestKeeper and RefFinder algorithms. The results showed that UBQ 10 and UBC were the most stable reference genes across all the tested samples. We measured the expression profiles of the heavy metal-inducible gene O. sativa METALLOTHIONEIN2b (OsMT2b) using the two most stable and one least stable reference genes in all samples. The relative expression of OsMT2b varied greatly according to the different reference genes. Our results may be beneficial for future studies involving the quantification of relative gene expression levels in rice plants.     

Genome-Wide Analysis of the ZRT,IRT-Like Protein (ZIP) Family and Their Responses to Metal Stress in <Emphasis Type="Italic">Populus trichocarpa</Emphasis>

The ZRT, IRT-like protein (ZIP) family plays an important role in the transport of zinc (Zn) and iron (Fe) across the cell membrane in many different species. However, studies on ZIP family are mainly limited in herbaceous species; hence, we investigated functional divergence of ZIP family in Populus trichocarpa. We identified 21 ZIP genes in P. trichocarpa and classified them into four groups based on phylogenetic analysis. Structural analyses revealed that most of the PtrZIP transporters have eight transmembrane domains (TMDs). PtrZIP members were unequally positioned in 19 P. trichocarpa linkage groups (LGs), with six tandem duplications and four segmental duplications. The promoter regions of PtrZIP genes contain Zn, Fe, copper (Cu), and other metal stress-related cis-elements. Additionally, tissue-specific expression of PtrZIP genes showed that most of them had relatively high expression levels in the root. Quantitative real-time RT-PCR (qRT-PCR) analysis revealed that the expression of PtrZIP genes were induced not only under deficiency or excess condition of Zn, Fe, Cu and manganese (Mn) but also under excess condition of cadmium (Cd) and lead (Pb) stress. These findings indicated that PtrZIP genes may have played potential roles in metal transporters. Genome-wide analysis of PtrZIP genes in P. trichocarpa provided more comprehensive insights on the structure and function of this gene family.     

Reference genes selection for quantitative gene expression studies in <Emphasis Type="Italic">Pinus massoniana</Emphasis> L.

Key message

Evaluation and selection of reference genes in Pinus massoniana L. (PM) for gene expression studies of various tissues, floral organ development, and abiotic stress.

Abstract

An important prerequisite for obtaining accurate gene expression results using quantitative real-time PCR is the selection of a reference gene or a group of genes having a highly stable level of expression. Pinus massoniana L. (PM) is the predominant fast-growing timber forest tree species in southern China. In this study of PM, we evaluated various tissues, flowers in different developmental phases, leaves from a cultivar with insect resistance, and leaves from plants under several types of abiotic stresses. Comprehensive Analysis was performed using BestKeeper, Normfinder, geNorm, and RefFinder software to select the most stable reference gene or gene group from among 25 candidate genes in these samples. The results showed that different experimental conditions require the use of different reference genes: ACT1 could be used as a reference gene for all samples in this study; UBI4 was the best gene for various tissues and zinc stress; CYP was the most stable gene for leaves from insect-resistant materials and Pb stress; Fbox and UBI11 were the best reference genes for salt stress; Fbox + RRP8, ARF + TUBA, and EF1B + IDH were the best reference groups for drought stress, low temperature stress, and flowers in different developmental phases, respectively. This study presents a reliable selection of reference genes for Masson pine, and the conclusions are meaningful for improving the accuracy of expression analyses in future molecular biology studies.

     

A panel of the most suitable reference genes for RT-qPCR expression studies of coffee: screening their stability under different conditions

The reliability of analyses using real-time quantitative polymerase chain reaction (RT-qPCR) depends on the selection of appropriate reference genes to correct for sample-to-sample and run-to-run variations. The aim of the present study was to select the most suitable reference genes for gene expression analyses in tissue samples from coffee, Coffea arabica L. (Arabica) grown under well-watered (WW) and water-deficit (WD) conditions and C. canephora Pierre ex A. Froehner (Robusta) grown under WW conditions. Expression profiles and stabilities were evaluated for 12 reference genes in different tissues from C. arabica and for 8 genes in tissues from C. canephora. The web-based RefFinder tool, which combines the geNorm, NormFinder, Bestkeeper, and Delta-Ct algorithms, was employed to assess the stability of the tested genes. The most stable reference genes identified for all tissues grouped (WW/WD) of C. arabica were clathrin adaptor protein medium subunit (AP47), ubiquitin (UBQ), 60S ribosomal protein L39 (RPL39), and elongation factor 1α (EF1α), while class III alcohol dehydrogenase (ADH2), β-actin (ACT), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and ubiquitin (UBQ) genes were the most stable for all tissues grouped (WW) of C. canephora tissues. Validation by the expression level analysis of CaACO-like demonstrated that the use of the best and the worst set of reference genes produced different expression results. The results reinforce the general assumption that there is no universal reference gene and that it is essential to select the most appropriate gene for each individual experiment to apply adequate normalization procedures of RT-qPCR data.     

The use of miRNAs as reference genes for miRNA expression normalization during <Emphasis Type="Italic">Lilium</Emphasis> somatic embryogenesis by real-time reverse transcription PCR analysis

Expression profiling of miRNAs has the ability to reveal the essence of somatic embryogenesis (SE). qRT-PCR is one of the most commonly used techniques for dynamic miRNA detection but requires optimal reference genes for data reliability. This is the first report on reference gene validation for miRNA expression normalization in Lilium (Lilium pumilum DC. Fisch. and Lilium davidii var. unicolor). In this study, seventeen miRNAs together with two snRNAs (U4, U6), one rRNA (5S rRNA) and three protein-coding genes (FP, ACT, GAPDH) were selected as reference candidates, and their expression stability was validated by qRT-PCR among eleven developing SE cultures in two lilies. Four normalization algorithms, including geNorm, BestKeeper, NormFinder and RefFinder, were also used to evaluate the stability of the reference candidates. For Lilium pumilum DC. Fisch., lpu-miR159a was the optimal reference gene during SE, followed by lpu-miR408b, while U6 was the least stable reference candidate. For Lilium davidii var. unicolor, FP presented greater stability than did half of the miRNA candidates, but the best reference gene was lda-miR162, followed by lda-miR159a. Further analysis of the expression level of miR156 and miR529 was used to evaluate the validity of the reference genes in both lilies. In general, miRNAs are superior to common protein-coding genes and snRNAs / rRNAs as reference genes for miRNA expression normalization during Lilium SE, and the most suitable reference miRNA is different between two species in the same Lilium genus. This is a pioneer study using suitable miRNAs as reference genes in Lilium and constitutes a small but essential step for the further exploration of miRNA function in Lilium, thus offering valuable references for other plants.     

Identification of genes showing differential expression profile associated with growth rate in skeletal muscle tissue of Landrace weanling pig

Suppression subtractive hybridization was used to identify genes showing differential expression profile associated with growth rate in skeletal muscle tissue of Landrace weanling pig. Two subtracted cDNA populations were generated from musculus longissimus muscle tissues of selected pigs with extreme expected breeding values at the age of 100 kg. Three upregulated genes (EEF1A2, TSG101 and TTN) and six downregulated genes (ATP5B, ATP5C1, COQ3, HADHA, MYH1 and MYH7) in pig with genetic propensity for higher growth rate were identified by sequence analysis of 12 differentially expressed clones selected by differential screening following the generation of the subtracted cDNA population. Real-time PCR analysis confirmed difference in expression profiles of the identified genes in musculus longissimus muscle tissues between the two Landrace weanling pig groups with divergent genetic propensity for growth rate. Further, differential expression of the identified genes except for the TTN was validated by Western blot analysis. Additionally, the eight genes other than the ATP5C1 co-localized with the same chromosomal positions as QTLs that have been previously identified for growth rate traits. Finally, the changes of expression predicted from gene function suggested association of upregulation of expression of the EEF1A2, TSG101 and TTN genes and downregulation of the ATP5B, ATP5C1, COQ3, HADHA, MYH1 and MYH7 gene expression with increased growth rate. The identified genes will provide an important insight in understanding of the molecular mechanism underlying growth rate in Landrace pig breed.     

Selection of reference genes for quantitative real-time PCR in <Emphasis Type="Italic">Casuarina equisetifolia</Emphasis> under salt stress

Real time quantitative PCR (qPCR) is widely used in gene expression analysis for its accuracy and sensitivity. Reference genes serving as endogenous controls are necessary for gene normalization. In order to select an appropriate reference gene to normalize gene expression in Casuarina equisetifolia under salt stress, 10 potential reference genes were evaluated using real time qPCR in the leaves and roots of plants grown under different NaCl concentrations and treatment durations. GeNorm, NormFinder, and BestKeeper analyses reveal that elongation factor 1-alpha (EF1α) and ubiquitin-conjugating enzyme E2 (UBC) were the most appropriate reference genes for real time qPCR under salt stress. However, β-tubulin (βTUB) and actin 7, which were widely used as reference genes in other plant species, were not always stably expressed. The combination of EF1α, UBC, uncharacterized protein 2, DNAJ homolog subfamily A member 2, and glyceraldehyde-3-phosphate dehydrogenase should be ideal reference genes for normalizing gene expression data in all samples under salt stress. It indicates the need for reference gene selection for normalizing gene expression in C. equisetifolia. In addition, the suitability of reference genes selected was confirmed by validating the expression of WRKY29-like and expansin-like B1. The results enable analysis of salt response mechanism and gene expression in C. equisetifolia.     

Validation of reference genes for gene expression analysis of response to anthocyanin induction in cell cultures of <Emphasis Type="Italic">Vitis davidii</Emphasis> (Rom. Caill.) Foëx

Real-time quantitative polymerase chain reaction (RT-qPCR) is an effective method for detecting changes of gene expression in plant cell metabolic regulation. A set of 15 reference gene candidates were selected for the present study of anthocyanin biosynthesis regulation, and stability. The suitability of their expression was evaluated in eight different experimental treatments in spine grape (Vitis davidii [Rom. Caill.] Foëx.) cell cultures. The results indicated that SAND family protein (SAND) and V-type proton ATPase subunit G (VAG) were the most stable reference genes for culture duration, tubulin alpha-3/alpha-5 chain (α-tubulin) and tubulin beta-1 chain (β-tubulin) for illumination conditions, ubiquitin-conjugating enzyme E2-17 kDa (UBQ) and VAG for UVB treatment, VAG and 60S ribosomal protein L18-2 (60SRP) for temperature treatment, AP47, clathrin adaptor complex subunit mu (AP-2) and 60SRP for cinnamic acid treatment, α-tubulin and UBQ for chitosan treatment, actin and alcohol dehydrogenase 2 (ADH2) for kinetin treatment, and β-tubulin and elongation factor 1-α (EF1-α) for cell line. Finally, the reliability of the selected reference genes was confirmed by investigating the expression profiles of the target gene dihydroflavonol 4-reductase (DFR) in spine grape cell cultures. The results of the present study offer the most robust platform for the most precise and broad application of RT-qPCR to investigate gene expression associated with anthocyanin biosynthesis in spine grape cell cultures.     

Plant reference genes for development and stress response studies

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

Validation of reference genes for RT-qPCR studies in <Emphasis Type="Italic">Stevia rebaudiana</Emphasis> in response to elicitor agents

Stevia rebaudiana is an important source of natural steviol glycosides and is of increasing interest in various fields of study. Therefore, understanding the molecular processes regulating its metabolism is of great importance. In this study, the stability of seven reference genes (18S ribosomal RNA, Actin, Aquaporin, Calmodulin, Eukaryote elongation factor 1-α, Malate dehydrogenase, and Ubiquitin) under the effect of three stress-related elicitors (methyl jasmonate, salicylic acid, and spermidine) was evaluated in stevia plants. We used RefFinder software, which makes use of the four main currently available algorithms for reference gene selection: geNorm, NormFinder, BestKeeper, and the Comparative ?Ct method. The results indicated that Ubiquitin and Actin can be used as reference genes under all tested experimental conditions. The genes, 18S ribosomal RNA, traditionally used as a reference gene, along with Calmodulin, showed the lowest stability. The expression of Deoxyxylulose-5-phosphate synthase and Kaurenoic acid hydroxylase genes was used to confirm the validated reference genes, showing that inadequacy of the reference gene may lead to erroneous results. This is the first study on the stability of reference genes in Stevia rebaudiana plants, and is of great relevance for further analysis of the gene expression of the steviol glycoside biosynthetic pathway.     

Glutathione reductase gene expression depends on chloroplast signals in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Glutathione reductase (EC 1.6.4.2) is one of the main antioxidant enzymes of the plant cell. In Arabidopsis thaliana, glutathione reductase is encoded by two genes: the gr1 gene encodes the cytosolic-peroxisomal form, and the gr2 gene encodes the chloroplast-mitochondrial form. Little is known about the regulation of expression of plant glutathione reductase genes. In the present work, we have demonstrated that gr2 (but not gr1) gene expression in Arabidopsis leaves changes depending on changes in redox state of the photosynthetic electron transport chain. Expression of both the gr1 and gr2 genes was induced by reactive oxygen species. In heterotrophic suspension cell culture of Arabidopsis, expression of both studied genes did not depend on H₂O₂ level or on changes in the redox state of the mitochondrial electron transport chain. Our data indicate that chloroplasts are involved in the regulation of the glutathione reductase gene expression in Arabidopsis.