Reliable reference genes for normalization of gene expression data in tea plants (Camellia sinensis) exposed to metal stresses

Metal tolerance protein CsMTP4 has dual functions in maintaining zinc homeostasis in tea plant

Plants have evolved a series of zinc (Zn) homeostasis mechanisms to cope with the fluctuating Zn in the environment. How Zn is taken up, translocated and tolerate by tea plant remains unknown. In this study, on the basis of RNA-Sequencing, we isolated a plasma membrane-localized Metal Tolerance Protein (MTP) family member CsMTP4 from Zn-deficient tea plant roots and investigated its role in regulation of Zn homeostasis in tea plant. Heterologous expression of CsMTP4 specifically enhanced the tolerance of transgenic yeast to Zn excess. Moreover, overexpression of CsMTP4 in tea plant hairy roots stimulated Zn uptake under Zn deficiency. In addition, CsMTP4 promoted the growth of transgenic Arabidopsis plants by translocating Zn from roots to shoots under Zn deficiency and conferred the tolerance to Zn excess by enhancing the efflux of Zn from root cells. Transcriptome analysis of the CsMTP4 transgenic Arabidopsis found that the expression of Zn metabolism-related genes were differentially regulated compared with wild-type plants when exposed to Zn deficiency and excess conditions. This study provides a mechanistic understanding of Zn uptake and translocation in plants and a new strategy to improve phytoremediation efficiency.

Transcriptome Analysis Reveals Differentially Expressed Genes Involved in Aluminum, Copper and Cadmium Accumulation in Tea 'Qianmei 419' and 'Qianfu 4'

Tea, as a global nonalcoholic beverage, is widely consumed due to its economic, health and cultural importance. Polyploids have the ability to solve the problems of low yield, cold resistance and insect resistance in tea tree varieties. However, the response mechanism to aluminum and heavy metal remains unclear. In this study, the content of Al, Cu and Cd were measured in the leaves and roots of 'Qianmei 419' and 'Qianfu 4', respectively. The content of Al, Cd and Cu in the roots of the 'Qianmei 419' tea variety were significantly higher than in 'Qianfu 4' roots. Only the content of Cu in the leaves of the 'Qianmei 419' tea variety was significantly higher than that in the roots of the 'Qianfu 4' tea variety. Moreover, we found that the content of Al, Cu and Cd in the soil around the root of 'Qianfu 4' were higher than in the soil around the root of 'Qianmei 419'. RNA-seq was performed to identify the DEGs involved in the accumulation of Al, Cu and Cd between 'Qianmei 419' and 'Qianfu 4'. A total of 23,813 DEGs were identified in the triploid tea variety, including 16,459 upregulated DEGs and 7354 downregulated DEGs. Among them, by analyzing the expression levels of some metal transporter genes, it was found that most of the metal transporter genes were downregulated in the triploid tea plants. In short, through the analysis of transcriptome data and metal content, it was found that changes in metal transporter gene expression affect the accumulation of metals in tea plants. These results provide candidate genes to enhance multi-metal tolerance through genetic engineering technology.

Evaluation and validation of reference genes in Cymbidium faberi for real-time quantitative PCR

For all organs at all Cymbidium faberi stages, ACT, UBQ3 and GAPDH can be selected as reference genes. For organs of the vegetative stage, UBQ2 and UBQ3 can be chosen for analysis of normalized gene expression. For the bud stage, ACT and UBQ3 can be used for analysis of gene expression. For the full blossom stage, ACT, UBQ3 and UBQ2 can be introduced into relative gene expression analysis. For vegetative organs, UBQ2 and ACT can be used as reference genes. For reproductive organs, ACT, UBQ3 and UBQ2 can be used as a reference for data processing. CfAG1 gene expression is more consistent when UBQ3, GAPDH and ACT are used as reference genes.

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

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

Identification of suitable reference genes for expression studies in pomegranate under different biotic and abiotic stress conditions

Pomegranate (Punica granatum L.) is an important economic fruit crop, facing many biotic and abiotic challenges during cultivation. Several research programs are in progress to understand both biotic and abiotic stress factors and mitigate these challenges using gene expression studies based on the qPCR approach. However, research publications are not available yet to select the standard reference gene for normalizing target gene expression values in pomegranate. The most suitable candidate reference gene is required to ensure precise and reliable results for qPCR analysis. Eight candidate reference genes' stability was evaluated under different stress conditions using different algorithms such as ∆Ct, geNorm, BestKeeper, NormFinder, and RefFinder. The various algorithms revealed that EFA1 and 18S rRNA were common and most stable reference genes (RGs) under abiotic and wilt stress. Whereas comprehensive ranking by RefFinder showed GAPDH and CYPF were the most stable RGs under combined biotic (pooled samples of all biotic stress) and bacterial blight samples. For normalizing target gene expression under wilt, nematode, bacterial blight, and abiotic stress conditions both GAPDH and CYPFreference genes are adequate for qPCR. The above data provide comprehensive details for the selection of a candidate reference gene in various stresses in pomegranate.

Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight

Gray blight (GB) is one of the most destructive diseases of tea plants, causing considerable damage and productivity losses; however, the dynamic roles of defense genes during pathogen infection remain largely unclear. To explore the numerous molecular interactions associated with GB stress in tea plants, we employed transcriptome, sRNAome and degradome sequencing from 1 to 13 days post-inoculation (dpi) at 3-day intervals. The transcriptomics results showed that differentially expressed genes (DEGs) related to flavonoid synthesis, such as chalcone synthase (CHS) and phenylalanine ammonia-lyase (PAL), were particularly induced at 4 dpi. Consistent with this, the contents of catechins (especially gallocatechin), which are the dominant flavonoids in tea plants, also increased in the leaves of tea plants infected with GB. Combined analysis of the sRNAome and degradome revealed that microRNAs could mediate tea plant immunity by regulating DEG expression at the post-transcriptional level. Co-expression network analysis demonstrated that miR530b-ethylene responsive factor 96 (ERF96) and miRn211-thaumatin-like protein (TLP) play crucial roles in the response to GB. Accordingly, gene-specific antisense oligonucleotide assays suggested that suppressing ERF96 decreased the levels of reactive oxygen species (ROS), whereas suppressing TLP increased the levels of ROS. Furthermore, ERF96 was induced, but TLP was suppressed, in susceptible tea cultivars. Our results collectively demonstrate that ERF96 is a negative regulator and TLP is a positive regulator in the response of tea plants to GB. Taken together, our comprehensive integrated analysis reveals a dynamic regulatory network linked to GB stress in tea plants and provides candidate genes for improvement of tea plants.

Indole primes defence signalling and increases herbivore resistance in tea plants

Upon herbivore attack, plants emit herbivore-induced plant volatiles (HIPVs). HIPVs can prime defences and resistance of intact plants. However, how HIPVs are decoded and translated into functional defence responses is not well understood, especially in long-lived woody plants. Here, we investigated the impact of the aromatic HIPV indole on defence-related early signalling, phytohormone accumulation, secondary metabolite biosynthesis and herbivore resistance in tea plants. We find that tea plants infested with tea geometrid caterpillars release indole at concentrations >450 ng*hr^-1 . Exposure to corresponding doses of synthetic indole primes the expression of early defence genes involved in calcium (Ca²⁺ ) signalling, MPK signalling and jasmonate biosynthesis. Indole exposure also primes the production of jasmonates and defence-related secondary metabolites. These changes are associated with higher herbivore resistance of indole-exposed tea plants. Chemical inhibition of Ca²⁺ and jasmonate signalling provides evidence that both are required for indole-mediated defence priming and herbivore resistance. Our systematic assessment of the impact of indole on defence signalling and deployment shows that indole acts by boosting Ca²⁺ signalling, resulting in enhanced jasmonate-dependent defence and resistance in a woody plant. Our work extends the molecular basis of HIPV-induced defence priming from annual plants to an economically important tree species.

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

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

Regulation of antioxidant enzymes and osmo-protectant molecules by salt and drought responsive genes in Bambusa balcooa

Bio-energy crops need to be grown on marginal salt and drought lands in India as per policy. Understanding environmental stress response in bio-energy crops might help in promoting cultivation of bio-energy feedstock on marginal salty and drought land. This is one of the first report for vegetative propagation of Bamboo (Bambusa balcooa) under salt and drought stress to understand antioxidant enzymes' gene regulations to combat stress through activation of antioxidant enzymes and osmo-protectant molecules to scavenge reactive oxygen species as measured by physiological changes. Morphological, physiological, and biochemical traits were noted as indicators of plant health upon different sodium chloride (NaCl) salt-stress while various drought conditions with correlation analysis. A significant up-regulation of genes related to most of the antioxidant enzymes was observed up to salinity of 14 mS cm^- 1 electric conductivity (EC) at 150 mM NaCl experimental salt stress which declined with higher salt-stress. While in the case of drought-stress, all genes remained up-regulated while proline dehydrogenase (PDH) remained down-regulated up-to 100% drought-stress having 4% soil moisture. The gene expressions of antioxidant enzymes were significantly correlated with their corresponding gene-products namely super-oxide dismutase (SOD), catalase (CAT), glutathione reductase (GR) and ascorbate peroxidase (APX) activities. Biochemical parameters such as, soluble sugar, proline, malondialdehyde (MDA), total amino acids, hydrogen peroxide and electrolyte leakage ratio also showed positive correlation (p = 0.001) with salt condition. Genetic and biochemical test parameters were significantly correlated with physiological attributes of plant health at soil EC of 14 mS cm^- 1 shown as 150 mM NaCl salt stress and 60% drought-stress having 17% soil moisture content, were the optimum stress tolerance limits observed. Application of these data would be useful to cultivate 0.63 million ha of salinity affected land and 10.05 million ha of drought affected land among wastelands in India to meet biofuel need.

Nanoparticles Affect the Expression Stability of Housekeeping Genes in Plant Cells

Purpose

We report on the expression stability of several housekeeping/reference genes that can be used in the normalization of target gene expression in quantitative real-time PCR (qRT-PCR) analysis of plant cells challenged with metal nanoparticles (NPs).

Materials and Methods

Uniform cell suspension cultures of Hypericum perforatum were treated with 25 mg/l silver and gold NPs (14-15 nm in diameter). Cells were collected after 0.5, 4.0, and 12 h. The total RNA isolated from the cells was analyzed for the stability of ACT2, ACT3, ACT7, EF1-α, GAPDH, H2A, TUB-α, TUB-β, and 18S rRNA genes using qRT-PCR. The cycle threshold (Ct) values of the genes were analyzed using the geNorm, NormFinder, BestKeeper, and RefFinder statistical algorithms to rank gene stability. The stability of the top-ranked genes was validated by normalizing the expression of HYP1.

Results

The expression of the tested housekeeping genes varied with treatment duration and NP types. EF1-α in gold NP treatment and TUB-α and EF1-α in silver NP treatment ranked among the top three positions. However, none of the genes retained their top ranking with time and across NP types.

Conclusion

EF1-α can be used as a reference for treatment involving both silver and gold NPs in H. perforatum cells. TUB-α can be used only for silver NP-treated cells. The expression instability of most of the housekeeping genes highlights the importance of systematic standardization of reference genes for NP treatment conditions to draw proper conclusions on the target gene expression.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance

BACKGROUND

Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.).

RESULTS

In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca²⁺. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response.

CONCLUSIONS

A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

JA-Ile-macrolactone 5b Induces Tea Plant (Camellia sinensis) Resistance to Both Herbivore Ectropis obliqua and Pathogen Colletotrichum camelliae

Jasmonates (JAs), the group of lipid-derived hormones, were found to control the defense responses in a myriad of plants. Meaningfully, the macrolactones of 12-hydroxy jasmonate isoleucine (12OH-JA-Ile) were reported to induce the defensive response of wild tobacco. However, little to nothing has been known about the elicitation effect of JA-Ile-macrolactones on woody plants to harmful organisms, let alone its underlying mechanisms. Here, we first optimized the synthetic routine using mild toxic reagent isobutyl chloroformate instead of ethyl chloroformate for conjugation, and we used acetonitrile (MeCN) instead of ethyl alcohol for the better dissolution of p-toluenesulfonic acid (p-TsOH) to gain JA-Ile-macrolactones. JA-Ile-macrolactone 5b-treated tea plants significantly inhibited the larvae weight gain of Ectropis obliqua larvae and the lesions caused by Colletotrichum camelliae. Furthermore, the expression level of CsOPR3 was significantly upregulated in 5b-treated leaves. Meanwhile, 5b reduced the accumulation of eriodictyol 7-O-glucuronide (EDG) in tea plants, which was confirmed to promote the growth rate of E. obliqua larvae by artificial diet assay. In conclusion, our study proved that the exogenous application of 5b could induce the tea plant resistance both to herbivore E. obliqua and pathogen C. camelliae, and EDG was identified as one of the secondary metabolites that could influence the growth rate of E. obliqua, but it did not directly influence the infection of C. camelliae in vitro. Further research should be carried out to clarify the mechanism through which 5b induces tea plant resistance to C. camelliae.

Identification and validation of superior housekeeping gene(s) for qRT-PCR data normalization in Agave sisalana (a CAM-plant) under abiotic stresses

The adaptive mechanisms in Agave species enable them to survive and exhibit remarkable tolerance to abiotic stresses. Quantitative real-time PCR is a highly reliable approach for validation of targeted differential gene expression. However, stable housekeeping gene(s) is prerequisite for accurate normalization of expression data by qRT-PCR. Till date, no systematic validation study for candidate housekeeping gene identification or evaluation has been carried-out in Agave species. A total of 17 candidate housekeeping genes were identified from the de novo assembled transcriptomic data of A. sisalana and rigorously analyzed for expression stability assessment under drought, heat, cold and NaCl stress. Different statistical algorithms like geNorm, BestKeeper, NormFinder, and RefFinder on expression data determined the superior housekeeping gene(s) for accurate normalization of the gene of interest (GOI). The comprehensive evaluation revealed the β-Tub 4, WIN-1 and CYC-A as the most stable, while EEF1α, GAPDH, and UBE2 were ranked as the least stable genes in pooled samples. Pairwise combination by geNorm showed that up to two housekeeping genes would be adequate to normalize the GOI expression data precisely. Validation of identified most and least stable housekeeping genes was carried-out by normalizing the expression data of AsHSP20 under abiotic stress conditions. Copy number of AsHSP20 gene supports the reliability of the genes used for normalization. This is the first report on the screening and validation of the housekeeping genes under abiotic stress condition in A. sisalana that would assist to understand the stress tolerance mechanisms by novel gene identification and accurate validation.

Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plants under differential biotic stresses

The selection of reliable reference genes (RGs) for normalization under given experimental conditions is necessary to develop an accurate qRT-PCR assay. To the best of our knowledge, only a small number of RGs have been rigorously identified and used in tea plants (Camellia sinensis (L.) O. Kuntze) under abiotic stresses, but no critical RG identification has been performed for tea plants under any biotic stresses till now. In the present study, we measured the mRNA transcriptional levels of ten candidate RGs under five experimental conditions; these genes have been identified as stable RGs in tea plants. By using the ΔCt method, geNorm, NormFinder and BestKeeper, CLATHRIN1 and UBC1, TUA1 and SAND1, or SAND1 and UBC1 were identified as the best combination for normalizing diurnal gene expression in leaves, stems and roots individually; CLATHRIN1 and GAPDH1 were identified as the best combination for jasmonic acid treatment; ACTIN1 and UBC1 were identified as the best combination for Toxoptera aurantii-infested leaves; UBC1 and GAPDH1 were identified as the best combination for Empoasca onukii-infested leaves; and SAND1 and TBP1 were identified as the best combination for Ectropis obliqua regurgitant-treated leaves. Furthermore, our results suggest that if the processing time of the treatment was long, the best RGs for normalization should be recommended according to the stability of the proposed RGs in different time intervals when intragroup differences were compared, which would strongly increase the accuracy and sensitivity of target gene expression in tea plants under biotic stresses. However, when the differences of intergroup were compared, the RGs for normalization should keep consistent across different time points. The results of this study provide a technical guidance for further study of the molecular mechanisms of tea plants under different biotic stresses.

miR477 targets the phenylalanine ammonia-lyase gene and enhances the susceptibility of the tea plant (Camellia sinensis) to disease during Pseudopestalotiopsis species infection

MAIN CONCLUSION: miR477 acts as a negative regulator in tea plant immunity against Pseudopestalotiopsis infection by repressing the expression of its target gene PAL. MicroRNA (miRNA)-mediated post-transcriptional regulation plays a fundamental role in various plant physiological processes, including responses to pathogens. Our previous research revealed that miR477 might be involved in the tea plant-Pseudopestalotiopsis interaction (data not shown). In the present study, the accumulation of miR477 significantly decreased in tea plants during Pseudopestalotiopsis species infection. Using miRNA and degradome data sets, the targeting of phenylalanine ammonia-lyase (PAL) by miR477 was validated by 5' RLM-RACE. GUS assay showed that the expression of PAL was post-transcriptionally regulated by miR477 and silenced by mRNA cleavage. A negative correlation between the expression of miR477 and PAL was found in tea plants infected by the pathogen. The transgenic lines overexpressing Csn-miR477 exhibited increased susceptibility to Pseudopestalotiopsis species, which was associated with reduced expression of PAL during infection. The degree of severity of the leaf lesions and the results of trypan blue staining showed that the plants overexpressing Csn-miR477 exhibited more severe damage upon pathogen infection than wild-type plants. In addition, more H2O2 and O2-, higher malondialdehyde (MDA) contents and less superoxide dismutase (SOD) and peroxidase (POD) activities were detected in the transgenic plants than in the wild-type plants after inoculation with Pseudopestalotiopsis species. Taken together, our results implied that Csn-miR477 might act as a negative regulator in pathogen-infected tea plants by inhibiting the expression of its target, PAL, and that Csn-miR477 is a candidate miRNA for improving the adaptation of tea plant to disease.

Validation of Reliable Reference Genes for RT-qPCR Studies of Target Gene Expression in Colletotrichum camelliae During Spore Germination and Mycelial Growth and Interaction With Host Plants

The tea plant [Camellia sinensis (L.) O. Kuntze] is one of the most important leaf crops, and it is widely used for the production of non-alcoholic beverages worldwide. Tea also has a long history of medicinal use. Colletotrichum camelliae Massee is one of the dominant fungal pathogens that infects tea leaves and causes severe tea anthracnose disease. To analyze the molecular biology of C. camelliae, the quantification of pathogen gene expression by the RT-qPCR method is necessary. Reliable RT-qPCR results require the use of stable reference genes for data normalization. However, suitable reference genes have not been reported in C. camelliae thus far. In this study, 12 candidate genes (i.e., CcSPAC6B12.04c, CcWDR83, Cchp11, Ccnew1, CcHplo, CcRNF5, CcHpcob, CcfaeB-2, CcYER010C, CcRNM1, CcUP18, and CcACT) were isolated from C. camelliae and assessed as potential reference genes. The expression stability of these genes in C. camelliae during spore germination and mycelial growth and interaction with host plants was first evaluated using several statistical algorithms, such as geNorm, NormFinder, and Bestkeeper. A web-based analysis program, Refinder, was then used to find the most suitable reference genes. Our results indicated that Cenew1, CcHplo, and CcSPAC6B12.04c were the most stable reference genes in C. camelliae under all conditions. Our work provided the most suitable reference genes for future studies performed to quantify the target gene expression levels of C. camelliae.

Selection of reliable reference genes for gene expression analysis in seeds at different developmental stages and across various tissues in Paeonia ostii

Paeonia ostii seeds have recently been identified as a new source of α-linolenic acid in China. Studying the gene expression patterns of unsaturated fatty acid-related genes would be helpful for understanding the mechanism of α-linolenic acid accumulation. Quantitative real-time polymerase chain reaction (qRT-PCR) is a useful method for reliably evaluating gene expression, and it is necessary to select reliable reference genes for data normalization in qRT-PCR analysis. In this study, we evaluated the expression stability of 12 candidate reference genes using four mathematical algorithms (∆C_t, BestKeeper, NormFinder, and geNorm). The web-based tool RefFinder was used to integrate the results and to provide a comprehensive ranking order. The expression stability ranking orders of reference genes were different caculated by these four algorithms, and the ranking order analyzed by the RefFinder was UBQ > Tip41 > UCE > EF-1α > α-TUB > PP2A > ACT > GAPDH > SAM > CYP > β-TUB > 18S at the different seed development stages, and UBQ > Tip41 > EF-1α > α-TUB > PP2A > UCE > GAPDH > SAM > ACT > CYP > 18S > β-TUB in P. ostii tissues. UBQ and Tip41 are the two most stable whereas 18S and β-TUB are the two least stable reference genes for gene expression in various tissues and seeds at different developmental stages in P. ostii.

Validation of reference genes for gene expression studies in post-harvest leaves of tea plant (Camellia sinensis)

Tea is one of three major non-alcoholic beverages that are popular all around the world. The economic value of tea product largely depends on the post-harvest physiology of tea leaves. The utilization of quantitative reverse transcription polymerase chain reaction is a widely accepted and precise approach to determine the target gene expression of tea plants, and the reliability of results hinges on the selection of suitable reference genes. A few reliable reference genes have been documented using various treatments and different tissues of tea plants, but none has been done on post-harvest leaves during the tea manufacturing process. The present study selected and analyzed 15 candidate reference genes: Cs18SrRNA, CsGADPH, CsACT, CsEF-1α, CsUbi, CsTUA, Cs26SrRNA, CsRuBP, CsCYP, CselF-4α, CsMON1, CsPCS1, CsSAND, CsPPA2, CsTBP. This study made an assessment on the expression stability under two kinds of post-harvest treatment, turn over and withering, using three algorithms-GeNorm, Normfinder, and Bestkeeper. The results indicated that the three commonly used reference genes, CsTUA, Cs18SrRNA, CsRuBP, together with Cs26SrRNA, were the most unstable genes in both the turn over and withering treatments. CsACT, CsEF-1α, CsPPA2, and CsTBP were the top four reference genes in the turn over treatment, while CsTBP, CsPCS1, CsPPA2, CselF-4α, and CsACT were the five best reference genes in the withering group. The expression level of lipoxygenase genes, which were involved in a number of diverse aspects of plant physiology, including wounding, was evaluated to validate the findings. To conclude, we found a basis for the selection of reference genes for accurate transcription normalization in post-harvest leaves of tea plants.

Selection and validation of reference genes for target gene analysis with quantitative real-time PCR in the leaves and roots of Carex rigescens under abiotic stress

Carex rigescens is an ornamental turfgrass in northern China which has a relatively low maintenance cost and robust tolerance to many adverse environmental conditions, so it could be considered a new material for researching into plant stress resistance. However, suitable reference genes are vacant for obtaining reliable results in quantitative real-time PCR (qRT-PCR) analysis of C. rigescens in adversity research. In this study, we tested the expression stability of nine potential reference genes in leaves and roots under five different abiotic stress conditions, including cold, salt, heat, osmotic and cadmium (Cd). We then selected the best reference genes according to the analysis results calculated by three algorithmic programs (geNorm, NormFinder and BestKeeper) and used the RankAggreg package to merge the outputted data. The results showed that combinations of at least two reference genes should be used for reliable normalization except in heat-treated root samples, which require three reference genes. eIF-4α, GADPH, SAND and PEPKR1 and their combination were found to be the most stably expressed reference genes, while SAM, TUA-α and UPL7 were the three least stable reference genes among most of experimental samples. In addition, five stress-induced genes (Cu-Zn SOD, P5CS, LEA, GST, and APX) were chosen to verify the stability of the selected reference genes in various tissues and under various stress conditions. The results of this study will provide an important fundamental basis both for gene expression verification for transcriptomic and proteomic analyses and for gene expression analysis for future gene function research in C. rigescens.

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

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

Selection of reference genes for tissue/organ samples on day 3 fifth-instar larvae in silkworm, Bombyx mori

The silkworm, Bombyx mori, is one of the world's most economically important insect. Surveying variations in gene expression among multiple tissue/organ samples will provide clues for gene function assignments and will be helpful for identifying genes related to economic traits or specific cellular processes. To ensure their accuracy, commonly used gene expression quantification methods require a set of stable reference genes for data normalization. In this study, 24 candidate reference genes were assessed in 10 tissue/organ samples of day 3 fifth-instar B. mori larvae using geNorm and NormFinder. The results revealed that, using the combination of the expression of BGIBMGA003186 and BGIBMGA008209 was the optimum choice for normalizing the expression data of the B. mori tissue/organ samples. The most stable gene, BGIBMGA003186, is recommended if just one reference gene is used. Moreover, the commonly used reference gene encoding cytoplasmic actin was the least appropriate reference gene of the samples investigated. The reliability of the selected reference genes was further confirmed by evaluating the expression profiles of two cathepsin genes. Our results may be useful for future studies involving the quantification of relative gene expression levels of different tissue/organ samples in B. mori.

Reference gene selection for qRT-PCR assays in Stellera chamaejasme subjected to abiotic stresses and hormone treatments based on transcriptome datasets

Background

Stellera chamaejasme Linn, an important poisonous plant of the China grassland, is toxic to humans and livestock. The rapid expansion of S. chamaejasme has greatly damaged the grassland ecology and, consequently, seriously endangered the development of animal husbandry. To draft efficient prevention and control measures, it has become more urgent to carry out research on its adaptive and expansion mechanisms in different unfavorable habitats at the genetic level. Quantitative real-time polymerase chain reaction (qRT-PCR) is a widely used technique for studying gene expression at the transcript level; however, qRT-PCR requires reference genes (RGs) as endogenous controls for data normalization and only through appropriate RG selection and qRT-PCR can we guarantee the reliability and robustness of expression studies and RNA-seq data analysis. Unfortunately, little research on the selection of RGs for gene expression data normalization in S. chamaejasme has been reported.

Method

In this study, 10 candidate RGs namely, 18S, 60S, CYP, GAPCP1, GAPDH2, EF1B, MDH, SAND, TUA1, and TUA6, were singled out from the transcriptome database of S. chamaejasme, and their expression stability under three abiotic stresses (drought, cold, and salt) and three hormone treatments (abscisic acid, ABA; gibberellin, GA; ethephon, ETH) were estimated with the programs geNorm, NormFinder, and BestKeeper.

Result

Our results showed that GAPCP1 and EF1B were the best combination for the three abiotic stresses, whereas TUA6 and SAND, TUA1 and CYP, GAPDH2 and 60S were the best choices for ABA, GA, and ETH treatment, respectively. Moreover, GAPCP1 and 60S were assessed to be the best combination for all samples, and 18S was the least stable RG for use as an internal control in all of the experimental subsets. The expression patterns of two target genes (P5CS2 and GI) further verified that the RGs that we selected were suitable for gene expression normalization.

Discussion

This work is the first attempt to comprehensively estimate the stability of RGs in S. chamaejasme. Our results provide suitable RGs for high-precision normalization in qRT-PCR analysis, thereby making it more convenient to analyze gene expression under these experimental conditions.

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.