Reference genes for the normalization of gene expression in eucalyptus species

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

Exogenous applications of brassinosteroids promote secondary xylem differentiation in Eucalyptus grandis

Brassinosteroids (BRs) play many pivotal roles in plant growth and development, especially in cell elongation and vascular development. Although its biosynthetic and signal transduction pathway have been well characterized in model plants, their biological roles in Eucalyptus grandis, a major hardwood tree providing fiber and energy worldwide, remain unclear. Here, we treated E. grandis plantlets with 24-epibrassinolide (EBL), the most active BR and/or BR biosynthesis inhibitor brassinazole. We recorded the plant growth and analyzed the cell structure of the root and stem with histochemical methods; then, we performed a secondary growth, BR synthesis, and signaling-related gene expression analysis. The results showed that the BRs dramatically increased the shoot length and diameter, and the exogenous BR increased the xylem area of the stem and root. In this process, EgrBRI1, EgrBZR1, and EgrBZR2 expression were induced by the BR treatment, and the expressions of HD-ZIPIII and cellulose synthase genes were also altered. To further verify the effect of BRs in secondary xylem development in Eucalyptus, we used six-month-old plants as the material and directly applied EBL to the xylem and cambium of the vertical stems. The xylem area, fiber cell length, and cell numbers showed considerable increases. Several key BR-signaling genes, secondary xylem development-related transcription factor genes, and cellulose and lignin biosynthetic genes were also considerably altered. Thus, BR had regulatory roles in secondary xylem development and differentiation via the BR-signaling pathway in this woody plant.

Development of reference genes for RT-qPCR analysis of gene expression in Pleurotus pulmonarius for biotechnological applications

Jatropha curcas is an oilseed crop with biorefinery applications. Whilst cake generated following oil extraction offers potential as a protein source for animal feed, inactivation of toxic phorbol esters present in the material is necessary. Pleurotus pulmonarius is a detoxifying agent for jatropha cake with additional potential as animal feed, edible mushroom and for enzyme production. For the characterization of fungal genes involved in phorbol ester degradation, together with other industrial applications, reverse transcription-quantitative PCR (RT-qPCR) is a tool that enables accurate quantification of gene expression. For this, reliable analysis requires reference genes for normalization of mRNA levels validated under conditions employed for target genes. The stability of potential reference genes β-TUB, ACTIN, GAPDH, PHOS, EF1α, TRPHO, LAC, MNP3, MYP and VP were evaluated following growth of P. pulmonarius on toxic, non-toxic jatropha cake and a combined treatment, respectively. NormFinder and geNorm algorithms for expression stability analysis identified PHOS, EF1α and MNP3 as appropriate for normalizing gene expression. Reference gene combinations contrasting in ranking were compared following normalization of relative expression of the CHU_2040 gene, encoding an esterase enzyme potentially involved in phorbol ester degradation. The reference genes for P. pulmonarius will facilitate the elucidation of mechanisms involved in detoxification of phorbol esters as well as analysis of target genes for application in biorefinery models.

Selection and verification of reliable internal reference genes in stem development of herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) has emerged in the cut flower market due to its beautiful appearance. The bending flower stems caused by a lack of mechanical strength is the main problem restricting the development of the cut P. lactiflora industry. So it is of great worth to reveal the basis of stem development changes in P. lactiflora to improve its cut flower quality. Quantitative research on gene expression characteristics can provide clues for understanding their biological functions, and the screening of relatively stable expression genes is a prerequisite for the quantitative study of gene expression characteristics. Thus, it is necessary to find appropriate genes during stem development so as to analyze the qRT‒PCR results. In this study, 10 genes were screened, and these genes expressed stably in stems of different stem strengths at three different developmental stages. Then, their expressions were evaluated by RefFinder, BestKeeper, NormFinder, and GeNorm programs. The results demonstrated that γ-tubulin (γ-TUB) was the most suitable gene, followed by α-tubulin (α-TUB) and β-D-glucosidase (β-GUS), whereas histone H3 (His) was the least suitable gene. Besides, the temporal and spatial expression characteristics of PlCOMT1, the key gene concerned with the synthesis of cell wall fillers in P. lactiflora, were also used to evaluate the suitability of genes. Consequently, γ-TUB and α-TUB are the two best combinations during stem development, and their combination can be used for the stem development of P. lactiflora. These findings will provide a reference for the selection of genes related to stem development and the study of molecular mechanisms related to stem development in P. lactiflora.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01325-5.

Expression Activity of Artificial Promoters for Disease Resistance in Transgenic Eucalyptus urophylla

The transcriptional properties of artificial promoters are closely related to the type and arrangement position of cis-elements. GWSF (374-bp) was an effective SPIP with four cis-element dimers. There were four pathogen-inducible cis-elements in the GWSF promoter (GST1-boxes, W-boxes, S-boxes, and F-boxes) and a minimal cauliflower mosaic virus 35S promoter. V-element dimers were inserted into the upstream (VGWSF), midstream (GWVSF), and downstream (GWSFV) regions of the original GWSF promoter sequence to examine their affect on the position. The expression activity of promoters was analyzed and estimated using the histochemical staining of leaf discs of eucalyptus with transient expression, an image digitization method to extract the color features, and the induction treatment by a plant pathogenic microorganism/inducer and qPCR assays. The histochemical staining results of the adventitious buds indicated that the promoters had been successfully integrated into the E. urophylla genome and that they drove the expression of the gus gene. There was a noticeable difference in the intensity of color between the adventitious buds on the same callus block, as well as the intensity of color within the same adventitious bud. According to the established two-factor model of blue value, there was a greater difference between the levels of the genotype factor than the promoter factor in eucalyptus leaf discs. Further, the basal and inducible transcriptional levels of the three improved promoters were investigated by qPCR. With the basal transcriptional level of the GWSF promoter normalized to one, the relative basal levels of VGWSF, GWVSF, and GWSFV were 1.40, 1.45, and 4.15, respectively. The qPCR results were consistent with the staining results of GUS histochemical staining. The three improved promoters all had the properties of being induced by salicylic acid, Ralstonia solanacearum, and Phytophthora capsici. The three improved promoters demonstrated a significantly higher TMV induction activity: their induction activity from high to low was GWSFV > GWVSF > VGWSF. The findings will be beneficial to the construction and optimization of artificial promoters for transgenic plants.

Identification and Validation of Reliable Reference Genes for Gene Expression Studies in Koelreuteria paniculata

RT-qPCR is considered a rapid and reliable technique for analyzing gene expression. This technique is commonly used to analyze the expression of various genes at diverse transcriptional levels in different samples. However, few studies have characterized ornamental Koelreuteria species for reliable reference genes. In this study, eight reference genes were evaluated as controls in RT-qPCR with SYBR green to quantify gene expression in different Koelreuteria paniculata samples. All selected reference genes showed a broad range of Ct values in all samples, which was supportive of their variable expression. Our results showed significant variation in the stable expression of K. paniculata genes. Sample data, analyzed using geNorm, NormFinder, and BestKeeper, showed that phospholipase (PLA2) and β-actin (ACT) were the most suitable and statistically reliable reference genes, whereas ribosomal protein L13 (RPL13) and elongation factor 1-α (EF1α) were less stable and unsuitable for use as internal controls. To compare gene expression levels, two or more reference genes should be used for data normalization. Thus, the stability and expression of both PLA2 and ACT were believed to provide better normalization and quantification of the transcript levels for gene expression studies in K. paniculata.

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

BACKGROUND

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

METHODS AND RESULTS

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

CONCLUSIONS

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

Screening of stable internal reference gene of Quinoa under hormone treatment and abiotic stress

Real-time quantitative polymerase chain reaction is the most commonly used method to accurately detect gene expression patterns. The method requires stable internal reference genes to standardize the data. However, studies have shown that there is no stable expression of internal reference genes in different tissues and under different treatments. Therefore, in order to study the optimal reference genes of quinoa under different hormones and abiotic stress, leaves and stems from quinoa seedlings treated with low temperature (4 °C), salt (200 mmol/L) and abscisic acid (200 mmol/L) were used as experimental materials. Using ACT-1, eIF, EF1α, GAPDH, TUA, TUB-9, TUB-1, H2A and L8-1 as candidate reference genes, the expression stability of these 9 quinoa candidate reference genes under different hormone treatment and abiotic stress was evaluated by using geNorm, NormFinder and BestKeeper software. The results showed that TUB-1 gene under salt stress, L8-1 gene under low temperature stress, EF-1α gene induced by ABA. PLIM2c WLIM1and WLIM2b were selected to verify the candidate internal reference genes, and finally the expression of GAPDH was most unstable under the three treatments, which was not suitable to be the internal reference gene of quinoa under specific conditions, while EF1α showed good stability under the three different treatments and was suitable to be used as the internal reference gene. In conclusion, the results of this study could provide an important reference for quantifying the expression level of reference genes in quinoa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01094-z.

Selection and validation of endogenous reference genes for RT-qPCR normalization in different stresses and tissues of the tiger milk mushroom, Pleurotus tuber-regium

Pleurotus tuber-regium is a sclerotium-producing mushroom with pharmacological and nutritional value. Many researches applied molecular biological strategies that allow analyzing numerous genes expressed in response to different growth stresses or tissues. In order to capture accurate and reliable results in gene expression studies, it is necessary to select appropriate internal reference genes. They must have relatively uniform expression in the majority of tested samples, thereby reducing the impact of the sample itself on the RT-qPCR results. The selection and validation of 14 candidate reference genes, including α-tub, β-tub, γ-tub, GAPDH, Tif-5a, Tef-1α, ATPase, acyl-CoA, U2ribo, S/G, Wlp, VtpAsE, E3upl, and Sdh, were carried out for gene expression analysis in P. tuber-regium mycelia in response to different temperatures, drought levels, and salinity shifts, fruitbody, and sclerotium. Four statistical algorithms, NormFinder, geNorm, BestKeeper, and RefFinder, were recruited to evaluate the transcription stability of candidate reference genes. The RT-qPCR expression stability analysis indicated that the E3upl and Tif-5a were the most stable expressed genes among all experimental samples, so the combination of these two reference genes is suitable for the analysis of gene expression patterns in P. tuber-regium.

Determination of the panel of reference genes for quantitative real-time PCR in fetal and adult rat intestines

In quantitative real-time PCR (qRT-PCR) detection, the stability of reference genes varies with different organs, tissue locations, sex and developmental stages. This study aimed to screen out and determine the optimal panel of reference genes of the intestine in pre- and post-natal rats of different sex. We used qRT-PCR to detect the mRNA expression of six commonly used reference genes (ACTB, GAPDH, HPRT1, B2M, RPLPO and SDHA) in rat intestines at gestational day 21 (GD21) and postnatal week 12 (PW12). Using GeNorm, BestKeeper and NormFinder software comprehensively analyzed the stability of candidate reference genes and screened out stable reference genes. Further, we used the pathological model of prenatal dexamethasone exposure (PDE) to verify the stability of the selected panel of reference genes. Based on the results of the software analysis, the optimal panel of reference genes in the fetal rat intestine was SDHA + ACTB, and the adult rat small intestine and colon were ACTB + HPRT1 and RPLP0 + GAPDH, respectively. There was no significant sex difference in the above results. Besides, in the PDE model, the results were consistent with those under physiological conditions. Therefore, the stability of intestinal reference genes in fetal rats and adult rats was different, and the intestinal reference genes of adult rats were intestinal segments-specific. The selected panel of reference genes was still stable under pathological conditions. This study determined the optimal panel of reference genes of pre- and post-natal rat intestines and provided reliable reference genes for the qRT-PCR analysis of rat intestines.

Monoterpene synthases responsible for the terpene profile of anther glands in Eucalyptus polybractea R.T. Baker (Myrtaceae)

Research on terpene biosynthesis in the genus Eucalyptus (Myrtaceae) is poorly developed, but recently large numbers of terpene synthase (TPS) genes have been identified. Few of these have been characterized or their expression localized to specific tissues. A prime candidate for detailed examination of TPS gene expression is the bisexual eucalypt flower-composed of male and female reproductive organs, and vegetative tissues that may express different TPS genes. We aimed to characterize and compare the terpene profile and TPS genes expressed in anthers and gynoecia in the high oil-yielding Eucalyptus polybractea R.T. Baker. We hypothesized that gynoecia will produce greater amounts of defensive terpenes, whereas anthers will have a terpene profile that is biased towards a role in pollination. Microscopy of isolated anthers showed them to possess a single, prominent oil gland. Chemical analysis of whole floral structures at different stages of development showed total oil per unit dry mass increased as flower buds expanded, with highest concentrations in mature flower buds just prior to flower opening. The oil profile of gynoecia was dominated by the monoterpene 1,8-cineole, whereas that of isolated anthers were enriched with the monoterpene α-pinene. Through transcriptomic analysis and recombinant protein expression, we were able to identify monoterpene synthases responsible for the different profiles. Synthases for α-pinene and 1,8-cineole were expressed in each tissue type, but the relative expression of the former was higher in anthers. Sequence comparison and site-directed mutagenesis of the α-pinene synthase allowed us to identify amino acids that influence the α-pinene to β-pinene ratio of the product profile. We suggest the terpene constituents of anthers may have multiple roles including attracting pollinators through emission of volatile α-pinene, deterrence of palynivores through emission of volatile 1,8-cineole and adhesion of pollen to pollinators via the release of sticky α-pinene onto the anther surface.

Panel of suitable reference genes and its gender differences of fetal rat liver under physiological conditions and exposure to dexamethasone during pregnancy

The panel of suitable reference genes in the fetal liver have not been reported. In this study, five commonly used reference genes (GAPDH, β-actin, Rn18 s, Rpl13a, and Rps29) were firstly selected as candidates. Bestkeeper, GeNorm, and NormFinder software were then used to screen out the panel of suitable reference genes of male and female fetal rat liver under physiological and prenatal dexamethasone exposure (PDE) conditions. Finally, we verified the reliability of the screened panel of reference genes by standardizing sterol regulatory element binding protein 1c (SREBP1c) expression with different reference genes. The results showed that GAPDH + Rn18 s and GAPDH + Rpl13a were respectively the panel of suitable reference genes in male and female rat fetal liver under the physiological model, while Rn18 s + Rps29 and GAPDH + Rn18 s were respectively under the PDE model. The results showed that different reference genes affected the statistical results of SREBP1c expression, and the screened panel of suitable reference genes under the PDE model had smaller intragroup differences, when compared with other reference genes under physiological and PDE models. In conclusion, we screened and determined that the panel of suitable reference genes were GAPDH + Rn18 s and Rn18 s + Rps29 in the male rat fetal liver under physiological and PDE models, while they were GAPDH + Rpl13a and GAPDH + Rn18 s in the females, and confirmed that the selection of the panel of suitable reference genes in the fetal liver had gender differences and pathological model specificity.

Selection and Validation of the Optimal Panel of Reference Genes for RT-qPCR Analysis in the Developing Rat Cartilage

Real-time fluorescence quantitative PCR (RT-qPCR) is widely used to detect gene expression levels, and selection of reference genes is crucial to the accuracy of RT-qPCR results. Minimum Information for Publication of RT-qPCR Experiments (MIQE) proposes that using the panel of reference genes for RT-qPCR is conducive to obtaining accurate experimental results. However, the selection of the panel of reference genes for RT-qPCR in rat developing cartilage has not been well documented. In this study, we selected eight reference genes commonly used in rat cartilage from literature (GAPDH, ACTB, 18S, GUSB, HPRT1, RPL4, RPL5, and SDHA) as candidates. Then, we screened out the optimal panel of reference genes in female and male rat cartilage of fetus (GD20), juvenile (PW6), and puberty (PW12) in physiology with stability analysis software of genes expression. Finally, we verified the reliability of the selected panel of reference genes with the rat model of intrauterine growth retardation (IUGR) induced by prenatal dexamethasone exposure (PDE). The results showed that the optimal panel of reference genes in cartilage at GD20, PW6, and PW12 in physiology was RPL4 + RPL5, which was consistent with the IUGR model, and there was no significant gender difference. Further, the results of standardizing the target genes showed that RPL4 + RPL5 performed smaller intragroup differences than other panels of reference genes or single reference genes. In conclusion, we found that the optimal panel of reference genes in female and male rat developing cartilage was RPL4 + RPL5, and there was no noticeable difference before and after birth.

EgPHI-1, a PHOSPHATE-INDUCED-1 gene from Eucalyptus globulus, is involved in shoot growth, xylem fiber length and secondary cell wall properties

MAIN CONCLUSION

EgPHI-1 is a member of PHI-1/EXO/EXL protein family. Its overexpression in tobacco resulted in changes in biomass partitioning, xylem fiber length, secondary cell wall thickening and composition, and lignification. Here, we report the functional characterization of a PHOSPHATE-INDUCED PROTEIN 1 homologue showing differential expression in xylem cells from Eucalyptus species of contrasting phenotypes for wood quality and growth traits. Our results indicated that this gene is a member of the PHI-1/EXO/EXL family. Analysis of the promoter cis-acting regulatory elements and expression responses to different treatments revealed that the Eucalyptus globulus PHI-1 (EgPHI-1) is transcriptionally regulated by auxin, cytokinin, wounding and drought. EgPHI-1 overexpression in transgenic tobacco changed the partitioning of biomass, favoring its allocation to shoots in detriment of roots. The stem of the transgenic plants showed longer xylem fibers and reduced cellulose content, while the leaf xylem had enhanced secondary cell wall thickness. UV microspectrophotometry of individual cell wall layers of fibers and vessels has shown that the transgenic plants exhibit differences in the lignification of S2 layer in both cell types. Taken together, the results suggest that EgPHI-1 mediates the elongation of secondary xylem fibers, secondary cell wall thickening and composition, and lignification, making it an attractive target for biotechnological applications in forestry and biofuel crops.

Reference gene selection for real-time quantitative PCR normalization in Hemarthria compressa and Hemarthria altissima leaf tissue

Hemarthria compressa and Hemarthria altissima are widely used as livestock forage and play important roles in tropical and subtropical grassland agricultural systems promoting healthy ecological environment and the development of animal husbandry. Leaf tissue of "Yaan" limpograss (H. compressa) and "H255" whip grass (H. altissima) were used to test the mRNA expression levels of 12 reference genes using RT-qPCR. The Delta-Ct method, BestKeeper (ver. 1.0), Genorm (ver. 3.5), Normfinder (ver. 0.953) and RefFinder were used to analyze the expression stability of the 12 reference genes under drought, salt, acid-aluminum and cold stresses to provide significant technical support for the study of gene expression under various abiotic stresses in Hemarthria. The results showed that the candidate reference genes showed divergent expression levels under various abiotic stresses. Among the genes that were selected, CL18892 showed the highest expression stability under salt stress in the leaf tissue. eEF-1α was the most stable gene under cold and acid-aluminum stresses and CL16384 was comparatively the most suitable genes under drought stress. As a whole, according to RefFinder analysis, CYP5, BMK.74327 and CL21527 were the most suitable reference genes for studying the effects of abiotic stress in Hemarthria. In general, CL16812 and CL18038 were not suitable reference genes under abiotic stress conditions that were examined in this study.

Genome-Wide Characterization, Evolution, and Expression Profiling of VQ Gene Family in Response to Phytohormone Treatments and Abiotic Stress in Eucalyptus grandis

VQ genes play important roles in plant development, growth, and stress responses. However, little information regarding the functions of VQ genes is available for Eucalyptus grandis. In our study, genome-wide characterization and identification of VQ genes were performed in E. grandis. Results showed that 27 VQ genes, which divided into seven sub-families (I-VII), were found, and all but two VQ genes showed no intron by gene structure and conserved motif analysis. To further identify the function of EgrVQ proteins, gene expression analyses were also developed under hormone treatments (brassinosteroids, methyl jasmonate, salicylic acid, and abscisic acid) and abiotic conditions (salt stress, cold 4 °C, and heat 42 °C). The results of a quantitative real-time PCR analysis indicated that the EgrVQs were variously expressed under different hormone treatments and abiotic stressors. Our study provides a comprehensive overview of VQ genes in E. grandis, which will be beneficial in the molecular breeding of E. grandis to promote its resistance to abiotic stressors; the results also provide a basis from which to conduct further investigation into the functions of VQ genes in E. grandis.

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

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

Temporal regulation of terpene synthase gene expression in Eucalyptus globulus leaves upon ozone and wounding stresses: relationships with stomatal ozone uptake and emission responses

Ozone and wounding are key abiotic factors but, their interactive effects on temporal changes in terpene synthase gene expression and emission responses are poorly understood. Here, we applied combined acute ozone and wounding stresses to the constitutive isoprenoid-emitter Eucalyptus globulus and studied how isoprene, 1,8-cineole, and isoledene synthase genes were regulated, and how the gene expression was associated with temporal changes in photosynthetic characteristics, product emission rates, and stomatal ozone uptake through recovery phase. Photosynthetic characteristics and emission rate of isoprene, 1,8-cineole, and isoledene were synergistically altered, while three TPS gene expressions were antagonistically altered by combined stress applications. A time-delay analysis indicated that the best correspondences between gene expression and product emission rates were observed for 0 h time-shift for wounding and 0-2 h time-shifts for separate ozone, and combined ozone and wounding treatments. The best correspondence between ozone uptake and gene expression was observed for 0-4 h time-shifts for separate ozone and combined ozone and wounding treatments. Overall, this study demonstrated that expression profiles of isoprene, the monoterpene 1,8-cineole, and the sesquiterpene isoledene synthase genes differentially influenced their corresponding product emissions for separate and combined ozone and wounding treatments through recovery.

Characterization of Brassinazole resistant (BZR) gene family and stress induced expression in Eucalyptus grandis

Brassinosteroids (BRs) are a group of plant hormones which play a pivotal role in modulating cell elongation, stress responses, vascular differentiation and senescence. In response to BRs, BRASSINAZOLE-RESISTANT (BZR) transcription factors (TFs) accumulate in the nucleus, where they modulate thousands of target genes and coordinate many biological processes, especially in regulating defense against biotic and abiotic stresses. In this study, 6 BZR TFs of Eucalyptus grandis (EgrBZR) from a genome-wide survey were characterized by sequence analysis and expression profiling against several abiotic stresses. The results showed that BZR gene family in Eucalyptus was slightly smaller compared to Populus and Arabidopsis, but all phylogenetic groups were represented. Various systematic in silico analysis of these TFs validated the basic properties of BZRs, whereas comparative studies showed a high degree of similarity with recognized BZRs of other plant species. In the organ-specific expression analyses, 4 EgrBZRs were expressed in vascular tissue indicating their possible functions in wood formation. Meanwhile, almost all EgrBZR genes showed differential transcript abundance levels in response to exogenously applied BR, MeJA, and SA, and salt and cold stresses. Besides, protein interaction analysis showed that all EgrBZR genes were associated with BR signaling directly or indirectly. These TFs were proposed as transcriptional activators or repressors of abiotic stress response and growth and development pathways of E. grandis by participating in BR signaling processes. These findings would be helpful in resolving the regulatory mechanism of EgrBZRs in stress resistance conditions but require further functional study of these potential TFs in Eucalyptus.

Exogenous GA3 application altered morphology, anatomic and transcriptional regulatory networks of hormones in Eucalyptus grandis

Gibberellins (GAs) play a key role in plant growth and development including cell elongation, cell expansion, and xylem differentiation. Eucalyptus are the world's most widely planted hardwood trees providing fiber and energy. However, the roles of GAs in Eucalyptus remain unclear and their effects on xylem development remain to be determined. In this study, E. grandis plants were treated with 0.10 mg L^-1 GA₃ and/or paclobutrazol (PAC, a GA inhibitor). The growth of shoot and root were recorded, transverse sections of roots and stems were stained using toluidine blue, and expression levels of genes related to hormone response and secondary cell wall biosynthesis were analyzed by quantitative real-time PCR. The results showed that GA₃ dramatically promoted the length of shoot and root, but decreased the diameter of root and stem. Exogenous GA₃ application also significantly promoted xylem development in both stem and root. Expression analysis revealed that exogenous GA₃ application altered the transcript levels of genes related to the GA biosynthetic pathway and GA signaling, as well as genes related to auxin, cytokinin, and secondary cell wall. These findings suggest that GAs may interact with other hormones (such as auxin and cytokinin) to regulate the expression of secondary cell wall biosynthesis genes and trigger xylogenesis in Eucalyptus plants.

A Whole-Transcriptome Approach to Evaluating Reference Genes for Quantitative Gene Expression Studies: A Case Study in Mimulus

While quantitative PCR (qPCR) is widely recognized as being among the most accurate methods for quantifying gene expression, it is highly dependent on the use of reliable, stably expressed reference genes. With the increased availability of high-throughput methods for measuring gene expression, whole-transcriptome approaches may be increasingly utilized for reference gene selection and validation. In this study, RNA-seq was used to identify a set of novel qPCR reference genes and evaluate a panel of traditional “housekeeping” reference genes in two species of the evolutionary model plant genus Mimulus. More broadly, the methods proposed in this study can be used to harness the power of transcriptomes to identify appropriate reference genes for qPCR in any study organism, including emerging and nonmodel systems. We find that RNA-seq accurately estimates gene expression means in comparison to qPCR, and that expression means are robust to moderate environmental and genetic variation. However, measures of expression variability were only in agreement with qPCR for samples obtained from a shared environment. This result, along with transcriptome-wide comparisons, suggests that environmental changes have greater impacts on expression variability than on expression means. We discuss how this issue can be addressed through experimental design, and suggest that the ever-expanding pool of published transcriptomes represents a rich and low-cost resource for developing better reference genes for qPCR.

The Eucalyptus Tonoplast Intrinsic Protein (TIP) Gene Subfamily: Genomic Organization, Structural Features, and Expression Profiles

Plant aquaporins are water channels implicated in various physiological processes, including growth, development and adaptation to stress. In this study, the Tonoplast Intrinsic Protein (TIP) gene subfamily of Eucalyptus, an economically important woody species, was investigated and characterized. A genome-wide survey of the Eucalyptus grandis genome revealed the presence of eleven putative TIP genes (referred as EgTIP), which were individually assigned by phylogeny to each of the classical TIP1-5 groups. Homology modeling confirmed the presence of the two highly conserved NPA (Asn-Pro-Ala) motifs in the identified EgTIPs. Residue variations in the corresponding selectivity filters, that might reflect differences in EgTIP substrate specificity, were observed. All EgTIP genes, except EgTIP5.1, were transcribed and the majority of them showed organ/tissue-enriched expression. Inspection of the EgTIP promoters revealed the presence of common cis-regulatory elements implicated in abiotic stress and hormone responses pointing to an involvement of the identified genes in abiotic stress responses. In line with these observations, additional gene expression profiling demonstrated increased expression under polyethylene glycol-imposed osmotic stress. Overall, the results obtained suggest that these novel EgTIPs might be functionally implicated in eucalyptus adaptation to stress.

Evaluation of stability and validation of reference genes for RT-qPCR expression studies in rice plants under water deficit

Many studies use strategies that allow for the identification of a large number of genes expressed in response to different stress conditions to which the plant is subjected throughout its cycle. In order to obtain accurate and reliable results in gene expression studies, it is necessary to use reference genes, which must have uniform expression in the majority of cells in the organism studied. RNA isolation of leaves and expression analysis in real-time quantitative polymerase chain reaction (RT-qPCR) were carried out. In this study, nine candidate reference genes were tested, actin 11 (ACT11), ubiquitin conjugated to E2 enzyme (UBC-E2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), beta tubulin (β-tubulin), eukaryotic initiation factor 4α (eIF-4α), ubiquitin 10 (UBQ10), ubiquitin 5 (UBQ5), aquaporin TIP41 (TIP41-Like) and cyclophilin, in two genotypes of rice, AN Cambará and BRS Querência, with different levels of soil moisture (20%, 10% and recovery) in the vegetative (V5) and reproductive stages (period preceding flowering). Currently, there are different softwares that perform stability analyses and define the most suitable reference genes for a particular study. In this study, we used five different methods: geNorm, BestKeeper, ΔCt method, NormFinder and RefFinder. The results indicate that UBC-E2 and UBQ5 can be used as reference genes in all samples and softwares evaluated. The genes β-tubulin and eIF-4α, traditionally used as reference genes, along with GAPDH, presented lower stability values. The gene expression of basic leucine zipper (bZIP23 and bZIP72) was used to validate the selected reference genes, demonstrating that the use of an inappropriate reference can induce erroneous results.

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

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

Selection of Reliable Reference Genes for Gene Expression Analysis under Abiotic Stresses in the Desert Biomass Willow, Salix psammophila

Salix psammophila is a desert shrub willow that has extraordinary adaptation to abiotic stresses and plays an important role in maintaining local ecosystems. Moreover, S. psammophila is regarded as a promising biomass feedstock because of its high biomass yields and short rotation coppice cycle. However, few suitable reference genes (RGs) for quantitative real-time polymerase chain reaction (qRT-PCR) constrain the study on normalization of gene expression in S. psammophila until now. Here, we investigated the expression stabilities of 14 candidate RGs across tissue types and under four abiotic stress treatments, including heat, cold, salt, and drought treatments. After calculation of PCR efficiencies, three different software, NormFinder, geNorm, and BestKeeper were employed to analyze systematically the qRT-PCR data, and the outputs were merged by RankAggreg software. The optimal RGs selected for gene expression analysis were EF1α (Elongation factor-1 alpha) and OTU (OTU-like cysteine protease family protein) for different tissue types, UBC (Ubiquitin-conjugating enzyme E2) and LTA4H (Leukotriene A-4 hydrolase homolog) for heat treatment, HIS (Histone superfamily protein H3) and ARF2 (ADP-ribosylation factor 2) for cold treatment, OTU and ACT7 (Actin 7) for salt treatment, UBC and LTA4H for drought treatment. The expression of UBC, ARF2, and VHAC (V-type proton ATPase subunit C) varied the least across tissue types and under abiotic stresses. Furthermore, the relative genes expression profiles of one tissue-specific gene WOX1a (WUSCHEL-related homeobox 1a), and four stress-inducible genes, including Hsf-A2 (Heat shock transcription factors A2), CBF3 (C-repeat binding factor 3), HKT1 (High-Affinity K⁺ Transporter 1), and GST (Glutathione S-transferase), were conducted to confirm the validity of the RGs in this study. These results provided an important RGs application guideline for gene expression characterization in S. psammophila.

Seasonal Variation of Carbon Metabolism in the Cambial Zone of Eucalyptus grandis

Eucalyptus species are the most widely hardwood planted in the world. It is one of the successful examples of commercial forestry plantation in Brazil and other tropical and subtropical countries. The tree is valued for its rapid growth, adaptability and wood quality. Wood formation is the result of cumulative annual activity of the vascular cambium. This cambial activity is generally related to the alternation of cold and warm, and/or dry and rainy seasons. Efforts have focused on analysis of cambial zone in response to seasonal variations in trees from temperate zones. However, little is known about the molecular changes triggered by seasonal variations in trees from tropical countries. In this work we attempted to establish a global view of seasonal alterations in the cambial zone of Eucalyptus grandis Hill ex Maiden, emphasizing changes occurring in the carbon metabolism. Using transcripts, proteomics and metabolomics we analyzed the tissues harvested in summer-wet and winter-dry seasons. Based on proteomics analysis, 70 proteins that changed in abundance were successfully identified. Transcripts for some of these proteins were analyzed and similar expression patterns were observed. We identified 19 metabolites differentially abundant. Our results suggest a differential reconfiguration of carbon partioning in E. grandis cambial zone. During summer, pyruvate is primarily metabolized via ethanolic fermentation, possibly to regenerate NAD(+) for glycolytic ATP production and cellular maintenance. However, in winter there seems to be a metabolic change and we found that some sugars were highly abundant. Our results revealed a dynamic change in E. grandis cambial zone due to seasonality and highlight the importance of glycolysis and ethanolic fermentation for energy generation and maintenance in Eucalyptus, a fast growing tree.

Concerted transcription of auxin and carbohydrate homeostasis-related genes underlies improved adventitious rooting of microcuttings derived from far-red treated Eucalyptus globulus Labill mother plants

Economically important plant species, such as Eucalyptus globulus, are often rooting recalcitrant. We have previously shown that far-red light enrichment applied to E. globulus donor-plants improved microcutting rooting competence and increased rooting zone/shoot carbohydrate ratio. To better understand this developmental response, the relative expression profiles of genes involved in auxin signaling (ARF6, ARF8, AGO1), biosynthesis (YUC3) and transport (AUX1, PIN1, PIN2); sucrose cleavage (SUS1, CWINV1), transport (SUC5), hexose phosphorylation (HXK1, FLN1) and starch biosynthesis (SS3) were quantified during adventitious rooting of E. globulus microcuttings derived from donor plants exposed to far-red or white light. Expression of auxin transport-related genes increased in the first days of root induction. Far-red enrichment of donor plants induced ARF6, ARF8 and AGO1 in microcuttings. The first two gene products could activate GH3 and other rooting related genes, whereas AGO1 deregulation of the repressor ARF17 may relief rooting inhibition. Increased sink strength at the basal stem with sucrose unloading in root tissue mediated by SUC and subsequent hydrolysis by SUS1 were also supported by gene expression profile. Fructose phosphorylation and starch biosynthesis could also contribute to proper carbon allocation at the site of rooting, as evidenced by increased expression of related genes. These data are in good agreement with increased contents of hexoses and starch at the cutting base severed from far-red exposed donor plants. To sum up, pathways integrating auxin and carbohydrate metabolism were activated in microcuttings derived from donor plants exposed to far red light enrichment, thereby improving rooting response in E. globulus.

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

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

Genome-Wide Identification and Validation of Reference Genes in Infected Tomato Leaves for Quantitative RT-PCR Analyses

The Gram-negative bacterium Xanthomonas campestris pv. vesicatoria (Xcv) causes bacterial spot disease of pepper and tomato by direct translocation of type III effector proteins into the plant cell cytosol. Once in the plant cell the effectors interfere with host cell processes and manipulate the plant transcriptome. Quantitative RT-PCR (qRT-PCR) is usually the method of choice to analyze transcriptional changes of selected plant genes. Reliable results depend, however, on measuring stably expressed reference genes that serve as internal normalization controls. We identified the most stably expressed tomato genes based on microarray analyses of Xcv-infected tomato leaves and evaluated the reliability of 11 genes for qRT-PCR studies in comparison to four traditionally employed reference genes. Three different statistical algorithms, geNorm, NormFinder and BestKeeper, concordantly determined the superiority of the newly identified reference genes. The most suitable reference genes encode proteins with homology to PHD finger family proteins and the U6 snRNA-associated protein LSm7. In addition, we identified pepper orthologs and validated several genes as reliable normalization controls for qRT-PCR analysis of Xcv-infected pepper plants. The newly identified reference genes will be beneficial for future qRT-PCR studies of the Xcv-tomato and Xcv-pepper pathosystems, as well as for the identification of suitable normalization controls for qRT-PCR studies of other plant-pathogen interactions, especially, if related plant species are used in combination with bacterial pathogens.

Selection of reference genes for transcriptional analysis of edible tubers of potato (Solanum tuberosum L.)

Potato (Solanum tuberosum) yield has increased dramatically over the last 50 years and this has been achieved by a combination of improved agronomy and biotechnology efforts. Gene studies are taking place to improve new qualities and develop new cultivars. Reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) is a bench-marking analytical tool for gene expression analysis, but its accuracy is highly dependent on a reliable normalization strategy of an invariant reference genes. For this reason, the goal of this work was to select and validate reference genes for transcriptional analysis of edible tubers of potato. To do so, RT-qPCR primers were designed for ten genes with relatively stable expression in potato tubers as observed in RNA-Seq experiments. Primers were designed across exon boundaries to avoid genomic DNA contamination. Differences were observed in the ranking of candidate genes identified by geNorm, NormFinder and BestKeeper algorithms. The ranks determined by geNorm and NormFinder were very similar and for all samples the most stable candidates were C2, exocyst complex component sec3 (SEC3) and ATCUL3/ATCUL3A/CUL3/CUL3A (CUL3A). According to BestKeeper, the importin alpha and ubiquitin-associated/ts-n genes were the most stable. Three genes were selected as reference genes for potato edible tubers in RT-qPCR studies. The first one, called C2, was selected in common by NormFinder and geNorm, the second one is SEC3, selected by NormFinder, and the third one is CUL3A, selected by geNorm. Appropriate reference genes identified in this work will help to improve the accuracy of gene expression quantification analyses by taking into account differences that may be observed in RNA quality or reverse transcription efficiency across the samples.

Identification of suitable reference genes for gene expression studies using quantitative polymerase chain reaction in lung cancer in vitro

The present study aimed to examine 10 housekeeping genes (HKGs), including 18s ribosomal RNA (18S), glyceraldehyde‑3‑phosphate dehydrogenase (GAPDH), ribosomal protein large P0 (RPLP0), β‑actin (ACTB), peptidylprolyl isomerase A (PPIA), phosphoglycerate kinase‑1 (PGK1), β‑2‑microglobulin (B2M), ribosomal protein LI3a (RPL13A), hypoxanthine phosphoribosyl transferase‑1 (HPRT1) and TATA box binding protein (TBP) in order to identify the most stable and suitable reference genes for use in expression studies in non‑small cell lung cancer. The mRNA expression encoding the panel of the 10 HKGs was determined using reverse transcription‑quantitative PCR (RT‑qPCR) in human lung cancer cell lines. Three software programs, BestKeeper, NormFinder and geNorm, were used to ascertain the most suitable reference genes to normalize the RNA input. The present study examined three lung cancer cell lines (A549, NCI‑H446 and NCI‑H460). The analysis of the experimental data using BestKeeper software revealed that all 10 HKGs were stable, with GADPH, followed by 18S being the most stable genes and PPIA and HPRT1 being the least stable genes. The NormFinder software results demonstrated that PPIA followed by ACTB were the most stable and B2M and RPLP0 were the least stable. The geNorm software results revealed that ACTB and PGK1, followed by PPIA were the most stable genes and B2M and RPLP0 were identified as the least stable genes. Due to discrepancies in the ranking orders of the reference genes obtained by different analyzing software programs, it was not possible to determine a single universal reference gene. The suitability of selected reference genes requires unconditional validation prior to each study. Based on the three analyzing programs, ACTB, PPIA and PGK1 were the most stable reference genes in lung cancer cell lines.

Identification and validation of reference genes for Populus euphratica gene expression analysis during abiotic stresses by quantitative real-time PCR

Populus euphratica is the only arboreal species that is established in the world's largest shifting-sand desert in China and is well-adapted to the extreme desert environment, so it is widely considered a model system for researching into abiotic stress resistance of woody plants. However, few P. euphratica reference genes (RGs) have been identified for quantitative real-time polymerase chain reaction (qRT-PCR) until now. Validation of suitable RGs is essential for gene expression normalization research. In this study, we screened 16 endogenous candidate RGs in P. euphratica leaves in six abiotic stress treatments, including abscisic acid (ABA), cold, dehydration, drought, short-duration salt (SS) and long-duration salt (LS) treatments, each with 6 treatment gradients. After calculation of PCR efficiencies, three different software tools, NormFinder, geNorm and BestKeeper, were employed to analyze the qRT-PCR data systematically, and the outputs were merged by means of a non-weighted unsupervised rank aggregation method. The genes selected as optimal for gene expression analysis of the six treatments were RPL17 (ribosomal protein L17) in ABA, EF1α (elongation factor-1 alpha) in cold, HIS (histone superfamily protein H3) in dehydration, GIIα in drought and SS, and TUB (tubulin) in LS. The expression of 60S (the 60S ribosomal protein) varied the least during all treatments. To illustrate the suitability of these RGs, the relative quantifications of three stress-inducible genes, PePYL1, PeSCOF-1 and PeSCL7 were investigated with different RGs. The results, calculated using qBasePlus software, showed that compared with the least-appropriate RGs, the expression profiles normalized by the recommended RGs were closer to expectations. Our study provided an important RG application guideline for P. euphratica gene expression characterization.

Comparative transcriptomics of the model mushroom Coprinopsis cinerea reveals tissue-specific armories and a conserved circuitry for sexual development

Background

It is well known that mushrooms produce defense proteins and secondary metabolites against predators and competitors; however, less is known about the correlation between the tissue-specific expression and the target organism (antagonist) specificity of these molecules. In addition, conserved transcriptional circuitries involved in developing sexual organs in fungi are not characterized, despite the growing number of gene expression datasets available from reproductive and vegetative tissue. The aims of this study were: first, to evaluate the tissue specificity of defense gene expression in the model mushroom Coprinopsis cinerea and, second, to assess the degree of conservation in transcriptional regulation during sexual development in basidiomycetes.

Results

In order to characterize the regulation in the expression of defense loci and the transcriptional circuitries controlling sexual reproduction in basidiomycetes, we sequenced the poly (A)-positive transcriptome of stage 1 primordia and vegetative mycelium of C. cinerea A43mutB43mut. Our data show that many genes encoding predicted and already characterized defense proteins are differentially expressed in these tissues. The predicted specificity of these proteins with regard to target organisms suggests that their expression pattern correlates with the type of antagonists these tissues are confronted with. Accordingly, we show that the stage 1 primordium-specific protein CC1G_11805 is toxic to insects and nematodes. Comparison of our data to analogous data from Laccaria bicolor and Schizophyllum commune revealed that the transcriptional regulation of nearly 70 loci is conserved and probably subjected to stabilizing selection. A Velvet domain-containing protein was found to be up-regulated in all three fungi, providing preliminary evidence of a possible role of the Velvet protein family in sexual development of basidiomycetes. The PBS-soluble proteome of C. cinerea primordia and mycelium was analyzed by shotgun LC-MS. This proteome data confirmed the presence of intracellular defense proteins in primordia.

Conclusions

This study shows that the exposure of different tissues in fungi to different types of antagonists shapes the expression pattern of defense loci in a tissue-specific manner. Furthermore, we identify a transcriptional circuitry conserved among basidiomycetes during fruiting body formation that involves, amongst other transcription factors, the up-regulation of a Velvet domain-containing protein.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-492) contains supplementary material, which is available to authorized users.

Identification of candidate reference genes in perennial ryegrass for quantitative RT-PCR under various abiotic stress conditions

Background

Quantitative real-time reverse-transcriptase PCR (qRT-PCR) is an important technique for analyzing differences in gene expression due to its sensitivity, accuracy and specificity. However, the stability of the expression of reference genes is necessary to ensure accurate qRT-PCR assessment of expression in genes of interest. Perennial ryegrass (Lolium perenne L.) is important forage and turf grass species in temperate regions, but the expression stability of its reference genes under various stresses has not been well-studied.

Methodology/Principal Findings

In this study, 11 candidate reference genes were evaluated for use as controls in qRT-PCR to quantify gene expression in perennial ryegrass under drought, high salinity, heat, waterlogging, and ABA (abscisic acid) treatments. Four approaches – Delta C_T, geNorm, BestKeeper and Normfinder were used to determine the stability of expression in these reference genes. The results are consistent with the idea that the best reference genes depend on the stress treatment under investigation. Eukaryotic initiation factor 4 alpha (eIF4A), Transcription elongation factor 1 (TEF1) and Tat binding protein-1 (TBP-1) were the three most stably expressed genes under drought stress and were also the three best genes for studying salt stress. eIF4A, TBP-1, and Ubiquitin-conjugating enzyme (E2) were the most suitable reference genes to study heat stress, while eIF4A, TEF1, and E2 were the three best reference genes for studying the effects of ABA. Finally, Ubiquitin (UBQ), TEF1, and eIF4A were the three best reference genes for waterlogging treatments.

Conclusions/Significance

These results will be helpful in choosing the best reference genes for use in studies related to various abiotic stresses in perennial ryegrass. The stability of expression in these reference genes will enable better normalization and quantification of the transcript levels for studies of gene expression in such studies.

Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses

Strawberry (Fragaria spp) is an emerging model for the development of basic genomics and recombinant DNA studies among rosaceous crops. Functional genomic and molecular studies involve relative quantification of gene expression under experimental conditions of interest. Accuracy and reliability are dependent upon the choice of an optimal reference control transcript. There is no information available on validated endogenous reference genes for use in studies testing strawberry-pathogen interactions. Thirteen potential pre-selected strawberry reference genes were tested against different tissues, strawberry cultivars, biotic stresses, ripening and senescent conditions, and SA/JA treatments. Evaluation of reference candidate's suitability was analyzed by five different methodologies, and information was merged to identify best reference transcripts. A combination of all five methods was used for selective classification of reference genes. The resulting superior reference genes, FaRIB413, FaACTIN, FaEF1α and FaGAPDH2 are strongly recommended as control genes for relative quantification of gene expression in strawberry. This report constitutes the first systematic study to identify and validate optimal reference genes for accurate normalization of gene expression in strawberry plant defense response studies.

Xylem transcription profiles indicate potential metabolic responses for economically relevant characteristics of Eucalyptus species

Background

Eucalyptus is one of the most important sources of industrial cellulose. Three species of this botanical group are intensively used in breeding programs: E. globulus, E. grandis and E. urophylla. E. globulus is adapted to subtropical/temperate areas and is considered a source of high-quality cellulose; E. grandis grows rapidly and is adapted to tropical/subtropical climates; and E. urophylla, though less productive, is considered a source of genes related to robustness. Wood, or secondary xylem, results from cambium vascular differentiation and is mostly composed of cellulose, lignin and hemicelluloses. In this study, the xylem transcriptomes of the three Eucalyptus species were investigated in order to provide insights on the particularities presented by each of these species.

Results

Data analysis showed that (1) most Eucalyptus genes are expressed in xylem; (2) most genes expressed in species-specific way constitutes genes with unknown functions and are interesting targets for future studies; (3) relevant differences were observed in the phenylpropanoid pathway: E. grandis xylem presents higher expression of genes involved in lignin formation whereas E. urophylla seems to deviates the pathway towards flavonoid formation; (4) stress-related genes are considerably more expressed in E. urophylla, suggesting that these genes may contribute to its robustness.

Conclusions

The comparison of these three transcriptomes indicates the molecular signatures underlying some of their distinct wood characteristics. This information may contribute to the understanding of xylogenesis, thus increasing the potential of genetic engineering approaches aiming at the improvement of Eucalyptus forest plantations productivity.

Reference genes for high-throughput quantitative reverse transcription-PCR analysis of gene expression in organs and tissues of Eucalyptus grown in various environmental conditions

Interest in the genomics of Eucalyptus has skyrocketed thanks to the recent sequencing of the genome of Eucalyptus grandis and to a growing number of large-scale transcriptomic studies. Quantitative reverse transcription-PCR (RT-PCR) is the method of choice for gene expression analysis and can now also be used as a high-throughput method. The selection of appropriate internal controls is becoming of utmost importance to ensure accurate expression results in Eucalyptus. To this end, we selected 21 candidate reference genes and used high-throughput microfluidic dynamic arrays to assess their expression among a large panel of developmental and environmental conditions with a special focus on wood-forming tissues. We analyzed the expression stability of these genes by using three distinct statistical algorithms (geNorm, NormFinder and ΔCt), and used principal component analysis to compare methods and rankings. We showed that the most stable genes identified depended not only on the panel of biological samples considered but also on the statistical method used. We then developed a comprehensive integration of the rankings generated by the three methods and identified the optimal reference genes for 17 distinct experimental sets covering 13 organs and tissues, as well as various developmental and environmental conditions. The expression patterns of Eucalyptus master genes EgMYB1 and EgMYB2 experimentally validated our selection. Our findings provide an important resource for the selection of appropriate reference genes for accurate and reliable normalization of gene expression data in the organs and tissues of Eucalyptus trees grown in a range of conditions including abiotic stresses.

Validation of reference genes from Eucalyptus spp. under different stress conditions

Background

The genus Eucalyptus consists of approximately 600 species and subspecies and has a physiological plasticity that allows some species to propagate in different regions of the world. Eucalyptus is a major source of cellulose for paper manufacturing, and its cultivation is limited by weather conditions, particularly water stress and low temperatures. Gene expression studies using quantitative reverse transcription polymerase chain reaction (qPCR) require reference genes, which must have stable expression to facilitate the comparison of the results from analyses using different species, tissues, and treatments. Such studies have been limited in eucalyptus.

Results

Eucalyptus globulus Labill, Eucalyptus urograndis (hybrid from Eucalyptus urophylla S.T. Blake X Eucalyptus grandis Hill ex-Maiden) and E. uroglobulus (hybrid from E. urograndis X E. globulus) were subjected to different treatments, including water deficiency and stress recovery, low temperatures, presence or absence of light, and their respective controls. Except for treatment with light, which examined the seedling hypocotyl or apical portion of the stem, the expression analyses were conducted in the apical and basal parts of the stem. To select the best pair of genes, the bioinformatics tools GeNorm and NormFinder were compared. Comprehensive analyses that did not differentiate between species, treatments, or tissue types, showed that IDH (isocitrate dehydrogenase), SAND (SAND protein), ACT (actin), and A-Tub (α-tubulin) genes were the most stable. IDH was the most stable gene in all of the treatments.

Conclusion

Comparing these results with those of other studies on eucalyptus, we concluded that five genes are stable in different species and experimental conditions: IDH, SAND, ACT, A-Tub, and UBQ (ubiquitin). It is usually recommended a minimum of two reference genes is expression analysis; therefore, we propose that IDH and two others genes among the five identified genes in this study should be used as reference genes for a wide range of conditions in eucalyptus.

Identification and validation of reference genes for expression studies in human keratinocyte cell lines treated with and without interferon-γ - a method for qRT-PCR reference gene determination

Based on the exquisite sensitivity, reproducibility and wide dynamic range of quantitative reverse-transcription real-time polymerase chain reaction (qRT-PCR), it is currently the gold standard for gene expression studies. Target gene expression is calculated relative to a stably expressed reference gene. An ideal reference should be uniformly expressed during all experimental conditions within the given experimental system. However, no commonly applicable 'best' reference gene has been identified. Thus, endogenous controls must be determined for every experimental system. As no appropriate reference genes have been reported for immunological studies in keratinocytes, we aimed at identifying and validating a set of endogenous controls for these settings. An extensive validation of sixteen possible endogenous controls in a panel of 8 normal and transformed keratinocyte cell lines in experimental conditions with and without interferon-γ was performed. RNA and cDNA quality was stringently controlled. Candidate reference genes were assessed by TaqMan(®) qRT-PCR. Two different statistical algorithms were used to determine the most stably and reproducibly expressed housekeeping genes. mRNA abundance was compared and reference genes with widely different ranges of expression than possible target genes were excluded. Subsequent geNorm and NormFinder analyses identified GAPDH, PGK1, IPO8 and PPIA as the most stably expressed genes in the keratinocyte panel under the given experimental conditions. We conclude that the geometric means of expression values of these four genes represents a robust normalization factor for qRT-PCR analyses in interferon-γ-dependent gene expression studies in keratinocytes. The methodology and results herein may help other researchers by facilitating their choice of reference genes.