Identification of Reference Genes for Quantitative Gene Expression Studies in a Non-Model Tree Pistachio (Pistacia vera L.)

Temporal transcriptome and metabolite analyses provide insights into the biochemical and physiological processes underlying endodormancy release in pistachio (Pistacia vera L.) flower buds

Pistachio (Pistacia vera L.), an economically and nutritionally important tree crop, relies on winter chill for bud endodormancy break and subsequent blooming and nut production. However, insufficient winter chill poses an increasing challenge in pistachio growing regions. To gain a better understanding of the physiological and biochemical responses of endodormant pistachio buds to chilling accumulation, we investigated the global gene expression changes in flower buds of pistachio cv. Kerman that were cultivated at three different orchard locations and exposed to increasing durations of winter chill. The expression of genes encoding β-1,3-glucanase and β-amylase, enzymes responsible for breaking down callose (β-1,3-glucan) and starch (α-1,4-glucan), respectively, increased during the endodormancy break of pistachio buds. This result suggested that the breakdown of callose obstructing stomata as well as the release of glucose from starch enables symplasmic trafficking and provides energy for bud endodormancy break and growth. Interestingly, as chilling accumulation increased, there was a decrease in the expression of nine-cis-epoxycarotenoid dioxygenase (NCED), encoding an enzyme that uses carotenoids as substrates and catalyzes the rate-limiting step in abscisic acid (ABA) biosynthesis. The decrease in NCED expression suggests ABA biosynthesis is suppressed, thus reducing inhibition of endodormancy break. The higher levels of carotenoid precursors and a decrease in ABA content in buds undergoing endodormancy break supports this suggestion. Collectively, the temporal transcriptome and biochemical analyses revealed that the degradation of structural (callose) and non-structural (starch) carbohydrates, along with the attenuation of ABA biosynthesis, are critical processes driving endodormancy break in pistachio buds.

Proteomics approach to investigating osmotic stress effects on pistachio

Osmotic stress can occur due to some stresses such as salinity and drought, threatening plant survival. To investigate the mechanism governing the pistachio response to this stress, the biochemical alterations and protein profile of PEG-treated plants was monitored. Also, we selected two differentially abundant proteins to validate via Real-Time PCR. Biochemical results displayed that in treated plants, proline and phenolic content was elevated, photosynthetic pigments except carotenoid decreased and MDA concentration were not altered. Our findings identified a number of proteins using 2DE-MS, involved in mitigating osmotic stress in pistachio. A total of 180 protein spots were identified, of which 25 spots were altered in response to osmotic stress. Four spots that had photosynthetic activities were down-regulated, and the remaining spots were up-regulated. The biological functional analysis of protein spots exhibited that most of them are associated with the photosynthesis and metabolism (36%) followed by stress response (24%). Results of Real-Time PCR indicated that two of the representative genes illustrated a positive correlation among transcript level and protein expression and had a similar trend in regulation of gene and protein. Osmotic stress set changes in the proteins associated with photosynthesis and stress tolerance, proteins associated with the cell wall, changes in the expression of proteins involved in DNA and RNA processing occur. Findings of this research will introduce possible proteins and pathways that contribute to osmotic stress and can be considered for improving osmotic tolerance in pistachio.

Selection of the Reference Gene for Expression Normalization in Salsola ferganica under Abiotic Stress

Salsola ferganica is a natural desert herbaceous plant in the arid area of western and northwestern China. Because of its salt tolerance and drought resistance, it is of great significance in desert afforestation and sand-fixing capacity. There has been much research on the genes involved in plants under desert stresses in recent years. The application of the best internal reference genes for standardization was a critical procedure in analyzing the gene expression under different types. Even so, the reference gene has not been reported in the application of gene expression normalization of S. ferganica. In this study, nine reference genes (TUA-1726, TUA-1760, TUB, GAPDH, ACT, 50S, HSC70, APT, and U-box) in S. ferganica were adopted and analyzed under six different treatments (ABA, heat, cold, NaCl, methyl viologen (MV), and PEG). The applicability of candidate genes was evaluated by statistical software, including geNorm, NormFinder, BestKeeper, and RefFinder, based on their stability values in all the treatments. These results indicated that the simultaneous selection of two stable reference genes would fully standardize the optimization of the normalization research. To verify the feasibility of the above internal reference genes, the CT values of AP2/ERF transcription factor family genes were standardized using the most (ACT) and least (GAPDH) stable reference genes in S. ferganica seedlings under six abiotic stresses. The research showed that HSC70 and U-box were the most appropriate reference genes in ABA stressed samples, and ACT and U-box genes were the optimal references for heat-stressed samples. TUA-1726 and U-box showed the smallest value in gene expression levels of cold treatment. The internal reference groups of the best applicability for the other samples were U-box and ACT under NaCl treatment, ACT and TUA-1726 under MV stress, HSC70 and TUB under PEG treatment, and ACT in all samples. ACT and U-box showed higher stability than the other genes based on the comprehensive stability ranking of RefFinder, as determined by the geometric mean in this study. These results will contribute to later gene expression studies in other closely related species and provide an important foundation for gene expression analysis in S. ferganica.

Root vacuolar sequestration and suberization are prominent responses of Pistacia spp. rootstocks during salinity stress

Understanding the mechanisms of stress tolerance in diverse species is needed to enhance crop performance under conditions such as high salinity. Plant roots, in particular in grafted agricultural crops, can function as a boundary against external stresses in order to maintain plant fitness. However, limited information exists for salinity stress responses of woody species and their rootstocks. Pistachio (Pistacia spp.) is a tree nut crop with relatively high salinity tolerance as well as high genetic heterogeneity. In this study, we used a microscopy-based approach to investigate the cellular and structural responses to salinity stress in the roots of two pistachio rootstocks, Pistacia integerrima (PGI) and a hybrid, P. atlantica x P. integerrima (UCB1). We analyzed root sections via fluorescence microscopy across a developmental gradient, defined by xylem development, for sodium localization and for cellular barrier differentiation via suberin deposition. Our cumulative data suggest that the salinity response in pistachio rootstock species is associated with both vacuolar sodium ion (Na+) sequestration in the root cortex and increased suberin deposition at apoplastic barriers. Furthermore, both vacuolar sequestration and suberin deposition correlate with the root developmental gradient. We observed a higher rate of Na+ vacuolar sequestration and reduced salt-induced leaf damage in UCB1 when compared to P. integerrima. In addition, UCB1 displayed higher basal levels of suberization, in both the exodermis and endodermis, compared to P. integerrima. This difference was enhanced after salinity stress. These cellular characteristics are phenotypes that can be taken into account during screening for sodium-mediated salinity tolerance in woody plant species.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress

Quantitative real-time reverse-transcriptase PCR (qRT-PCR) has been frequently used for detecting gene expression. To obtain reliable results, selection of suitable reference genes is a fundamental and necessary step. Garlic (Allium sativum), a member from Alliaceae family, has been used both as a food flavoring and as a traditional medicine. In the present study, garlic plants were exposed to salt stress (200 mM NaCl) for 0, 1, 4 and 12 h, and garlic roots, bulbs, and leaves were harvested for subsequent analysis. The expression stability of eight candidate reference genes, eukaryotic translation initiation factor 4α (eIF-4α), actin (ACTIN), tubulin β-7 (TUB7), TAP42-interacting protein of 41 kDa (TIP41), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), SAND family protein (SAND), elongation factor 1 alpha (EF-1α), and protein phosphatase 2A (PP2A) were evaluated by geNorm, NormFinder, and BestKeeper. All genes tested displayed variable expression profiles under salt stress. In the leaf and root group, ACTIN was the best reference gene for normalizing gene expression. In garlic clove, ACTIN and SAND were the least variable, and were suitable for gene expression studies under salt stress; these two genes also performed well in all samples tested. Based on our results, we recommend that it is essential to use specific reference genes in different situations to obtain accurate results. Using a combination of multiple stable reference genes, such as ACTIN and SAND, to normalize gene expression is encouraged. The results from the study will be beneficial for accurate determination of gene expression in garlic and other plants.

Reference gene selection for qRT-PCR analysis of flower development in Lagerstroemia indica and L. speciosa

Quantitative real-time polymerase chain reaction (qRT-PCR) is a prevalent method for gene expression analysis, depending on the stability of the reference genes for data normalization. Lagerstroemia indica and L. speciosa are popular ornamental plants which are famous for the long flowering period. However, no systematic studies on reference genes in Lagerstroemia have yet been conducted. In the present study, we selected nine candidate reference genes (GAPDH, TUA, TUB, 18S, RPII, EF-1α, ATC, EIF5A and CYP) and evaluated their expression stability in different tissues during floral development of L. indica and L. speciosa using four algorithms (geNorm, NormFinder, BestKeeper and, RefFinder). Results showed that RPII and EF-1α were the most stably expressed and suitable reference genes for both of Lagerstroemia species. Moreover, ACT exhibited high expression stability in L. indica and GAPDH was a suitable reference gene for L. speciosa in different flower development stages. TUB was an unsuitable reference gene for gene expression normalization due to significant variations in expression across all samples. Finally, we verified the reliability of the selected candidate reference genes by amplifying an AGAMOUS homolog (LsAG1) of Arabidopsis thaliana. This study provides a list of suitable reference genes, thereby broadening the genetic basis of the gene expression patterns in Lagerstroemia species.

A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery

Background

Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress.

Results

RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars.

Conclusion

This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance.

Electronic supplementary material

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