Selection of reference genes for gene expression studies in Siberian Apricot (Prunus sibirica L.) Germplasm using quantitative real-time PCR

Advancements in Reference Gene Selection for Fruit Trees: A Comprehensive Review

Real-time quantitative polymerase chain reaction (qRT-PCR) has been widely used in gene expression analyses due to its advantages of sensitivity, accuracy and high throughput. The stability of internal reference genes has progressively emerged as a major factor affecting the precision of qRT-PCR results. However, the stability of the expression of the reference genes needs to be determined further in different cells or organs, physiological and experimental conditions. Methods for evaluating these candidate internal reference genes have also evolved from simple single software evaluation to more reliable and accurate internal reference gene evaluation by combining different software tools in a comprehensive analysis. This study intends to provide a definitive reference for upcoming research that will be conducted on fruit trees. The primary focus of this review is to summarize the research progress in recent years regarding the selection and stability analysis of candidate reference genes for different fruit trees.

15-cis-Phytoene Desaturase and 15-cis-Phytoene Synthase Can Catalyze the Synthesis of β-Carotene and Influence the Color of Apricot Pulp

Fruit color affects its commercial value. β-carotene is the pigment that provides color for many fruits and vegetables. However, the molecular mechanism of β-carotene metabolism during apricot ripening is largely unknown. Here, we investigated whether β-carotene content affects apricot fruit color. First, the differences in β-carotene content between orange apricot 'JTY' and white apricot 'X15' during nine developmental stages (S1-S9) were compared. β-carotene contents highly significantly differed between 'JTY' and 'X15' from S5 (color transition stage) onwards. Whole-transcriptome analysis showed that the β-carotene synthesis genes 15-cis-phytoene desaturase (PaPDS) and 15-cis-phytoene synthase (PaPSY) significantly differed between the two cultivars during the color transition stage. There was a 5 bp deletion in exon 11 of PaPDS in 'X15', which led to early termination of amino acid translation. Gene overexpression and virus-induced silencing analysis showed that truncated PaPDS disrupted the β-carotene biosynthesis pathway in apricot pulp, resulting in decreased β-carotene content and a white phenotype. Furthermore, virus-induced silencing analysis showed that PaPSY was also a key gene in β-carotene biosynthesis. These findings provide new insights into the molecular regulation of apricot carotenoids and provide a theoretical reference for breeding new cultivars of apricot.

Metabolomic and transcriptomice analyses of flavonoid biosynthesis in apricot fruits

Introduction

Flavonoids, as secondary metabolites in plants, play important roles in many biological processes and responses to environmental factors.

Methods

Apricot fruits are rich in flavonoid compounds, and in this study, we performed a combined metabolomic and transcriptomic analysis of orange flesh (JN) and white flesh (ZS) apricot fruits.

Results and discussion

A total of 222 differentially accumulated flavonoids (DAFs) and 15855 differentially expressed genes (DEGs) involved in flavonoid biosynthesis were identified. The biosynthesis of flavonoids in apricot fruit may be regulated by 17 enzyme-encoding genes, namely PAL (2), 4CL (9), C4H (1), HCT (15), C3'H (4), CHS (2), CHI (3), F3H (1), F3'H (CYP75B1) (2), F3'5'H (4), DFR (4), LAR (1), FLS (3), ANS (9), ANR (2), UGT79B1 (6) and CYP81E (2). A structural gene-transcription factor (TF) correlation analysis yielded 3 TFs (2 bHLH, 1 MYB) highly correlated with 2 structural genes. In addition, we obtained 26 candidate genes involved in the biosynthesis of 8 differentially accumulated flavonoids metabolites in ZS by weighted gene coexpression network analysis. The candidate genes and transcription factors identified in this study will provide a highly valuable molecular basis for the in-depth study of flavonoid biosynthesis in apricot fruits.

Genome-wide identification and analysis of the WRKY gene family and low-temperature stress response in Prunus sibirica

BACKGROUND

WRKY transcription factors are a prominent gene family in plants, playing a crucial role in various biological processes including development, metabolism, defense, differentiation, and stress response. Although the WRKY gene family has been extensively studied and analysed in numerous plant species, research on Prunus sibirica's WRKY genes (PsWRKY) remains lacking.

RESULTS

This study analysed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, and Gene ontology (GO) annotation of PsWRKY gene family members using bioinformatics methods based on the whole-genome data of P. sibirica. In total, 55 WRKYs were identified in P. sibirica and were heterogeneously distributed on eight chromosomes. Based on the phylogenetic analysis, these WRKYs were classified into three major groups: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III. Members of different subfamilies have different cis-acting elements, conserved motifs, and intron-exon structures, indicating functional heterogeneity of the WRKY family. Prediction of subcellular localisation indicated that PsWRKYs were mainly located in the nucleus. Twenty pairs of duplicated genes were identified, and segmental duplication events may play an important role in PsWRKY gene family expansion. Analysis of the Ka/Ks ratio showed that the PsWRKY family's homologous genes were primarily purified by selection. Additionally, GO annotation analysis showed that the WRKY gene family was mainly involved in responses to stimuli, immune system processes, and reproductive processes. Furthermore, quantitative real-time PCR (qRT-PCR) analysis showed that 23 PsWRKYs were highly expressed in one or more tissues (pistils and roots) and PsWRKYs showed specific expression patterns under different low-temperature stress conditions.

CONCLUSIONS

Our results provide a scientific basis for the further exploration and functional validation of WRKYs in P. sibirica.

Comparative Anatomical and Transcriptomics Reveal the Larger Cell Size as a Major Contributor to Larger Fruit Size in Apricot

Fruit size is one of the essential quality traits and influences the economic value of apricots. To explore the underlying mechanisms of the formation of differences in fruit size in apricots, we performed a comparative analysis of anatomical and transcriptomics dynamics during fruit growth and development in two apricot cultivars with contrasting fruit sizes (large-fruit Prunus armeniaca 'Sungold' and small-fruit P. sibirica 'F43'). Our analysis identified that the difference in fruit size was mainly caused by the difference in cell size between the two apricot cultivars. Compared with 'F43', the transcriptional programs exhibited significant differences in 'Sungold', mainly in the cell expansion period. After analysis, key differentially expressed genes (DEGs) most likely to influence cell size were screened out, including genes involved in auxin signal transduction and cell wall loosening mechanisms. Furthermore, weighted gene co-expression network analysis (WGCNA) revealed that PRE6/bHLH was identified as a hub gene, which interacted with 1 TIR1, 3 AUX/IAAs, 4 SAURs, 3 EXPs, and 1 CEL. Hence, a total of 13 key candidate genes were identified as positive regulators of fruit size in apricots. The results provide new insights into the molecular basis of fruit size control and lay a foundation for future breeding and cultivation of larger fruits in apricot.

Mining for genes related to pistil abortion in Prunus sibirica L

In Prunus sibirica, the phenomenon of pistil abortion is very common and seriously affects its fruit quality and yield; however, the molecular mechanisms of pistil abortion remains unclear. In this study, we identified differentially expressed genes (DEGs) and pathways associated with pistil abortion using transcriptome sequencing. After comparative analysis, a total of 1,950 DEGs were identified, of which 1,000 were upregulated, and 950 were downregulated. Gene Ontology (GO) functional enrichment analysis of DEGs showed that metabolic process, cellular process, single-organism process, membrane, membrane part, cell, binding, catalytic activity, and transporter activity contained the largest number of DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the plant-pathogen interaction, starch and sucrose metabolism, and plant hormone signal transduction pathways contained the largest number of DEGs. The NAC, bHLH, and B3 transcription factor families contained the largest number of DEGs. qRT-PCR detection confirmed that the gene expression levels were consistent with the transcriptome sequencing results. This study provides a theoretical basis and scientific basis for further research on the molecular mechanisms of P. sibirica pistil abortion.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Tibetan Medicinal Plant Saussurea Laniceps Callus under Abiotic Stresses and Hormone Treatments

Real-time quantitative PCR (RT-qPCR) is an important technique for studying gene expression analysis, but accurate and reliable results depend on the use of a stable reference gene. This study proposes to test the expression stability of candidate reference genes in the callus of Saussurea laniceps, a unique Tibetan medicinal plant. Based on the S. laniceps callus transcriptome, eleven candidate reference genes, including TUA2, TUB3, TUB8, TIF3B1, TIF3H1, ELF5A, PP2AA2, UEV1D, UBL5, UBC36, and SKIP1), were validated for RT-qPCR normalization in the callus under abiotic stress (salt, cold, and UV) and hormone treatments (abscisic acid, MeJA, and salicylic acid). The stability of the candidate genes was evaluated in all the samples of S. laniceps. Comprehensive analysis of all samples showed that the best reference genes were UBC36 and UBL5. ELF5A and TIF3B1 were ranked as the most stable genes in the sample sets under abiotic stress. For hormone stimulation, UBC36 and TIF3H1 genes had the best stability. This study provides useful guidelines and a starting point for reference gene selection for expression analysis using RT-qPCR techniques in S. laniceps.

Effects of exogenous abscisic acid on oil content, fatty acid composition, biodiesel properties and lipid components in developing Siberian apricot (Prunus sibirica) seeds

Previous studies in Siberian apricot (Prunus sibirica) seed kernel (SASK) have suggested the involvement of abscisic acid (ABA) signaling pathway in oil accumulation. However, there are few reports on the effects of ABA on the metabolism of fatty acids (FA) in seed development. Here, we first evaluated the response of developing SASK to ABA treatment, with a focus on oil content, FA composition, biodiesel properties, lipid compounds and gene expressions. Compared with control samples, the application of exogenous ABA increased the total oil content by 6.55% in mature SASK. The C18:1 content markedly increased in ABA treatment, and conversely C16:0 decreased. Exogenous ABA also improved the biodiesel properties of SASK oil, making it better suited to the specifications of biodiesel standards. Furthermore, the molecular species of phosphatidylcholine (PC), phosphatidic acid (PA), diacylglycerol (DAG) and triacylglycerol (TAG) were detected using lipidomics analysis. The 18:1/18:1 was the main component in PA, PC and DAG, while the main components of 18:1/18:1/18:2, 18:1/18:1/18:3, 18:2/18:2/18:2 and 18:1/18:1/18:1 in TAG. Most lipid species gradually increased with SASK maturity. In addition, the relative contents of TAG-18:1/18:1/18:2 and TAG-18:1/18:1/18:1 in developing SASK increased with the application of exogenous ABA. We also detected elevated gene expression of key genes involved in ABA chemical pathway, which likely affected FA biosynthesis and accumulation. Our results provide insight into the effects of ABA on the oil accumulation in developing SASK, which has direct applications to improving the quality of SASK-derived biodiesel.

Analysis of Metabolites and Gene Expression Changes Relative to Apricot (Prunus armeniaca L.) Fruit Quality During Development and Ripening

Apricot (Prunus armeniaca L.) is a valuable worldwide agronomical crop, with a delicious fruit highlighted as a functional food with both nutritional and bioactive properties, remarkably beneficial to human health. Apricot fruit ripening is a coordinated developmental process which requires change in the expression of hundreds to thousands of genes to modify many biochemical and physiological processes arising from quality characteristics in ripe fruit. In addition, enhancing fruit and nutraceutical quality is one of the central objectives to be improved in the new varieties developed by breeding programs. In this study we analyzed the contents of main metabolites linked to the nutraceutical value of apricot fruits, together with the most important pomological characteristics and biochemical contents of fruit during the ripening process in two contrasted apricot genotypes. Additionally, the gene expression changes were analyzed using RNA-Seq and real time qPCR. Results showed that genes with differential expression in the biosynthetic pathways, such as phenylpropanoids, flavonoids, starch and sucrose and carotenoid metabolism, could be possible candidates as molecular markers of fruit quality characteristics for fruit color and soluble solid content. The gene involves in carotenoid metabolism carotenoid cleavage dioxygenase 4, and the gene sucrose synthase in starch and sucrose metabolism were identified as candidate genes in the ripening process for white skin ground color and flesh color and high soluble sugar content. The application of these candidate genes on marker-assisted selection in apricot breeding programs may contribute to the early selection of high-quality fruit genotypes with suitable nutraceutical values.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

WRINKLED1 transcription factor orchestrates the regulation of carbon partitioning for C18:1 (oleic acid) accumulation in Siberian apricot kernel

WRINKLED1 (WRI1), an APETALA2 (AP2)-type transcription factor, has been shown to be required for the regulation of carbon partitioning into fatty acid (FA) synthesis in plant seeds. To our knowledge, the regulatory network of WRI1 remains unknown in Prunus sibirica kernel (PSK), a novel woody biodiesel feedstock in China. In this study, based on the transcriptional data from developing oilseeds of multiple plant species, we identified 161 WRI1-coexpressed genes using weighted gene co-expression network analysis (WGCNA). The major portion of WRI1-coexpressed genes was characterized to be involved in carbon partitioning and FA biosynthesis. Additionally, we detected the temporal patterns for oil content and FA compositions in developing PSK from two different germplasms (AS-85 and AS-86). The major differences between the two germplasms are higher contents of oil and C18:1 in AS-85 than in AS-86 at a mature stage. Thus, AS-85 and AS-86 are desirable materials to explore the molecular and metabolic mechanisms of oil accumulation in Siberian apricot. Expression analysis in developing PSK of AS-85 and AS-86 indicated that the expression level of P. sibirica WRI1 (PsWRI1) was closely correlated to accumulative rate of oil. Also, the comparison of expression profiles in developing PSK of AS-85 and AS-86 displayed that the pPK, E1-α, E2, TAL, BC, MCMT, BS, SAD and FAD2 have a high correlation with PsWRI1. Transient expression showed that Pro_SAD- and Pro_BS-driving GUS expression showed no substantial difference between AS-85 and AS-86, while the expression level of Pro_PEPCK-AS-85 driving GUS was significantly higher than that of Pro_PEPCK-AS-86 driving GUS. Additionally, transient co-transformation with PsWRI1 revealed that Pro_SAD, Pro_PEPCK and Pro_BS activity could be specifically up-regulated by PsWRI1. This regulatory mechanism of PsWRI1 may create a steep concentration difference, thereby facilitating carbon flux into C18:1 accumulation in developing PSK. Overall, all our findings imply a versatile mechanism of WRI1 to optimize carbon allocation for oil accumulation, which can provide reference for researching the woody biodiesel plants.

Comprehensive evaluation of fuel properties and complex regulation of intracellular transporters for high oil production in developing seeds of Prunus sibirica for woody biodiesel

BACKGROUND

Based on our previous studies of 17 Prunus sibirica germplasms, one plus tree with high quality and quantity of seed oils has emerged as novel potential source of biodiesel. To better develop P. sibirica seed oils as woody biodiesel, a concurrent exploration of oil content, FA composition, biodiesel yield and fuel properties as well as prediction model construction for fuel properties was conducted on developing seeds to determine the optimal seed harvest time for producing high-quality biodiesel. Oil synthesis required supply of carbon source, energy and FA, but their transport mechanisms still remains enigmatic. Our recent 454 sequencing of P. sibirica could provide long-read sequences to identify membrane transporters for a better understanding of regulatory mechanism for high oil production in developing seeds.

RESULTS

To better develop the seed oils of P. sibirica as woody biodiesel, we firstly focused on a temporal and comparative evaluation of growth tendency, oil content, FA composition, biodiesel yield and fuel properties as well as model construction for biodiesel property prediction in different developing seeds from P. sibirica plus tree (accession AS-80), revealing that the oils from developing seeds harvested after 60 days after flowering (DAF) could be as novel potential feedstock for producing biodiesel with ideal fuel property. To gain new insight into membrane transport mechanism for high oil yield in developing seeds of P. sibirica, we presented a global analysis of transporter based on our recent 454 sequencing data of P. sibirica. We annotated a total of 116 genes for membrane-localized transporters at different organelles (plastid, endoplasmatic reticulum, tonoplast, mitochondria and peroxisome), of which some specific transporters were identified to be involved in carbon allocation, metabolite transport and energy supply for oil synthesis by both RT-PCR and qRT-PCR. Importantly, the transporter-mediated model was well established for high oil synthesis in developing P. sibirica seeds. Our findings could help to reveal molecular mechanism of increased oil production and may also present strategies for engineering oil accumulation in oilseed plants.

CONCLUSIONS

This study presents a temporal and comparative evaluation of developing P. sibirica seed oils as a potential feedstock for producing high-quality biodiesel and a global identification for membrane transporters was to gain better insights into regulatory mechanism of high oil production in developing seeds of P. sibirica. Our findings may present strategies for developing woody biodiesel resources and engineering oil accumulation.

Identification of AUXIN RESPONSE FACTOR gene family from Prunus sibirica and its expression analysis during mesocarp and kernel development

BACKGROUND

Auxin response factors (ARFs) in auxin signaling pathway are an important component that can regulate the transcription of auxin-responsive genes involved in almost all aspects of plant growth and development. To our knowledge, the comprehensive and systematic characterization of ARF genes has never been reported in Prunus sibirica, a novel woody biodiesel feedstock in China.

RESULTS

In this study, we identified 14 PsARF genes with a perfect open reading frame (ORF) in P. sibirica by using its previous transcriptomic data. Conserved motif analysis showed that all identified PsARF proteins had typical DNA-binding and ARF domain, but 5 members (PsARF3, 8 10, 16 and 17) lacked the dimerization domain. Phylogenetic analysis of the ARF proteins generated from various plant species indicated that ARFs could be categorized into 4 major groups (Class I, II, III and IV), in which all identified ARFs from P. sibirica showed a closest relationship with those from P. mume. Comparison of the expression profiles of 14 PsARF genes in different developmental stages of Siberian apricot mesocarp (SAM) and kernel (SAK) reflected distinct temporal or spatial expression patterns for PsARF genes. Additionally, based on the expressed data from fruit and seed development of multiple plant species, we identified 1514 ARF-correlated genes using weighted gene co-expression network analysis (WGCNA). And the major portion of ARF-correlated gene was characterized to be involved in protein, nucleic acid and carbohydrate metabolic, transport and regulatory processes.

CONCLUSIONS

In summary, we systematically and comprehensively analyzed the structure, expression pattern and co-expression network of ARF gene family in P. sibirica. All our findings provide theoretical foundation for the PsARF gene family and will pave the way for elucidating the precise role of PsARF genes in SAM and SAK development.

Expression Stabilities of Candidate Reference Genes for RT-qPCR in Chinese Jujube (Ziziphus jujuba Mill.) under a Variety of Conditions

Reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) is a powerful method for evaluating patterns of gene expression. Jujube whole-genome sequencing has been completed, and analysis of gene function, an important part of any follow-up study, requires the appropriate selection of reference genes. Indeed, suitable reference gene selection for RT-qPCR is critical for accurate normalization of target gene expression. In this study, the software packages geNorm and NormFinder were employed to examine the expression stabilities of nine candidate reference genes under a variety of conditions. Actin-depolymerizing factor 1 (ACT1), Histone-H3 (His3), and Polyadenylate-binding protein-interacting protein (PAIP) were determined to be the most stably expressed genes during five stages of fruit development and ACT1, SiR-Fd, BTF3, and Tubulin alpha chain (TUA) across different tissues/organs. Whereas ACT1, Basic Transcription factor 3 (BTF3), Glyceraldehyde-3-phosphate dehydrogenase (GADPH), and PAIP were the most stable under dark conditions. ACT1, PAIP, BTF3, and Elongation factor 1- gamma (EF1γ) were the most stably expressed genes under phytoplasma infection. Among these genes, SiR-Fd and PAIP are here first reported as stable reference genes. When normalized using these most stable reference genes, the expression patterns of four target genes were found to be in accordance with physiological data, indicating that the reference genes selected in our study are suitable for use in such analyses. This study provides appropriate reference genes and corresponding primers for further RT-qPCR studies in Chinese jujube and emphasizes the importance of validating reference genes for gene expression analysis under variable experimental conditions.

Evaluation and selection of candidate reference genes for normalization of quantitative RT-PCR in Withania somnifera (L.) Dunal

Quantitative real-time PCR (qRT-PCR) is now globally used for accurate analysis of transcripts levels in plants. For reliable quantification of transcripts, identification of the best reference genes is a prerequisite in qRT-PCR analysis. Recently, Withania somnifera has attracted lot of attention due to its immense therapeutic potential. At present, biotechnological intervention for the improvement of this plant is being seriously pursued. In this background, it is important to have comprehensive studies on finding suitable reference genes for this high valued medicinal plant. In the present study, 11 candidate genes were evaluated for their expression stability under biotic (fungal disease), abiotic (wounding, salt, drought, heat and cold) stresses, in different plant tissues and in response to various plant growth regulators (methyl jasmonate, salicylic acid, abscisic acid). The data as analyzed by various software packages (geNorm, NormFinder, Bestkeeper and ΔCt method) suggested that cyclophilin (CYP) is a most stable gene under wounding, heat, methyl jasmonate, different tissues and all stress conditions. T-SAND was found to be a best reference gene for salt and salicylic acid (SA) treated samples, while 26S ribosomal RNA (26S), ubiquitin (UBQ) and beta-tubulin (TUB) were the most stably expressed genes under drought, biotic and cold treatment respectively. For abscisic acid (ABA) treated samples 18S-rRNA was found to stably expressed gene. Finally, the relative expression level of the three genes involved in the withanolide biosynthetic pathway was detected to validate the selection of reliable reference genes. The present work will significantly contribute to gene analysis studies in W. somnifera and facilitate in improving the quality of gene expression data in this plant as well as and other related plant species.

Transcriptomic analysis revealed the mechanism of oil dynamic accumulation during developing Siberian apricot (Prunus sibirica L.) seed kernels for the development of woody biodiesel

BACKGROUND

Siberian apricot (Prunus sibirica L.) has emerged as a novel potential source of biodiesel in China, but the molecular regulatory mechanism of oil accumulation in Siberian apricot seed kernels (SASK) is still unknown at present. To better develop SASK oil as woody biodiesel, it is essential to profile transcriptome and to identify the full repertoire of potential unigenes involved in the formation and accumulation of oil SASK during the different developing stages.

RESULTS

We firstly detected the temporal patterns for oil content and fatty acid (FA) compositions of SASK in 7 different developing stages. The best time for obtaining the high quality and quantity of SASK oil was characterized at 60 days after flowering (DAF), and the representative periods (10, 30, 50, 60, and 70 DAF) were selected for transcriptomic analysis. By Illumina/Solexa sequencings, approximately 65 million short reads (average length = 96 bp) were obtained, and then assembled into 124,070 unigenes by Trinity strategy (mean size = 829.62 bp). A total of 3,000, 2,781, 2,620, and 2,675 differentially expressed unigenes were identified at 30, 50, 60, and 70 DAF (10 DAF as the control) by DESeq method, respectively. The relationship between the unigene transcriptional profiles and the oil dynamic patterns in developing SASK was comparatively analyzed, and the specific unigenes encoding some known enzymes and transcription factors involved in acetyl-coenzyme A (acetyl-CoA) formation and oil accumulation were determined. Additionally, 5 key metabolic genes implicated in SASK oil accumulation were experimentally validated by quantitative real-time PCR (qRT-PCR). Our findings could help to construction of oil accumulated pathway and to elucidate the molecular regulatory mechanism of increased oil production in developing SASK.

CONCLUSIONS

This is the first study of oil temporal patterns, transcriptome sequencings, and differential profiles in developing SASK. All our results will serve as the important foundation to further deeply explore the regulatory mechanism of SASK high-quality oil accumulation, and may also provide some reference for researching the woody biodiesel plants.

Selection of reliable reference genes for quantitative real-time PCR gene expression analysis in Jute (Corchorus capsularis) under stress treatments

To accurately measure gene expression using quantitative reverse transcription PCR (qRT-PCR), reliable reference gene(s) are required for data normalization. Corchorus capsularis, an annual herbaceous fiber crop with predominant biodegradability and renewability, has not been investigated for the stability of reference genes with qRT-PCR. In this study, 11 candidate reference genes were selected and their expression levels were assessed using qRT-PCR. To account for the influence of experimental approach and tissue type, 22 different jute samples were selected from abiotic and biotic stress conditions as well as three different tissue types. The stability of the candidate reference genes was evaluated using geNorm, NormFinder, and BestKeeper programs, and the comprehensive rankings of gene stability were generated by aggregate analysis. For the biotic stress and NaCl stress subsets, ACT7 and RAN were suitable as stable reference genes for gene expression normalization. For the PEG stress subset, UBC, and DnaJ were sufficient for accurate normalization. For the tissues subset, four reference genes TUBβ, UBI, EF1α, and RAN were sufficient for accurate normalization. The selected genes were further validated by comparing expression profiles of WRKY15 in various samples, and two stable reference genes were recommended for accurate normalization of qRT-PCR data. Our results provide researchers with appropriate reference genes for qRT-PCR in C. capsularis, and will facilitate gene expression study under these conditions.