Nitrogen starvation, salt and heat stress in coffee (Coffea arabica L.): identification and validation of new genes for qPCR normalization

RGeasy: a reference gene analysis tool for gene expression studies via RT-qPCR

Gene expression through RT-qPCR can be performed by the relative quantification method, which requires the expression normalization through reference genes. Therefore, it is essential to validate, experimentally, the candidate reference genes. Thus, although there are several studies that are performed to identify the most stable reference genes, most them validate genes for very specific conditions, not exploring the whole potential of the research since not all possible combinations of treatments and/or conditions of the study are explored. For this reason, new experiments must be conducted by researchers that have interest in analyzing gene expression of treatments and/or conditions present, but not explored, in these studies. Here, we present the RGeasy tool, which aims to facilitate the selection of reference genes, allowing the user to choose genes for a greater number of combinations of treatments/conditions, compared to the ones present in the original articles, through just a few clicks. RGeasy was validated with RT-qPCR data from gene expression studies performed in two coffee species, Coffea arabica and Coffea canephora, and it can be used for any animal, plant or microorganism species. In addition to displaying a rank of the most stable reference genes for each condition or treatment, the user also has access to the primer pairs for the selected reference genes.

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

Crosstalk Between Ethylene and Abscisic Acid During Changes in Soil Water Content Reveals a New Role for 1-Aminocyclopropane-1- Carboxylate in Coffee Anthesis Regulation

Coffee (Coffea arabica L.) presents an asynchronous flowering regulated by an endogenous and environmental stimulus, and anthesis occurs once plants are rehydrated after a period of water deficit. We evaluated the evolution of Abscisic Acid (ABA), ethylene, 1-aminocyclopropane-1-carboxylate (ACC) content, ACC oxidase (ACO) activity, and expression analysis of the Lysine Histidine Transporter 1 (LHT1) transporter, in the roots, leaves, and flower buds from three coffee genotypes (C. arabica L. cv Oeiras, Acauã, and Semperflorens) cultivated under field conditions with two experiments. In a third field experiment, the effect of the exogenous supply of ACC in coffee anthesis was evaluated. We found an increased ACC level, low ACO activity, decreased level of ethylene, and a decreased level of ABA in all tissues from the three coffee genotypes in the re-watering period just before anthesis, and a high expression of the LHT1 in flower buds and leaves. The ethylene content and ACO activity decreased from rainy to dry period whereas the ABA content increased. A higher number of opened and G6 stage flower buds were observed in the treatment with exogenous ACC. The results showed that the interaction of ABA-ACO-ethylene and intercellular ACC transport among the leaves, buds, and roots in coffee favors an increased level of ACC that is most likely, involved as a modulator in coffee anthesis. This study provides evidence that ACC can play an important role independently of ethylene in the anthesis process in a perennial crop.

LysM receptors in Coffea arabica: Identification, characterization, and gene expression in response to Hemileia vastatrix

Pathogen‐associated molecular patterns (PAMPs) are recognized by pattern recognition receptors (PRRs) localized on the host plasma membrane. These receptors activate a broad-spectrum and durable defense, which are desired characteristics for disease resistance in plant breeding programs. In this study, candidate sequences for PRRs with lysin motifs (LysM) were investigated in the Coffea arabica genome. For this, approaches based on the principle of sequence similarity, conservation of motifs and domains, phylogenetic analysis, and modulation of gene expression in response to Hemileia vastatrix were used. The candidate sequences for PRRs in C. arabica (Ca1-LYP, Ca2-LYP, Ca1-CERK1, Ca2-CERK1, Ca-LYK4, Ca1-LYK5 and Ca2-LYK5) showed high similarity with the reference PRRs used: Os-CEBiP, At-CERK1, At-LYK4 and At-LYK5. Moreover, the ectodomains of these sequences showed high identity or similarity with the reference sequences, indicating structural and functional conservation. The studied sequences are also phylogenetically related to the reference PRRs described in Arabidopsis, rice, and other plant species. All candidates for receptors had their expression induced after the inoculation with H. vastatrix, since the first time of sampling at 6 hours post‐inoculation (hpi). At 24 hpi, there was a significant increase in expression, for most of the receptors evaluated, and at 48 hpi, a suppression. The results showed that the candidate sequences for PRRs in the C. arabica genome display high homology with fungal PRRs already described in the literature. Besides, they respond to pathogen inoculation and seem to be involved in the perception or signaling of fungal chitin, acting as receptors or co-receptors of this molecule. These findings represent an advance in the understanding of the basal immunity of this species.

Constitutive Defense Strategy of Coffee Under Field Conditions: A Comparative Assessment of Resistant and Susceptible Cultivars to Rust

Coffea arabica is the most economically important coffee species worldwide. However, its production is severely limited by diseases such as rust. The mechanisms underlying constitutive defense responses in coffee are still poorly understood, compared with induced defense mechanisms. We aimed to characterize constitutive defense responses of thirteen cultivars of C. arabica. Cultivars were classified under field conditions according to the level of resistance to rust: resistant (R), moderately resistant (MR), and susceptible (S). Based on this classification, the stability of eight reference genes (RGs) was evaluated. The most stable RGs were EF1α, APT1, and 24S. We also evaluated the expression of CaWRKY1, CaPAL1, CaCAD1, and CaPOX1, and activities of PAL, CAD, and POX, which are involved in lignin biosynthesis, and leaf content of total phenolic compounds and lignin. Gene expression and enzymatic activity were not correlated with defense metabolites in the R cultivar group but showed a negative correlation with phenolic compounds in MR cultivars. Cultivar S showed positive correlations of gene expression and enzyme activity with phenolic compounds. These results may assist coffee breeding programs regarding selection of genotypes and in optimization of rust resistance.

The urea transporter DUR3 is differentially regulated by abiotic and biotic stresses in coffee plants

The high costs of N fertilizers in the coffee production emphasizes the need to optimize fertilization practices and improve nitrogen use efficiency. Urea is widespread in nature, characterizing itself as a significant source of nitrogen for the growth and development of several organisms. Thus, the characterization of genes involved in urea transport in coffee plants is an important research topic for the sustainable production of this valuable cash crop. In the current study, we evaluated the expression of the DUR3 gene under abiotic and biotic stresses in coffee plants. Here, we show that the expression of a high-affinity urea transporter gene (CaDUR3) was up-regulated by N starvation in leaves and roots of two out of three C. arabica cultivars examined. Moreover, the CaDUR3 gene was differentially expressed in coffee plants under different abiotic and biotic stresses. In plants of cv. IAPAR59, CaDUR3 showed an increased expression in leaves after exposure to water deficit and heat stress, while it was downregulated in plants under salinity. Upon infection with H. vastatrix (coffee rust), the CaDUR3 was markedly up-regulated at the beginning of the infection process in the disease susceptible Catuaí Vermelho 99 in comparison with the resistant cultivar. These results indicate that besides urea acquisition and N-remobilization, CaDUR3 gene may be closely involved in the response to various stresses.

Low-Copy Genes in Terpenoid Metabolism: The Evolution and Expression of MVK and DXR Genes in Angiosperms

Terpenoids are a diverse class of metabolites that impact plant metabolism in response to environmental cues. They are synthesized either via a predominantly cytosolic (MVA) pathway or a plastidic pathway (MEP). In Arabidopsis, several enzymes from the MVA and MEP pathways are encoded by gene families, excluding MVK and DXR, which are single-copy genes. In this study, we assess the diversity, evolution and expression of DXR and MVK genes in selected angiosperms and Coffea arabica in particular. Evolutionary analysis revealed that DXR and MVK underwent purifying selection, but the selection effect for DXR was stronger than it was for MVK. Digital gene expression (DGE) profile analysis of six species revealed that expression levels of MVK in flowers and roots were high, whereas for DXR peak values were observed in leaves. In C. arabica, both genes were highly expressed in flowers, and CaDXR was upregulated in response to methyl jasmonate. C. arabica DGE data were validated by assessing gene expression in selected organs, and by plants treated with hexanoic acid (Hx) using RT-qPCR. MVK expression was upregulated in roots treated with Hx. CaDXR was downregulated in leaves by Hx treatment in a genotype-specific manner, indicating a differential response to priming.

Identification and evaluation of reference genes for quantitative real-time PCR analysis in Passiflora edulis under stem rot condition

Passion fruit (Passiflora edulis), an important tropical and subtropical fruit, has a high edible and medicinal value. Stem rot disease is one of the most important diseases of passion fruit. An effective way for control and prevention of this disease is to identify the genes associated with resistance to this disease. Quantitative real-time PCR (RT-qPCR) has mainly been widely applied to detect gene expression because of its simplicity, fastness, low cost and high sensitivity. One of the requirements for RT-qPCR is the availability of suitable reference genes for normalization of gene expression. However, currently, no Passiflora edulis reference genes have been identified andthus it has hindered the gene expression studies in this plant. The present study aimed to address this issue. We analyzed sixteen candidate reference genes, including nine common (GAPDH, UBQ, ACT1, ACT2, EF-1α-1, EF-1α-2, TUA, NADP, and GBP) and seven novel genes (C13615, C24590, C27182, C10445, C21209, C22199, and C22526), in different tissues (stem, leaf, flower and fruit) of two accessions under stem rot condition. We calculated the expression stability in twenty-four samples using the ΔCt, GeNorm, NormFinder, BestKeeper and RefFinder. The results showed that both C21209 and EF-1α-2 were sufficient to normalize gene expression under stem rot, whereas the commonly used reference genes, GAPDH and UBQ, were the least stable ones. The expression patterns of PeUFC under stem rot condition normalized by stable and unstable reference genes indicated the suitability of using the optimal reference genes. To our knowledge, this is the first systematic study of reference genes in Passiflora edulis, which identified a number of reliable reference genes suitable for gene expression studies in Passiflora edulis by RT-qPCR.

Primer design and amplification efficiencies are crucial for reliability of quantitative PCR studies of caffeine biosynthetic N-methyltransferases in coffee

Primers having suboptimal amplification efficiencies were shown to falsely represent fold change expression of the N-methyltransferases gene family involved in caffeine biosynthesis in Coffea canephora. To study this phenomenon, the role of stability of the internal reference gene, as well as the amplification efficiency correction of the primers was investigated. GAPDH and Ubiquitin exhibited a good stability for studying the ontogeny of endosperm tissue, as well as the leaf transcriptome during stress from salicylic acid, methyl jasmonate, PEG-mediated drought and sudden exposure to light. Ubiquitin manifested low variation in Cq under all these stress regimes and in endosperm ontogeny with 30.1-30.9 in the best dataset and 28.8-30.9 in the most deviating dataset. It was observed that problems arising due to improper amplification efficiency of the target or reference genes or both could lead to misinterpretation of gene expression levels. Quantitative RT-PCR performed at a sub-optimal efficiency of GAPDH reference gene at 1.68 led to the faulty interpretation of 2.007 folds upregulation by the 2^-ΔΔCt method and 1.705 folds upregulation by Efficiency method for the first NMT (Xanthosine methyltransferase), which actually is repressed during dark acclimatization of coffee plants. Efficiency correction improved the reliability of the expression data and also indicated a downregulation of this gene by 0.485 folds and 0.474 folds using 2^-ΔΔCt and E method, respectively, in concordance to earlier reports. Hence, efficiency correction of the primers having suboptimal efficiencies is an absolute prerequisite for the accurate calculation of fold change using quantitative RT-PCR.

An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

Cowpea and abiotic stresses: identification of reference genes for transcriptional profiling by qPCR

BACKGROUND

Due to cowpea ability to fix nitrogen in poor soils and relative tolerance to drought and salt stresses, efforts have been directed to identifying genes and pathways that confer stress tolerance in this species. Real-time quantitative PCR (qPCR) has been widely used as the most reliable method to measure gene expression, due to its high accuracy and specificity. In the present study, nine candidate reference genes were rigorously tested for their application in normalization of qPCR data onto roots of four distinct cowpea accessions under two abiotic stresses: root dehydration and salt (NaCl, 100 mM). In addition, the regulation of four target transcripts, under the same referred conditions was also scrutinized.

RESULTS

geNorm, NormFinder, BestKeeper, and ΔCt method results indicated a set of three statistically validated RGs for each stress condition: (I) root dehydration (actin, ubiquitin-conjugating enzyme E2 variant 1D, and a Phaseolus vulgaris unknown gene-UNK), and (II) salt (ubiquitin-conjugating enzyme E2 variant 1D, F-box protein, and UNK). The expression profile of the target transcripts suggests that flavonoids are important players in the cowpea response to the abiotic stresses analyzed, since chalcone isomerase and chalcone synthase were up-regulated in the tolerant and sensitive accessions. A lipid transfer protein also participates in the cowpea tolerance mechanisms to root dehydration and salt stress. The referred transcript was up-regulated in the two tolerant accessions and presented no differential expression in the sensitive counterparts. Chitinase B, in turn, generally related to plant defense, was an important target transcript under salt stress, being up-regulated at the tolerant, and down-regulated in the sensitive accession.

CONCLUSIONS

Reference genes suitable for qPCR analyses in cowpea under root dehydration and salt stress were identified. This action will lead to a more accurate and reliable analysis of gene expression on this species. Additionally, the results obtained in this study may guide future research on gene expression in cowpea under other abiotic stress types that impose osmotic imbalance. The target genes analyzed, in turn, deserve functional evaluation due to their transcriptional regulation under stresses and biotechnological potential.

Identification of Endogenous Reference Genes for RT-qPCR Expression Analysis in Urochloa brizantha Under Abiotic Stresses

Urochloa brizantha is one of the most important warm season forage grasses in tropical countries. Despite its importance, there are few studies on gene expression in this species under stressful conditions. Real-time (RT-qPCR) is an accurate technique for gene quantification analysis, but reference genes must be validated under the same conditions used to assess the expression of the target genes. Here, we evaluated the stability of nine reference genes: Actin 12, Eukaryotic initiation factor 4 A, Elongation factor-1 alpha, FTSH protease 4, U2 auxiliary fator, Succinol Co-enzyme A, Tubulin alfa-5, Tubulin beta-6, Ubiquitin conjugating enzyme. Total RNA was extract from leaf tissues of U. brizantha subjected to 6, 12 and 24 h of cold and heat stresses (10 and 45 °C, respectively), and drought, including moderate (−0.5 to −0.7 MPa), severe (−1.1 to −1.8 MPa) and recovery after re-watering. The RefFinder web-based tool was used to rank the most stable reference genes for each stress. Elongation factor-1 alpha, Elongation factor-1 alpha or Ubiquitin conjugating enzyme, and Eukaryotic initiation factor 4 A were the most stable genes for heat, cold and drought stress, respectively. The expression of Rubisco large subunit gene was normalized against the most stable gene selected by ReFfinder for each stress.

Genome-wide identification, classification and transcriptional analysis of nitrate and ammonium transporters in Coffea

Nitrogen (N) is quantitatively the main nutrient required by coffee plants, with acquisition mainly by the roots and mostly exported to coffee beans. Nitrate (NO3-) and ammonium (NH4+) are the most important inorganic sources for N uptake. Several N transporters encoded by different gene families mediate the uptake of these compounds. They have an important role in source preference for N uptake in the root system. In this study, we performed a genome-wide analysis, including in silico expression and phylogenetic analyses of AMT1, AMT2, NRT1/PTR, and NRT2 transporters in the recently sequenced Coffea canephora genome. We analyzed the expression of six selected transporters in Coffea arabica roots submitted to N deficiency. N source preference was also analyzed in C. arabica using isotopes. C. canephora N transporters follow the patterns observed for most eudicots, where each member of the AMT and NRT families has a particular role in N mobilization, and where some of these are modulated by N deficiency. Despite the prevalence of putative nitrate transporters in the Coffea genome, ammonium was the preferential inorganic N source for N-starved C. arabica roots. This data provides an important basis for fundamental and applied studies to depict molecular mechanisms involved in N uptake in coffee trees.

Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis

Coffea arabica L. is an important crop in several developing countries. Despite its economic importance, minimal transcriptome data are available for fruit tissues, especially during fruit development where several compounds related to coffee quality are produced. To understand the molecular aspects related to coffee fruit and grain development, we report a large-scale transcriptome analysis of leaf, flower and perisperm fruit tissue development. Illumina sequencing yielded 41,881,572 high-quality filtered reads. De novo assembly generated 65,364 unigenes with an average length of 1,264 bp. A total of 24,548 unigenes were annotated as protein coding genes, including 12,560 full-length sequences. In the annotation process, we identified nine candidate genes related to the biosynthesis of raffinose family oligossacarides (RFOs). These sugars confer osmoprotection and are accumulated during initial fruit development. Four genes from this pathway had their transcriptional pattern validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Furthermore, we identified ~24,000 putative target sites for microRNAs (miRNAs) and 134 putative transcriptionally active transposable elements (TE) sequences in our dataset. This C. arabica transcriptomic atlas provides an important step for identifying candidate genes related to several coffee metabolic pathways, especially those related to fruit chemical composition and therefore beverage quality. Our results are the starting point for enhancing our knowledge about the coffee genes that are transcribed during the flowering and initial fruit development stages.

Validation of Suitable Reference Genes for RT-qPCR Data in Achyranthes bidentata Blume under Different Experimental Conditions

Real-time quantitative polymerase chain reaction (RT-qPCR) is a sensitive technique for gene expression studies. However, choosing the appropriate reference gene is essential to obtain reliable results for RT-qPCR assays. In the present work, the expression of eight candidate reference genes, EF1-α (elongation factor 1-α), GAPDH (glyceraldehyde 3-phosphate dehydrogenase), UBC (ubiquitin-conjugating enzyme), UBQ (polyubiquitin), ACT (actin), β-TUB (β-tubulin), APT1 (adenine phosphoribosyltransferase 1), and 18S rRNA (18S ribosomal RNA), was evaluated in Achyranthes bidentata samples using two algorithms, geNorm and NormFinder. The samples were classified into groups according to developmental stages, various tissues, stresses (cold, heat, drought, NaCl), and hormone treatments (MeJA, IBA, SA). Suitable combination of reference genes for RT-qPCR normalization should be applied according to different experimental conditions. In this study, EF1-α, UBC, and ACT genes were verified as the suitable reference genes across all tested samples. To validate the suitability of the reference genes, we evaluated the relative expression of CAS, which is a gene that may be involved in phytosterol synthesis. Our results provide the foundation for gene expression analysis in A. bidentata and other species of Amaranthaceae.

Selection and Validation of Reference Genes for Accurate RT-qPCR Data Normalization in Coffea spp. under a Climate Changes Context of Interacting Elevated [CO2] and Temperature

World coffee production has faced increasing challenges associated with ongoing climatic changes. Several studies, which have been almost exclusively based on temperature increase, have predicted extensive reductions (higher than half by 2,050) of actual coffee cropped areas. However, recent studies showed that elevated [CO₂] can strongly mitigate the negative impacts of heat stress at the physiological and biochemical levels in coffee leaves. In addition, it has also been shown that coffee genotypes can successfully cope with temperatures above what has been traditionally accepted. Altogether, this information suggests that the real impact of climate changes on coffee growth and production could be significantly lower than previously estimated. Gene expression studies are an important tool to unravel crop acclimation ability, demanding the use of adequate reference genes. We have examined the transcript stability of 10 candidate reference genes to normalize RT-qPCR expression studies using a set of 24 cDNAs from leaves of three coffee genotypes (CL153, Icatu, and IPR108), grown under 380 or 700 μL CO₂ L^-1, and submitted to increasing temperatures from 25/20°C (day/night) to 42/34°C. Samples were analyzed according to genotype, [CO₂], temperature, multiple stress interaction ([CO₂], temperature) and total stress interaction (genotype, [CO₂], and temperature). The transcript stability of each gene was assessed through a multiple analytical approach combining the Coeficient of Variation method and three algorithms (geNorm, BestKeeper, NormFinder). The transcript stability varied according to the type of stress for most genes, but the consensus ranking obtained with RefFinder, classified MDH as the gene with the highest mRNA stability to a global use, followed by ACT and S15, whereas α-TUB and CYCL showed the least stable mRNA contents. Using the coffee expression profiles of the gene encoding the large-subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (RLS), results from the in silico aggregation and experimental validation of the best number of reference genes showed that two reference genes are adequate to normalize RT-qPCR data. Altogether, this work highlights the importance of an adequate selection of reference genes for each single or combined experimental condition and constitutes the basis to accurately study molecular responses of Coffea spp. in a context of climate changes and global warming.

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.

Selection and Validation of Appropriate Reference Genes for Quantitative Real-Time PCR Analysis of Gene Expression in Lycoris aurea

Lycoris aurea (L' Hér.) Herb, a perennial grass species, produces a unique variety of pharmacologically active Amaryllidaceae alkaloids. However, the key enzymes and their expression pattern involved in the biosynthesis of Amaryllidaceae alkaloids (especially for galanthamine) are far from being fully understood. Quantitative real-time polymerase chain reaction (qRT-PCR), a commonly used method for quantifying gene expression, requires stable reference genes to normalize its data. In this study, to choose the appropriate reference genes under different experimental conditions, 14 genes including YLS8 (mitosis protein YLS8), CYP2 (Cyclophilin 2), CYP 1 (Cyclophilin 1), TIP41 (TIP41-like protein), EXP2 (Expressed protein 2), PTBP1 (Polypyrimidine tract-binding protein 1), EXP1 (Expressed protein 1), PP2A (Serine/threonine-protein phosphatase 2A), β-TUB (β-tubulin), α-TUB (α-tubulin), EF1-α (Elongation factor 1-α), UBC (Ubiquitin-conjugating enzyme), ACT (Actin) and GAPDH (Glyceraldehyde 3-phosphate dehydrogenase) were selected from the transcriptome datasets of L. aurea. And then, expressions of these genes were assessed by qRT-PCR in various tissues and the roots under different treatments. The expression stability of the 14 candidates was analyzed by three commonly used software programs (geNorm, NormFinder, and BestKeeper), and their results were further integrated into a comprehensive ranking based on the geometric mean. The results show the relatively stable genes for each subset as follows: (1) EXP1 and TIP41 for all samples; (2) UBC and EXP1 for NaCl stress; (3) PTBP1 and EXP1 for heat stress, polyethylene glycol (PEG) stress and ABA treatment; (4) UBC and CYP2 for cold stress; (5) PTBP1 and PP2A for sodium nitroprusside (SNP) treatment; (6) CYP1 and TIP41 for methyl jasmonate (MeJA) treatment; and (7) EXP1 and TIP41 for various tissues. The reliability of these results was further enhanced through comparison between part qRT-PCR result and RNA sequencing (RNA-seq) data. In summary, our results identified appropriate reference genes for qRT-PCR in L. aurea, and will facilitate gene expression studies under these conditions.

Transcriptome analysis in Coffea eugenioides, an Arabica coffee ancestor, reveals differentially expressed genes in leaves and fruits

Studies in diploid parental species of polyploid plants are important to understand their contributions to the formation of plant and species evolution. Coffea eugenioides is a diploid species that is considered to be an ancestor of allopolyploid Coffea arabica together with Coffea canephora. Despite its importance in the evolutionary history of the main economic species of coffee, no study has focused on C. eugenioides molecular genetics. RNA-seq creates the possibility to generate reference transcriptomes and identify coding genes and potential candidates related to important agronomic traits. Therefore, the main objectives were to obtain a global overview of transcriptionally active genes in this species using next-generation sequencing and to analyze specific genes that were highly expressed in leaves and fruits with potential exploratory characteristics for breeding and understanding the evolutionary biology of coffee. A de novo assembly generated 36,935 contigs that were annotated using eight databases. We observed a total of ~5000 differentially expressed genes between leaves and fruits. Several genes exclusively expressed in fruits did not exhibit similarities with sequences in any database. We selected ten differentially expressed unigenes in leaves and fruits to evaluate transcriptional profiles using qPCR. Our study provides the first gene catalog for C. eugenioides and enhances the knowledge concerning the mechanisms involved in the C. arabica homeologous. Furthermore, this work will open new avenues for studies into specific genes and pathways in this species, especially related to fruit, and our data have potential value in assisted breeding applications.

Biochemical traits and proteomic changes in postharvest flowers of medicinal chrysanthemum exposed to enhanced UV-B radiation

The article studied UV-B effects on biochemical traits and proteomic changes in postharvest flowers of medicinal chrysanthemum. The experiment about UV-B effects on biochemical traits in flowers included six levels of UV-B treatments (0 (UV0), 50 (UV50), 200 (UV200), 400 (UV400), 600 (UV600) and 800 (UV800) μWcm(-2)). UV400, UV600 and UV800 treatments significantly increased the contents of hydrogen peroxide, malondialdehyde and UV-B absorbing compounds, and the activity of phenylalanine ammonia lyase enzyme over the control. The contents of chlorogenic acid and flavone in flowers were significantly increased by UV-B treatments (except for UV50 and UV800). Two-dimensional gel electrophoresis was utilized to analyze proteomic changes in flowers with or without UV-B radiation. Results indicated that 43 protein spots (>1.5-fold difference in volume) were detected, including 19 spots with a decreasing trend and 24 spots with an increasing trend, and 19 differentially expressed protein spots were successfully indentified by MALDI-TOF MS. The indentified proteins were classified based on functions, the most of which were involved in photosynthesis, respiration, protein biosynthesis and degradation and defence. An overall assessment using biochemical and differential proteomic data revealed that UV-B radiation could affect biochemical reaction and promote secondary metabolism processes in postharvest flowers.

Galactinol synthase transcriptional profile in two genotypes of Coffea canephora with contrasting tolerance to drought

Increased synthesis of galactinol and raffinose family oligosaccharides (RFOs) has been reported in vegetative tissues in response to a range of abiotic stresses. In this work, we evaluated the transcriptional profile of a Coffea canephora galactinol synthase gene (CcGolS1) in two clones that differed in tolerance to water deficit in order to assess the contribution of this gene to drought tolerance. The expression of CcGolS1 in leaves was differentially regulated by water deficit, depending on the intensity of stress and the genotype. In clone 109A (drought-susceptible), the abundance of CcGolS1 transcripts decreased upon exposure to drought, reaching minimum values during recovery from severe water deficit and stress. In contrast, CcGolS1 gene expression in clone 14 (drought-tolerant) was stimulated by water deficit. Changes in galactinol and RFO content did not correlate with variation in the steady-state transcript level. However, the magnitude of increase in RFO accumulation was higher in the tolerant cultivar, mainly under severe water deficit. The finding that the drought-tolerant coffee clone showed enhanced accumulation of CcGolS1 transcripts and RFOs under water deficit suggests the possibility of using this gene to improve drought tolerance in this important crop.

Homeologous genes involved in mannitol synthesis reveal unequal contributions in response to abiotic stress in Coffea arabica

Polyploid plants can exhibit transcriptional modulation in homeologous genes in response to abiotic stresses. Coffea arabica, an allotetraploid, accounts for 75% of the world's coffee production. Extreme temperatures, salinity and drought limit crop productivity, which includes coffee plants. Mannitol is known to be involved in abiotic stress tolerance in higher plants. This study aimed to investigate the transcriptional responses of genes involved in mannitol biosynthesis and catabolism in C. arabica leaves under water deficit, salt stress and high temperature. Mannitol concentration was significantly increased in leaves of plants under drought and salinity, but reduced by heat stress. Fructose content followed the level of mannitol only in heat-stressed plants, suggesting the partitioning of the former into other metabolites during drought and salt stress conditions. Transcripts of the key enzymes involved in mannitol biosynthesis, CaM6PR, CaPMI and CaMTD, were modulated in distinct ways depending on the abiotic stress. Our data suggest that changes in mannitol accumulation during drought and salt stress in leaves of C. arabica are due, at least in part, to the increased expression of the key genes involved in mannitol biosynthesis. In addition, the homeologs of the Coffea canephora subgenome did not present the same pattern of overall transcriptional response, indicating differential regulation of these genes by the same stimulus. In this way, this study adds new information on the differential expression of C. arabica homeologous genes under adverse environmental conditions showing that abiotic stresses can influence the homeologous gene regulation pattern, in this case, mainly on those involved in mannitol pathway.

Reference gene selection for quantitative real-time RT-PCR normalization in Iris. lactea var. chinensis roots under cadmium, lead, and salt stress conditions

Quantitative real time PCR (RT-qPCR) has emerged as an accurate and sensitive method to measure the gene expression. However, obtaining reliable result depends on the selection of reference genes which normalize differences among samples. In this study, we assessed the expression stability of seven reference genes, namely, ubiquitin-protein ligase UBC9 (UBC), tubulin alpha-5 (TUBLIN), eukaryotic translation initiation factor (EIF-5A), translation elongation factor EF1A (EF1 α ), translation elongation factor EF1B (EF1b), actin11 (ACTIN), and histone H3 (HIS), in Iris. lactea var. chinensis (I. lactea var. chinensis) root when the plants were subjected to cadmium (Cd), lead (Pb), and salt stress conditions. All seven reference genes showed a relatively wide range of threshold cycles (C t ) values in different samples. GeNorm and NormFinder algorithms were used to assess the suitable reference genes. The results from the two software units showed that EIF-5A and UBC were the most stable reference genes across all of the tested samples, while TUBLIN was unsuitable as internal controls. I. lactea var. chinensis is tolerant to Cd, Pb, and salt. Our results will benefit future research on gene expression in response to the three abiotic stresses.

Selection of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in Brassica napus under various stress conditions

Data normalization is essential for reliable output of quantitative real-time reverse-transcription polymerase chain reaction (qRT-PCR) assays, as the unsuitable choice of reference gene(s), whose expression might be influenced by exogenous treatments in plant tissues, could cause misinterpretation of results. To date, no systematic studies on reference genes have been performed in stressed Brassica napus. In this study, we investigated the expression variations of nine candidate reference genes in 40 samples of B. napus leaves subjected to various exogenous treatments. Parallel analyses by geNorm and NormFinder revealed that optimal reference genes differed across the different sets of samples. The best-ranked reference genes were PP2A and TIP41 for salt stress, TIP41 and ACT7 for heavy metal (Cr(6+)) stress, PP2A and UBC21 for drought stress, F-box and SAND for cold stress, F-box and ZNF for salicylic acid stress, TIP41, ACT7, and PP2A for methyl jasmonate stress, TIP41 and ACT7 for abscisic acid stress, and TIP41, UBC21, and PP2A for Sclerotinia sclerotiorum stress. Two newly employed reference genes, TIP41 and PP2A, showed better performances, suggesting their suitability in multiple conditions. To further validate the suitability of the reference genes, the expression patterns of BnWRKY40 and BnMKS1 were studied in parallel. This study is the first systematic analysis of reference gene selection for qRT-PCR normalization in B. napus, an agriculturally important crop, under different stress conditions. The results will contribute toward more accurate and widespread use of qRT-PCR in gene analysis of the genus Brassica.

Selection of suitable reference genes for qPCR normalization under abiotic stresses in Oenanthe javanica (BI.) DC

Accurate normalization of gene expression data is an absolute prerequisite to obtain reliable results in qPCR analysis. Oenanthe javanica, an aquatic perennial herb, belongs to the Oenanthe genus in Apiaceae family, with known medicinal properties. In the current study, O. javanica was subjected to hormone stimuli (gibberellin, salicylic acid, methyl jasmonate, and abscisic acid) and abiotic stresses (heat, cold, salt, and drought), and the expression of nine candidate reference genes (eIF-4α, ACT7, TIP41, GAPDH, SAND, EF-1α, PP2A, TBP, and TUB) was evaluated. Stability of the genes was assessed using geNorm, NormFinder and BestKeeper. All the genes presented distinct expression profiles under the experimental conditions analyzed. Under abiotic stress conditions, ACT7 and PP2A genes displayed the maximum stability; PP2A and SAND were the most stable genes under hormone stimuli. Even though PP2A gene was most stable across all the samples, individual analysis revealed changes in expression profile. To further validate the suitability of the reference genes identified in this study, the expression level of M6PR gene under salt treatment was studied. Based on our data, we propose that it is essential to normalize the target gene expression with specific reference genes under different experimental conditions for most accurate results. To our knowledge, this is the first systematic analysis for reference genes under abiotic stress and hormone stimuli conditions in O. javanica. This will be beneficial for future studies on O. javanica and other plants in Apiaceae family at molecular level.

Validation of suitable reference genes for gene expression analysis in the halophyte Salicornia europaea by real-time quantitative PCR

Real-time quantitative polymerase chain reaction (RT-qPCR), a reliable technique for quantifying gene expression, requires stable reference genes to normalize its data. Salicornia europaea, a stem succulent halophyte with remarkable salt resistance and high capacity for ion accumulation, has not been investigated with regards to the selection of appropriate reference genes for RT-qPCR. In this study, the expression of 11 candidate reference genes, GAPDH (Glyceraldehyde 3-phosphate dehydrogenase), Actin, α-Tub (α-tubulin), β-Tub (β-tubulin), EF1-α (Elongation factor 1-α), UBC (Ubiquitin-conjugating enzyme), UBQ (Polyubiquitin), CYP (Cyclophilin), TIP41 (TIP41-like protein), CAC (Clathrin adaptor complexes), and DNAJ (DnaJ-like protein), was analyzed in S. europaea samples, which were classified into groups according to various abiotic stresses (NaCl, nitrogen, drought, cold and heat), tissues and ages. Three commonly used software programs (geNorm, NormFinder and BestKeeper) were applied to evaluate the stability of gene expression, and comprehensive ranks of stability were generated by aggregate analysis. The results show that the relatively stable genes for each group are the following: (1) CAC and UBC for whole samples; (2) CAC and UBC for NaCl stress; (3) Actin and α-Tub for nitrogen treatment; (4) Actin and GAPDH for drought stress; (5) α-Tub and UBC for cold stress; (6) TIP41 and DNAJ for heat stress; (7) UBC and UBQ for different tissues; and (8) UBC and Actin for various developmental stages. These genes were validated by comparing transcriptome profiles. Using two stable reference genes was recommended in the normalization of RT-qPCR data. This study identifies optimal reference genes for RT-qPCR in S. europaea, which will benefit gene expression analysis under these conditions.