Identification and validation of quantitative real-time reverse transcription PCR reference genes for gene expression analysis in teak (Tectona grandis L.f.)

Selection and validation of appropriate reference genes for RT-qPCR analysis of Nitraria sibirica under various abiotic stresses

BACKGROUND

Nitraria sibirica Pall. is a halophytic shrub with strong environmental adaptability that can survive in extremely saline-alkali and drought-impacted environments. Gene expression analysis aids in the exploration of the molecular mechanisms of plant responses to abiotic stresses. RT-qPCR is the most common technique for studying gene expression. Stable reference genes are a prerequisite for obtaining accurate target gene expression results in RT-qPCR analysis.

RESULTS

In this study, a total of 10 candidate reference genes were selected from the transcriptome of N. sibirica, and their expression stability in leaves and roots under different treatment conditions (salt, alkali, drought, cold, heat and ABA) was evaluated with the geNorm, NormFinder, BestKeeper, comparative ΔCt and RefFinder programs. The results showed that the expression stability of the candidate reference genes was dependent on the tissue and experimental conditions tested. ACT7 combined with R3H, GAPDH, TUB or His were the most stable reference genes in the salt- or alkali-treated leaves, salt-treated roots and drought-treated roots, respectively; R3H and GAPDH were the most suitable combination for drought-treated leaves, heat-treated root samples and ABA-treated leaves; DIM1 and His maintained stable expression in roots under alkali stress; and TUB combined with R3H was stable in ABA-treated roots. TBCB and GAPDH exhibited stable expression in heat-treated leaves; TBCB, R3H, and ERF3A were stable in cold-treated leaves; and the three most stable reference genes for cold-treated roots were TBCB, ACT11 and DIM1. The reliability of the selected reference genes was further confirmed by evaluating the expression patterns of the NsP5CS gene under the six treatment conditions.

CONCLUSION

This study provides a theoretical reference for N. sibirica gene expression standardization and quantification under various abiotic stress conditions and will help to reveal the molecular mechanisms that confer stress tolerance to N. sibirica.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Overexpression of the Tectona grandis TgNAC01 regulates growth, leaf senescence and confer salt stress tolerance in transgenic tobacco plants

NAC transcription factors play critical roles in xylem secondary development and in regulation of stress response in plants. NAC proteins related to secondary cell wall development were recently identified and characterized in Tectona grandis (teak), one of the hardwood trees of highest economic importance in the world. In this work, we characterized the novel TgNAC01 gene, which is involved in signaling pathways that mediate teak response to stress. Abscisic acid (ABA) increases TgNAC01 expression in teak plants. Therefore, this gene may have a role in signaling events that mediate ABA-dependent osmotic stress responsive in this plant species. Stable expression in tobacco plants showed that the TgNAC01 protein is localized in the cell nucleus. Overexpression of TgNAC01 in two out three independent transgenic tobacco lines resulted in increased growth, leaf senescence and salt tolerance compared to wild type (WT) plants. Moreover, the stress tolerance of transgenic plants was affected by levels of TgNAC01 gene expression. Water potential, gas exchange and chlorophyll fluorescence were used to determine salt stress tolerance. The 35S:TgNAC01-6 line under 300 mM NaCl stress responded with a significant increase in photosynthesis rate, stomatal conductance, transpiration and carboxylation efficiency, but lower water potential compared to WT plants. The data indicate that the TgNAC01 transcription factor acts as a transcriptional activator of the ABA-mediated regulation and induces leaf senescence.

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

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

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

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

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.

Analysis of NAC Domain Transcription Factor Genes of Tectona grandis L.f. Involved in Secondary Cell Wall Deposition

NAC proteins are one of the largest families of plant-specific transcription factors (TFs). They regulate diverse complex biological processes, including secondary xylem differentiation and wood formation. Recent genomic and transcriptomic studies of Tectona grandis L.f. (teak), one of the most valuable hardwood trees in the world, have allowed identification and analysis of developmental genes. In the present work, T. grandis NAC genes were identified and analyzed regarding to their evolution and expression profile during wood formation. We analyzed the recently published T. grandis genome, and identified 130 NAC proteins that are coded by 107 gene loci. These proteins were classified into 23 clades of the NAC family, together with Populus, Eucalyptus, and Arabidopsis. Data on transcript expression revealed specific temporal and spatial expression patterns for the majority of teak NAC genes. RT-PCR indicated expression of VND genes (Tg11g04450-VND2 and Tg15g08390-VND4) related to secondary cell wall formation in xylem vessels of 16-year-old juvenile trees. Our findings open a way to further understanding of NAC transcription factor genes in T. grandis wood biosynthesis, while they are potentially useful for future studies aiming to improve biomass and wood quality using biotechnological approaches.

Physiological and molecular responses to drought stress in teak (Tectona grandis L.f.)

Drought stress is an increasingly common and worrying phenomenon because it causes a loss of production in both agriculture and forestry. Teak is a tropical tree which needs alternating rainy and dry seasons to produce high-quality wood. However, a robust understanding about the physiological characteristics and genes related to drought stress in this species is lacking. Consequently, after applying moderate and severe drought stress to teak seedlings, an infrared gas analyzer (IRGA) was used to measure different parameters in the leaves. Additionally, using the root transcriptome allowed finding and analyzing the expression of several drought-related genes. As a result, in both water deficit treatments a reduction in photosynthesis, transpiration, stomatal conductance and leaf relative water content was found. As well, an increase in free proline levels and intrinsic water use efficiency was found when compared to the control treatment. Furthermore, 977 transcripts from the root contigs showed functional annotation related to drought stress, and of these, TgTPS1, TgDREB1, TgAREB1 and TgPIP1 were selected. The expression analysis of those genes along with TgHSP1, TgHSP2, TgHSP3 and TgBI (other stress-related genes) showed that with moderate treatment, TgTPS1, TgDREB1, TgAREB1, TgPIP1, TgHSP3 and TgBI genes had higher expression than the control treatment, but with severe treatment only TgTPS1 and TgDREB1 showed higher expression than the control treatment. At the end, a schematic model for the physiological and molecular strategies under drought stress in teak from this study is provided. In conclusion, these physiological and biochemical adjustments in leaves and genetic changes in roots under severe and prolonged water shortage situations can be a limiting factor for teak plantlets' growth. Further studies of those genes under different biotic and abiotic stress treatments are needed.

Genetically Modified Heat Shock Protein90s and Polyamine Oxidases in Arabidopsis Reveal Their Interaction under Heat Stress Affecting Polyamine Acetylation, Oxidation and Homeostasis of Reactive Oxygen Species

One Sentence Summary

Heat shock proteins90 (HSP90s) induce acetylation of polyamines (PAs) and interact with polyamine oxidases (PAOs) affecting oxidation of PAs and hydrogen peroxide (H₂O₂) homeostasis in Arabidopsis thaliana.

Abstract

The chaperones, heat shock proteins (HSPs), stabilize proteins to minimize proteotoxic stress, especially during heat stress (HS) and polyamine (PA) oxidases (PAOs) participate in the modulation of the cellular homeostasis of PAs and reactive oxygen species (ROS). An interesting interaction of HSP90s and PAOs was revealed in Arabidopsis thaliana by using the pLFY:HSP90RNAi line against the four AtHSP90 genes encoding cytosolic proteins, the T-DNA Athsp90-1 and Athsp90-4 insertional mutants, the Atpao3 mutant and pharmacological inhibitors of HSP90s and PAOs. Silencing of all cytosolic HSP90 genes resulted in several-fold higher levels of soluble spermidine (S-Spd), acetylated Spd (N⁸-acetyl-Spd) and acetylated spermine (N¹-acetyl-Spm) in the transgenic Arabidopsis thaliana leaves. Heat shock induced increase of soluble-PAs (S-PAs) and soluble hydrolyzed-PAs (SH-PAs), especially of SH-Spm, and more importantly of acetylated Spd and Spm. The silencing of HSP90 genes or pharmacological inhibition of the HSP90 proteins by the specific inhibitor radicicol, under HS stimulatory conditions, resulted in a further increase of PA titers, N⁸-acetyl-Spd and N¹-acetyl-Spm, and also stimulated the expression of PAO genes. The increased PA titers and PAO enzymatic activity resulted in a profound increase of PAO-derived hydrogen peroxide (H₂O₂) levels, which was terminated by the addition of the PAO-specific inhibitor guazatine. Interestingly, the loss-of-function Atpao3 mutant exhibited increased mRNA levels of selected AtHSP90 genes. Taken together, the results herein reveal a novel function of HSP90 and suggest that HSP90s and PAOs cross-talk to orchestrate PA acetylation, oxidation, and PA/H₂O₂ homeostasis.

Expression profiling of CTR1-like and EIN2-like genes in buds and leaves of Populus tremula, and in vitro study of the interaction between their polypeptides

In Arabidopsis, the serine/threonine protein kinase Constitutive Triple Response 1 (CTR1) and Ethylene Insensitive 2 polypeptide (EIN2) functions are key negative and positive components, respectively, in the ethylene signalling route. Here, we report on an in silico study of members of the CTR1-like and EIN2-like polypeptide families from poplars. The expression of CTR1-like and EIN2-like genes such as Ptre-CTR1, Ptre-CTR3 and Ptre-EIN2a was studied in Populus tremula buds and leaves in response to dehydration, various light conditions and under senescence. In buds under dehydration, the maximal fold-change of the Ptre-CTR1, Ptre-CTR3 and Ptre-EIN2a expression level recorded almost identical values. This suggests that maintenance of a constant ratio between the transcript levels of genes encoding positive and negative ethylene signalling components is required under stress. The expression of the studied genes was 1.4-to 3-fold higher in response to darkness, but 4.5- to 51.2-fold and 21.6- to 51.2-fold higher under the early and moderate leaf senescence, respectively. It is worth noting that the senescence-related Ptre-EIN2a and Ptre-CTR3a expression profiles were very similar. Using in vitro assays, we demonstrated the ability of the catalytic domain of Ptre-CTR1 to phosphorylate the Ptre-EIN2a-like polypeptide, which is similar to that in Arabidopsis. The target substrate, the Ptre-CEND2a polypeptide (C-terminal part of Ptre-EIN2a), was only phosphorylated by the protein kinase Ptre-CTR1 and not by Ptre-CTR3. Moreover, the addition of Ptre-CTR3 polypeptides (-CTR3a or -CTR3b forms) to the reaction mixture had an inhibitory effect on Ptre-CTR1 auto- and trans-phosphorylation. In contrast to Ptre-CTR1, Ptre-CTR3 may act as a positive regulator in ethylene signalling in poplar; however, this hypothesis requires in vivo confirmation. Thus, the ethylene signalling route in poplar seems to be under the control of certain additional mechanisms which have not been reported in Arabidopsis.

Selection and evaluation of housekeeping genes as endogenous controls for quantification of mRNA transcripts in Theileria parva using quantitative real-time polymerase chain reaction (qPCR)

The reliability of any quantitative real-time polymerase chain reaction (qPCR) experiment can be seriously compromised by variations between samples as well as between PCR runs. This usually result from errors in sample quantification, especially with samples that are obtained from different individuals and tissues and have been collected at various time intervals. Errors also arise from differences in qPCR efficiency between assays performed simultaneously to target multiple genes on the same plate. Consequently, the derived quantitative data for the target genes become distorted. To avoid this grievous error, an endogenous control, with relatively constant transcription levels in the target individual or tissue, is included in the qPCR assay to normalize target gene expression levels in the analysis. Several housekeeping genes (HKGs) have been used as endogenous controls in quantification studies of mRNA transcripts; however, there is no record in the literature of the evaluation of these genes for the tick-borne protozoan parasite, Theileria parva. Importantly, the expression of these genes should be invariable between different T. parva stocks, ideally under different experimental conditions, to gain extensive application in gene expression studies of this parasite. Thus, the expression of several widely used HKGs was evaluated in this study, including the genes encoding β-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 28S rRNA, cytochrome b and fructose-2.6-biphosphate aldolase (F6P) proteins. The qPCR analysis revealed that the expression of genes encoding cytochrome b, F6P and GAPDH varied considerably between the two T. parva stocks investigated, the cattle-derived T. parva Muguga and the buffalo-derived T. parva 7014. 28S rRNA and β-actin gene expression was the most stable; thus, these genes were considered suitable candidates to be used as endogenous control genes for mRNA quantification studies in T. parva.

Selection and validation of reference genes for quantitative gene expression analyses in black locust (Robinia pseudoacacia L.) using real-time quantitative PCR

Black locust (Robinia pseudoacacia L.) is an easy to raise, fast growing, medium-sized deciduous tree species highly tolerant to harsh eco-conditions, i.e., drought and harsh winters, and it is widely adaptable to sandy, loamy, and marshy soils. The basis for this adaptability remains to be investigated at the transcriptomic level using real-time quantitative PCR (qPCR). Selection of a reliable gene for the normalization of qPCR data is important for obtaining accurate results in gene expression. The goal of this study was to identify an appropriate reference gene from 12 candidate genes for gene expression analysis in black locust exposed to various stressors such as abscisic acid (ABA), NaCl, polyethylene glycol (PEG) and varying temperatures. In GeNorm and NormFinder analyses, ACT (actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase) gene expression were the most stable in all conditions except heat stress, but in BestKeeper analysis, GAPDH and helicase gene expression were the most stable under NaCl and heat stress. In contrast, ACT and GAPDH were highest under abscisic acid (ABA), GAPDH and βTUB (beta tubulin) under cold stress, and helicase and EF1α (elongation factor 1 alpha) under PEG stress. We found that the most stable reference gene combination for all conditions was ACT and GAPDH. Additionally, the expression pattern of NAC2 (a transcription factor) and BGL2 in different tissues and under different stress conditions was analyzed relative to ACT and GAPDH and UBQ (ubiquitin) the least stably expressed gene. NAC2 and BGL2 both had highest expression in flowers and pods under ABA stress at 48h. This study provides useful reference genes for future gene expression studies in black locust.

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.

Análisis in silico y expresión génica del factor de transcripción TgNAC01 implicado en xilogénesis y estrés abiótico en Tectona grandis

El xilema secundario es el componente más abundante de la biomasa vegetal. Por tanto, conocer los genes que regulan su formación ayudaría a diseñar estrategias para el mejoramiento genético de la madera. Así, el objetivo de este trabajo fue realizar el análisis computacional de la estructura primaria y secundaria del factor de transcripción (FT) TgNAC01 de Tectona grandis, además de evaluar su historia evolutiva, dominios conservados y expresión génica en tejidos lignificados de árboles de 12 y 60 años. Para ello, se realizó una evaluación del potencial de interacción ion-electrón (PIIE), mediante el método del espectro de la información (MEI) utilizando la librería SFAPS de R-Project, seguido del modelamiento estructural utilizando el software MODELLER y visualizado mediante PyMol. Además, el análisis de alineamiento de secuencia múltiple y filogenia fue mediante el software Bioedit y MrBayes respectivamente. También se evaluó los niveles de síntesis del FT TgNAC01 mediante qRT-PCR. Como resultados, se evidencio que el FT mantiene una estructura β-hoja antiparalela retorcida, que se compacta contra una α-hélice en la región N-terminal, teniendo así tres dominios α hélice y siete dominios β plegada. Asimismo, mediante el MEI se demostró que tiene alrededor de cinco funciones biológicas y mutaciones sobre los aminoácidos con mayor PIIE, lo que conlleva a evoluciones sobre las redes de regulación genética. Finalmente, el FT TgNAC01 juega un papel fundamental en la organización y desarrollo de las partes que componen la albura, como las células radiales de la zona cambial, los vasos, fibras y los anillos de crecimiento.

Evaluation of Candidate Reference Genes for Quantitative Gene Expression Analysis in Spodoptera exigu a after Long-time Exposure to Cadmium

Studies on the transcriptional control of gene expression play an important role in many areas of biology. Reference genes, which are often referred to as housekeeping genes, such as GAPDH, G3PDH, EF2, RpL7A, RpL10, TUBα and Actin, have traditionally been assumed to be stably expressed in all conditions, and they are frequently used to normalize mRNA levels between different samples in qPCR analysis. However, it is known that the expression of these genes is influenced by numerous factors, such as experimental conditions. The difference in gene expression underlies a range of biological processes, including development, reproduction and behavior. The aim of this study was to show the problems associated with using reference genes in the qPCR technique, in a study on inbred strains of Spodoptera exigua selected toward cadmium resistance. We present and discuss our results and observations, and give some recommendations concerning the use and limitations of housekeeping genes as internal standards, especially in research on insects. Our results suggest that holometabolism and poikilothermia, as well as time since metamorphosis and the level of exposure to the selective factor (cadmium in this case), have a significant effect on the expression of reference genes.

Selection of reliable reference genes for gene expression studies in Trichoderma afroharzianum LTR-2 under oxalic acid stress

An appropriate reference gene is required to get reliable results from gene expression analysis by quantitative real-time reverse transcription PCR (qRT-PCR). In order to identify stable and reliable reference genes in Trichoderma afroharzianum under oxalic acid (OA) stress, six commonly used housekeeping genes, i.e., elongation factor 1, ubiquitin, ubiquitin-conjugating enzyme, glyceraldehyde-3-phosphate dehydrogenase, α-tubulin, actin, from the effective biocontrol isolate T. afroharzianum strain LTR-2 were tested for their expression during growth in liquid culture amended with OA. Four in silico programs (comparative ΔCt, NormFinder, geNorm and BestKeeper) were used to evaluate the expression stabilities of six candidate reference genes. The elongation factor 1 gene EF-1 was identified as the most stably expressed reference gene, and was used as the normalizer to quantify the expression level of the oxalate decarboxylase coding gene OXDC in T. afroharzianum strain LTR-2 under OA stress. The result showed that the expression of OXDC was significantly up-regulated as expected. This study provides an effective method to quantify expression changes of target genes in T. afroharzianum under OA stress.

Identification of suitable reference genes for quantitative real-time PCR normalization in blotched snakehead Channa maculata

A systematic study was conducted to identify reliable reference genes for normalization of gene expression analysis in the blotched snakehead Channa maculata under normal physiological conditions. Firstly, the partial complementary (c)DNA of nine candidate reference genes (actb, tmem104, ube2l3, ef1α, churc1, tmem256, rpl13a, sep15 and g6pd) were cloned from C. maculata. The expression levels of these genes were then assessed in embryos of different developmental stages and various tissue types of adult fish using quantitative real-time (qrt-)PCR. RefFinder algorithm was used to evaluate the expression stability of these genes based on their cycle-threshold (C_t ) values in the qrt-PCR analysis. Results showed that there was no single best reference gene for all stages of embryos and adult tissues tested. Furthermore, it was found that, among the nine candidate genes tested, actb and tmem104 were the most stable reference genes across adult tissue types, while sep15 and tmem256 were the most stable ones across developmental stages of embryos. These stable reference genes are recommended for normalization of gene expression analysis in C. maculata.

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 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 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.

Identification of Novel Reference Genes Suitable for qRT-PCR Normalization with Respect to the Zebrafish Developmental Stage

Reference genes used in normalizing qRT-PCR data are critical for the accuracy of gene expression analysis. However, many traditional reference genes used in zebrafish early development are not appropriate because of their variable expression levels during embryogenesis. In the present study, we used our previous RNA-Seq dataset to identify novel reference genes suitable for gene expression analysis during zebrafish early developmental stages. We first selected 197 most stably expressed genes from an RNA-Seq dataset (29,291 genes in total), according to the ratio of their maximum to minimum RPKM values. Among the 197 genes, 4 genes with moderate expression levels and the least variation throughout 9 developmental stages were identified as candidate reference genes. Using four independent statistical algorithms (delta-CT, geNorm, BestKeeper and NormFinder), the stability of qRT-PCR expression of these candidates was then evaluated and compared to that of actb1 and actb2, two commonly used zebrafish reference genes. Stability rankings showed that two genes, namely mobk13 (mob4) and lsm12b, were more stable than actb1 and actb2 in most cases. To further test the suitability of mobk13 and lsm12b as novel reference genes, they were used to normalize three well-studied target genes. The results showed that mobk13 and lsm12b were more suitable than actb1 and actb2 with respect to zebrafish early development. We recommend mobk13 and lsm12b as new optimal reference genes for zebrafish qRT-PCR analysis during embryogenesis and early larval stages.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BACKGROUND

Currently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.

RESULTS

We present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.

CONCLUSIONS

Expression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

Selection of reliable reference genes for gene expression studies in Clonostachys rosea 67-1 under sclerotial induction

Reference genes are important to precisely quantify gene expression by real-time PCR. In order to identify stable and reliable expressed genes in mycoparasite Clonostachys rosea in different modes of nutrition, seven commonly used housekeeping genes, 18S rRNA, actin, β-tubulin, elongation factor 1, ubiquitin, ubiquitin-conjugating enzyme and glyceraldehyde-3-phosphate dehydrogenase, from the effective biocontrol isolate C. rosea 67-1 were tested for their expression under sclerotial induction and during vegetative growth on PDA medium. Analysis by three software programs showed that differences existed among the candidates. Elongation factor 1 was most stable; the M value in geNorm, SD value in Bestkeeper and stability value in Normfinder analysis were 0.405, 0.450 and 0.442, respectively, indicating that the gene elongation factor 1 could be used to normalize gene expression in C. rosea in both vegetative growth and parasitic process. By using elongation factor 1, the expression of a serine protease gene, sep, in different conditions was assessed, which was consistent with the transcriptomic data. This research provides an effective method to quantitate expression changes of target genes in C. rosea, and will assist in further investigation of parasitism-related genes of this fungus.

Transcriptome and expression profiling analysis link patterns of gene expression to antennal responses in Spodoptera litura

Background

The study of olfaction is key to understanding the interaction of insects with their environment and provides opportunities to develop novel tactics for control of pest species. Recent developments in transcriptomic approaches enable the molecular basis of olfaction to be studied even in species with limited genomic information. Here we use transcriptome and expression profiling analysis to characterize the antennal transcriptome of the noctuid moth and polyphagous pest Spodoptera litura.

Results

We identify 74 candidate genes involved in odor detection and recognition, encoding 26 ORs, 21 OBPs, 18 CSPs and 9 IRs. We examine their expression levels in both sexes and seek evidence for their function by relating their expression with levels of EAG response in male and female antennae to 58 host and non-host plant volatiles and sex pheromone components. The majority of olfactory genes showed sex-biased expression, usually male-biased in ORs. A link between OR gene expression and antennal responses to odors was evident, a third of the compounds tested evoking a sex-biased response, in every case also male-biased. Two candidate pheromone receptors, OR14 and OR23 were especially strongly expressed and male-biased and we suggest that these may respond to the two female sex pheromone components of S. litura, Z9E11-14:OAc and Z9E12-14:OAc, which evoked strongly male-biased EAG responses.

Conclusions

Our results provide the molecular basis for elucidating the olfactory profile of moths and the sexual divergence of their behavior and could enable the targeting of particular genes, and behaviors for pest management.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1375-x) contains supplementary material, which is available to authorized users.

Validation of reference genes for quantitative real-time PCR during latex regeneration in rubber tree

In rubber tree, latex regeneration is one of the decisive factors influencing the rubber yield, although its molecular regulation is not well known. Quantitative real-time PCR (qPCR) is a popular and powerful tool used to understand the molecular mechanisms of latex regeneration. However, the suitable reference genes required for qPCR are not available to investigate the expressions of target genes during latex regeneration. In this study, 20 candidate reference genes were selected and evaluated for their expression stability across the samples during the process of latex regeneration. All reference genes showed a relatively wide range of the threshold cycle values, and their stability was validated by four different algorithms (comparative delta Ct method, Bestkeeper, NormFinder and GeNorm). Three softwares (comparative delta Ct method, NormFinder and GeNorm) exported similar results that identify UBC4, ADF, UBC2a, eIF2 and ADF4 as the top five suitable references, and 18S as the least suitable one. The application of the screened references would improve accuracy and reliability of gene expression analysis in latex regeneration experiments.

Selection of appropriate reference genes for quantitative real-time PCR in Oxytropis ochrocephala Bunge using transcriptome datasets under abiotic stress treatments

BACKGROUND

Oxytropis ochrocephala Bunge, an indigenous locoweed species in China, poses great threats to livestock on grasslands. There is a need for further genetic study in the plants per se, for understanding the basis of its acclimation mechanism in various unfavorable environmental conditions and to implement effective control measures. Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) is the most commonly used method for gene expression analysis. To facilitate gene expression studies and obtain more accurate qRT-PCR data, normalization relative to stable reference genes is required. The aim of this study was to select the most stable reference genes for transcriptional analysis in O. ochrocephala.

RESULTS

We selected 12 candidate reference genes, 18S ribosomal RNA (18S RNA), actin2/7 (ACT7), β-actin (ACTB), actin101 (ACT101), actin11 (ACT11), β-tubulin (TUB), α-tubulin (TUA), glyceraldehyde-3-phosphate dehydrogenase-1 (GAPDH1), GAPDH2, metallothionein-like protein (MET), fructose-bisphosphate aldolase (FBA) and histone H3 (HIS), from the transcriptome datasets of O. ochrocephala and determined the suitability by analyzing their expression levels when exposed to a range of abiotic stress conditions. By employing software packages including geNorm, NormFinder and BestKeeper, HIS, ACT7, and ACT101 were assessed as the most suitable set for normalization in all samples. When normalized with the most stable reference genes, the expression patterns of the three target genes were in accordance with those in the transcriptome data, indicating that the reference genes selected in this study are suitable.

CONCLUSIONS

The study provided appropriate reference genes for accurate normalization in qRT-PCR analysis in O. ochrocephala and emphasized the importance of validating reference genes for gene expression analysis under specific experimental condition. The usage of inappropriate reference gene would cause misinterpretation.

Transcriptional changes induced by in vivo exposure to pentachlorophenol (PCP) in Chironomus riparius (Diptera) aquatic larvae

Pentachlorophenol (PCP) has been extensively used worldwide as a pesticide and biocide and is frequently detected in the aquatic environment. In the present work, the toxicity of PCP was investigated in Chironomus riparius aquatic larvae. The effects following short- and long-term exposures were evaluated at the molecular level by analyzing changes in the transcriptional profile of different endocrine genes, as well as in genes involved in the stress response and detoxification. Interestingly, although no differences were found after 12- and 24-h treatments, at 96-h exposures PCP was able to induce significant increases in transcripts from the ecdysone receptor gene (EcR), the early ecdysone-inducible E74 gene, the estrogen-related receptor gene (ERR), the Hsp70 gene and the CYP4G gene. In contrast, the Hsp27 gene appeared to be downregulated, while the ultraspiracle gene (usp) (insect ortholog of the retinoid X receptor) was not altered in any of the conditions assayed. Moreover, Glutathione-S-Transferase (GST) activity was not affected. The results obtained show the ability of PCP to modulate transcription of different biomarker genes from important cellular metabolic activities, which could be useful in genomic approaches to monitoring. In particular, the significant upregulation of hormonal genes represents the first evidence at the genomic level of the potential endocrine disruptive effects of PCP on aquatic invertebrates.