Identification of two putative reference genes from grapevine suitable for gene expression analysis in berry and related tissues derived from RNA-Seq data

Systematic identification of reference genes for qRT-PCR of Ardisia kteniophylla A. DC under different experimental conditions and for anthocyanin-related genes studies

Ardisia kteniophylla A. DC, widely known as folk medicinal herb and ornamental plant, has been extensively investigated due to its unique leaf color, anti-cancer and other pharmacological activities. The quantitative real-time PCR (qRT-PCR) was an excellent tool for the analysis of gene expression with its high sensitivity and quantitative properties. Normalizing gene expression with stable reference genes was essential for qRT-PCR accuracy. In addition, no studies have yet been performed on the selection, verification and stability of internal reference genes suitable for A. kteniophylla, which has greatly hindered the biomolecular researches of this species. In this study, 29 candidate genes were successfully screened according to stable expression patterns of large-scale RNA seq data that from a variety of tissues and the roots of different growth stages in A. kteniophylla. The candidates were then further determined via qRT-PCR in various experimental samples, including MeJA, ABA, SA, NaCl, CuSO4, AgNO3, MnSO4, CoCl2, drought, low temperature, heat, waterlogging, wounding and oxidative stress. To assess the stability of the candidates, five commonly used strategies were employed: delta-CT, geNorm, BestKeeper, NormFinder, and the comprehensive tool RefFinder. In summary, UBC2 and UBA1 were found to be effective in accurately normalizing target gene expression in A. kteniophella regardless of experimental conditions, while PP2A-2 had the lowest stability. Additionally, to verify the reliability of the recommended reference genes under different colored leaf samples, we examined the expression patterns of six genes associated with anthocyanin synthesis and regulation. Our findings suggested that PAP1 and ANS3 may be involved in leaf color change in A. kteniphella. This study successfully identified the ideal reference gene for qRT-PCR analysis in A. kteniphella, providing a foundation for future research on gene function, particularly in the biosynthesis of anthocyanins.

Metabolic and Molecular Rearrangements of Sauvignon Blanc (Vitis vinifera L.) Berries in Response to Foliar Applications of Specific Dry Yeast

Dry yeast extracts (DYE) are applied to vineyards to improve aromatic and secondary metabolic compound content and wine quality; however, systematic information on the underpinning molecular mechanisms is lacking. This work aimed to unravel, through a systematic approach, the metabolic and molecular responses of Sauvignon Blanc berries to DYE treatments. To accomplish this, DYE spraying was performed in a commercial vineyard for two consecutive years. Berries were sampled at several time points after the treatment, and grapes were analyzed for sugars, acidity, free and bound aroma precursors, amino acids, and targeted and untargeted RNA-Seq transcriptional profiles. The results obtained indicated that the DYE treatment did not interfere with the technological ripening parameters of sugars and acidity. Some aroma precursors, including cys-3MH and GSH-3MH, responsible for the typical aromatic nuances of Sauvignon Blanc, were stimulated by the treatment during both vintages. The levels of amino acids and the global RNA-seq transcriptional profiles indicated that DYE spraying upregulated ROS homeostatic and thermotolerance genes, as well as ethylene and jasmonic acid biosynthetic genes, and activated abiotic and biotic stress responses. Overall, the data suggested that the DYE reduced berry oxidative stress through the regulation of specific subsets of metabolic and hormonal pathways.

Phenolic, Polysaccharides Composition, and Texture Properties during Ripening and Storage Time of New Table Grape Cultivars in Chile

The aim of this study is to determine the phenolic and polysaccharidic composition, texture properties, and gene expression of new seedless table grape cultivars Timco™ and Krissy™ and compare them to the traditional table grape variety Crimson Seedless (Vitis vinifera L.), during ripening and in commercial postharvest conditions. According to the results, phenolic compounds were present in very different proportions. The total anthocyanins responsible for skin color increased during maturation and the majority anthocyanin in the three cultivars was peonidin-3-glucoside, followed by malvidin-3-glucoside. The phenolic compounds presented a different behavior (decreasing or increasing) during postharvest. The total skin soluble polysaccharides decreased during ripening and postharvest in Crimson Seedless and Krissy™ and remained constant from technological maturity to postharvest storage in Timco™. In all cultivars, the majority soluble polysaccharide fraction was that with a molecular mass between 500 and 35 KDa. The skin mechanical properties of table grapes were good parameters for differentiating varieties, with better results for the new cultivars, compared to the traditional Crimson Seedless, especially in postharvest. Genes involved in the flavonoid pathway and cell wall metabolism in skins exhibited an increase in expression from veraison to remaining constant at the end of the berry ripening.

Selection and validation of novel stable reference genes for qPCR analysis in EMT and MET

Quantitative real-time polymerase chain reaction is a powerful tool for quantifying gene expression. The relative quantification relies on normalizing the data to reference genes or internal controls not modulated by the experimental conditions. The most widely used internal controls occasionally show changed expression patterns in different experimental settings, such as the mesenchymal to epithelial transition. Thus, identifying appropriate internal controls is of utmost importance. We analyzed multiple RNA-Seq datasets using a combination of statistical approaches such as percent relative range and coefficient of variance to define a list of candidate internal control genes, which was then validated experimentally and by using in silico analyses as well. We identified a group of genes as strong internal control candidates with high stability compared to the classical ones. We also presented evidence for the superiority of the percent relative range method for calculating expression stability in data sets with larger sample sizes. We used multiple methods to analyze data collected from several RNA-Seq datasets; we identified Rbm17 and Katna1 as the most stable reference genes in EMT/MET studies. The percent relative range approach surpasses other methods when analyzing datasets of larger sample sizes.

Transcriptomic Approach for Global Distribution of SNP/Indel and Plant Genotyping

Single Nucleotide Polymorphisms (SNPs) are the most common structural variants found in any genome. They have been used for different genetic studies, from the understanding of genetic structure of populations to the development of breeding selection markers. In this chapter we present the use of transcriptomic data obtained from contrasting phenotypes for a target trait, in searching of SNPs and insertions/deletions (InDels). This approach has the advantage that the identified markers are in or close to differentially expressed genes, and so they have higher chances to tag the genes underlying the phenotypic expression of a particular trait.

Validation of Reference Genes in a Population of Blueberry (Vaccinium corymbosum) Plants Regenerated in Colchicine

For the first time we report the validation of reference genes in plants from a population of blueberry (Vaccinium corymbosum) clones cultured in vitro on a colchicine-supplemented medium. Nodal segment explants of the cultivar Duke were regenerated by organogenesis under different periods of colchicine 1 mg/L exposure (1, 2, 3, 5, 30 days). The clones selected for the study showed variability for phenotypic traits after 2 years of adaptation to field conditions, compared to plants of the donor genotype that were regenerated on a medium without colchicine. Vaccinium myrtillus (GAPDH) and Vaccinium macrocarpon (ATP1, NADH, RPOB and COX2) were used as reference genomes for primer design. The results show that colchicine treatments can cause genomic changes in blueberry plants. At the molecular level, exposure of plants to colchicine in early periods could promote an increase in gene expression of specific genes such as ATP1, COX2, GAPDH, MATK, NADH and RPOB. However, prolonged exposure (30 days) could decrease gene expression of the genes studied. For qPCR assays, the primers designed for ATP1, COX2, GAPDH and MATK genes showed high efficiency. In addition, the GAPDH, ATP1, NADH and COX2 genes showed high stability and could be recommended as potential reference genes for gene expression assays in Vaccinium.

Investigation of chicken housekeeping genes using next-generation sequencing data

Accurate normalization of the gene expression assays, using housekeeping genes (HKGs), is critically necessary. To do so, selection of a proper set of HKGs for a specific experiment is of great importance. Despite many studies, there is no consensus about the suitable set of HKGs for implementing in the quantitative real-time PCR analyses of chicken tissues. A limited number of HKGs have been widely used. However, wide utilization of a little number of HKGs for all tissues is challenging. The emergence of high-throughput gene expression RNA-seq data has enabled the simultaneous comparison of the stability of multiple HKGs. Therefore, employing the average coefficient of variations of at least three datasets per tissue, we sorted all reliably expressed genes (REGs; with FPKM ≥ 1 in at least one sample) and introduced the top 10 most suitable and stable reference genes for each of the 16 chicken tissues. We evaluated the consistency of the results of five tissues using the same methodology on other datasets. Furthermore, we assessed 96 previously widely used HKGs (WU-HKGs) in order to challenge the accuracy of the previous studies. The New Tuxedo software suite was used for the main analyses. The results revealed novel, different sets of reference genes for each of the tissues with 17 common genes among the top 10 genes lists of 16 tissues. The results did disprove the suitability of WU-HKGs such as Actb, Ldha, Scd, B2m, and Hprt1 for any of the tissues examined. On the contrary, a total of 6, 13, 14, 23, and 32 validated housekeeping genes (V-HKGs) were discovered as the most stable and suitable reference genes for muscle, spleen, liver, heart, and kidney tissues, respectively. Although we identified a few new HKGs usable for multiple tissues, the selection of suitable HKGs is required to be tissue specific. The newly introduced reference genes from the present study, despite lacking experimental validation, will be able to contribute to the more accurate normalization for future expression analysis of chicken genes.

Identification of Reference Genes for Reverse Transcription-Quantitative PCR Analysis of Ginger Under Abiotic Stress and for Postharvest Biology Studies

Gene expression analysis largely improves our understanding of the molecular basis underpinning various plant biological processes. Stable reference genes play a foundational role during the normalization of gene expression levels. However, until now, there have been few reference genes suitable for ginger reverse transcription-quantitative PCR (RT-qPCR) research. In this study, 29 candidate reference genes with stable expression patterns across multiple ginger tissues and 13 commonly used reference genes were selected to design RT-qPCR primers. After amplification specificity validation, 32 candidates were selected and further evaluated by RT-qPCR using samples from various organs subjected to NaCl, drought, heat, waterlogging, and chilling stress. Four strategies, including delta-CT, BestKeeper, geNorm, and NormFinder, were used to rank the stability of reference genes, and the ranks produced by these four strategies were comprehensively evaluated by RefFinder to determine the final rank. Overall, the top three stability reference genes indicated by RefFinder were RBP > ATPase > 40S_S3. Their expression pattern correlation analysis showed that the coefficients among each pair of RBP, ATPase, and 40S_S3 were larger than 0.96, revealing consistent and stable expression patterns under various treatments. Then, the expression of three pathogenesis-related (PR) genes and seven MYB genes in rhizomes during postharvest storage and subjected to pathogen infection was normalized by RBP, ATPase, 40S_S3, RBP and ATPase, ATPase and 40S-S3, and RBP and 40S-S3. The results showed that PR and MYB genes were induced by postharvest deterioration and pathogen infection. The correlation coefficients of RBP/ATPase, RBP/40S_S3, ATPase/40S_S3, RBP and ATPase/ATPase and 40S-S3, RBP and ATPase/RBP and 40S-S3, and ATPase and 40S-S3/RBP and 40S-S3 were 0.99, 0.96, 0.99, 0.99, 1.00, and 1.00, respectively, which confirmed the stability of these three reference genes in postharvest biology studies of ginger. In summary, this study identified appropriate reference genes for RT-qPCR in ginger and facilitated gene expression studies under biotic and abiotic stress conditions.

Genome-wide screening of novel RT-qPCR reference genes for study of GLRaV-3 infection in wine grapes and refinement of an RNA isolation protocol for grape berries

BACKGROUND

Grapevine, as an essential fruit crop with high economic values, has been the focus of molecular studies in diverse areas. Two challenges exist in the grapevine research field: (i) the lack of a rapid, user-friendly and effective RNA isolation protocol for mature dark-skinned berries and, (ii) the lack of validated reference genes that are stable for quantification of gene expression across desired experimental conditions. Successful isolation of RNA with sufficient yield and quality is essential for downstream analyses involving nucleic acids. However, ripe berries of dark-skinned grape cultivars are notoriously challenging in RNA isolation due to high contents of polyphenolics, polysaccharides, RNase and water.

RESULTS

We have optimized an RNA isolation protocol through modulating two factors at the lysis step that could impact results of RNA isolation - 2-ME concentration and berry mass. By finding the optimal combination among the two factors, our refined protocol was highly effective in isolating total RNA with high yield and quality from whole mature berries of an array of dark-skinned wine grape cultivars. Our protocol takes a much shorter time to complete, is highly effective, and eliminates the requirement for hazardous organic solvents. We have also shown that the resulting RNA preps were suitable for multiple downstream analyses, including the detection of viruses and amplification of grapevine genes using reverse transcription-polymerase chain reaction (RT-PCR), gene expression analysis via quantitative reverse transcription PCR (RT-qPCR), and RNA Sequencing (RNA-Seq). By using RNA-Seq data derived from Cabernet Franc, we have identified seven novel reference gene candidates (CYSP, NDUFS8, YLS8, EIF5A2, Gluc, GDT1, and EF-Hand) with stable expression across two tissue types, three developmental stages and status of infection with grapevine leafroll-associated virus 3 (GLRaV-3). We evaluated the stability of these candidate genes together with two conventional reference genes (actin and NAD5) using geNorm, NormFinder and BestKeeper. We found that the novel reference gene candidates outperformed both actin and NAD5. The three most stable reference genes were CYSP, NDUFS8 and YSL8, whereas actin and NAD5 were among the least stable. We further tested if there would be a difference in RT-qPCR quantification results when the most stable (CYSP) and the least stable (actin and NAD5) genes were used for normalization. We concluded that both actin and NAD5 led to erroneous RT-qPCR results in determining the statistical significance and fold-change values of gene expressional change.

CONCLUSIONS

We have formulated a rapid, safe and highly effective protocol for isolating RNA from recalcitrant berry tissue of wine grapes. The resulting RNA is of high quality and suitable for RT-qPCR and RNA-Seq. We have identified and validated a set of novel reference genes based on RNA-Seq dataset. We have shown that these new reference genes are superior over actin and NAD5, two of the conventional reference genes commonly used in early studies.

The grapevine (Vitis vinifera L.) floral transcriptome in Pinot noir variety: identification of tissue-related gene networks and whorl-specific markers in pre- and post-anthesis phases

The comprehension of molecular processes underlying the development and progression of flowering in plants is a hot topic, not only because that often the products of interest for human and animal nutrition are linked to the development of fruits or seeds, but also because the processes of gametes formation occurring in sexual organs are at the basis of recombination and genetic variability which constitutes the matter on which evolution acts, whether understood as natural or human driven. In the present study, we used an NGS approach to produce a grapevine flower transcriptome snapshot in different whorls and tissues including calyx, calyptra, filament, anther, stigma, ovary, and embryo in both pre- and post-anthesis phases. Our investigation aimed at identifying hub genes that unequivocally distinguish the different tissues providing insights into the molecular mechanisms that are at the basis of floral whorls and tissue development. To this end we have used different analytical approaches, some now consolidated in transcriptomic studies on plants, such as pairwise comparison and weighted-gene coexpression network analysis, others used mainly in studies on animals or human's genomics, such as the tau (τ) analysis aimed at isolating highly and absolutely tissue-specific genes. The intersection of data obtained by these analyses allowed us to gradually narrow the field, providing evidence about the molecular mechanisms occurring in those whorls directly involved in reproductive processes, such as anther and stigma, and giving insights into the role of other whorls not directly related to reproduction, such as calyptra and calyx. We believe this work could represent an important genomic resource for functional analyses of grapevine floral organ growth and fruit development shading light on molecular networks underlying grapevine reproductive organ determination.

Probing into the Effects of Grapevine Leafroll-Associated Viruses on the Physiology, Fruit Quality and Gene Expression of Grapes

Grapevine leafroll is one of the most widespread and highly destructive grapevine diseases that is responsible for great economic losses to the grape and wine industries throughout the world. Six distinct viruses have been implicated in this disease complex. They belong to three genera, all in the family Closteroviridae. For the sake of convenience, these viruses are named as grapevine leafroll-associated viruses (GLRaV-1, -2, -3, -4, -7, and -13). However, their etiological role in the disease has yet to be established. Furthermore, how infections with each GLRaV induce the characteristic disease symptoms remains unresolved. Here, we first provide a brief overview on each of these GLRaVs with a focus on genome structure, expression strategies and gene functions, where available. We then provide a review on the effects of GLRaV infection on the physiology, fruit quality, fruit chemical composition, and gene expression of grapevine based on the limited information so far reported in the literature. We outline key methodologies that have been used to study how GLRaV infections alter gene expression in the grapevine host at the transcriptomic level. Finally, we present a working model as an initial attempt to explain how infections with GLRaVs lead to the characteristic symptoms of grapevine leafroll disease: leaf discoloration and downward rolling. It is our hope that this review will serve as a starting point for grapevine virology and the related research community to tackle this vastly important and yet virtually uncharted territory in virus-host interactions involving woody and perennial fruit crops.

Identification of SNPs and InDels associated with berry size in table grapes integrating genetic and transcriptomic approaches

BACKGROUND

Berry size is considered as one of the main selection criteria in table grapes breeding programs, due to the consumer preferences. However, berry size is a complex quantitive trait under polygenic control, and its genetic determination of berry weight is not yet fully understood. The aim of this work was to perform marker discovery using a transcriptomic approach, in order to identify and characterize SNP and InDel markers associated with berry size in table grapes. We used an integrative analysis based on RNA-Seq, SNP/InDel search and validation on table grape segregants and varieties with different genetic backgrounds.

RESULTS

Thirty SNPs and eight InDels were identified using a transcriptomic approach (RNA-Seq). These markers were selected from SNP/InDel found among segregants from a Ruby x Sultanina population with contrasting phenotypes for berry size. The set of 38 SNP and InDel markers was distributed in eight chromosomes. Genotype-phenotype association analyses were performed using a set of 13 RxS segregants and 41 table grapes varieties with different genetic backgrounds during three seasons. The results showed several degrees of association of these markers with berry size (10.2 to 30.7%) as other berry-related traits such as length and width. The co-localization of SNP and /or InDel markers and previously reported QTLs and candidate genes associated with berry size were analysed.

CONCLUSIONS

We identified a set of informative and transferable SNP and InDel markers associated with berry size. Our results suggest the suitability of SNPs and InDels as candidate markers for berry weight in seedless table grape breeding. The identification of genomic regions associated with berry weight in chromosomes 8, 15 and 17 was achieved with supporting evidence derived from a transcriptome experiment focused on SNP/InDel search, as well as from a QTL-linkage mapping approach. New regions possibly associated with berry weight in chromosomes 3, 6, 9 and 14 were identified.

Identification of Reference Genes in Chicken Intraepithelial Lymphocyte Natural Killer Cells Infected with Very-virulent Infectious Bursal Disease Virus

Due to the limitations in the range of antibodies recognising avian viruses, quantitative real-time PCR (RT-qPCR) is still the most widely used method to evaluate the expression of immunologically related genes in avian viruses. The objective of this study was to identify suitable reference genes for mRNA expression analysis in chicken intraepithelial lymphocyte natural killer (IEL-NK) cells after infection with very-virulent infectious bursal disease virus (vvIBDV). Fifteen potential reference genes were selected based on the references available. The coefficient of variation percentage (CV%) and average count of these 15 genes were determined by NanoString technology for control and infected samples. The M and V values for shortlisted reference genes (ACTB, GAPDH, HMBS, HPRT1, SDHA, TUBB1 and YWHAZ) were calculated using geNorm and NormFinder. GAPDH, YWHAZ and HMBS were the most stably expressed genes. The expression levels of three innate immune response related target genes, CASP8, IL22 and TLR3, agreed in the NanoString and RNA sequencing (RNA-Seq) results using one or two reference genes for normalisation (not HMBS). In conclusion, GAPDH and YWHAZ could be used as reference genes for the normalisation of chicken IEL-NK cell gene responses to infection with vvIBDV.

Type II Ice-Binding Proteins Isolated from an Arctic Microalga Are Similar to Adhesin-Like Proteins and Increase Freezing Tolerance in Transgenic Plants

Microalgal ice-binding proteins (IBPs) in the polar region are poorly understood at the genome-wide level, although they are important for cold adaptation. Through the transcriptome study with the Arctic green alga Chloromonas sp. KNF0032, we identified six Chloromonas IBP genes (CmIBPs), homologous with the previously reported IBPs from Antarctic snow alga CCMP681 and Antarctic Chloromonas sp. They were organized with multiple exon/intron structures and low-temperature-responsive cis-elements in their promoters and abundantly expressed at low temperature. The biological functions of three representative CmIBPs (CmIBP1, CmIBP2 and CmIBP3) were tested using in vitro analysis and transgenic plant system. CmIBP1 had the most effective ice recrystallization inhibition (IRI) activities in both in vitro and transgenic plants, and CmIBP2 and CmIBP3 had followed. All transgenic plants grown under nonacclimated condition were freezing tolerant, and especially 35S::CmIBP1 plants were most effective. After cold acclimation, only 35S::CmIBP2 plants showed slightly increased freezing tolerance. Structurally, the CmIBPs were predicted to have β-solenoid forms with parallel β-sheets and repeated TXT motifs. The repeated TXT structure of CmIBPs appears similar to the AidA domain-containing adhesin-like proteins from methanogens. We have shown that the AidA domain has IRI activity as CmIBPs and phylogenetic analysis also supported that the AidA domains are monophyletic with ice-binding domain of CmIBPs, and these results suggest that CmIBPs are a type of modified adhesins.

Genome-Wide Identification and Evaluation of New Reference Genes for Gene Expression Analysis Under Temperature and Salinity Stresses in Ciona savignyi

Rapid adaptation/accommodation to changing environments largely contributes to maximal survival of invaders during biological invasions, usually leading to success in crossing multiple barriers and finally in varied environments in recipient habitats. Gene expression is one of the most important and rapid ways during responses to environmental stresses. Selection of proper reference genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). Here we identified eight candidate novel reference genes from the RNA-Seq data in an invasive model ascidian Ciona savignyi under temperature and salinity stresses. Subsequently, the expression stability of these eight novel reference genes, as well as other six traditionally used reference genes, was evaluated using RT-qPCR and comprehensive tool RefFinder. Under the temperature stress, two traditional reference genes, ribosomal proteins S15 and L17 (RPS15, RPL17), and one novel gene Ras homolog A (RhoA), were recommended as the top three stable genes, which can be used to normalize target genes with a high and moderate expression level, respectively. Under the salinity stress, transmembrane 9 superfamily member (TMN), MOB kinase activator 1A-like gene (MOB) and ubiquitin-conjugating enzyme (UBQ2) were suggested as the top three stable genes. On the other hand, several commonly used reference genes such as α-tubulin (TubA), β-tubulin (TubB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) showed unstable expressions, thus these genes should not be used as internal controls for gene expression analysis. We also tested the expression level of an important stress response gene, large proline-rich protein bag6-like gene (BAG) using different reference genes. As expected, we observed different results and conclusions when using different normalization methods, thus suggesting the importance of selection of proper reference genes and associated normalization methods. Our results provide a valuable reference gene resource for the normalization of gene expression in the study of environmental adaptation/accommodation during biological invasions using C. savignyi as a model.

Transcriptome-based identification and validation of reference genes for plant-bacteria interaction studies using Nicotiana benthamiana

RT-qPCR is a widely used technique for the analysis of gene expression. Accurate estimation of transcript abundance relies strongly on a normalization that requires the use of reference genes that are stably expressed in the conditions analyzed. Initially, they were adopted from those used in Northern blot experiments, but an increasing number of publications highlight the need to find and validate alternative reference genes for the particular system under study. The development of high-throughput sequencing techniques has facilitated the identification of such stably expressed genes. Nicotiana benthamiana has been extensively used as a model in the plant research field. In spite of this, there is scarce information regarding suitable RT-qPCR reference genes for this species. Employing RNA-seq data previously generated from tomato plants, combined with newly generated data from N. benthamiana leaves infiltrated with Pseudomonas fluorescens, we identified and tested a set of 9 candidate reference genes. Using three different algorithms, we found that NbUbe35, NbNQO and NbErpA exhibit less variable gene expression in our pathosystem than previously used genes. Furthermore, the combined use of the first two is sufficient for robust gene expression analysis. We encourage employing these novel reference genes in future RT-qPCR experiments involving N. benthamiana and Pseudomonas spp.

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

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

Transcriptome analysis of grapevine under salinity and identification of key genes responsible for salt tolerance

The negative effects of soil salinity towards grape yield depend upon salt concentration, cultivar type, developmental stage, and rootstock. Thompson Seedless variety of grape plant is considered moderately sensitive to salinity when grown upon its own root stock. In recent epoch, identification of key genes responsive to salinity offers hope to generate salinity-tolerant crop plants by their overexpression through genetic manipulation. In the present report, salt responsive transcriptome analysis of Thompson Seedless grape variety was done to identify vital genes involved in salinity tolerance which could be used further to generate salt liberal grape plant or other crop plants. Transcriptome libraries for control and 150-mM-NaCl-treated grape leaves were sequenced on Illumina platform where 714 genes were found to be differentially expressed. Gene ontology analysis indicated that under salinity conditions, the genes involved in metabolic process were highly enriched. Keto Encyclopedia of Genes and Genomes analysis revealed that, among the top 22 enriched pathways for the salt stress upregulated genes, the carbohydrate metabolism, signal transduction, energy metabolism, amino acid metabolism, biosynthesis of secondary metabolite, and lipid metabolism pathways possessed the largest number of transcripts. Key salinity-induced genes were selected and validated through qRT-PCR analysis which was comparable to RNA-seq results. Real-time PCR analysis also revealed that after 24 days of salinity, the expression of most of the selected key genes was highest. These salinity-induced genes will be characterized further in a model plant and also in Vitis vinifera through transgenic approach to disclose their role towards salt tolerance.

Identification of reference genes for RT-qPCR in the Antarctic moss Sanionia uncinata under abiotic stress conditions

Sanionia uncinata is a dominant moss species in the maritime Antarctic. Due to its high adaptability to harsh environments, this extremophile plant has been considered a good target for studying the molecular adaptation mechanisms of plants to a variety of environmental stresses. Despite the importance of S. uncinata as a representative Antarctic plant species for the identification and characterization of genes associated with abiotic stress tolerance, suitable reference genes, which are critical for RT-qPCR analyses, have not yet been identified. In this report, 11 traditionally used and 6 novel candidate reference genes were selected from transcriptome data of S. uncinata and the expression stability of these genes was evaluated under various abiotic stress conditions using three statistical algorithms; geNorm, NormFinder, and BestKeeper. The stability ranking analysis selected the best reference genes depending on the stress conditions. Among the 17 candidates, the most stable references were POB1 and UFD2 for cold stress, POB1 and AKB for drought treatment, and UFD2 and AKB for the field samples from a different water contents in Antarctica. Overall, novel genes POB1 and AKB were the most reliable references across all samples, irrespective of experimental conditions. In addition, 6 novel candidate genes including AKB, POB1 and UFD2, were more stable than the housekeeping genes traditionally used for internal controls, indicating that transcriptome data can be useful for identifying novel robust normalizers. The reference genes validated in this study will be useful for improving the accuracy of RT-qPCR analysis for gene expression studies of S. uncinata in Antarctica and for further functional genomic analysis of bryophytes.

Exogenous gibberellic acid application induces the overexpression of key genes for pedicel lignification and an increase in berry drop in table grape

Most table grape (Vitis vinifera L.) varieties require gibberellic acid (GA₃) applications to obtain an adequate berry size in order to satisfy market requirements. However, GA₃ treatments also produce severe berry drop in some cultivars, which occurs mainly after a cold storage period during post-harvest. Berry drop in bunches treated with GA₃ has been related to the hardening and thickening of the pedicel produced by the over-accumulation of cellulose and its lignification. The main goal of this study was to compare the morphology and gene expression in pedicel samples of genotypes contrasting for berry drop susceptibility. These genotypes are Thompson Seedless, which exhibits a low incidence of berry drop, and a genetic line (Line #23) of INIA's breeding program that is very susceptible to berry drop at harvest and after storage in bunches sprayed with GA₃. The parameters measured to study this phenomenon during fruit growth and post-harvest storage included fruit detachment force (FDF), hardness and thickness of the pedicel and berry drop frequency. Histological analyses of pedicel structures at harvest showed an increase in cell size and deposition of lignin in the cortex zone in both contrasting genotypes treated with GA₃. The expression profile in both genotypes of the key lignin biosynthesis genes Vv4CL4, VvCCR1L and VvCAD1 analyzed by quantitative real time PCR (qPCR) revealed evident changes in response to GA₃ treatments. In particular, gene VvCAD1 is overexpressed (100X) in pedicels of line #23 treated with GA₃ after 30 and 45 days in cold storage compared to control. Moreover, the frequency of berry drop was higher for Line #23 treated with GA₃ than for the control (23% vs. 1%). Our results suggest that gibberellic acid regulates the expression of the biosynthesis of lignin genes, generating changes in cell wall composition and pedicel structure that result in an increase in berry drop.

Plant reference genes for development and stress response studies

Many reference genes are used by different laboratories for gene expression analyses to indicate the relative amount of input RNA/DNA in the experiment. These reference genes are supposed to show least variation among the treatments and with the control sets in a given experiment. However, expression of reference genes varies significantly from one set of experiment to the other. Thus, selection of reference genes depends on the experimental conditions. Sometimes the average expression of two or three reference genes is taken as standard. This review consolidated the details of about 120 genes attempted for normalization during comparative expression analysis in 16 different plants. Plant species included in this review are Arabidopsis thaliana, cotton (Gossypium hirsutum), tobacco (Nicotiana benthamiana and N. tabacum), soybean (Glycine max), rice (Oryza sativa), blueberry (Vaccinium corymbosum), tomato (Solanum lycopersicum), wheat (Triticum aestivum), potato (Solanum tuberosum), sugar cane (Saccharum sp.), carrot (Daucus carota), coffee (Coffea arabica), cucumber (Cucumis sativus), kiwi (Actinidia deliciosa) and grape (Vitis vinifera). The list includes model and cultivated crop plants from both monocot and dicot classes. We have categorized plant-wise the reference genes that have been used for expression analyses in any or all of the four different conditions such as biotic stress, abiotic stress, developmental stages and various organs and tissues, reported till date. This review serves as a guide during the reference gene hunt for gene expression analysis studies.

Selection and evaluation of reference genes for qRT-PCR analysis in Euscaphis konishii Hayata based on transcriptome data

BACKGROUND

Quantitative real-time reverse transcription-polymerase chain reaction has been widely used in gene expression analysis, however, to have reliable and accurate results, reference genes are necessary to normalize gene expression under different experimental conditions. Several reliable reference genes have been reported in plants of Traditional Chinese Medicine, but none have been identified for Euscaphis konishii Hayata.

RESULTS

In this study, 12 candidate reference genes, including 3 common housekeeping genes and 9 novel genes based on E. konishii Hayata transcriptome data were selected and analyzed in different tissues (root, branch, leaf, capsule and seed), capsule and seed development stages. Expression stability was calculated using geNorm and NormFinder, the minimal number of reference genes required for accurate normalization was calculated by Vn/Vn + 1 using geNorm. EkEEF-5A-1 and EkADF2 were the two most stable reference genes for all samples, while EkGSTU1 and EkGAPDH were the most stable reference genes for tissue samples. For seed development stages, EkGAPDH and EkEEF-5A-1 were the most stable genes, whereas EkGSTU1 and EkGAPDH were identified as the two most stable genes in the capsule development stages. Two reference genes were sufficient to normalize gene expression across all sample sets.

CONCLUSION

Results of this study revealed that suitable reference genes should be selected for different experimental samples, and not all the common reference genes are suitable for different tissue samples and/or experimental conditions. In this study, we present the first data of reference genes selection for E. konishii Hayata based on transcriptome data, our data will facilitate further studies in molecular biology and gene function on E. konishii Hayata and other closely related species.

Relationship among color development, anthocyanin and pigment-related gene expression in 'Crimson Seedless' grapes treated with abscisic acid and sucrose

'Crimson Seedless' is one of the most important table grape varieties in Chile, but under certain environmental conditions, the fruit exhibits inadequate red color development, causing economic losses due to lower product quality. The use of plant growth regulators, such as abscisic acid (ABA) and ethylene, during development increases the anthocyanin content of the skin, improving the color of the berry. Recently, sucrose has been identified as a signaling molecule capable of regulating the expression of genes of the anthocyanin biosynthesis pathway. The aim of this study was to analyze the effect of application of ABA and/or sucrose on color development and their relationship with anthocyanin metabolism. Applications of ABA (400 ppm or 200 ppm) and/or sucrose (90 mM) were performed close to the véraison stage. During development and at harvest, quality attributes such as berry firmness, total soluble solids and titratable acidity were not affected by these treatments. Increased red color development was observed in fruits treated with ABA and/or sucrose, due to accumulation of anthocyanins. Fruits subjected to sucrose treatment showed higher levels of anthocyanins than untreated fruits but lower levels than fruits treated with ABA. Increased expression of genes involved in anthocyanin biosynthesis was observed in ABA- and sucrose-treated fruits compared to untreated fruits. Based on these findings, we demonstrated that sucrose improved fruit color development by increasing synthesis and accumulation of anthocyanins, thus allowing earlier harvests and improving table grape quality.

Use of RNA-seq data to identify and validate RT-qPCR reference genes for studying the tomato-Pseudomonas pathosystem

The agronomical relevant tomato-Pseudomonas syringae pv. tomato pathosystem is widely used to explore and understand the underlying mechanisms of the plant immune response. Transcript abundance estimation, mainly through reverse transcription-quantitative PCR (RT-qPCR), is a common approach employed to investigate the possible role of a candidate gene in certain biological process under study. The accuracy of this technique relies heavily on the selection of adequate reference genes. Initially, genes derived from other techniques (such as Northern blots) were used as reference genes in RT-qPCR experiments, but recent studies in different systems suggest that many of these genes are not stably expressed. The development of high throughput transcriptomic techniques, such as RNA-seq, provides an opportunity for the identification of transcriptionally stable genes that can be adopted as novel and robust reference genes. Here we take advantage of a large set of RNA-seq data originating from tomato leaves infiltrated with different immunity inducers and bacterial strains. We assessed and validated 9 genes that are much more stable than two traditional reference genes. Specifically, ARD2 and VIN3 were the most stably expressed genes and consequently we propose they be adopted for RT-qPCR experiments involving this pathosystem.

Selection and Validation of Appropriate Reference Genes for Quantitative Real-Time PCR Normalization in Staminate and Perfect Flowers of Andromonoecious Taihangia rupestris

Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) is the most commonly used and powerful method for gene expression analysis due to its high sensitivity, specificity, and high throughput, and the accuracy of this approach depends on the stability of reference genes used for normalization. Taihangia rupestris Yu and Li (Rosaceae), an andromonoecious plant, produces both bisexual flowers and unisexual male flowers within the same individual. Using qRT-PCR technique, investigation of the gene expression profiling in staminate and perfect flowers would improve our understanding of the molecular mechanism in regulation of flower formation and sex differentiation in andromonoecious T. rupestris. To accurate normalize the gene expression level in Taihangia flower, 16 candidate reference genes, including 10 traditional housekeeping genes, and 6 newly stable genes, were selected based on transcriptome sequence data and previous studies. The expressions of these genes were assessed by qRT-PCR analysis in 51 samples, including 30 staminate and perfect flower samples across developmental stages and 21 different floral tissue samples from mature flowers. By using geNorm, NormFinder, BestKeeper, and comprehensive RefFinder algorithms, ADF3 combined with UFD1 were identified as the optimal reference genes for staminate flowers, while the combination of HIS3/ADF3 was the most accurate reference genes for perfect floral samples. For floral tissues, HIS3, UFD1, and TMP50 were the most suitable reference genes. Furthermore, two target genes, TruPI, and TruFBP24, involved in floral organ identity were selected to validate the most and least stable reference genes in staminate flowers, perfect flowers, and different floral tissues, indicating that the use of inappropriate reference genes for normalization will lead to the adverse results. The reference genes identified in this study will improve the accuracy of qRT-PCR quantification of target gene expression in andromonoecious T. rupestris flowers, and will facilitate the functional genomics studies on flower development and sex differentiation in the future.

LbCML38 and LbRH52, two reference genes derived from RNA-Seq data suitable for assessing gene expression in Lycium barbarum L

For quantitative real-time PCR (qRT-PCR) analysis, the key prerequisite that determines result accuracy is the selection of appropriate reference gene(s). Goji (Lycium barbarum L.) is a multi-branched shrub belonging to the Solanaceae family. To date, no systematic screening or evaluation of reference gene(s) in Goji has been performed. In this work, we identified 18 candidate reference genes from the transcriptomic sequencing data of 14 samples of Goji at different developmental stages and under drought stress condition. The expression stability of these candidate genes was rigorously analyzed using qRT-PCR and four different statistical algorithms: geNorm, BestKeeper, NormFinder and RefFinder. Two novel reference genes LbCML38 and LbRH52 showed the most stable expression, whereas the traditionally used reference genes such as LbGAPDH, LbHSP90 and LbTUB showed unstable expression in the tested samples. Expression of a target gene LbMYB1 was also tested and compared using optimal reference genes LbCML38 and LbRH52, mediocre reference gene LbActin7, and poor reference gene LbHSP90 as normalization standards, respectively. As expected, calculation of the target gene expression by normalization against LbCML38, LbActin7 or LbHSP90 showed significant differences. Our findings suggest that LbCML38 and LbRH52 can be used as reference genes for gene expression analysis in Goji.

Transcriptome profiling of grapevine seedless segregants during berry development reveals candidate genes associated with berry weight

BACKGROUND

Berry size is considered as one of the main selection criteria in table grape breeding programs. However, this is a quantitative and polygenic trait, and its genetic determination is still poorly understood. Considering its economic importance, it is relevant to determine its genetic architecture and elucidate the mechanisms involved in its expression. To approach this issue, an RNA-Seq experiment based on Illumina platform was performed (14 libraries), including seedless segregants with contrasting phenotypes for berry weight at fruit setting (FST) and 6-8 mm berries (B68) phenological stages.

RESULTS

A group of 526 differentially expressed (DE) genes were identified, by comparing seedless segregants with contrasting phenotypes for berry weight: 101 genes from the FST stage and 463 from the B68 stage. Also, we integrated differential expression, principal components analysis (PCA), correlations and network co-expression analyses to characterize the transcriptome profiling observed in segregants with contrasting phenotypes for berry weight. After this, 68 DE genes were selected as candidate genes, and seven candidate genes were validated by real time-PCR, confirming their expression profiles.

CONCLUSIONS

We have carried out the first transcriptome analysis focused on table grape seedless segregants with contrasting phenotypes for berry weight. Our findings contributed to the understanding of the mechanisms involved in berry weight determination. Also, this comparative transcriptome profiling revealed candidate genes for berry weight which could be evaluated as selection tools in table grape breeding programs.

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

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

Validation of Reference Genes for RT-qPCR Studies of Gene Expression in Preharvest and Postharvest Longan Fruits under Different Experimental Conditions

Reverse transcription quantitative PCR (RT-qPCR) as the accurate and sensitive method is use for gene expression analysis, but the veracity and reliability result depends on whether select appropriate reference gene or not. To date, several reliable reference gene validations have been reported in fruits trees, but none have been done on preharvest and postharvest longan fruits. In this study, 12 candidate reference genes, namely, CYP, RPL, GAPDH, TUA, TUB, Fe-SOD, Mn-SOD, Cu/Zn-SOD, 18SrRNA, Actin, Histone H3, and EF-1a, were selected. Expression stability of these genes in 150 longan samples was evaluated and analyzed using geNorm and NormFinder algorithms. Preharvest samples consisted of seven experimental sets, including different developmental stages, organs, hormone stimuli (NAA, 2,4-D, and ethephon) and abiotic stresses (bagging and girdling with defoliation). Postharvest samples consisted of different temperature treatments (4 and 22°C) and varieties. Our findings indicate that appropriate reference gene(s) should be picked for each experimental condition. Our data further showed that the commonly used reference gene Actin does not exhibit stable expression across experimental conditions in longan. Expression levels of the DlACO gene, which is a key gene involved in regulating fruit abscission under girdling with defoliation treatment, was evaluated to validate our findings. In conclusion, our data provide a useful framework for choice of suitable reference genes across different experimental conditions for RT-qPCR analysis of preharvest and postharvest longan fruits.

Suppression of the D-class MADS-box AGL11 gene triggers seedlessness in fleshy fruits

Seedlessness, one of the most desired traits in fleshy fruits, can be obtained altering solely AGL11 gene, a D -class MADS-box. Opposite to overlapping functions described for ovule identity. AGAMOUS like-11 (AGL11) is a D-class MADS-box gene that determines ovule identity in model species. In grapevine, VviAGL11 has been proposed as the main candidate gene responsible for seedlessness because ovules develop into seeds after fertilization. Here, we demonstrate that AGL11 has a direct role in the determination of the seedless phenotype. In grapevine, broad expression analysis revealed very low expression levels of the seedless allele compared to the seeded allele at the pea-size berry stage. Heterozygous genotypes have lower transcript accumulation than expected considering the diploid nature of grapevine, thereby revealing that the dominant phenotype previously described for seedlessness is based on its expression level. In a seeded somatic variant of Sultanina (Thompson Seedless) that has well-developed seeds, Sultanine Monococco, structural differences were identified in the regulatory region of VviAGL11. These differences affect transcript accumulation levels and explain the phenotypic differences between the two varieties. Functional experiments in tomato demonstrated that SlyAGL11 gene silencing produces seedless fruits and that the degree of seed development is proportional to transcript accumulation levels. Furthermore, the genes involved in seed coat development, SlyVPE1 and SlyVPE2 in tomato and VviVPE in grapevine, that are putatively controlled by SlyAGL11 and VviAGL11, respectively, are expressed at lower levels in silenced tomato lines and in seedless grapevine genotypes. In conclusion, this work provides evidence that the D-class MADS-box AGL11 plays a major and direct role in seed development in fleshy fruits, providing a valuable tool for further analysis of fruit development.

Reference gene validation for quantitative RT-PCR during biotic and abiotic stresses in Vitis vinifera

Grapevine is one of the most cultivated fruit crop worldwide with Vitis vinifera being the species with the highest economical importance. Being highly susceptible to fungal pathogens and increasingly affected by environmental factors, it has become an important agricultural research area, where gene expression analysis plays a fundamental role. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is currently amongst the most powerful techniques to perform gene expression studies. Nevertheless, accurate gene expression quantification strongly relies on appropriate reference gene selection for sample normalization. Concerning V. vinifera, limited information still exists as for which genes are the most suitable to be used as reference under particular experimental conditions. In this work, seven candidate genes were investigated for their stability in grapevine samples referring to four distinct stresses (Erysiphe necator, wounding and UV-C irradiation in leaves and Phaeomoniella chlamydospora colonization in wood). The expression stability was evaluated using geNorm, NormFinder and BestKeeper. In all cases, full agreement was not observed for the three methods. To provide comprehensive rankings integrating the three different programs, for each treatment, a consensus ranking was created using a non-weighted unsupervised rank aggregation method. According to the last, the three most suitable reference genes to be used in grapevine leaves, regardless of the stress, are UBC, VAG and PEP. For the P. chlamydospora treatment, EF1, CYP and UBC were the best scoring genes. Acquaintance of the most suitable reference genes to be used in grapevine samples can contribute for accurate gene expression quantification in forthcoming studies.