ISSN: 2167-0870
+44 1478 350008
Research Article - (2022)Volume 12, Issue 3
Background: Vaccinium bracteatum Thunb. (VBT) is widely distributed in the mountainous areas in eastern and southern China. VBT leaves have great medical value and can be used to stain rice to produce “Wumifan”. Its fruits also contain rich nutrients. However, there has been limited attention to exploring the molecular content of VBT. Previously, we performed RNA-seq on three typical VBT fruits that were at various stages of ripening, although a reliable reference gene was lost in validation.
Results: In this study, we selected ten candidate reference genes based on previous studies and transcriptomics analyses. Subsequently, these genes were evaluated using a combination of methods, including geNorm, NormFinder, and Bestkeeper, with a comprehensive ranking assessment. As a result, we found that the actin2, NADH, and ADK genes have high reliability for analysing the expression levels of genes involved in fruit development. Furthermore, the transcript levels of 15 DEGs from transcriptomic analysis were assessed using NADH as a reference gene, and RT-qPCR data were highly consistent with the transcriptomic data.
Conclusions: These results provide a reliable reference gene for further studying gene expression, which will be beneficial for comprehensively exploring VBT.
Vaccinium bracteatum; Reference genes; Gene expression; RT-qPCR
Vaccinium bracteatum Thunb. (VBT), which belongs to the Ericaceous family and is also called wild blueberry, ”Nanzhu”, and “Wufanzi”, is diffusely distributed in mountainous areas of eastern and southern China [1,2]. The plant is famous for producing “Wumifan”, as its leaves contain pigments that can be used to stain the rice. Although VBT fruit is like blueberries, its phenological period is different from that of blueberries. The flowering phase of VBT is June, while fruit ripen in November and December, and its fruit is smaller than blueberries. The species is also used extensively as an herbal medicine in China, the original description as an herbal medicine was in the Tang Dynasty, recorded by the book “Ben Cao Shi Yi”[3,4]. Recently, the species was demonstrated to confer prominent benefits for hypoglaemia, sedation and hypnosis, antineoplastic, antimicrobial and antioxidant activity [3-9]. Additionally, VBT contains high levels of nutrients, such as amino acids, minerals and unsaturated fatty acids [3,5,10]. Therefore, VBT is an important resource for the traditional diet as well as for its pharmacological effects, for which it has considerable prospects for development.
Recently, an increasing number of studies have focused on the herbal effects of VBT and its potential prospects in food production. Polysaccharides [11] and flavonoids [12] were reported as major components of VBT leaves. In addition, several reports focused on the formation of “Wumi” [2,6,13,14]. The findings showed that iridoid glycosides may be the precursor compound, among a profile of pigments, to bind to rice protein molecules [15]. These studies provided novel insights for understanding VBT and identified prospects for exploration of VBT products. However, most of these studies focused on chemical compounds [4,7,11,16,17], and very few studies have concentrated on VBT’s genetic and biochemically active molecules, which may promote its application in improving the breeding of crops.
Since the leaves and fruit of VBT have great benefits to human health, the nutrient and key component synthetic pathways should be further elucidated. To discover the molecular mechanism of VBT in pigment synthesis, nutrient enrichment, and the production of other bioactive components, genomics and transcriptomics analyses offer effective methods. Transcritpomics analysis is needed to confirm the quantitative, real- time, Reverse Transcriptase Polymerase Chain Reaction (RT-qPCR), which is always used to test gene expression levels in biological samples and tissues [18-20]. RT-qPCR was calculated based on reference genes, such as actin [21,22], tubulin [22], EF1 [20,22], GAPDH [21], ubiquitin [23], and others. Actin, GAPDH, tubulin, and ubiquitin are used commonly in analysing gene expression levels in plants [23-26]. However, there are no standards to use as reference genes in various species or tissues [24]. The best reference genes were generally confirmed by analytical tools, such as GeNorm, NormFinder, and BestKeeper [25]. However, the scores for reference genes evaluated sing these three methods are often not uniform. Therefore, researchers have selected genes that have high scores in GeNorm and NormFinder or used a combination of the three methods [24]. As a result, the best reference gene remains to be identified.
VBT is a unique species in China. Because of its high medicinal value and rich nutrient qualities, the molecular mechanisms involved in nutrient and active biochemical component synthesis in VBT should be better understood. Thus, our group sequenced the genome and transcriptome of VBT. However, validation of the differentially expressed genes in the transcriptome involved in various fruit- development stages using the actin gene as a normalizer was not stable and the results showed pronounced variability across samples of fruits at the same developmental stage. Moreover, the best reference gene in the development of fruits is not actin which, therefore, implies that it is not adaptable for validation of genes involved in fruit development. Consequently we selected ten candidate reference genes from the VBT transcriptomic data that had similar gene expression levels in different fruit development stages. Additionally, we evaluated the ten canidates using GeNorm, NormFinder, and BestKeeper. The results will benefit the future analysis of gene expression in VBT.
Sample collection, RNA extraction and cDNA synthesis
The fruits were collected from a blue berry orchard (30°51′N, 118°23′E) in Wuhu city. According to monitoring fruit ripeness, we went to the company at different times. Then, three green (S1), red (S2), and blue (S3) fruits were collected in 50 mL tubes, and three tubes were prepared for each fruit. At the same time, they were placed in liquid nitrogen and later kept in an ultracold storage freezer for future study.
Total RNA of fruits from the three ripening stages was extracted using a mirVana miRNA Isolation kit (Ambion). The total RNA was later evaluated using gel electrophoresis, and integrity was assessed using Agilent Technologies (Santa Clara, CA, USA). Once the RNA samples were of good quality, cDNA was synthesized using a PrimeScript RT reagent kit with gDNA Eraser (Perfect Real Time, Takara, Japan) following the manufacturer’s instructions.
Selection of candidate reference genes
According to a previous study and the transcriptomics data, ten candidate reference genes were selected [18]. The genes consist of RP (50S ribosomal protein, CL31684Contig1), actin5 (actin related protein-5, CL33133Contig1), UBQ (Ubiquitin domain-containing protein, CL33177Contig1), actin1 (actin related protein-1,CL34455Contig1), actin4 (actin related protein-4, CL8799Contig1), Ubiquitin- Conjugating Enzyme (UBE, CL31240Contig1), actin3 (actin related protein-3, CL40530Contig1), actin2 (actin related protein-2, CL11169Contig1), NADH (NADH dehydrogenase, CL31981Contig1), and Adenylate Kinase (ADK, CL27968Contig1) (Table 1). These genes had similar FPKM values in the three ripening stages according to the de novo transcriptomics analysis.
id | Sequence | foldChange | log2FoldChange | pval | padj | fpkm_S1-1 | fpkm_S1-2 | fpkm_S1-3 | fpkm_S3-1 | fpkm_S3-2 | fpkm_S3-3 | length | swiss.Des | named in this study |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CL34455Contig1 | AGACTCGTCTTGTTTATGTACGTTTGGTTGGTAAGCAAATCAAGCTCAAACACCTTGATAATAGATGAGCTGAGTTTGAAACAGACACGAGCCGAGTTCGAGCTCTACTAATTTTGGACCAGCCGAGCTTGAATAATAGTACAAAATTCAGCTAGGTTCCTCTTTAGCAGCCCAACTATAAGCCCAAAAAATTGTTGCTTTTGTCTATTTCACATGACAATAATGCAAAGCTCCATGCGTAAAATGAGTAAAGGCACTCCAAGCATGGAAAAGCTACAATCTGTAACTCTCTCTCTAGAAAAGAGTTCAAGCTTCCGGAATTAAGAAAAGCATGGAATTTCGATCCGATTTCATCTTCCTCAGTCTGTAAAGCAAGAAAAATCCACGCTTTGGAACCACTTCAATTCAAAAAAAACCCAAAGCAAAAAGAGAAAAAGACGCACACAAACCCACACCCACAAAGAAACTCTTCAGAATTTGTAAAGAAAACTTGAGAATTTCTAGAGAGAGAGAGAGAGTGATGATACTTTTGCAGTCTCACTCGAGGTATCTTCTACAGATCTTGTCAAATCGGGTTCAGAATGTTGAGAAAGGGGTTGAATTGGATAGCTACTGGGTTGAATTTGGTGATATACGTTACCATATTCAGGCATCAATGAAGAATCCCCATCTCCTGTTGTTATCTGTCTCATTACCCACCCCGCCTCCAGAAACTGTTTTCTTTGGGGGACTTCCTCTTGGAGCCATAGAAGCTGTAAAAGCTGCTTATGGTCTGGTCATACAAATTCTTGATCCTCCAAGAGATGGCTTTAATCTCACCGTGAAGTTGAATTTAACTAAACTTCCTCCAGATGAAGAGAACAAAAACAATCTTTTGGTGAAGATTGCATCTCTGAGGGAAGTAGTAATGGGTGCCCCATTAAGAGTAATGCTAAAGCATCTTGCTTCAAGGACTGTTGCTCCGGATATAGATAGGCTTGTTGCCCTTGTGCATCGGCCTACAGAGTCTTTTTTCCTTGCTCCGCAGGCAGAAAAGGTTACTGTGGTTTTTCCTATGAGATTCAAGGATTCCATAGATACTGCTCTTGCAACTTCATTCCTGCAGGAATTTGTGGAAGCAAGACGTACAGCTGGACTTTTCAATGCCCCTCCTTGTTCATGGTCTCCTTCTCCTCCCTCAGAACTGAAGGGTGCGCCTGTTGAGGCGTTGTCTGCAAATGCCGGATTTGTCTCTTTTGTTCTTTTCCCACGTCACGTAGAAGGCAATAAACTCGACAAAACTGTTTGGAATTTATCGACGTTTCATGCGTATGTCAATTATCATGTAAAGTGTTCAGAAGGATTCATGCACACTCGAATGAGGCGGCGTGTGGAGTCCCTGATTCAGGCTCTGGATCGGGCCAAACCAGATAAAGAGAAAGCAAAGAAAACTGGACAAAGTAGATCTTTCAAGCGGCTGAGCCTCAAGGAAGCCAGGAGCGGGTTGAAGTCTTGAAGTTGCGGATGCTGTAATCAAGGGGAGAACATGTGAATTGGCTGGGTTTGCCCAGATCTTTGAAAGAAATGGTTCTGTTTTTGCTTGGAAGATCACTGCAGGCACTTGAGGCCAAGAAGACGTCTAGAGTTTAGAGACTATAAAATGTGATTTTTGTCGGGCTGTACTGGTGCCACTTAACACATTGTAAAATTTCCCCATTAGTAGGTGTTACATCGCTGGTTAACTCATGTTATTTGGGTCATGACCCTGTAAATCTAGCTGCGATCGCTTGAGCATGTGTAACGTCTTGAACAATCTTATCTAGCCCTAATAATAATCCGACACATGCTTCAAAG> | 0.842626927 | -0.247034076 | 0.142005735 | 0.329670456 | 15.30243 | 14.53465 | 12.84122 | 12.37452 | 12.51225 | 12.06213 | 1833 | Actin-related protein 2/3 complex subunit 2A | actin1 |
CL11169Contig1 | CTCGAAATGATTCCCCACTCTCTCTCTTCAAAATTTAAAAATCGGAATTTAGTCCTCATTTCATTTGCATTTCAAATCGCCATTATACCTCACTATATTATCTCCCTATCTACTTCCTTTCCTCTCCTCATTTCGTTTTGGCCGTTCGACGGAGGAGAGAGAAAACAGTGAGGTCTAGTGAGAGAAAGCTTGCGTGATTTTCGATCTACAATCCAAACTACGTTATAAAGAATGGCTGATGCAGAGGATATTCAACCCCTTGTCTGTGACAATGGAACTGGAATGGTCAAGGCTGGTTTTGCCGGTGATGATGCTCCCAGGGCAGTTTTCCCCAGTATTGTAGGTCGACCCAGGCATACTGGTGTCATGGTTGGGATGGGCCAGAAAGATGCTTATGTAGGTGACGAAGCCCAATCCAAAAGGGGTATTCTTACCTTGAAGTATCCCATCGAACACGGTATTGTCAGCAACTGGGATGACATGGAAAAAATTTGGCATCACACATTCTACAATGAGCTTCGTGTTGCTCCAGAAGAGCACCCAGTTCTTCTTACAGAGGCACCACTCAACCCTAAGGCCAACAGAGAAAAGATGACCCAAATTATGTTTGAGACTTTCAATGTTCCTGCCATGTATGTTGCCATACAGGCCGTACTCTCTCTATATGCCAGTGGTCGTACAACTGGTATTGTGCTGGATTCTGGTGATGGTGTGAGTCACACTGTGCCTATCTATGAGGGATATGCTCTTCCTCATGCGATCCTCCGTCTTGATCTTGCTGGTCGCGATCTTACAGATGCCCTCATGAAGATCCTTACTGAGAGAGGTTACATGTTCACCACCACAGCCGAACGGGAAATTGTCCGTGACATGAAGGAAAAGCTTGCATACGTGGCTCTTGACTACGAGCAAGAGCTGGAGACTGCAAAGAGCAGCTCCTCAGTTGAAAAGAACTACGAGCTTCCTGATGGTCAGGTCATTACCATTGGGGCTGAGAGATTCAGATGCCCAGAAGTCCTCTTCCAGCCATCTTTGATTGGAATGGAAGCTGCGGGAATCCACGAGACTACGTACAACTCAATCATGAAGTGCGATGTGGATATTAGGAAGGATCTTTACGGAAACATTGTGCTCAGTGGTGGTTCAACTATGTTCCCTGGTATTGCAGACCGTATGAGCAAGGAAATTACTGCTCTTGCTCCAAGCAGCATGAAGATCAAGGTTGTGGCACCACCAGAGAGAAAATACAGTGTCTGGATTGGAGGATCTATCCTTGCATCCCTCAGCACCTTCCAACAGATGTGGATTTCAAAGGGTGAATATGACGAATCTGGCCCATCCATTGTCCACAGGAAGTGCTTCTGATTTCTCAATTGCTTTGGTGGTGAGTTTCTGTTTCGTTAATTTGTTGGCTTTTCATGTTAATTGTCATGAAGTCAGTCTGGTTGACATGGAATTGCTGAGGTGGGGCTTATGAATGGAGGGATTAATTTGATTGGGATGGATTGGTTTCTCTTCATCATCAGATGCGGTTCATTTTTAGACGCTGGCCTTGCTTCATCTGTGATTCAACCGTCTTTTTAGTAGGATGCTTGTAGGTTGAGAGTGGTTGTGGTGGATTTTTCTTTTTTTCTTTCTTTTTTCCCTGTCTTTTCTTTCAATTTTCAGATTTGTTTCTTTCCCCTGTTATTATGCTAGTATTTACTGTATGAGAATTTTTTATTAGCGTCTGTTTGCGGTCAGAAATTATTATGGAGATCATATTTTGATACTGGGTTTAAACATTCCAAGTATTATGCATGGGTAATAAAGTACTTTGATAACCTTTGTGTCTGCTGGTCAATCCAAATCTAAAGCAATGAAACACAGTTGTGGGCACGGATGTTAGGCACATGCCGAACACGGATGCTTTCGGGATGTGTTTGAGCACCTTTTTGTTGTGGGGTTCACGGTTCTGGACTCATCTAAGTATGCATCTGTGTATGTATATCTTTGTGATCAAACTTCTCCAATATGTTGCGGAGGATCTTCATTTAAGTTATGTTTTCAGCAGATGCCTCCCATCTTGTTTCTTGCAATTTGCCAGGGAATGCTGGTGTCTTGTCCTGTCATATTCATGTAATTTGTAAAGCATGCGAACTCATCTTATAGAGGGGTGGGGTCACTTGAATCTTGGTCAGATTTGCCGAATGTTATCAGAATGGCC> | 0.935542942 | -0.096124219 | 0.162582637 | 0.363935442 | 713.6686 | 712.0665 | 763.8563 | 708.4312 | 692.0251 | 704.3399 | 2206 | Actin-7 | actin2 |
CL40530Contig1 | GAAGGCAAAAAGAACTTACCGCTTTATTGTATTCAAGATTGAGGAGAAACAAAAGCAGGTGGTTGTGGAAAAGGTTGGTGAGCCCACTGAGAGCTATGATGACTTCACTGCGAGCCTTCCTGCTGATGAGTGTCGATATGCTGTCTACGATTATGACTTTGTGACAAATGAGAATTGCCAAAAGAGCAGGATTTTCTTCATTGCATGGTCCCCCGATACATCACGGGTGAGGAGCAAAATGATTTATGCAAGCTCCAAGGACAGATTCAAGAGGGAGCTAGATGGTATCCAGGTAGAACTACAAGCAACAGATCCAACCGAGATAGGCCTTGATGTTATCAAATGCCGCGCTGAAAAATAATTGATTCTCTTGTCGTGCTGCTCGTGCATTGGAGTACTTTGCTCGTAGCTTCACTGAAGAGAAATGTTAATAGTTTCGGGACAGTGGAACCTGTTTAAGATATGATATTGGTTTGCTGTTCAGAGCTCGACATTATAATATTTGTCCTCATCTAGAACGAGTTCAGTAATTTGGAGCTTATTTGATGACTTCATCTTGTGTGGAACTTAATATACTTGTGTCTAATTCGCAAGACCAATGTATGTCAATATCTGCTATCATCTTTTGGTTTGTGTTTCCCGACCTTGCGTGAGACTTGTTCTGATTTGCAAGACCAGCGTCTGCTATCACCTTTTGGTTTGGTTTTAAGTGGGTAGTGCTTGAAAACAACTGTACATCACGGGGCCGAAGTAAGAAATTTTTTACCCCAGTATGTTTGAATTTGTTCACCATCTGGTATTCTGGTGAGGAAGTAAGCTCGTAGAGCTGGGAATAGTGTCCTTTGTTTTGAATATTATGATCGATTAGCTGAT> | 0.990478096 | -0.013803024 | 0.974913988 | 1 | 116.8511 | 143.9433 | 126.9255 | 135.1737 | 133.0853 | 127.0351 | 873 | Actin-depolymerizing factor 2 | actin3 |
CL8799Contig1 | AATGGCCCTTTTCCCCCAAGAACGTGTCAACCGGGTCCGTACATCCAGACTGTAATCCGATCCAACGGTCCGGATCAAATCTCTCCCCCCTCCTTTTTCTATAACCCACAAACTCCATTCTCCCATTTCGCGTCCATTCATCTCCGCATACTTGACTTACCCTCCTCTTACTTTCTCTCATCTCTCTCTCTCTAGGATTTCCAGGAGGTAGCATAAGATGGCAGAAGCTGAGGATATTCAGCCACTCGTGTGCGACAATGGAACAGGAATGGTTAAGGCTGGATTTGCTGGAGATGATGCCCCGAGGGCTGTGTTCCCTAGCATTGTAGGTCGACCTCGTCACACTGGTGTGATGGTTGGCATGGGCCAGAAAGATGCTTATGTTGGAGATGAAGCCCAATCCAAGCGTGGTATCTTAACATTGAAATACCCAATCGAGCATGGTATTGTGAGCAACTGGGATGATATGGAGAAGATCTGGCATCACACGTTCTACAACGAGCTCCGAGTTGCCCCAGAGGAGCACCCTGTCCTCCTTACTGAAGCACCTCTGAACCCCAAGGCTAACCGTGAGAAGATGACCCAAATCATGTTTGAAACCTTCAACACCCCTGCTATGTATGTTGCTATTCAGGCTGTGCTTTCCCTATATGCTAGTGGTCGTACAACTGGAATCGTGCTGGACTCTGGTGATGGTGTGAGCCATACGGTCCCTATATACGAAGGATACGCACTCCCACACGCCATCCTTCGTCTTGACTTGGCGGGCCGCGATCTCACCGATGCTTTGATGAAAATCCTTACTGAGCGCGGCTACTCCTTCACTACTACAGCAGAGCGGGAAATCGTGAGGGACATGAAGGAGAAGCTCGCTTATATTGCCCTAGACTATGAACAAGAGCTAGAGACTTCCAAAACAAGCTCCTCTGTTGAGAAAAGCTACGAGCTTCCCGATGGCCAAGTGATCACCATTGGCGCGGAGCGTTTCCGGTGCCCCGAGGTTCTTTTCCAGCCGTCGATGATTGGAATGGAAGCGGCGGGAATTCACGAGACCACTTATAATTCTATCATGAAGTGTGATGTGGATATCAGAAAAGATCTCTACGGAAACATTGTGCTCAGTGGT> | 1.193350158 | 0.255017426 | 0.019414568 | 0.069357816 | 55.1562 | 56.18276 | 61.27122 | 68.94726 | 72.84101 | 69.88845 | 1126 | Actin-97 | actin4 |
CL33133Contig1 | AGAGAGAGAGGAGAGAGAGAGAGTCGTTGGGTAGTCATTACCAAGGAAAGAATCCGATACACAGAGAGAGAGAGAAGATAATTCAATGATTACAATTCATATGCTCGTCTCTGAAATATGAAATGGATGGAGTGAAAACCATCTTCTCCGACCCACCACCACTTCACTTTCCGTCGTTGAGCCGCCCTTTTCGCTTCCAATCAACGTCGTCTTCTCCACTACCATGATCACCCCCTTCCTTAAACCCCTTCAAGGATCTCAGCAGCAGATCCACCACCTCCTCGCTCTTTCTCTCTCCTCTCTCCTCTCTCATGGACCCCGCCACCTCTCGCCCCGCCGCTGTCATCGACAACGGCACCGGGTATACTAAGATGGGATTTGCTGGTAATGTAGAGCCGTGTTTTATACTTCCTACAGTCGTCGCCGTCAACGATTCATTCGTAAACCAAACTCGAAGCTCTTCAAAGGGAAGCAATTGGCTAGCGCAGCACAGCGCAGGCGTGATGGCCGATCTCGATTTCTTCATAGGAGACGAGGGTTTGTCTAGATCCAAATCGAATAGTACGTATACTCTAACGCATCCGATTAGGAATGGCCAGGTTGATAATTGGGATGCCATGGAGAGGTTTTGGCAGCAGTGTATATTCAATTACTTGCGGTGTGATCCAGAGGATCATTATTTCTTGTTGACTGAGAGTCCGCTGACTCCGCCGGAGAGCCGGGAGTACACAGGCGAGATTATGTTCGAGACTTTCAATGTGCCGGGGCTTTATATCGCGGTGCAGCCTGTGCTTGCACTTGCAGCCGGGTACACGACGTCTAAGTGTGAGATGACAGGGGTTGTAGTGGATATTGGAGATGGGGCTACTCATGTTGTCCCTGTTGCAGATGGTTACGTTATTGGGAGCAGCATTAAGTCGATTCCAATTGCTGGAAAAGATGTAACTCTCTTCATCCAACAGCTTATGCGGGAAAGAGGAGAGCATGTTCCTCCTGAAGATTCATTTGAAGTAGCTCGAAAGGTGAAGGAAACATATTGCTACACTTGTTCTGATATTGTCAAGGAGTATAATAAGCATGACAAAGAACCAGCTAAGTACATCAAGCATTACAGAGGTATTAAACCTAAAACGGGGGCACCATATTCTTGTGATGTTGGTTATGAACGTTTTCTTGGCCCTGAGGTTTTCTTTAATCCTGAGATTTACAGCAACGATTTTGCAATGCCCCTACCAGCTGTAATTGACAAGTGTATTCAGTCTGCACCAATTGACACAAGAAGGGCTTTGTATAAGAATATTGTATTGTCTGGTGGGTCAACAATGTTCAAAGACTTCAACCGAAGATTACAACGGGATCTAAAGAAGATTGTTGACGCTCGTGTACTTGCATCTGATGCCCGGCTCGGTGGAGAGTTAAAGTCACAACCAGTGGAGGTCAATGTTGTCAGCCATCCCATCCAGCGGTATGCTGTCTGGTTCGGAGGCTCTGTCCTTGCATCAACACCAGAATTTTTTGCGGCTTGTCACACCAAAGCAGAATATGAAGAATATGGGGCAAGCATTTGCCGTACAAATCCTGTTTTCAAGGGTATGTATTGACGCAAATAACTATGTCAGAAGATGTTGTATGGGCATTCAACGATTGTAATGCCTCATCATTTCCACTTCCAACGTGCTCCAGATTTGGAGATCCAGTCTACAGGCGGAAATCATGTGCAATTTAGAGGCTTAGCTTGCAGGAGGTGGTGTACAATATATCTTTGGTTTGATTTGTACGAAAGTGAAGTTTTTGTATTCTTTCATTAGCGTAGGGGGCAGCAAGTAGTCAAAAGACTGACCATGTCTAAAGGCATATCTACATTGGTAACACTGATATAAGAAGGTTTTTTTTTTTTTCACTGTTGTTAGGGAGTTGATATTCCCTTATAGGACATTCTCAATTAGGTGGTTTGATTCATTCGCATTGGATTGAGGAGGAATAGGGTTGGCGTAATAAGTAGAAAAGTTCCCTTTATTTGTAATCAAAATCTGAACTTTTCATTAGTGTCTGTACATTGTTAATCTTTACATATTGGTGAACTGCTTGCAGACTTATTAGTTGAGGTTGTTGATAGAAGCAGCGAGAT> | 1.046441256 | 0.065491325 | 0.617640862 | 0.828321542 | 26.51006 | 22.48087 | 23.33269 | 25.11693 | 25.60212 | 26.98382 | 2128 | Actin-related protein 3 | actin5 |
CL31684Contig1 | GATGATGTCATCAAAACTAAAACAAAAACTAAATTGGGAAACAAACAGGGCCAATATAATTTGGAGATAAAAATTTGCGTTAATCGGTCGGTCAATAACATAATTTTGAGTCTGTTGGGGCTAATCGGTGATATTCAGAAGTTGTAGGGCGAAATTTGCGGAACTAAGAAAATAAACATGTATCAGCCCTAAATTACAATTCTATTTACCGCTCCAAAACCTTCATCGTAAATCGAACAAAGCTGGCGGAAGCACCGGCGACTATCGGCGGCTTCGGCGAAGAGTAGCACAGGAAGCAGGAAGATGAATTTTGTGACGTCACTTTGTAGAAGATTGAATATCAATGAGCTGGTTACAAATGCTCCTGTATACAGCTCTGCTAGTGATATAACTGGGGGAGGATTAAGTTTGACATTCAGACGCTGGGCTACAAAAAAGACGGCGGGATCCACAAAGAATGGTCGAGACTCCAATCCAAAGTATCTTGGAGTAAAAAAGTTCGGTGGAGAGAGAGTGATTCCTGGGAATATCATTGTTCGTCAAAGAGGAACCCGATTCCATCCTGGAGACTATGTTGGAATGGGAAGGGATCACACCCTTTATGCTTTGAAAGCAGGTTGTGTCAAATTTGAACGCCATAAGCTGAGCGGAAGGAAGTGGGTGCATGTTGAGCCGAAGGAAGGCCATGTCCTCCACCCTGTCTATTCAGATGCTTCATCAGCTCCTAAGATGCAGACTGCTTTCTAGACTTCATCCTTTCATTTTCATCGAGTGTTCAATATCTGATGCAACTTCAAGTGACTGTTCAGAAACAGTTCTGTTTTCAGTTCTCTCTCTTTTTTTTTTTTTGCTTTTGTTTCTTCTGTTTTTAGAAACAGGCAGCCCTATTACTTGATTGGAAACAAGCTGCATCTGAAATCTTGTAGCCAAAAAAATACTGCCTTTGATTCATGAATTGGACTCAACATGCGTCTTGTGCTGTTGTCACGCCAAGATTCAGAGCTAAATAAAGTGTATAAGTACTTTTATGTGTTTTTCGTTTTCAATTACCGGAGATGAGTATATCGTTTGGGTTTGTCTTTTGTTTTTCATTCTCCTCCAGTCGATGAGGGCAGGAAGCTTGGTTTGGTTTTGTTTCTTGTTTTGAAGTGTGGTGTTCAATAAGGCCACACCAAATTCAATTGCAG> | 0.870193741 | -0.200591455 | 0.116041215 | 0.342337153 | 58.09411 | 62.86819 | 52.09703 | 48.09272 | 48.79577 | 42.17338 | 1187 | 50S ribosomal protein L27 {ECO:0000255|HAMAP-Rule:MF_00539} | RP |
CL31240Contig1 | CTCTCTCTCTCTCTCTCTCATCGATCGGAGACTTGTTGGTTCACCGTAACTACTACGCCGGCGATCTATTGAATCGTACAGTTTCTCATTCGACGTCAGTAGATAATTTGTACTCAAGATTACCAAGATGTCATCTCCAAGCAAGAGGAGAGATATGGATGTCATGAAATTGATGATGAGTGACTATGCTGTGGAGCCTATAAATGATGGAATTAGCGAGTTCAATGTGGAATTCCATGGTCCAAAAGAAAGCCTTTATGAAGGTGGTGTATGGAAAGTCCGTGTTGAGCTACCGGACGCTTACCCTTACAAGTCTCCTTCCATTGGGTTTCTGAACAAAATTTTCCACCCAAATGTTGACGAGCTGTCTGGTTCTGTATGTTTGGACGTAATTAATCAGTCATGGAGTCCAATGTTTGATCTGCTCAATGTGTTTGAGGTTTTTCTTCCACAGCTTTTGCTTTACCCAAATCCTTCAGACCCATTAAATGGTGATGCAGCATCGCTAATGATGAAGGATCAGAAGCAATATGAACAGAAAGTGAAAGAGTATTGTGAACGTTATGCAAAGAAAGAGAATATCACCAGTACTTTGTCAGACGAGAGTGATGAGGAGATAAGTGACGAGGATGTCAGTGATGGCCGAAGTTTGTCTAGTGAAGATGATGTTGCTGGAGATGCAGACCCTTAAAAGTGATGAAAAGCCCTGGACTTCAATGTGCCAAATCTATTTCATGTTTGGTTTTTGACTTTTACTTAAATGTATGACAAAATGCTATTTTCTTGTAAATTTATTTATGGAGCAGTAGTTAGAGACTGGATACAAGGTACCATGCTTGGTATCTAATCCATTATATATTCTGGCACATTTGTTCAAAGAAGAATAGAAATGAAGTTGTCCTTTTCATTCAATCTCTTGTTACTGTTCCTGTAATTATTTGAGCATGTAACAAATTCTCGGATTAGGCGTTTATATTAAGTTTGTCGTACTTTAGTGATTTGTTATGCAAACCAATAAAGCAGGGCTCATCCTTATCACTCATCAAAGTCGTGTCTAGCTGATCCAGATCACTGCCTGGCTCTCACATGAGCTTGGTGCGAGTCTTTCTAACTTTTTCAAGCAT> | 1.027003472 | 0.038441059 | 0.832424243 | 0.957814758 | 18.12283 | 15.90809 | 14.38538 | 16.19619 | 16.93001 | 17.95536 | 1124 | Ubiquitin-conjugating enzyme E2 5 | UBE |
CL33177Contig1 | GTTTTCAGCTACACTCAGCAAGGAAATCCGTCCTGTTTGCAAGAAATTTATGCAAGATGTAATGTCCCACAGCCAGCTCTACTTGAATTTAGAAGAACCCTTTTCCGTAAACGCCCTGCTTCTGTAGGCGATGACTAGATAGTTGGGGTTGTGGTGATTGCACCGTTATACCTTCTTTCTTTCTTTTTTTTTTTTTTAAATTGATTTGGGAATGGGTGCGGAAGGGGAGGCAGCTGAGTCGAGCTTGGGCGTGAGCGGTGGTGAGGAAGAATTGGTAGCTCTCAACATTCGATGCTCCAACGGATCGAAATTCTCGGTGCGGATGAGCCTCGAGGCGACCGTCGAGGCGTTCAAAACTGTTCTTGCGCAGAACTGTGATGTGCCGGCTGATCAGCAGAGGTTGATTTACAAGGGTCGGATCTTGAAGAATGACCAAACCCTAGCTAGCTATGGTTTACAGGCAGATCACACTGTCCATATGGTACGCGGTTTTGCTCCAGCAGCAGCATCAACGCCTGCTGCTGCTGCTGCTACTGGCAATCCTACTCCTACCTCTGGCAATGCACAGGGTGTTGGTTTGAATGAGGGTGCGGGCGTTGGAGGTGTTGGTCTTGGTCCGTCTCTATTTCCTGGGCTTGGTCTCAACTCCCTGGGAGGCTCTGGAGCATCTGGATTATTTGGAGCGGGGCTGCCGGAATTTGAACAAGTGCAGCAACAATTGACTCAGAATCCCAACATGATGAGAGAAATAATGAACATGCCAGCCATTCAGGGCCTAATGAATAACCCTGAGCTGATGCGCAGCATGATCATGAACAACCCTCAAATGCGCGAGATCATTGACCGTAATCCAGAGCTTGCTCATGTTCTCAATGATCCTGGTATCCTTCGACAGACATTAGAAGCAGCACGAAACCCTGAGCTTATGCGGGAGATGATGCGAAACACTGACAGGGCAATGAGTAACATTGAATCTTCTCCCGAGGGATTCAACATGCTCAGGAGAATGTATGAAAATGTTCAGGAACCATTTCTCAATGCTACAACTCTGGCTGGAGGCACTGGAAATGATCCAGGGTCAAACCCGTTTGCAGCTCTTTTGGGAAATCAAGGTGGTGTTTCCCAAACCAGAGATGGATCTAACAATACCTCAACCACTGGGTCTCAGACAACTGCAGGGTCTACTGTTCCAAACTCCAATCCCCTTCCTAACCCTTGGGGTGGAACTGGTGGAGGTCCCCAAGCAAACACTACTGCAAGGTCAAATCCTTCTGGGGATGCCAGGGCACCAACCATTGGTGGGTTAGGAGGGCTTGGTGTTCCAGAGTTGGAGCGCATGTTTGGCATGCCTGACCCCTCTTCAATGAATCAGCTCTTGCAGAATCCAGCTGTGTCACAAATGATGCAGAGCTTATTATCCAACCCCCAATACATGAATCAGTTTCTTGGTCTGAATCCGCAACTCCGGAATATGCTCGACTCGAACCCTCAGTTAAGGGAAATGATGCAAAACCCAGAGATGCTTAGGCAATTAACCTCGCCTGAGACTATGCAGCAAATGTTGACTGTACAGCAGGGTCTTTTGTCTCAGCTTAACAGGCAGCAATCGACTCAGGAAACAACTCAGACTGGTGGGACTCAAGGAACGCAAAACAATATGGGACTTGAGTTGTTGATGAACATGTTAGGTGGACTTGGAGGTGGTGGCTTGGCTGTTCCAAACACGCCTAATGTGCCTCCAGAAGAGCTTTATGCAACTCAACTGCAACAGCTTCAAGAAATGGGTTTCTTTGATACTGCAGAGAATATTAGGGCACTTCGTGCCACCTCAGGAAATGTTCATGCTGCTGTGGAGCGACTTTTGAACAACCCTGGCCTGTAGTTAGTGCACCTGTATTGATGTTCAGTCAAAAATACAGTTGGGCGTGTAAAAGTGGACCATTTTCTTGGAAAGGATGAAAAACGTGTTGGCACTCTCATTTTGATTTCGGGCCTCGTGAATAAGTACACTTTCAATCAATGTTTAGAGGTGAACCTGGCTGACAAAGTATTTTGCATTACCACTGTTACATGCCCTCTCTCTTTCTCTTGTTCTCTTTCCCTCCCTCTCTGTGACCGTGCGCTGATGTACATACACAATTAAGTATTCTTCCGTGGAACATTTCAAAATTGTATGGATGTTTGAGGATGTGCAATTCCTGTCTTCTGTTTGTGTTTTAGTTGAATTGCCGTTATGCTGGGAAATTTGCTGATCTTTGTCTAGGGTTGACCGGAGCCCAAGTCTCAATTCTTATCAGCGGATAGAATTATGTTTATCTTCTTCTTCTTCT> | 1.084343028 | 0.116821221 | 0.242721939 | 0.478244521 | 55.6577 | 46.33849 | 50.76082 | 56.82728 | 56.23069 | 56.86632 | 2330 | Ubiquitin domain-containing protein DSK2b | UBQ |
CL31981Contig1 | CACACAGAGAGAGAGAGAGAGAGTGAGAGAGAGAGGAAAAGCAGCGATGGGGAGGAAGAAGGGAGTGGTGTTCGACGAGGGGGCGCCGGATGATTTCGATCCAGAGAACCCGTACAAGGACCCAGTGGCGATGCTGGAGATGAGGGAGCATATAGTGAGAGAGAAGTGGATCCAAATCGAGACGGCCAAGATCATCCGTGAGAAGCTCAAGTGGTGTTACCGCATCGAGGGCGTCAACCACCACCAGAAGTGCCGACACCTCGTCCAGCAATATCTCGATTCCACTCGCGGCATAGGCTGGGGCAAGGACCAGCGTCCCCTTTCTCTTCACGGTCCTAAGGTTGAGGCGGTTGAGGCTGAGTGATTTCTTCCAGGTCACTGAAATTAGTTGTCTCGTTGACAAATAAAGTGACGATTGCGAGGGTTACCTGTTTAGCATCTGGTAAATTGTTCTGCTGTTGGTTGTTTTCTTTCTTGTTGAACTAGAGTCAAATTATGCGGGCTTTGGCTTTATGTTTAGTCTTTTGATTTTTGGCATTGTAAACTGATAAGATTTGTATCAAAGCAACAATTCAAAGGACATCAGAAGATGGAATAAGTTTCTTTTCTTTTTCGGTGGAAAAATGCAGTGTTTGCCCGTATTTAGCTTATTTGTGTTTGAAAGTGCTATGTTGCAGACTTGCAGTATCTTGAACGAGGACTGTCTTATTCTAGCAACAACTATAGGATGAAGTTGTAATGTCTCATAATGCGTGGACCTTGTAAGGGACAAATGACAGATTGGTCCCCCTGACTTCTGGGAAATTGTATTTTTGCCTATAACCCCCTTACGAAGGGCAAAATCATTCTTTCCTCCCATACCACCGCCACATACTCATTCAAATTTTTGTGCGGGATAGAAGGCCCCTGGAGCACTATTTGGCTGTGCTTCAAGAGCATTGAATAATCATTTCCCATTTTGCCCCTCAAATCCTATGTCCAAAAAAATGGATGTACCGGTGTTTCGATGCATATCTCGTAACGTGAAACGATGAGAAACAAATGCACCAAAATACCGGTAGGAACTTATTTTTTTGACATAGAGGATTTGGGCTCCCCATTTTGCTAATTTGACCCATTTTTCACAACCTTCTTTGAAGGCTCAAGTGTAATGGATCCCCCTCTAGTTCCTAGAAATCCAATTCTAAGTGACCCTATCAGAGGATGAATGATGTACAGTGTTAGTGTTATGTATTAATAATTGAATTTAGTATGATTTATGTCGTGATCCTCTGATCTTTTTTATTTAT> | 1.000958034 | 0.00138149 | 0.887321765 | 0.981271639 | 103.2273 | 110.473 | 95.14808 | 108.8576 | 108.6172 | 100.3717 | 1289 | NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 10-B | NADH |
CL27968Contig1 | CCATCTCTCTCTCTCTCTCTCTCAAACTGCACTAAGCTTCAGATCGAAAACTTAGAGAGAGAAAACTCCCGGAAGGTCAACAATGGCGAGCCCTTCGTCCTCGTC | 1.014181162 | 0.020315383 | 0.856072983 | 0.988815686 | 27.026 | 28.76875 | 27.58738 | 25.21378 | 27.02683 | 25.80532 | 1219 | Adenylate kinase 4 | ADK |
Table 1: Candidate reference were selected from the RNA-seq data.
Primers for the candidate reference genes and DEGs were designed using Premier Primer 6 software based on the following criteria: primer length of ~18–22 bp, melting temperature 59-61°C, GC content 40%-60%, and amplicon length of 80-200 b Later, the primers were synthesized by Sangon Biotech Company (Shanghai, China) (Table 1).
RT-qPCR and data analyses
RT-qPCR was performed using Tiangen SuperReal PreMix Plus (SYBR Green) kit (Tiangen, Beijing, China) on Bio-Rad CFX Manager (Bio-Rad, CA, USA) following the manufacturer’s instructions. When RT-qPCR was completed, the melting curve and Cycle threshold (Ct) values were promptly exported to confirm primer specificity and expression levels of the reference gene, respectively. The stability of the candidate reference genes was evaluated using geNorm, NormFinder, and BestKeeper, and comprehensive rankings were assessed by RefFinder (http://blooge. cn/RefFinder/) based on the Ct values. The genes ACO, PG, and other DEGs were acquired from transcriptomics, and their primers were designed for RT-qPCR (Table 1). The expression levels of the genes were normalized to the transcript levels of the genes NADH, ADK, actin5, and NADH+ADK to confirm the reference genes. The DEGs were further analysed using the NADH, and the Pearson correlation value (R) between the RT-qPCR data and RNA-seq data was analysed using the formulas in Excel.
Prediction of candidate reference genes in VBT transcriptomics and primer design. The fruit transcriptome includes green fruits, red fruits and blue fruits, which represent the progressive stages of fruit ripeness. According to previous studies and the transcriptomic analysis, we identified ten reference genes. There were five actinrelated genes (actin-related protein 2/3 complex subunit 2A (actin 1); actin-7 (actin 2); actin- depolymerizing factor 2 (actin 3); actin-97 (actin 4); actin-related protein 3 (actin 5)), 50S Ribosomal Protein L27 (RP), ubiquitin-conjugating enzyme E25 (UBE), ubiquitin domain-containing protein DSK2b (UBQ), NADH dehydrogenase (NADH), and adenylate kinase 4 (ADK), all of which have been widely used for gene validation in plants and fruits (Table 2). Simultaneously, RNA from the green, red, and blue fruits was extracted, and pure total RNA was quantified. Furthermore, the primers for RT– qPCR were designed (Table 1), and later these primers were confirmed using melting curves. The results showed that a single peak for each primer pair was present, which indicates exceptionally high primer specificity (Figure 1).
Gene ID | Gene name | Forward primer | Reverse primer |
---|---|---|---|
5’ to 3’ | 5’ to 3’ | ||
CL12314Contig1 | HCT | AAGCCGCATCTACCACATCT | GCTCTCCTTCGCTCAATCTC |
CL32463Contig1 | FLS | CGGTAATCTCCCTCTCCAAG | AACTCCTCTCCAACCTCCTT |
CL5691Contig1 | DFR | GATGATGGCGATGGAGGAAG | AGGAGGTGAGAGGCGATAAG |
CL11378Contig1 | LDOX | AAGAAGTTGGCGGAAAGGAC | GTTGTGGAGGATGAAGGTGAG |
CL11556Contig1 | F35H2 | GCGTGGATGGACTTACAAGG | CCGACTCACTGTGTTCCTTC |
CL41813Contig1 | UGT | CGACCAACTGACTGATGCTAAG | TCTGATCTCAAGCCTTCTCTCC |
CL36852Contig1 | At4g | GGAACACGCAACATTGAAGG | GGCGATTAACTGACCTCCAT |
CL28668Contig1 | MYB113 | GAACCCATCCCAAACACAAC | GCTTCTTCACCGTTATCCAC |
CL20586Contig1 | S7AT | TCTTTCATCCGAGCCGTATC | AGTGAGGTGGAGTTGGAGAT |
CL17343Contig1 | Ace | CATCAAGCCAAAGACCTCTCTC | CCTATCAGCACCGAAACTACTG |
CL13518Contig1 | ODD | CAGAGCCAATAGTCCGAGTC | TCCAGTGGTGAAGGAGAAGG |
CL14366Contig1 | F3MO | GCTCACTCTTGTAGTCTCCTC | CGGGCTCCAAAGAATGTCTC |
CL29736Contig1 | LAR2 | ATCTACGGTGATGGCTCTGT | TGCGTCCAATCTTCTTCTCC |
CL24482Contig1 | ANR2 | GGAATCGACCTCTAAGTGACTC | TCCTACTTGCCTGGTTCTGT |
CL26328Contig1 | CcoA | GGAACCGTATCAATGTCAGAGG | TAGCACTCCTCATGGCACTT |
CL34455Contig1 | Actin1 | GCTGCTTATGGTCTGGTCAT | CACCCATTACTACTTCCCTCAG |
CL11169Contig1 | Actin2 | AGTCCTCTTCCAGCCATCTT | GAACCACCACTGAGCACAAT |
CL40530Contig1 | Actin3 | CTCCAAGGACAGATTCAAGAGG | GAGCAGCACGACAAGAGAAT |
CL8799Contig1 | Actin4 | GGACTCTGGTGATGGTGTGAG | GCTCTGCTGTAGTAGTGAAGG |
CL33133Contig1 | Actin5 | GTCGTCTTCTCCACTACCAT | ACACGGCTCTACATTACCAG |
CL31684Contig1 | RP | GCACAGCGATTCAGTTCTAG | CAACTACGGACAGAGCACTT |
CL31240Contig1 | UBE | CCCAAATCCTTCAGACCCAT | CCTCGTCACTTATCTCCTCATC |
CL33177Contig1 | UBQ | CTACTGGCAATCCTACTCCTAC | GGAAATAGAGACGGACCAAGAC |
CL31981Contig1 | NADH | ACACCTCGTCCAGCAATATC | CTCAACCGCCTCAACCTTAG |
CL27968Contig1 | ADK | GTGTCACTTGGCTACTGGAG | GGCTTCTTCAGTGCTTCATC |
comp20029_c0_seq3_2 | ACO1 | GTCCATAGCCTCGTTCTACAAC | GGTACTTTCCACAGCCTTCATG |
CL5504Contig1 | ACO2 | CCCAAACATTCAGACACCAC | GCTCCACGCTTATGAACTTG |
CL30130Contig1 | PG | ATGGAACTACGGTGATGAGC | GGTACGAGAGGAAGGAGGAG |
Table 2: Primers used in this study.
Figure 1: Melting curves of the reference genes based on RT-qPCR data.
Note: a: actin1; b: actin2; c: actin3; d: actin4; e: actin5; f: R; g: UBE; h: UBQ; i: NADH; J: ADK, The melting curves for each gene was generated on Bio-Rad CFX Manager according to instructions
Expression levels of the candidate reference genes
To evaluate the expression levels of the candidate reference genes, the cycle threshold (Ct) values of the ten candidate reference genes in green, red, and blue fruits were acquired with three replicates. The expression levels of these ten reference genes were dramatically varied, with Ct values ranging from 17.61 to 29.73 cycles. According to the Ct values, actin 2 expression levels in the three fruit samples were the highest, and actin 5 was the lowest (Figure 2). In the green fruits, all the gene expression levels ranged from actin2>actin3>NADH>actin4>UBE>ADK>UBQ>actin1>RP>actin5 (Table 3). In the red and blue fruits, the arrangement changed in the order actin2>actin3>NADH>actin4>ADK>UBE>UBQ>act in1>RP>actin5 (Table 3). According to these data, sswe also found that the expression levels of actin3, NADH, ADK, and UBE had narrow variances across the three fruit development stages, which convincingly suggested that their expression levels were stable in the three fruits (Figure 2). These data demonstrated that these genes could be used as reference genes.
actin1 | actin2 | actin3 | actin4 | actin5 | RP | UBE | UBQ | NADH | ADK | |
---|---|---|---|---|---|---|---|---|---|---|
S1 | 25.26 | 18.34 | 18.63 | 22.32 | 27.91 | 26.56 | 22.77 | 24.18 | 20.52 | 23.06 |
25.11 | 18.41 | 18.58 | 22.24 | 28.13 | 26.49 | 22.86 | 24.1 | 20.38 | 23.01 | |
25.32 | 18.25 | 18.72 | 22.21 | 27.68 | 26.61 | 22.72 | 24.13 | 20.4 | 23.02 | |
S2 | 26.03 | 18.91 | 19.04 | 23.24 | 29.49 | 29.04 | 23.4 | 25.41 | 20.61 | 23.27 |
25.91 | 18.81 | 19.13 | 23.16 | 29.53 | 28.07 | 23.36 | 25.36 | 20.61 | 23.29 | |
25.96 | 18.82 | 19.07 | 23.35 | 29.73 | 28.21 | 23.32 | 25.14 | 20.71 | 23.41 | |
S3 | 24.21 | 17.65 | 18.22 | 21.86 | 27.09 | 25.39 | 22.59 | 23.14 | 19.71 | 22.12 |
24.42 | 17.61 | 18.3 | 21.89 | 27.33 | 25.39 | 22.46 | 23.11 | 19.75 | 22.29 | |
24.35 | 17.64 | 18.38 | 21.91 | 27.33 | 25.4 | 22.52 | 23.33 | 19.77 | 22.27 |
Table 3: Ct values of ten reference genes in three fruit ripen stage.
Figure 2: Ct values variance among three fruit development stage.
RT-qPCR was performed using Bio-Rad CFX Manager, and Ct values of each reference genes in the samples (S1, S2 and S3) were acquired later. The picture were drawn by R, using these Ct values.
Expression stability of candidate reference genes
geNorm analysis: The average expression stability (M value) of the ten genes was calculated by geNorm (version 3.5). The M values of all ten genes across green, red, and blue fruits were lower than the threshold of 1.5, which means that all the selected genes could be used as reference genes (Table 1). However, the M value of these genes showed a great variability; the M value of the RP gene was the highest, while those of the NADH and ADK genes were the smallest, suggesting that they had the highest expression stability (Table 1). These results showed that NADH and ADK are the best reference genes, followed by actin2. To further assess the optimal number of reference genes, pairwise variation (Vn/Vn + 1) was analysed using geNorm, and 0.15 was used as a cut-off value according to a previous study. The results showed that all the Vn/ Vn+ 1 values were less than 0.15. Moreover, the pairwise variation (V2/V3) value was significantly less than 0.15 (Figure 3), suggesting that only two reference genes were necessary for the gene profile; thus, we surmised that no additional reference genes were required.
Figure 3: Pairwise variation (Vn/Vn+1) values for three fruit ripen stage.
Vn/Vn+1 > 0.15 indicates another (n+1) reference was required, while Vn/Vn+1 ≤ 0.15 demonstrates only n reference was enough.
NormFinder analysis
In addition to the geNorm analysis, NormFinder analysis is also a necessary method for confirming the best reference gene. The gene with the smallest S-value was the most suitable reference gene. According to the NormFinder analysis, the S-value of all the candidate genes was less than 1.0, which indicates that all the reference genes we evaluated could be used for RT–qPCR analysis. Of them, the S-value of the actin1 gene (0.03) was the smallest and that of RP (0.785) was the largest (Table 1), which indicates that the actin1 gene is the most stable reference gene, and the results were not uniform with geNorm analysis.
BestKeeper analysis
Finally, BestKeeper was used to evaluate the stability of the reference genes through the standard deviation (SD) and the coefficient of variance (CV) of the Ct values. If the SD value of the Ct values were >1, the genecould not be set as a reference gene. According to BestKeeper analysis, both Ct values of all the genes were smaller than 1.0, which indicates that all the genes we selected met the standard of BestKeeper. Among them, however, actin3, UBE, and NADH are the three genes with lowest SD values. Like the SD value, these three genes had smaller CV values(Table 1). However, UBQ, actin5, and RP both had higher SD values and CV values (Table 1). Thus, while these results demonstrate that all the ten genes we selected could be used as reference genes, the actin3, UBE, and NADH genes are the best choices overall.
Recommended comprehensive ranking
Since the analysis by the three methods is not uniform, we compared and ranked the reference genes by RefFinder, which is a comprehensive, web-based tool combined with geNorm, Normfiner, BestKeeper, and Ct values. The overall final ranking based on the ranking from each program is shown in Table 2. These results showed that actin2, ADK, actin1, and NADH were the four genes with smaller values, and the overall ranking values (actin2, 2.59; ADK, 2.83; actin1, 3.15; NADH, 3.20) were very close (Table 2). However, partial Ct values of actin1 in red fruit were over 26, suggesting that its expression levels are slightly low and may not be suitable for RT–qPCR analysis of genes in all stages of fruit development. Additionally, both geNorm and BestKeeper analysis showed that ADK and NADH are better than actin2, which are normally used as the main methods for evaluating the stability of reference genes (Tables 1 and 2). Considering all these results, we conclude that the two reference genes,ADK and NADH, are the best reference genes for RT–qPCR analysis across all stages of fruit development.
Identification of reference genes using marker genes in fruit development
Since ACO and PG genes are generally upregulated in fruit ripening programs, assessing their expression levels should be the best way to confirm the reliability of reference genes. The expression levels of ACO1, ACO2, and PG were calculated using NADH, ADK, NADH+ADK, RP and actin5. We found that the transcript levels of ACO genes normalized by RP were not in accordance with their expression features (Figure 4). Moreover, the transcript levels of PG normalized by actin5 in red fruit did not agree with the PG expression feature and the transcriptomic analysis (Figure 4). However, the results calculated by NADH and ADK are like the expression features of these genes and the transcriptomic analysis (Figure 4). Furthermore, NADH+ADK as the coreference gene is not consistent with NADH or ADK alone (Figure 4). Taken together, either NADH or ADK are the most suitable reference genes for RT–qPCR analysis in VBT fruit development.
Figure 4: Evaluation of reference genes using the mark genes involved in fruit ripeningS1 indicates the green fruits, S2 indicates the red fruits, S3
indicates the blue fruits. Relative expression levels ot the three genes were normalized to reference genes NADH, ADK, NADH+ADK, RP, and actin5,
respectively. Labels for each reference genes as normalizer were shown in the picture. Pictures were drawn using GraghPad Prism 7.0. Note: NADK+ADK;
NADH;
ADK;
Acrion%5;
RP.
Differential expressed genes from transcriptomic analysis validation using NADH
To further evaluate the reliability of the reference genes, the expression levels of 15 flavonoid synthesis-related genes were normalized to the expression level of NADH. These genes encode shikimate O-hydroxycinnamoyltransferase (SHT), flavonol synthase (FLS), flavonoid 3'-monooxygenase (F3M1), flavonoid 3',5'-hydroxylase 2 (F35H2), dihydroflavonol 4-reductase (DFR), leucoanthocyanidin dioxygenase (LDOX), leucoanthocyanidin reductase (LAR), anthocyanidin reductase (ANR), UDPglycosyltransferase 74B1 (UGT), caffeoyl-CoA O-methyltransferase (CcoA), probable caffeoyl-CoA O-methyltransferase (At4g), salutaridinol 7-O-acetyltransferase (S7AT), transcription factor MYB113 (MYB113), acetylajmalan esterase (Ace), and 2-Oxoglutarate-Dependent Dioxygenase (ODD). The results showed that the Pearson correlation coefficients of the 12 genes are above 0.9 (Figure 5, Table 3), which indicates that their RTqPCR results are highly consistent with the FPKM values in the transcriptome. In addition, the Pearson correlation coefficients of FLS and LAR2, and ANR2 arranged from 0.81 to 0.89, indicating a lower correlation. However, their expression patterns were also highly consistent with the transcriptomic data. Taken together, the reference gene NADH is highly reliable for RT–qPCR analysis in the fruit ripening stage.
Figure 5: Validation of the DEGs in three ripening stage using NADH.
Note: S1 indicates the green fruits, S2 indicates the red fruits, S3 indicates the blue fruits. The blue line indicates RT-qPCR data normalized by NADH, and the red lines represent FPKM values. All the details showed in the picture. Pictures were drawn in the Microsoft excel based the RT-qPCR data and FPKM value.
The VBT plant is well-known for producing “Wumifan” and is also used widely as an herbal medicine in China. However, the nutrients in its fruit have been ignored in previously reported studies. As a “wild-blueberry,” there are diverse and abundant amino acids, minerals, and unsaturated fatty acids in fruits, which are beneficial for human health[1]. Unlike blueberry, VBT plants are universally grown in the mountains in southern China, which suggests that this species is adapted the environment in China. Additionally, the fruit ripening stage is in autumn. Therefore, these features make VBT definitively specific to the Ericaceous family of plants. Blueberry are always ripened in early summer, when the rainy season has come, which dramatically impacts their production. Discovery of the genes controlling the fruit ripen time of VBT would allow their use in blueberry breeding, which would considerably benefit blueberry cultivation and production.
To comprehensively explore the VBT genome, we sequenced its genome and transcriptome. We found that no standard reference gene is suitable for RT–qPCR analysis. Reference genes, such as actin, tubulin, tef1, NADP, and ubiquitin, are widely used in different species [18,20,22-24]. However, appropriate reference genes are not uniform across diverse species, and even in different treatments or tissues of the same species [18-23]. We also attempted to use actin as the reference gene for validation of the expression levels according to the transcriptomic data. However, the RT–qPCR results are not stable, especially the Ct values, which are dramatically varied in various tissues. Therefore, reference screening is necessary for RT–qPCR analysis in VBT. To screen for the best reference genes in analysing the genes involved in fruit ripening of VBT, we obtained ten genes that had similar FPKM values in various ripening stages. The genes were analysed using GeNorm, NormFinder, BestKeeper. Thus, because the scores of a former reference gene that was evaluated by the three methods are not uniform, many researchers then selected genes with high scores in GeNorm and NormFinder, or selected genes using a combination of the three methods [18,20,21,23]. As a result, the best reference gene remained to be identified.
VBT is a unique species in China. Because of its high medicinal value and rich nutrient qualities, the molecular mechanisms involvd in nutrient and active biochemical component synthesis in VBT should be better understood [1,2,15,27]. Our group sequenced the genome and transcriptome of VBT. However, the validation of the differentially expressed genes in the transcriptome involved in different fruit˗development stages using the actin gene as a normalizer was not stable [23], and the results showed great differences in the same samples. The best reference gene in studying the development of fruits is not actin, which implies that it is not adaptable for validation of the genes involved in fruit development. Therefore, we obtained ten candidate reference genes from the transcriptomic data, that have similar gene expression levels in different fruit development stages. Furthermore, they wereevaluated by GeNorm, NormFinder, BestKeeper and a comprehensive ranking. The analysis showed that the comprehensive ranking values of several genes were remarkably close, suggesting that actin2, ADK, actin1 and NADH can be used as reference genes for gene expression analysis (Table 2). However, the Ct values for the actin1 gene were too high in various ripen fruit stages of VBT, and both geNorm and Bestkeeper analyses showed that ADK and NADH are better than actin2 (Table 1), which should not be beneficial for further analysis. Therefore, two genes, ADK and NADH, were selected for further analysis.
Ethylene is a key chemical involve in fruit development. During VBT fruitripening ethylene is continuously synthetized, accumulated, and released to the environment. To synthetize ethylene, the marker gene ACOs is dramatically expressed at the early stages of ripening, at the red fruit stage in blueberry, which means that these genes could be used to evaluate the reliability of reference genes. Therefore, we evaluated the expression levels of ACO genes using ADK, NADH, and NADH+ADK as reference genes. As a result, we found that both ADK and NADH are dependable for RT– qPCR analysis. Furthermore, the relative transcript levels of 15 DEGs predicted by transcriptomic analysis were also normalized to NADH based on RT–qPCR data. Pearson correlation coefficients between transcriptomic data and RT–qPCR data of the 12 genes normalized by the expression level of the NADH gene were above 0.9 in VBT fruit, which demonstrates the high reliability of the reference genes for characterizing gene expression levels in the fruit ripening stage.
VBT plants have great value for producing “Wumifan,” enhancing herbal medicine, and advancing blueberry breeding. To explore its molecular content, further studies on the VBT genome and proteome should be conducted. RT–qPCR is a major method for evaluating gene expression levels. However, no candidate reference genes had been comprehensively screened in VBT fruits. In this study, we assessed ten reference genes and found that NADH and ADK genes are the most reliable reference genes for RT– qPCR analysis across all stages of VBT fruit ripening. These results will benefit further studies on the function of VBT genes involved in fruit ripening and in fruit nutrient and bioactive substance production.
F.H. and J.X. conceived and designed this study. L.G. performed related experiments and data analyses. F.H. wrote the manuscript, and J.X. revised. The others all participated in revision of text and pictures, or data reanalysis.
This work was supported by grants from Key research and development projects of Anhui Province (202204c06020029), the Anhui Provincial Natural Science Foundation (2008085MC77 and 2008085QC129), the Major Science and Technology Projects in Anhui Province (No. 2020003a06020009), the Program of Fujian Key Laboratory for Monitoring and Management of Crop Pests (026052020008).
The authors declare that they have no competing interests.
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Crossref] [Google Scholar] [PubMed]
[Google Scholar] [PubMed]
Citation: He F, Gui L, Zhang Y, Zhu B, Zhang X, Shen M et al. (2022) Validation of Reference Genes for Gene Expression Analysis in Fruit Development of Vaccinium bracteatum Thunb using Quantitative Real-Time PCR. J Clin Trials. 12:501.
Received: 28-Jun-2022, Manuscript No. JCTR-22-18129; Editor assigned: 30-Jun-2022, Pre QC No. JCTR-22-18129 (PQ); Reviewed: 14-Jul-2022, QC No. JCTR-22-18129; Revised: 21-Jul-2022, Manuscript No. JCTR-22-18129 (R); Published: 28-Jul-2022 , DOI: 10.35248/2167-0870.22.12.501
Copyright: © 2022 He F, et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.