Identification of best housekeeping gene for normalization in qRT-PCR data under different development stages and stress conditions in Vigna mungo

Identification of best housekeeping gene for normalization in qRT-PCR data under different development stages and stress conditions in Vigna mungo | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Identification of best housekeeping gene for normalization in qRT-PCR data under different development stages and stress conditions in Vigna mungo Sudhir Kumar Jha, Pritee Sagar, Prateek Kumar, Shivangi Gupta, and 7 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5819031/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 26 Sep, 2025 Read the published version in Scientific Reports → Version 1 posted 8 You are reading this latest preprint version Abstract The qRT-PCR has emerged as the most sought-after local gene expression profiling technique. The absolute quantification of the mRNA is highly affected by the handling artifacts. In this regard, the internal control genes serve a vital role in negating the handling errors. So far, the most stable housekeeping genes to be used for the normalization of qRT-PCR data are lacking in V. mungo , commonly known as urdbean or blackgram. Therefore, in the current study we intend to explore the stability of the UFO, TUB2, ACT2, RHA, ISTL6, CyP, DNAJ, RBCS3B, RPS34, USP21, USB1, EfTu, UBQ10, and UBQE9 genes of V. mungo in 17 different developmental stages and 4 different abiotic stresses. The well-acknowledged algorithms such as geNorm, NormFinder, BestKeeper, the comparative ΔCt method and comprehensive RefFinder were used to assess the expression stability. According to the comprehensive ranking, RPS34 and RHA were identified as the most appropriate genes for qRT-PCR data normalization throughout all developmental stages, whereas ACT2 and RPS34 were optimal under abiotic stress conditions. Further, the suitability of these three candidate housekeeping genes for normalization has been demonstrated in an independent experiment. Based on the findings obtained, the present research advocates the use of a combination of HKGs (RPS34 and RHA for developing tissues and ACT2 and RPS34 for stress-treated tissues) for the normalization of qRT-PCR data. The proposed reference genes will guarantee the reliability and consistency of qRT-PCR results. Biological sciences/Biological techniques Biological sciences/Biotechnology Biological sciences/Plant sciences Housekeeping gene Normalization qRT-PCR urdbean Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Figure 10 Key Message The gene stability ranking indicated that RPS34 and RHA were best for qRT-PCR data normalization in various developmental stages, whereas ACT2 and RPS34 under abiotic stress. Introduction With the advent of sequencing technologies, abundant information and genomics resources are being developed on a daily basis in several legumes, which were otherwise considered as “orphan crops.” To date, genome sequences for several legumes, including Cajanus cajan , Cicer arietinum , Vigna mungo , and Cyamopsis tetragonoloba have been published [ 1 – 4 ]. Parallelly, several of these legume crops have been supplemented with multiple transcriptomics experiments to mine genes related to any development stage and their biotic-abiotic response and have been reviewed [ 5 – 7 ]. Gene expression profiling, either ‘global’ or 'local', is becoming vital in any advanced molecular biology and molecular breeding experiments. Even the sole transcriptomics experiments require validation of key genes to make the results more reliable and statistically more robust. The expression of any gene is spatiotemporal and dynamic, as it is influenced by several internal as well as external factors. Quantification of mRNA is essential for understanding the gene expression levels in the biological samples. Various techniques, including quantitative reverse transcription polymerase chain reaction (qRT-PCR), northern blotting, RNA-seq, micro-array, in-situ hybridization, the Nanostring nCounter system, and the Serial Analysis of Gene Expression (SAGE), have been developed for this purpose, each with its own advantages and limitations. Among these techniques, qRT-PCR is the most widely used technique due to its robustness and accuracy. For accurate mRNA quantification, the qRT-PCR results are normalized using a reference gene/housekeeping gene (HKG) whose expression is stable across different experimental conditions, allowing for the quantification of changes that are due to biological differences rather than technical artifacts [ 8 – 10 ]. Urdbean, also known as blackgram ( V. mungo ), is the third most important grain legume after chickpea and pigeonpea in the Indian subcontinent. It is a nutrient-rich pulse possessing a rich blend of protein, fiber, and essential micronutrients like iron, zinc and potassium. Its adaptability to diverse agro-ecologies, shorter cropping duration, and ability to fix atmospheric nitrogen attract the small and marginal farmers [ 11 ]. As the world grapples with the burning issues like climate change, ever increasing population and nutritional security, the importance of the urdbean as a nutrient-rich and agriculturally important crop is being recognized recently. Despite this significance, it has long been overlooked by mainstream research and is treated as an “orphan crop”. Recently, efforts are being made to develop some genetic and genomic resources in urdbean. The availability of draft genome assembly, high density linkage map, along with the development of molecular markers such as SSRs and SNPs has enabled the identification of genes/QTLs linked with different traits. These genomic resources hold great promise for accelerating the urdbean improvement using the advanced genomics strategies [ 11 ]. Once the putative candidate genes are known their expression analysis needs to be performed, which requires a suitable HKG for the data normalization. Several HKGs have been explored for normalization of the qRT-PCR data. In general, the HKGs used across the studies belong to central processes that are critical for the survival of the cell including, genes essential for the cell morphology, cell division, growth and development. Because these aforementioned genes are tied to fundamental and non-specialized processes, they are less likely to fluctuate in response to experimental variables that might alter the expression of genes involved in specific pathways. However, their stability should be validated for specific experiments, as no single gene is universally optimal for all situations. Therefore, in this work, we evaluated 14 genes of different functional classes (significantly reducing the chance that the genes will be co-regulated) using the qRT-PCR techniques in multiple stages, starting from cotyledon (germination) to pod (maturity) stages. The most suitable HKG for qRT-PCR in V. mungo was evaluated employing different algorithms developed to assess the stability of a normalizing gene through NormFinder, geNorm, BestKeeper, and the ΔCt method implemented in RefFinder software [ 12 – 16 ]. The most potential HKGs identified in the current study were further validated to demonstrate their suitability in normalizing the qRT-PCR data. The most stable HKGs identified in the current study will be highly useful for the researcher working on the advanced genomics experiments in Vigna species, including V. mungo . Materials and methods Plant materials In the current work, an urdbean variety, “IPU13-1” released from the ICAR-Indian Institute of Pulses Research, Kanpur, was used as plant material. The variety has a potential yield of 10–12 q/ha with a maturation duration of 70–75 days. The seeds were sterilized for 30 seconds in 70% ethanol before sowing in clay pots filled with autoclaved soil. Four developmental stages, including germination [5 days after sowing (DAS)], the vegetative stage (30 DAS), the reproductive stage (45 DAS), and podding (60 DAS), were selected for the gene expression analysis. From the germination stage, plumule, radicle, and cotyledon tissues were harvested. While leaves, stems, and root tissue were harvested from all remaining three stages, including the vegetative stage, reproductive stage, and podding stage. Likewise, vegetative meristem, floral bud, and 15-day-old pods were harvested at the vegetative, reproductive, and podding stages, respectively. Further, the pods were also split into pod walls and seeds. All tissues were freshly harvested and snap-frozen into liquid nitrogen before storing them at -80°C until RNA isolation. All the samples were collected in three independent biological replicates. Treatments Four different treatments, including drought, salt, aluminium, and cold stress, were given to the sterilized seeds of IPU13-1 at the germination stage. The experiment for screening in salt, drought, and aluminium stress conditions was done in the customized growth chamber of REMI (CHM 12+) with a Floutronix grow light. The temperature was set at 35°C ± 2°C with relative humidity set to 65% and a light/dark period of 12 hours each. For salinity and aluminium stress, seeds were germinated in germination paper put in a petri-plate, followed by soaking with 120 mM NaCl and a 100 mM solution of Al 2 (SO 4 ) 3 , respectively [ 17 , 18 ]. Likewise, the sterilized seedlings were germinated in MS media supplemented with 1% w/v of PEG 6000 for imposing drought stress [ 19 ]. Similarly, to collect tissues under cold stress, seeds were germinated in germination paper with the temperature set to 18°C ± 2°C with relative humidity set to 65% and a light/dark period of 12 hours each [ 20 ]. Tissues of plumule, radicle, and cotyledon from each of the aforementioned treatments were collected at 15 DAS in three independent biological replicates for the downstream analysis. Total RNA extraction and cDNA synthesis RNA isolation for each sample was done using the RNeasy Plant Mini Kit, Qiagen, Germany (Catalogue no. 74904), following the manufacturer’s protocol with slight modifications [ 21 ]. The DNase treatment was given to the isolated RNAs to remove the genomic DNA contamination. Total RNA integrity was checked in 2% agarose gel electrophoresis and the amount was quantified using a spectrophotometer (NanoDrop 2000c, Thermo-Scientific, USA). From the total quantity of extracted RNA, 1 µg was used for cDNA synthesis using Maxima H Minus Double-Stranded cDNA Synthesis Kit, Thermo-Scientific, USA (Catalogue no. K1612), following the manufacturer’s instructions. Finally, the cDNA was diluted with molecular-grade water in a ratio of 1:10 and used for qRT-PCR analysis [ 21 ]. Selection of housekeeping genes and primer synthesis A total of 14 housekeeping genes were shortlisted based on the literature survey. The majority of these genes have not been tested in urdbean to date. The identifier for the selected genes was used to fetch the Arabidopsis protein sequences for the selected genes from TAIR ( https://www.arabidopsis.org/ ). The reference protein sequences for the Urdbean genome were downloaded from NCBI ( https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_036885675.1/ ). Further, the NCBI blastp suit ( https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins ) was used for finding the urdbean homologs of the selected Arabidopsis protein sequences based on the similarity. The Urdbean protein with the maximum similarity for each Arabidopsis protein sequence was shortlisted for the downstream analysis. The protein sequences were further confirmed for functional conservation by analysis of conserved domains in both Arabidopsis and Vigna mungo homologs via NCBI-Batch CDD ( https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi ; [ 22 ]). Primers from the selected genes were designed using the IDT's PrimerQuest tool ( https://www.idtdna.com/PrimerQuest ). PrimerQuest was chosen due to its robust primer design algorithms, which incorporate melting temperature, GC content, and secondary structure considerations. Further, the oligonucleotides were synthesized by Europhin Genomics, India. The primer sequences for these genes are provided in the supplementary file (Table S1 ). The qRT-PCR analysis of the selected HKGs Temperature gradient PCR was performed to optimize the specific Tm of each primer. Each 20 µL qRT-PCR reaction consisted of 10 µL master mix (Brilliant III Ultra-fast SYBR Green qRT-PCR Master Mix, Agilent Technologies), 0.3 µM forward and reverse primers, 7.4 µL PCR grade water and 50 ng cDNA as the template. The qRT-PCR reactions were carried out in an Agilent AriaMx Real-Time qRT-PCR system. The PCR conditions involved an initial denaturation at 94°C for 3 min, followed by 40 cycles of 95°C for 30 s, 60°C for 15 s, and 72°C for 30 s [ 21 ]. The amplicons were separated in 1.2% agarose gel using TBE buffer. Expression stability analysis For qRT-PCR data, raw Ct values (quantification cycle) were used for the analysis, and the mean Ct value of each biological replicate was considered to ensure the accuracy. For the expression stability analysis, the RefFinder web server ( https://www.ciidirsinaloa.com.mx/RefFinder-master/ ; [ 16 ]) was employed. The RefFinder tool comprehensively utilizes several other gene stability analysis algorithms such as geNorm [ 13 ], Normfinder [ 12 ], BestKeeper [ 14 ], and the comparative ΔCt method [ 15 ]. The RefFinder tool utilizes only the original Ct value from quantitative real-time PCR (qRT-PCR) for the downstream analysis [ 16 ]. Accordingly, the mean values of each replicate were used for the analysis using the RefFinder tool. The geNorm algorithms calculate the gene expression stability measure (M-value) based on pairwise variation analysis and also determine the optimal number of reference genes required for the normalization. Whereas, the NormFinder algorithm employs an ANOVA-based model to estimate intra- and inter-group variation, providing a direct stability value for the HKGs tested. Likewise, the BestKeeper algorithm assesses the stability based on the standard deviation (SD) and coefficient of variation (CV) of the Ct values. The HKGs with SD > 1 are considered unreliable for normalization. Similarly, the ΔCt method compares the relative expression stability of the HKGs by analyzing the variation in their Ct values. The RefFinder algorithm comprehends an overall ranking based on the geometric mean of rankings from these four methods. The statistical significance of candidate reference gene stability was determined based on rankings derived from the combined analysis. Further, to plot the mean stability value (geNorm), stability value (NormFinder), standard deviations (BestKeeper), arithmetic mean of the standard deviation (ΔCt method), and the geometric mean of the stability ranking (RefFinder) WPS office was used. The box-plot representing the distribution of the Ct values was also plotted using the GraphPad v10.4.1. Gene normalizing efficiency of the candidate HKGs Further, to validate the reliability of the most suitable HKGs (RPS34, RHA, and ACT2) identified in the current study, they were used as an internal control in the expression analysis of the 5 test genes previously reported by our group (Table S2 ) [ 23 ]. Parallelly, an external Actin (ACT) gene that has been previously reported and used in urdbean qRT-PCR was also used as an internal control. The control and infected leaf samples from Yellow Vein Mosaic (YMV) tolerant (DPU88-31) and susceptible (LBG17) genotypes were collected as previously reported by our group [ 23 ]. Further, the RNA isolation, cDNA library preparation, and qRT-PCR were performed using the aforementioned methodologies. The mean Ct values were calculated in Microsoft Excel, 2016 using all 4 candidate HKGs, and a one-way ANOVA followed by a Tukey's multiple comparisons test was performed at a significance level of p ≤ 0.05. Results Vigna mungo is emerging as an important legume crop after chickpea and pigeonpea due to its high nutritive value, short cropping duration, and multiple cuisines made from it. It is an important blend in the vegetarian diets for the majority of the population in South Asia due to its high nutritive value [11]. Besides, legume crops are being targeted to achieve the second sustainable development goal of “Zero Hunger” [24]; urdbean can be a vital component in this battle. The Vigna species are highly versatile and can be accommodated to different cropping systems and environments, making them a good choice for advanced genomics experiments. Looking into the wide adaptability as well as utility of the Vigna mungo , it has attracted the farmers as well as researchers. The recent enrichment in genetic and genomic resources of urdbean has paved the way for the functional characterization of many candidate genes. And gene expression profiling is a key experiment in this regard, which further depends on HKGs for the normalization to get statistically robust results, taking care of the technical artifacts. Therefore, to identify the most suitable HKGs for normalization of qRT-PCR data in urdbean, we have conducted a comprehensive expression profiling of 14 candidate housekeeping genes in 17 different tissues. Identification of the V. mungo orthologs of selected Arabidopsis HKG s A total of 14 V. mungo candidate genes corresponding to each selected Arabidopsis gene were identified through NCBI blastp with the “max hit = 1” parameter. These 14 genes were considered for local gene expression analysis using the qRT-PCR technique in the 17 selected tissues (Table 1). The percent similarity between the gene orthologs of urdbean and Arabidopsis ranges between 47.337 and 99.738, among which 12 out of 14 have a percent similarity of over 65 percent, reflecting the sequence conservation (Table 1). All of the 14 genes were found to be annotated as hypothetical proteins, as the gene prediction pipelines are based on homology-based annotations but lack experimental validation. Since the genome of V. mungo was recently sequenced [3], this could be explained. The conserved domain analysis using NCBI-Batch CDD confirms the presence of the same domain as it was present in Arabidopsis ortholog (Figure S1). The higher percent similarity between the V. mungo HKGs and that of the Arabidopsis orthologs, as well as conserved domain analysis, advocates the functional as well as sequence conservation between them (Figure S1). Sampling of diverse plant tissues A total of 17 plant tissues representing different stages from germination to pod development were selected for the expression analysis. The tissues used in the current study capture all developmental stages in urdbean. The diversity in the type and stages of the tissues reinforces the utility of the HKGs identified in the current study to the urdbean research community (Figure 1). RNA quality assessment and p rimer specificity The quality of the RNA used in this work was evaluated to guarantee a proper estimate of gene expression. The 260/280 absorbance ratio of RNA samples was evaluated using the spectrophotometer and only the samples with absorbance between 1.9 and 2.1 were considered. This ensures the quality as well as purity of the RNA isolated. Sharp and clear bands matching 18S rRNA and 28S rRNA on the denaturing gel revealed intact, undamaged RNA molecules, therefore indicating strong RNA integrity (Figure 2a). Further, the temperature gradient PCR for each primer set was conducted to identify the optimal annealing temperature. The length of the amplified PCR product was verified using 2% agarose gel electrophoresis using a DNA ladder (100 bp) as a reference. In our study, the amplicon size was found to range between 85 bp and 260 bp (Figure 2b). Gene expression stability analysis in different developmental stages The qRT-PCR analysis of all selected 14 genes was conducted across the 17 tissues, representing the germination, vegetative, reproductive and pod development stages. The mean Ct (quantification cycle) values of all 14 genes from three biological groups were utilized for the estimation of stability across the three datasets (Figure 3). The candidate HKGs showed a large variation in mean Ct ranging from 10.84 (RBCS3B) to 39.82 (ISTL6 and CyP) while considering the total dataset (Figure 3). The whiskers represent the quartile where the Ct values for most of the HKGs fall, while the middle line represents the middle (median) Ct value in the dataset. The broader whisker means more variability in Ct value, while the smaller whisker supports the stable gene expression across the tissues (figure 3). Further, the data generated was analysed through the RefFinder algorithm, which also has other gene stability analysis algorithms such as geNorm, Normfinder, BestKeeper, and the comparative ΔCt method implemented in it. The geNorm algorithm calculates the mean stability value of the housekeeping genes by comparing their average pairwise variation with all other housekeeping genes considered. The geNorm algorithm relies on the principle that the expression ratio of two ideal internal control genes remains consistent across all samples, irrespective of experimental condition or cell type [25]. Therefore, the genes with the lowest average pairwise variation are considered to be the most stable. Using the geNorm algorithm, RPS34 and RHA genes were identified to be the most stable, with a mean stability value of 1.238, whereas the UFO gene was identified to be the least stable, with a mean stability value of 3.374. In the geNorm algorithm, any “M” value (geometric mean) below 1.5 is considered indicative of stable gene expression [13]. In our analysis, the top three genes i.e. RPS34, RHA and CyP were having a geometric mean stability value below 1.5. The mean stability value of all 14 genes as estimated through the geNorm algorithm is provided (Figure 4). Similarly, in the comparative ΔCt method, the ΔCt values of genes are compared like geNorm. In the ΔCt comparison method, the arithmetic mean of the standard deviation in the pairwise comparison is considered. The lower the arithmetic mean, the higher will be the stability of the housekeeping gene. Similar to the geNorm method, in the comparative ΔCt method, RPS34 and RHA were found to be the most stable housekeeping genes, with the arithmetic mean and standard deviation of 2.51 and 2.73, respectively (Figure 5). Likewise, the UFO gene was found to be the least stable housekeeping gene with an arithmetic mean of standard deviation of 4.95 (Figure 5). The gene stability ranking of all 14 genes as estimated through the ΔCt method algorithm is provided (Figure 5). In the BestKeeper algorithm, an index value based on the correlation analysis between the Ct values of the genes is calculated. Further, standard deviation, percent covariance explained, and power are calculated. Thereafter, the genes are ranked based on their stability metrics, with the most stable gene having the lowest standard deviation and highest correlation with others. In our analysis, the BestKeeper algorithm identified the USB1 gene to be the most stable, while EfTu was identified to be the least stable (Figure 6). The gene stability ranking of all 14 genes as estimated through the BestKeeper method algorithm is provided (Figure 6). The NormFinder algorithm estimates not only the overall variation of the candidate reference genes but also the variation between subgroups of samples. Unlike the geNorm and ΔCt comparison methods, NormFinder assesses the expression stability of each candidate independently. In our analysis, the RPS34 and ACT2 were found to be most stable using the NormFinder algorithm (Figure 7). The gene stability ranking of all 14 genes as estimated through the ΔCt method algorithm is provided (Figure 7). The RefFinder algorithm compares the weights estimated for each reference gene through the aforementioned four algorithms to calculate the geometric mean of the weights and make a comprehensive ranking. Through the RefFinder, we identified RPS34 and RHA as the best and most stable housekeeping genes, while EfTu and RBCS3B were the least suitable (Figure 8). Among the four algorithms implemented in RefFinder, RPS34 was found to be the best housekeeping gene in the geNorm, ΔCt comparison method, as well as the NormFinder algorithm, except for the BestKeeper algorithm, where USB1 was found to be the best as depicted in the Venn diagram when the top three most stable genes were compared (Figure 9). Similarly, RHA was among the top three in geNorm, the ΔCt comparison method, and BestKeeper when the top three most stable genes from each of these algorithms are compared (Figure 9). Overall, the RPS34 and RHA genes were also found to be most stable and suitable to be used as control genes for normalization study through the RefFinder algorithm. This was also supported by the whisker plot of the Ct values (Figure 3). The gene stability ranking (Better--Good--Average) of the 14 genes in different developmental stages is provided (Table 2). Gene expression stability analysis in different abiotic stress The expression pattern of all of the 14 genes was further explored under four abiotic stresses viz. salt, drought, aluminium, and cold stress in the germination stage tissues (plumule, radicle and cotyledon). In the four abiotic stresses, RPS34 was found to be the most stable, followed by USB1, while RBCS3B was least stable using the geNorm algorithm (Figure S2). Likewise, using the ΔCt method, ACT2 and RPS34 genes were identified to be the most stable HKGs, while RBCS3B was found to be the least stable (Figure S3). In the four abiotic stresses, USB1 was found to be the most stable followed by DNAJ, while EfTu was least stable using the BestKeeper algorithm (Figure S4). Likewise, using ACT2 and DNAJ genes were identified to be the most stable HKGs, while RBCS3B was found to be the least stable using the NormFinder algorithm (Figure S5). Comprehensively, RefFinder identifies ACT2 and RPS34 as the most stable HKGs (Figure S6). The gene stability ranking (Better--Good--Average) of the 14 genes in different abiotic stresses in seedling stages is provided (Table S3). Validation of the candidate HKGs identified in the current study For validation of the best HKGs identified in the current study, we used the Ct value of RPS34, RHA and ACT2 HKGs to normalize the expression of 5 previously reported genes i.e., CAM, PR1, HSP70, NAC, and DEF genes, along with an external control gene, i.e., Actin (Table S2). For the experiment, the tissues and primers were generated and used as previously reported [23]. The gene expression pattern of all 5 test genes was in concurrence with the previous reports [23]. The larger deviation in the relative fold change value of the test gene when normalized with ACT as a control was clearly evident (Figure 10a-e). On the other hand, when the test genes were normalized with the HKGs identified in the current study, a substantial variation was observed. The normalized relative gene expression level clearly reflects the stability as well as the suitability of using RPS34, RHA, and ACT2 genes over ACT in the qRT-PCR experiment (Figure 10a-e). Discussion The gene expression is a highly spatiotemporal event influenced by both external and internal factors. The profiling of gene expression patterns is among the most sought-after preliminary experiments in both animal and plant sciences [ 26 ]. Majorly done by qRT-PCR, there are chances of error due to pipetting and RNA quality affecting the imprecise estimation of gene expression level [ 27 ]. These probable errors are taken care of by using a good reference gene for the normalization of the data generated in the qRT-PCR experiments [ 13 ]. Specifically, the expression of the reference gene needs to be consistent across the spatiotemporal conditions and the different development stages. Therefore, while doing a gene expression analysis, selecting the most trustworthy reference genes has become a crucial research endeavor [ 28 ]. Limited efforts have been made to identify the suitable genes in legume crops, including Lupin angustifolius , Vigna angularis , Lens culinaris , Vicia faba , and Cicer arietinum [ 29 – 33 ]. So far only a handful of studies involving expression profiling have been conducted in V. mungo for different biotic and abiotic stresses [ 23 , 24 , 34 ]. The current study holds the potential to facilitate future research aiming to accurately explore the expression profiling of any selected gene across the developmental stages. To ensure the preciseness of the study, due care has been taken right from RNA isolation to qRT-PCR. Further, the majority of the candidate genes selected for the present study were novel to V. mungo and have not been explored to the best of our knowledge. A clear consensus on which HKGs to be used for data normalization in V. mungo is lacking. Therefore, comparing the expression Ct values obtained from profiling multiples samples using different statistical algorithms will pave the way forward. Hence, the current study was designed to profile a set of 14 HKGs in 17 different tissues representing different stages right from seed germination to seed maturity. The genes used in the current study represent diverse metabolic and developmental pathways, including floral development, cellular cytoskeleton, RNA and protein processing, and protein trafficking. Diversity in the metabolic pathways and developmental stages used in the current study makes the study more reliable and of better utility for the V. mungo researcher community. Further, five different algorithms, including geNorm [ 13 ], NormFinder [ 12 ], BestKeeper [ 14 ], the comparative ΔCt method [ 15 ], and RefFinder were utilized in the current study [ 16 ]. Interestingly, two genes, RPS34 and RHA, were found to be the two most stable HKGs using the different statistical algorithms (Table 2 ), with the exception of the BestKeeper algorithm. Likewise, RBSC3, EfTu and UFO were found to be least stable across the different developmental stages in the current study (Table 2 ). The ribosomal proteins have been reported to be major HKGs in all domains of life [ 35 ]. Since RPS34 is an integral part of protein biosynthesis machinery, its stable expression could be explained. Surprisingly, any of the ribosomal proteins, although expressing stably, have not been explored much as a suitable HKGs for qRT-PCR normalization. The current study advocates to exploring the ribosomal proteins for gene normalization studies in other species. Likewise, RNA helicases are enzymes that are crucial for RNA metabolism. They unwind RNA secondary structures and/or displace RNA-binding proteins and are essential for various cellular processes involving RNA, including pre-mRNA splicing, stability and decay [ 36 ]. Since RNA helicase is ubiquitous and performs several functions vital for cell growth and development, their listing in the most stable HKGs is not surprising. Despite modest changes in the ranking of other genes, likely due to variation in the algorithms used, the findings clearly indicated that RPS34 and RHA consistently were most stable. Although ACT2 outperforms the RPS34 gene under all four abiotic stress conditions, RPS34 was still among the top stable HKGs. Kundu and his coworkers put an effort to identify the best suitable HKGs in V. mungo under biotic (YMV) and abiotic stress (drought) but with only 7 HKGs in two stress conditions [ 34 ]. Comparably, our study not only uses more HKGs but also represents an array of developmental stages (17 tissues) along with the 4 stress conditions (salt, drought, aluminium and heat). The candidate HKGs identified in the current study were further demonstrated to be stably expressed and normalize the Ct value of the 5 test genes under independent experiment conditions. The HKGs identified in the current study i.e. RPS34, RHA, and ACT2 genes hold a greater potential to be widely used as HKGs in any gene expression experiment. Since this is the most comprehensive report on the identification of HKGs in Vigna species, the HKGs identified in the current study might be useful to the other Vigna species. Conclusion The current study explores the stability pattern of 14 HKGs representing diverse metabolic and development pathways in 17 different developmental stages from germination to seed maturity in Vigna mungo . The Ct values obtained were implemented to get a gene-wise ranking based on geNorm, NormFinder, BestKeeper, comparative ΔCt method and the RefFinder algorithm. Based on the comprehensive ranking, RPS34 and RHA were found to be most suitable for qRT-PCR data normalization in different developmental stages, while ACT2 and RPS34 were suitable under abiotic stress conditions. Based on the results obtained, the current study advocates using a combination of HKGs (RPS34 and RHA for developmental tissues and ACT2 and RPS34 for stress-treated tissues) for normalization of qRT-PCR data. The reference genes suggested here will ensure reliability and uniformity in qRT-PCR findings. Further, the expression of any gene is spatiotemporal and is influenced by both internal and external factors. And the genotypic background also plays an important role in it. The key findings of the current study need to be reconfirmed in independent background, including germplasm lines, to get more comprehensive insight. Parallelly, these genes further may also be tested and used in other related Vigna species such as V. radiata , V. unguiculata and V. angularis , which share closer ancestry to V. mungo . Declarations Author contributions KK has conceptualized the whole study. PS, PK, SG, AAC, VK, KT performed the experiment and collected the plant tissues. PS, and AD performed the qRT-PCR. MR and KK performed statistical analysis. SKJ, KK, AD, MR, GPD wrote the draft manuscript. MR, GPD and KK edited the manuscript to prepare final draft. MR and GPD supervised and administered the whole study. Acknowledgements We acknowledge the support provided by the Director, ICAR-IIPR, Kanpur. The research work carried out in this manuscript is a part of Institute funded project. Conflict of Interest There is no conflict of interest, either financial or personal, with other people or organizations that could inappropriately influence or perceived to influence this work. Data availability statement The data that support the findings of this study are available in the supporting information of this article References Singh, N. K. et al. The first draft of the pigeonpea genome sequence. Journal of plant biochemistry and biotechnology. 21 :98-112 (2012). https://doi.org/10.1007/s13562-011-0088-8 Varshney, R. K. et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nature biotechnology. 31( 3 ):240-246 (2013). https://doi.org/10.1038/nbt.2491 Junaid, A. et al. Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo. Life Science Alliance. 7( 1 ), (2024). https://doi.org/10.26508/lsa.202302074 Gaikwad, K. et al. The chromosome-scale genome assembly of cluster bean provides molecular insight into edible gum (galactomannan) biosynthesis family genes. Scientific Reports. 13( 1 ):9941 (2023). https://doi.org/10.1038/s41598-023-33762-3 Ameline-Torregrosa, C. et al. Transcriptomic approaches to unravel plant–pathogen interactions in legumes. Euphytica. 147 :25-36 (2006). https://doi.org/10.1007/s10681-006-6767-1 Garg, R & Jain, M. Transcriptome analyses in legumes: A resource for functional genomics. The Plant Genome. 6( 3 ) (2013). https://doi.org/10.3835/plantgenome2013.04.0011 Abdelrahman, M. et al. Legume genetic resources and transcriptome dynamics under abiotic stress conditions. Plant, cell & environment. 41( 9 ):1972-1983 (2018). https://doi.org/10.1111/pce.13123 Nicot, N. et al. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. Journal of Experimental Botany. 56( 421 ):2907-2914 (2005). https://doi.org/10.1093/jxb/eri285 Udvardi, M. K. Eleven golden rules of quantitative RT-PCR. The Plant Cell. 20( 7 ):1736-1737 (2008). https://doi.org/10.1105/tpc.108.061143 Bustin, S. A. et al. MIQE precis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Molecular Biology. 11 :1-5 (2010). https://doi.org/10.1186/1471-2199-11-74 Manu, B. et al. Genomic designing for abiotic stress tolerance in mungbean and urdbean. In Genomic Designing for Abiotic Stress Resistant Pulse Crops (pp. 271-343). Cham: Springer International Publishing (2022). https://doi.org/10.1007/978-3-030-91039-6_8 Andersen, C. L. et al. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer research. 64( 15 ):5245-5250 (2004). https://doi.org/10.1158/0008-5472.CAN-04-0496 Vandesompele, J. et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3 :1-12 (2002). https://doi.org/10.1186/gb-2002-3-7-research0034 Pfaffl, M. W. et al. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnology letters. 26 :509-515 (2004). https://doi.org/10.1023/b:bile.0000019559.84305.47 Silver, N. et al. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Molecular Biology. 7 :1-9 (2006). https://doi.org/10.1186/1471-2199-7-33 Xie, F. et al. RefFinder: a web-based tool for comprehensively analyzing and identifying reference genes. Functional & Integrative Genomics. 23( 2 ):125 (2023). https://doi.org/10.1007/s10142-023-01055-7 Dantas, B. F. et al. Germination, initial growth and cotyledon protein content of bean cultivars under salinity stress. Revista Brasileira de Sementes. 29 :106-110 (2007). https://doi.org/10.1590/S0101-31222007000200014 Roso, R. et al. Effects of Aluminum and pH on Germination of Echium plantagineum L. Journal of Agricultural Studies . 8( 1 ):313-327 (2020). https://doi.org/10.5296/jas.v8i1.16179 Islam, M. R. Evaluation of Polyethylene Glycol (PEG-6000) Induced Drought Stress Tolerant Mungbean Genotypes by using Correlation, Principal Component, Hierarchical Clustering and Multi-Trait Genotype-Ideotypes Distance Index Analysis. Current Applied Science And Technology. e0260871-e0260871 (2023). https://doi.org/10.55003/cast.2024.260871 Manasa, L. S. et al. Cold tolerance mechanisms in mungbean (Vigna radiata L.) genotypes during germination. Agriculture. 13( 2 ):315 (2023). https://doi.org/10.3390/agriculture13020315 Kumar, K. et al. Identification and characterization of NHX gene family for their role under salt stress in Vigna mungo. Physiologia Plantarum. 176( 5 ):e14563 (2024). https://doi.org/10.1111/ppl.14563 Marchler-Bauer, A. et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Research. 45( D1 ), D200-D203 (2017). https://doi.org/10.1093/nar/gkw1129 Kumar, S. et al. Gene expression and biochemical profiling of contrasting Vigna mungo genotypes against Mungbean Yellow Mosaic India Virus (MYMIV). Journal of Food Legumes. 35( 2 ):107-116 (2022). Kumar, K. et al. Single trait versus principal component based association analysis for flowering related traits in pigeonpea. Scientific Reports. 12( 1 ):10453 (2022). https://doi.org/10.1038/s41598-022-14568-1 Vijayakumar S & Sakuntala M. Validation of reference gene stability for normalization of RT-qPCR in Phytophthora capsici Leonian during its interaction with Piper nigrum L. Scientific Reports. 14(1):7331 (2024). https://doi.org/10.1038/s41598-024-58139-y Nolan, T. Hands, R. E. & Bustin, S. A. Quantification of mRNA using real-time RT-PCR. Nature Protocols. 1( 3 ):1559-1582 (2006). https://doi.org/10.1038/nprot.2006.236 Hu, R. et al. Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR. BMC Molecular Biology. 10:1-12 (2009). https://doi.org/10.1186/1471-2199-10-93 de Vega-Bartol, J. J. et al. Normalizing gene expression by quantitative PCR during somatic embryogenesis in two representative conifer species: Pinus pinaster and Picea abies. Plant Cell Reports. 32:715-729 (2013). https://doi.org/10.1007/s00299-013-1407-4 Taylor, C. M. et al. Identifying stable reference genes for qRT-PCR normalisation in gene expression studies of narrow-leafed lupin ( Lupinus angustifolius L.). PLoS One. 11(2):e0148300 (2016). https://doi.org/10.1371/journal.pone.0148300 Chi, C. et al. Selection and validation of reference genes for gene expression analysis in Vigna angularis using quantitative real-time RT-PCR. PloS one. 11(12):e0168479 (2016). https://doi.org/10.1371/journal.pone.0168479 Sinha, R. Sharma, T. R., Singh, A. K. Validation of reference genes for qRT-PCR data normalisation in lentil (Lens culinaris) under leaf developmental stages and abiotic stresses. Physiology and Molecular Biology of Plants. 25( 1 ):123-134 (2019). https://doi.org/10.1007/s12298-018-0609-1 Gutierrez, N. et al. Assessment of candidate reference genes for expression studies in Vicia faba L. by real-time quantitative PCR. Molecular Breeding. 28 :13-24 (2011). https://doi.org/10.1007/s11032-010-9456-7 Reddy, D. S. et al. Identification and validation of reference genes and their impact on normalized gene expression studies across cultivated and wild Cicer species. PloS one. 11(2):e0148451 (2016). https://doi.org/10.1371/journal.pone.0148451 Kundu, A. et al. Complex molecular mechanisms underlying MYMIV-resistance in Vigna mungo revealed by comparative transcriptome profiling. Scientific Reports. 9( 1 ):8858 (2019). https://doi.org/10.1038/s41598-019-45383-w Kumar J & Kumar V. Ribosome proteins—Their balanced production. In: Emerging Concepts in Ribosome Structure, Biogenesis, and Function (pp. 47-87). Academic Press (2021). https://doi.org/10.1016/B978-0-12-816364-1.00003-2 Bourgeois, C. F., Mortreux, F. & Auboeuf, D. The multiple functions of RNA helicases as drivers and regulators of gene expression. Nature Reviews Molecular Cell Biology . 17( 7 ):426-438 (2016). https://doi.org/10.1038/nrm.2016.50 Tables Table 1: List of genes used in the current study and their function S.no. Arabidopsis gene Gene name Vigna mungo ortholog used in the current study Percent homology with Arabidopsis homolog Function AT1G30950.1 Unusual Floral Organ (UFO) WVZ23601.1 66.743 Encodes Unusual floral organ genes which is required for the proper identity of the floral meristem and is constitutively expressed in the floral tissues AT5G62690.1 TUBULIN BETA CHAIN 2 (TUB2) WVY95928.1 95.778 Encodes tubulin beta-2/beta-3 chain, which are components of microtubules in the cytoskeleton and play essential roles in various cellular processes AT3G18780.2 ACTIN 2 (ACT2) WVY99582.1 92.042 Encodes an actin that is constitutively expressed in vegetative structures AT1G58050.1 RNA helicase (RHA) WVZ12136.1 47.337 Encodes RNA helicase that unwind RNA secondary structures and/or displace RNA-binding proteins in an ATP-dependent manner, that is essential for various cellular processes involving RNA AT1G13340.1 ISTL6 WVZ09401.1 51.507 ISTL6 gene is crucial for the proper trafficking and degradation of membrane proteins within the cell. Proper functioning of this pathway is vital for maintaining cellular homeostasis and responding to environmental stressors 6. AT2G36130.1 Cyclophilin (CyP) WVZ05807.1 87.805 Cyclophilins are a family of highly conserved proteins with peptidyl-prolyl cis-trans isomerase (PPIase) activity, which catalyzes the isomerization of peptide bonds at proline residues that is crucial for protein folding, and function 7. AT2G20560.1 DNAJ WVZ01865.1 79.769 DNAJ heat shock family protein are molecular chaperones involved in various cellular processes and play essential roles in protein folding, stabilization, degradation, and trafficking 8. AT5G38410.3 Rubisco small subunit 3b (RBCS3B) WVZ13955.1 67.935 Encodes a member of the Rubisco small subunit (RBCS) multi-gene family 9. AT5G58990.1 28S ribosomal S34 protein (RPS34) WVZ02079.1 84.286 (RPS34) is a component of the small subunit of the ribosome in eukaryotes that are responsible for protein synthesis, and their subunits (small and large) work together to translate mRNA into polypeptides 10. AT4G27320.1 Universal Stress Protein 21 (USP21) WVZ04305.1 70.892 Protein that contains a universal stress protein domain. They are essential for plant stress tolerance as well as cellular homeostasis 11. AT5G51170.1 U6 snRNA phosphodiesterase like protein (USB1) WVZ06716.1 65.351 Encodes U6 snRNA phosphodiesterase-like protein, that are involved in RNA processing, specifically within the spliceosome machinery, where they contribute to pre-mRNA splicing 12. AT1G07920.1 Elongation Factor-TU (EfTu) WVY95914.1 96.197 GTP binding Elongation factor Tu family protein 13. AT4G05320.2 Ubiquitin 10 (UBQ10) WVZ24872.1 99.738 Encodes Ubiquitin 10 that is a member of the ubiquitin protein family, which plays a central role in post-translational modification of proteins through ubiquitination 14. AT4G27960.2 Ubiquitin Conjugating Enzyme 9 (UBQE9) WVY93904.1 94.667 Encodes UBIQUITIN CONJUGATING ENZYME 9 is a key enzyme in the ubiquitination pathway, primarily involved in SUMOylation Table 2: Gene stability ranking of the 14 HKGs used in the current study through different algorithms using 17 developmental tissues Ranking Order (Better--Good--Average) Method 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Delta CT RPS34 RHA CyP ACT2 DNAJ UBQC9 ISTL6 UBQ10 USB1 TUB2 USP21 EfTu RBCS3B UFO BestKeeper USB1 DNAJ RHA ACT2 UFO CyP RPS34 UBQC9 UBQ10 ISTL6 TUB2 USP21 RBCS3B EfTu NormFinder RPS34 ACT2 DNAJ CyP RHA UBQC9 ISTL6 UBQ10 USB1 TUB2 USP21 EfTu RBCS3B UFO geNorm RPS34 RHA CyP ACT2 DNAJ ISTL6 UBQC9 UBQ10 USB1 TUB2 USP21 EfTu RBCS3B UFO Comprehensive ranking RPS34 RHA ACT2 DNAJ CyP USB1 UBQC9 ISTL6 UBQ10 TUB2 UFO USP21 EfTu RBCS3B Additional Declarations No competing interests reported. Supplementary Files 6Supplementarytables.docx 5Supplementaryfigures.docx Cite Share Download PDF Status: Published Journal Publication published 26 Sep, 2025 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 19 May, 2025 Reviews received at journal 09 May, 2025 Reviews received at journal 01 May, 2025 Reviewers agreed at journal 18 Apr, 2025 Reviewers agreed at journal 16 Apr, 2025 Reviewers invited by journal 16 Apr, 2025 Submission checks completed at journal 10 Apr, 2025 First submitted to journal 29 Mar, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5819031","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":435671789,"identity":"0dc95477-2dcd-4f89-b2ee-be175066b407","order_by":0,"name":"Sudhir Kumar Jha","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Sudhir","middleName":"Kumar","lastName":"Jha","suffix":""},{"id":435671791,"identity":"c11ec863-ddde-4938-bece-2c3553ce5287","order_by":1,"name":"Pritee Sagar","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Pritee","middleName":"","lastName":"Sagar","suffix":""},{"id":435671792,"identity":"a1b039a5-b83c-436b-a02d-621f0a790414","order_by":2,"name":"Prateek Kumar","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Prateek","middleName":"","lastName":"Kumar","suffix":""},{"id":435671793,"identity":"408211e7-a727-432a-b2c5-0ef94d3030bf","order_by":3,"name":"Shivangi Gupta","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Shivangi","middleName":"","lastName":"Gupta","suffix":""},{"id":435671795,"identity":"bd857235-ab94-4ed2-91d5-ee9aa5e918bf","order_by":4,"name":"Ankita Anil Chakote","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Ankita","middleName":"Anil","lastName":"Chakote","suffix":""},{"id":435671796,"identity":"c12107ae-8818-44d5-8f96-e590666ffb20","order_by":5,"name":"Vaibhav Kumar","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Vaibhav","middleName":"","lastName":"Kumar","suffix":""},{"id":435671797,"identity":"01bff4ba-5b0d-449c-9e2d-00cfa85898e0","order_by":6,"name":"Kalpana Tewari","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Kalpana","middleName":"","lastName":"Tewari","suffix":""},{"id":435671800,"identity":"63625dda-6773-48a8-bdfe-42b9c42dee99","order_by":7,"name":"Antara Das","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Antara","middleName":"","lastName":"Das","suffix":""},{"id":435671801,"identity":"44558484-d80d-4770-8a7e-2d4bb30a4261","order_by":8,"name":"Meenal Rathore","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Meenal","middleName":"","lastName":"Rathore","suffix":""},{"id":435671802,"identity":"6a055a7f-1011-4f1f-acfe-a8ed8108872b","order_by":9,"name":"Girish Prasad Dixit","email":"","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":false,"prefix":"","firstName":"Girish","middleName":"Prasad","lastName":"Dixit","suffix":""},{"id":435671803,"identity":"ac99199f-d19f-452c-ac7a-5a173f23d1d1","order_by":10,"name":"Kuldeep Kumar","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA20lEQVRIiWNgGAWjYBACCQbmBgkgzbiBgfkAiC9DhBZGmBa2BBCfhxQtPAYgAcJaJGckNt74UHFHdjv7mc+vbtRY8DCwHz66AZ8WaYnEZssZZ54Z7+zJ3WadcwzoMJ60tBv4tMhJJLZJ87YdTtxwg3ebcQ4bUIsEjxlhLX//gbTwPDPO+UeEFmmQFsYGsBbmx7ltRGiR7HnYbNlzDOSXNDPm3D4JHjZCfpE4nnzwxo8aUIgdfvw551udHD/74WN4tUDBARDBJgEmiVAO18L8gUjVo2AUjIJRMMIAAIcsTRl7Y8CPAAAAAElFTkSuQmCC","orcid":"","institution":"ICAR-Indian Institute of Pulses Research","correspondingAuthor":true,"prefix":"","firstName":"Kuldeep","middleName":"","lastName":"Kumar","suffix":""}],"badges":[],"createdAt":"2025-01-13 10:23:23","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5819031/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5819031/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-025-10295-5","type":"published","date":"2025-09-26T15:56:52+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":80028851,"identity":"ab89616d-757c-4ce5-9925-b0abbc03941e","added_by":"auto","created_at":"2025-04-07 07:07:49","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":288841,"visible":true,"origin":"","legend":"\u003cp\u003eTissues used in the identification of best HKGs for normalization of qRT-PCR data in blackgram (a) germination stage from which plumule, radicle and cotyledon tissues were collected at 5 days after sowing; (b) Vegetative stage from which vegetative leaves, meristem, stem, and roots tissues were collected at 30 days after sowing; (c) Reproductive stage from which reproductive buds, reproductive leaves, roots, and stem tissues were collected at 45 days after sowing; (d) Podding stage from which leaves, whole pod, stem and root tissues were collected at 60 days after sowing and (e) pods were split into pod wall and seed.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/a4110935f6886c38274b32b8.png"},{"id":80028679,"identity":"3748eabf-83bb-4df9-9b53-c5ff6536deef","added_by":"auto","created_at":"2025-04-07 06:59:49","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":84139,"visible":true,"origin":"","legend":"\u003cp\u003eThe agarose gel displaying the quality of the RNA isolated from the selected tissues (a), and the specific amplicon of the 14 housekeeping genes confirmed by their respective band sizes. Where P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11, P12, P13, P14, and L stands for UBQE9, TUB2, USP21, RPS34, USB1, RHA, ISTL6, EfTu, ACT2, DNAJ, CyP, RBCS3B, UFO, UBQ10 and 100 base pair Ladder, respectively\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/50b9872fc9c2eec0dd876bf1.png"},{"id":80028682,"identity":"c3ff8c17-f54d-4fe9-bae1-f645cb09d9c2","added_by":"auto","created_at":"2025-04-07 06:59:49","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":119583,"visible":true,"origin":"","legend":"\u003cp\u003eVariation in Ct values of all 14 HKGs considered for the current study across all the tissues. The central line within the box represents the mean value. The gray and yellow boxes represents the quartile where the Ct values for most of the HKGs falls. The middle line represents the middle (median) Ct value in the dataset. The larger vertical lines above the yellow and below gray whiskers represent the presence of outliers in the Ct values in some tissues.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/af9712395a5b6ded5ba22f32.png"},{"id":80028707,"identity":"f14b11c4-33d7-4f8f-bee8-f7f8f28414a5","added_by":"auto","created_at":"2025-04-07 07:00:02","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":166025,"visible":true,"origin":"","legend":"\u003cp\u003eGene stability analysis using the geNorm algorithm identifies RPS34 and RHA as the most suitable while RBCS3B and UFO as the least suitable for normalization of qRT-PCR data across different tissues.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/5246f7370a2a0d57650fe415.png"},{"id":80029773,"identity":"9448d4c4-3120-4641-95c4-b9fb42c75a79","added_by":"auto","created_at":"2025-04-07 07:15:49","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":163927,"visible":true,"origin":"","legend":"\u003cp\u003eThe ΔCt comparison method for gene stability analysis suggests RPS34 and RHA as the most suitable while RBCS3B and UFO as the least suitable for normalization of qRT-PCR data across different tissues.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/63b8b5935c90f4d13f541e7e.png"},{"id":80028684,"identity":"79ebacff-0086-413b-939b-8e17f83478e8","added_by":"auto","created_at":"2025-04-07 06:59:49","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":155124,"visible":true,"origin":"","legend":"\u003cp\u003eThe BestKeeper algorithm for gene stability analysis identifies USB1 and DNAJ as the most suitable while RBCS3B and EfTu as least suitable for normalization of qRT-PCR data across different tissues.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/638018b7c98c452ea580fb85.png"},{"id":80028680,"identity":"4be88abd-5c84-40c9-8390-ffab6e4d626d","added_by":"auto","created_at":"2025-04-07 06:59:49","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":167185,"visible":true,"origin":"","legend":"\u003cp\u003eGene stability analysis using the Normfinder algorithm infers RPS34 and ACT2 to be most suitable while RBCS3B and UFO as least suitable for normalization of qRT-PCR data across different tissues\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/76a6c4f074da1a7dbd89f35b.png"},{"id":80028688,"identity":"8193b28e-ce90-4ef6-882d-7bfea13b8746","added_by":"auto","created_at":"2025-04-07 06:59:49","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":167210,"visible":true,"origin":"","legend":"\u003cp\u003eGene stability analysis using the RefFinder algorithm infers RPS34 and RHA to be most suitable while EfTu and RBCS3B as least suitable for normalization of qRT-PCR data across different tissues\u003c/p\u003e","description":"","filename":"8.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/a7667db6cfc7388ed306e109.png"},{"id":80028862,"identity":"de46ff71-7c45-46ff-a5c2-b920d36764f3","added_by":"auto","created_at":"2025-04-07 07:07:49","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":206747,"visible":true,"origin":"","legend":"\u003cp\u003eConsistency in the stability of the HKGs using multiple algorithms explained via Venn diagram. The RHA (ΔCt comparison, geNorm and BestKeeper) and RPS34 (BestKeeper, NormFinder and geNorm) genes were found to be most suitable via three algorithms whereas, CyP (ΔCt comparison and geNorm) and DNAJ (NormFinder and BestKeeper) were identified consistently via two algorithms.\u003c/p\u003e","description":"","filename":"9.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/250269cd844f5a757d389456.png"},{"id":80028704,"identity":"5aaa907a-09b0-4e20-a818-ce5272402507","added_by":"auto","created_at":"2025-04-07 06:59:50","extension":"png","order_by":10,"title":"Figure 10","display":"","copyAsset":false,"role":"figure","size":138225,"visible":true,"origin":"","legend":"\u003cp\u003eRelative expression level of 5 test gene i.e. CAM (a), PR1 (b), HSP70 (c), NAC (d) and DEF (e) when normalized with the RPS34, PHS and ACT2 along with the external control ACT. The relative expression level clearly reflect the consistency in the expression pattern of the test genes when normalized with the HKGS identified in the current study. A larger deviation can be seen in the expression pattern of the test genes when normalized with external control ACT gene.\u003c/p\u003e","description":"","filename":"10.png","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/e2dbe4ff5247010f6c42a75f.png"},{"id":92430414,"identity":"78344543-910c-466b-9786-2e1fb4833ac8","added_by":"auto","created_at":"2025-09-29 16:01:49","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2471579,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/105bc1f8-7304-494b-a609-b2b00f575cda.pdf"},{"id":80028706,"identity":"d71dab1f-1bb6-4e2c-842d-32a603152541","added_by":"auto","created_at":"2025-04-07 06:59:57","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":24285,"visible":true,"origin":"","legend":"","description":"","filename":"6Supplementarytables.docx","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/0a6175b0389dfcf4c0634d42.docx"},{"id":80028853,"identity":"21060509-c10c-4965-afbe-74531f482a51","added_by":"auto","created_at":"2025-04-07 07:07:49","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":291921,"visible":true,"origin":"","legend":"","description":"","filename":"5Supplementaryfigures.docx","url":"https://assets-eu.researchsquare.com/files/rs-5819031/v1/d0a60c3f15e1e5d2f090540a.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Identification of best housekeeping gene for normalization in qRT-PCR data under different development stages and stress conditions in Vigna mungo","fulltext":[{"header":"Key Message","content":"\u003cp\u003eThe gene stability ranking indicated that RPS34 and RHA were best for qRT-PCR data normalization in various developmental stages, whereas ACT2 and RPS34 under abiotic stress.\u003c/p\u003e"},{"header":"Introduction","content":"\u003cp\u003eWith the advent of sequencing technologies, abundant information and genomics resources are being developed on a daily basis in several legumes, which were otherwise considered as \u0026ldquo;orphan crops.\u0026rdquo; To date, genome sequences for several legumes, including \u003cem\u003eCajanus cajan\u003c/em\u003e, \u003cem\u003eCicer arietinum\u003c/em\u003e, \u003cem\u003eVigna mungo\u003c/em\u003e, and \u003cem\u003eCyamopsis tetragonoloba\u003c/em\u003e have been published [\u003cspan additionalcitationids=\"CR2 CR3\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. Parallelly, several of these legume crops have been supplemented with multiple transcriptomics experiments to mine genes related to any development stage and their biotic-abiotic response and have been reviewed [\u003cspan additionalcitationids=\"CR6\" citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. Gene expression profiling, either \u0026lsquo;global\u0026rsquo; or 'local', is becoming vital in any advanced molecular biology and molecular breeding experiments. Even the sole transcriptomics experiments require validation of key genes to make the results more reliable and statistically more robust.\u003c/p\u003e \u003cp\u003eThe expression of any gene is spatiotemporal and dynamic, as it is influenced by several internal as well as external factors. Quantification of mRNA is essential for understanding the gene expression levels in the biological samples. Various techniques, including quantitative reverse transcription polymerase chain reaction (qRT-PCR), northern blotting, RNA-seq, micro-array, \u003cem\u003ein-situ\u003c/em\u003e hybridization, the Nanostring nCounter system, and the Serial Analysis of Gene Expression (SAGE), have been developed for this purpose, each with its own advantages and limitations. Among these techniques, qRT-PCR is the most widely used technique due to its robustness and accuracy. For accurate mRNA quantification, the qRT-PCR results are normalized using a reference gene/housekeeping gene (HKG) whose expression is stable across different experimental conditions, allowing for the quantification of changes that are due to biological differences rather than technical artifacts [\u003cspan additionalcitationids=\"CR9\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eUrdbean, also known as blackgram (\u003cem\u003eV. mungo\u003c/em\u003e), is the third most important grain legume after chickpea and pigeonpea in the Indian subcontinent. It is a nutrient-rich pulse possessing a rich blend of protein, fiber, and essential micronutrients like iron, zinc and potassium. Its adaptability to diverse agro-ecologies, shorter cropping duration, and ability to fix atmospheric nitrogen attract the small and marginal farmers [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. As the world grapples with the burning issues like climate change, ever increasing population and nutritional security, the importance of the urdbean as a nutrient-rich and agriculturally important crop is being recognized recently. Despite this significance, it has long been overlooked by mainstream research and is treated as an \u0026ldquo;orphan crop\u0026rdquo;. Recently, efforts are being made to develop some genetic and genomic resources in urdbean. The availability of draft genome assembly, high density linkage map, along with the development of molecular markers such as SSRs and SNPs has enabled the identification of genes/QTLs linked with different traits. These genomic resources hold great promise for accelerating the urdbean improvement using the advanced genomics strategies [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Once the putative candidate genes are known their expression analysis needs to be performed, which requires a suitable HKG for the data normalization.\u003c/p\u003e \u003cp\u003eSeveral HKGs have been explored for normalization of the qRT-PCR data. In general, the HKGs used across the studies belong to central processes that are critical for the survival of the cell including, genes essential for the cell morphology, cell division, growth and development. Because these aforementioned genes are tied to fundamental and non-specialized processes, they are less likely to fluctuate in response to experimental variables that might alter the expression of genes involved in specific pathways. However, their stability should be validated for specific experiments, as no single gene is universally optimal for all situations. Therefore, in this work, we evaluated 14 genes of different functional classes (significantly reducing the chance that the genes will be co-regulated) using the qRT-PCR techniques in multiple stages, starting from cotyledon (germination) to pod (maturity) stages. The most suitable HKG for qRT-PCR in \u003cem\u003eV. mungo\u003c/em\u003e was evaluated employing different algorithms developed to assess the stability of a normalizing gene through NormFinder, geNorm, BestKeeper, and the ΔCt method implemented in RefFinder software [\u003cspan additionalcitationids=\"CR13 CR14 CR15\" citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. The most potential HKGs identified in the current study were further validated to demonstrate their suitability in normalizing the qRT-PCR data. The most stable HKGs identified in the current study will be highly useful for the researcher working on the advanced genomics experiments in Vigna species, including \u003cem\u003eV. mungo\u003c/em\u003e.\u003c/p\u003e"},{"header":"Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePlant materials\u003c/h2\u003e \u003cp\u003eIn the current work, an urdbean variety, \u0026ldquo;IPU13-1\u0026rdquo; released from the ICAR-Indian Institute of Pulses Research, Kanpur, was used as plant material. The variety has a potential yield of 10\u0026ndash;12 q/ha with a maturation duration of 70\u0026ndash;75 days. The seeds were sterilized for 30 seconds in 70% ethanol before sowing in clay pots filled with autoclaved soil. Four developmental stages, including germination [5 days after sowing (DAS)], the vegetative stage (30 DAS), the reproductive stage (45 DAS), and podding (60 DAS), were selected for the gene expression analysis. From the germination stage, plumule, radicle, and cotyledon tissues were harvested. While leaves, stems, and root tissue were harvested from all remaining three stages, including the vegetative stage, reproductive stage, and podding stage. Likewise, vegetative meristem, floral bud, and 15-day-old pods were harvested at the vegetative, reproductive, and podding stages, respectively. Further, the pods were also split into pod walls and seeds. All tissues were freshly harvested and snap-frozen into liquid nitrogen before storing them at -80\u0026deg;C until RNA isolation. All the samples were collected in three independent biological replicates.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eTreatments\u003c/h3\u003e\n\u003cp\u003eFour different treatments, including drought, salt, aluminium, and cold stress, were given to the sterilized seeds of IPU13-1 at the germination stage. The experiment for screening in salt, drought, and aluminium stress conditions was done in the customized growth chamber of REMI (CHM 12+) with a Floutronix grow light. The temperature was set at 35\u0026deg;C\u0026thinsp;\u0026plusmn;\u0026thinsp;2\u0026deg;C with relative humidity set to 65% and a light/dark period of 12 hours each. For salinity and aluminium stress, seeds were germinated in germination paper put in a petri-plate, followed by soaking with 120 mM NaCl and a 100 mM solution of Al\u003csub\u003e2\u003c/sub\u003e(SO\u003csub\u003e4\u003c/sub\u003e)\u003csub\u003e3\u003c/sub\u003e, respectively [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Likewise, the sterilized seedlings were germinated in MS media supplemented with 1% w/v of PEG 6000 for imposing drought stress [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. Similarly, to collect tissues under cold stress, seeds were germinated in germination paper with the temperature set to 18\u0026deg;C\u0026thinsp;\u0026plusmn;\u0026thinsp;2\u0026deg;C with relative humidity set to 65% and a light/dark period of 12 hours each [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. Tissues of plumule, radicle, and cotyledon from each of the aforementioned treatments were collected at 15 DAS in three independent biological replicates for the downstream analysis.\u003c/p\u003e\n\u003ch3\u003eTotal RNA extraction and cDNA synthesis\u003c/h3\u003e\n\u003cp\u003eRNA isolation for each sample was done using the RNeasy Plant Mini Kit, Qiagen, Germany (Catalogue no. 74904), following the manufacturer\u0026rsquo;s protocol with slight modifications [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. The DNase treatment was given to the isolated RNAs to remove the genomic DNA contamination. Total RNA integrity was checked in 2% agarose gel electrophoresis and the amount was quantified using a spectrophotometer (NanoDrop 2000c, Thermo-Scientific, USA). From the total quantity of extracted RNA, 1 \u0026micro;g was used for cDNA synthesis using Maxima H Minus Double-Stranded cDNA Synthesis Kit, Thermo-Scientific, USA (Catalogue no. K1612), following the manufacturer\u0026rsquo;s instructions. Finally, the cDNA was diluted with molecular-grade water in a ratio of 1:10 and used for qRT-PCR analysis [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e].\u003c/p\u003e\n\u003ch3\u003eSelection of housekeeping genes and primer synthesis\u003c/h3\u003e\n\u003cp\u003eA total of 14 housekeeping genes were shortlisted based on the literature survey. The majority of these genes have not been tested in urdbean to date. The identifier for the selected genes was used to fetch the Arabidopsis protein sequences for the selected genes from TAIR (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.arabidopsis.org/\u003c/span\u003e\u003cspan address=\"https://www.arabidopsis.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The reference protein sequences for the Urdbean genome were downloaded from NCBI (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/datasets/genome/GCA_036885675.1/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_036885675.1/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Further, the NCBI blastp suit (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins\u003c/span\u003e\u003cspan address=\"https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e)\u003c/span\u003e was used for finding the urdbean homologs of the selected Arabidopsis protein sequences based on the similarity. The Urdbean protein with the maximum similarity for each Arabidopsis protein sequence was shortlisted for the downstream analysis. The protein sequences were further confirmed for functional conservation by analysis of conserved domains in both Arabidopsis and \u003cem\u003eVigna mungo\u003c/em\u003e homologs via NCBI-Batch CDD (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e; [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]). Primers from the selected genes were designed using the IDT's PrimerQuest tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.idtdna.com/PrimerQuest\u003c/span\u003e\u003cspan address=\"https://www.idtdna.com/PrimerQuest\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e).\u003c/span\u003e PrimerQuest was chosen due to its robust primer design algorithms, which incorporate melting temperature, GC content, and secondary structure considerations. Further, the oligonucleotides were synthesized by Europhin Genomics, India. The primer sequences for these genes are provided in the supplementary file (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e).\u003c/p\u003e\n\u003ch3\u003eThe qRT-PCR analysis of the selected HKGs\u003c/h3\u003e\n\u003cp\u003eTemperature gradient PCR was performed to optimize the specific Tm of each primer. Each 20 \u0026micro;L qRT-PCR reaction consisted of 10 \u0026micro;L master mix (Brilliant III Ultra-fast SYBR Green qRT-PCR Master Mix, Agilent Technologies), 0.3 \u0026micro;M forward and reverse primers, 7.4 \u0026micro;L PCR grade water and 50 ng cDNA as the template. The qRT-PCR reactions were carried out in an Agilent AriaMx Real-Time qRT-PCR system. The PCR conditions involved an initial denaturation at 94\u0026deg;C for 3 min, followed by 40 cycles of 95\u0026deg;C for 30 s, 60\u0026deg;C for 15 s, and 72\u0026deg;C for 30 s [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. The amplicons were separated in 1.2% agarose gel using TBE buffer.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eExpression stability analysis\u003c/h2\u003e \u003cp\u003eFor qRT-PCR data, raw Ct values (quantification cycle) were used for the analysis, and the mean Ct value of each biological replicate was considered to ensure the accuracy. For the expression stability analysis, the RefFinder web server (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ciidirsinaloa.com.mx/RefFinder-master/\u003c/span\u003e\u003cspan address=\"https://www.ciidirsinaloa.com.mx/RefFinder-master/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e; [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]) was employed. The RefFinder tool comprehensively utilizes several other gene stability analysis algorithms such as geNorm [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e], Normfinder [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e], BestKeeper [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e], and the comparative ΔCt method [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. The RefFinder tool utilizes only the original Ct value from quantitative real-time PCR (qRT-PCR) for the downstream analysis [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. Accordingly, the mean values of each replicate were used for the analysis using the RefFinder tool.\u003c/p\u003e \u003cp\u003eThe geNorm algorithms calculate the gene expression stability measure (M-value) based on pairwise variation analysis and also determine the optimal number of reference genes required for the normalization. Whereas, the NormFinder algorithm employs an ANOVA-based model to estimate intra- and inter-group variation, providing a direct stability value for the HKGs tested. Likewise, the BestKeeper algorithm assesses the stability based on the standard deviation (SD) and coefficient of variation (CV) of the Ct values. The HKGs with SD\u0026thinsp;\u0026gt;\u0026thinsp;1 are considered unreliable for normalization. Similarly, the ΔCt method compares the relative expression stability of the HKGs by analyzing the variation in their Ct values. The RefFinder algorithm comprehends an overall ranking based on the geometric mean of rankings from these four methods.\u003c/p\u003e \u003cp\u003eThe statistical significance of candidate reference gene stability was determined based on rankings derived from the combined analysis. Further, to plot the mean stability value (geNorm), stability value (NormFinder), standard deviations (BestKeeper), arithmetic mean of the standard deviation (ΔCt method), and the geometric mean of the stability ranking (RefFinder) WPS office was used. The box-plot representing the distribution of the Ct values was also plotted using the GraphPad v10.4.1.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eGene normalizing efficiency of the candidate HKGs\u003c/h3\u003e\n\u003cp\u003eFurther, to validate the reliability of the most suitable HKGs (RPS34, RHA, and ACT2) identified in the current study, they were used as an internal control in the expression analysis of the 5 test genes previously reported by our group (Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e) [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. Parallelly, an external Actin (ACT) gene that has been previously reported and used in urdbean qRT-PCR was also used as an internal control. The control and infected leaf samples from Yellow Vein Mosaic (YMV) tolerant (DPU88-31) and susceptible (LBG17) genotypes were collected as previously reported by our group [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. Further, the RNA isolation, cDNA library preparation, and qRT-PCR were performed using the aforementioned methodologies. The mean Ct values were calculated in Microsoft Excel, 2016 using all 4 candidate HKGs, and a one-way ANOVA followed by a Tukey's multiple comparisons test was performed at a significance level of p\u0026thinsp;\u0026le;\u0026thinsp;0.05.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cem\u003eVigna mungo\u0026nbsp;\u003c/em\u003eis emerging as an important legume crop after chickpea and pigeonpea due to its high nutritive value, short cropping duration, and multiple cuisines made from it. It is an important blend in the vegetarian diets for the majority of the population in South Asia due to its high nutritive value [11]. Besides, legume crops are being targeted to achieve the second sustainable development goal of \u0026ldquo;Zero Hunger\u0026rdquo; [24]; urdbean can be a vital component in this battle. The Vigna species are highly versatile and can be accommodated to different cropping systems and environments, making them a good choice for advanced genomics experiments. Looking into the wide adaptability as well as utility of the \u003cem\u003eVigna mungo\u003c/em\u003e, it has attracted the farmers as well as researchers. The recent enrichment in genetic and genomic resources of urdbean has paved the way for the functional characterization of many candidate genes. And gene expression profiling is a key experiment in this regard, which further depends on HKGs for the normalization to get statistically robust results, taking care of the technical artifacts. Therefore, to identify the most suitable HKGs for normalization of qRT-PCR data in urdbean, we have conducted a comprehensive expression profiling of 14 candidate housekeeping genes in 17 different tissues.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of the \u003cem\u003eV. mungo\u0026nbsp;\u003c/em\u003eorthologs of selected Arabidopsis HKG\u003c/strong\u003e\u003cstrong\u003es\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA total of 14 \u003cem\u003eV. mungo\u003c/em\u003e candidate genes corresponding to each selected Arabidopsis gene were identified through NCBI blastp with the \u0026ldquo;max hit = 1\u0026rdquo; parameter. These 14 genes were considered for local gene expression analysis using the qRT-PCR technique in the 17 selected tissues (Table 1). The percent similarity between the gene orthologs of urdbean and Arabidopsis ranges between 47.337 and 99.738, among which 12 out of 14 have a percent similarity of over 65 percent, reflecting the sequence conservation (Table 1). All of the 14 genes were found to be annotated as hypothetical proteins, as the gene prediction pipelines are based on homology-based annotations but lack experimental validation. Since the genome of \u003cem\u003eV. mungo\u003c/em\u003e was recently sequenced [3], this could be explained. The conserved domain analysis using NCBI-Batch CDD confirms the presence of the same domain as it was present in Arabidopsis ortholog (Figure S1). The higher percent similarity between the \u003cem\u003eV. mungo\u0026nbsp;\u003c/em\u003eHKGs and that of the Arabidopsis orthologs, as well as conserved domain analysis, advocates the functional as well as sequence conservation between them \u0026nbsp;(Figure S1).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSampling of diverse plant tissues\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA total of 17 plant tissues representing different stages from germination to pod development were selected for the expression analysis. The tissues used in the current study capture all developmental stages in urdbean. The diversity in the type and stages of the tissues reinforces the utility of the HKGs identified in the current study to the urdbean research community (Figure 1).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRNA quality assessment and p\u003c/strong\u003e\u003cstrong\u003erimer specificity\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe quality of the RNA used in this work was evaluated to guarantee a proper estimate of gene expression. The 260/280 absorbance ratio of RNA samples was evaluated using the spectrophotometer and only the samples with absorbance between 1.9 and 2.1 were considered. This ensures the quality as well as purity of the RNA isolated. Sharp and clear bands matching 18S rRNA and 28S rRNA on the denaturing gel revealed intact, undamaged RNA molecules, therefore indicating strong RNA integrity (Figure 2a). Further, the temperature gradient PCR for each primer set was conducted to identify the optimal annealing temperature. The length of the amplified PCR product was verified using 2% agarose gel electrophoresis using a DNA ladder (100 bp) as a reference. In our study, the amplicon size was found to range between 85 bp and 260 bp (Figure 2b).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGene expression stability analysis\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;in different developmental stages\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe qRT-PCR analysis of all selected 14 genes was conducted across the 17 tissues, representing the germination, vegetative, reproductive and pod development stages. The mean Ct (quantification cycle) values of all 14 genes from three biological groups were utilized for the estimation of stability across the three datasets (Figure 3). The candidate HKGs showed a large variation in mean Ct ranging from 10.84 (RBCS3B) to 39.82 (ISTL6 and CyP) while considering the total dataset (Figure 3). The whiskers represent the quartile where the Ct values for most of the HKGs fall, while the middle line represents the middle (median) Ct value in the dataset. The broader whisker means more variability in Ct value, while the smaller whisker supports the stable gene expression across the tissues (figure 3). Further, the data generated was analysed through the RefFinder algorithm, which also has other gene stability analysis algorithms such as geNorm, Normfinder, BestKeeper, and the comparative \u0026Delta;Ct method implemented in it.\u003c/p\u003e\n\u003cp\u003eThe geNorm algorithm calculates the mean stability value of the housekeeping genes by comparing their average pairwise variation with all other housekeeping genes considered. The geNorm algorithm relies on the principle that the expression ratio of two ideal internal control genes remains consistent across all samples, irrespective of experimental condition or cell type [25]. Therefore, the genes with the lowest average pairwise variation are considered to be the most stable. Using the geNorm algorithm, RPS34 and RHA genes were identified to be the most stable, with a mean stability value of 1.238, whereas the UFO gene was identified to be the least stable, with a mean stability value of 3.374. In the geNorm algorithm, any \u0026ldquo;M\u0026rdquo; value (geometric mean) below 1.5 is considered indicative of stable gene expression [13]. In our analysis, the top three genes i.e. RPS34, RHA and CyP were having a geometric mean stability value below 1.5. The mean stability value of all 14 genes as estimated through the geNorm algorithm is provided (Figure 4).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSimilarly, in the comparative \u0026Delta;Ct method, the \u0026Delta;Ct values of genes are compared like geNorm. In the \u0026Delta;Ct comparison method, the arithmetic mean of the standard deviation in the pairwise comparison is considered. The lower the arithmetic mean, the higher will be the stability of the housekeeping gene. Similar to the geNorm method, in the comparative \u0026Delta;Ct method, RPS34 and RHA were found to be the most stable housekeeping genes, with the arithmetic mean and standard deviation of 2.51 and 2.73, respectively (Figure 5). Likewise, the UFO gene was found to be the least stable housekeeping gene with an arithmetic mean of standard deviation of 4.95 (Figure 5). The gene stability ranking of all 14 genes as estimated through the \u0026Delta;Ct method algorithm is provided (Figure 5).\u003c/p\u003e\n\u003cp\u003eIn the BestKeeper algorithm, an index value based on the correlation analysis between the Ct values of the genes is calculated. Further, standard deviation, percent covariance explained, and power are calculated. Thereafter, the genes are ranked based on their stability metrics, with the most stable gene having the lowest standard deviation and highest correlation with others. In our analysis, the BestKeeper algorithm identified the USB1 gene to be the most stable, while EfTu was identified to be the least stable (Figure 6). The gene stability ranking of all 14 genes as estimated through the BestKeeper method algorithm is provided (Figure 6). The NormFinder algorithm estimates not only the overall variation of the candidate reference genes but also the variation between subgroups of samples. Unlike the geNorm and \u0026Delta;Ct comparison methods, NormFinder assesses the expression stability of each candidate independently. In our analysis, the RPS34 and ACT2 were found to be most stable using the NormFinder algorithm (Figure 7). The gene stability ranking of all 14 genes as estimated through the \u0026Delta;Ct method algorithm is provided (Figure 7).\u003c/p\u003e\n\u003cp\u003eThe RefFinder algorithm compares the weights estimated for each reference gene through the aforementioned four algorithms to calculate the geometric mean of the weights and make a comprehensive ranking. Through the RefFinder, we identified RPS34 and RHA as the best and most stable housekeeping genes, while EfTu and RBCS3B were the least suitable (Figure 8). Among the four algorithms implemented in RefFinder, RPS34 was found to be the best housekeeping gene in the geNorm, \u0026Delta;Ct comparison method, as well as the NormFinder algorithm, except for the BestKeeper algorithm, where USB1 was found to be the best as depicted in the Venn diagram when the top three most stable genes were compared (Figure 9). Similarly, RHA was among the top three in geNorm, the \u0026Delta;Ct comparison method, and BestKeeper when the top three most stable genes from each of these algorithms are compared (Figure 9). Overall, the RPS34 and RHA genes were also found to be most stable and suitable to be used as control genes for normalization study through the RefFinder algorithm. This was also supported by the whisker plot of the Ct values (Figure 3). The gene stability ranking (Better--Good--Average) of the 14 genes in different developmental stages is provided (Table 2).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGene expression stability analysis\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;in different abiotic stress\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe expression pattern of all of the 14 genes was further explored under four abiotic stresses \u003cem\u003eviz.\u0026nbsp;\u003c/em\u003esalt, drought, aluminium, \u0026nbsp;and cold stress in the germination stage tissues (plumule, radicle and cotyledon). In the four abiotic stresses, RPS34 was found to be the most stable, followed by USB1, while RBCS3B was least stable using the geNorm algorithm (Figure S2). Likewise, using the \u0026Delta;Ct method, ACT2 and RPS34 genes were identified to be the most stable HKGs, while RBCS3B was found to be the least stable (Figure S3). In the four abiotic stresses, USB1 was found to be the most stable followed by DNAJ, while EfTu was least stable using the BestKeeper algorithm (Figure S4). Likewise, using ACT2 and DNAJ genes were identified to be the most stable HKGs, while RBCS3B was found to be the least stable using the NormFinder algorithm (Figure S5). Comprehensively, RefFinder identifies ACT2 and RPS34 \u0026nbsp;as the most stable HKGs (Figure S6). The gene stability ranking (Better--Good--Average) of the 14 genes in different abiotic stresses in seedling stages is provided (Table S3).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eValidation of the candidate HKGs identified in the current study\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor validation of the best HKGs identified in the current study, we used the Ct value of RPS34, RHA and ACT2 HKGs to normalize the expression of 5 previously reported genes i.e., CAM, PR1, HSP70, NAC, and DEF genes, along with an external control gene, i.e., Actin (Table S2). For the experiment, the tissues and primers were generated and used as previously reported [23]. The gene expression pattern of all 5 test genes was in concurrence with the previous reports [23]. The larger deviation in the relative fold change value of the test gene when normalized with ACT as a control was clearly evident (Figure 10a-e). On the other hand, when the test genes were normalized with the HKGs identified in the current study, a substantial variation was observed. The normalized relative gene expression level clearly reflects the stability as well as the suitability of using RPS34, RHA, and ACT2 genes over ACT in the qRT-PCR experiment (Figure 10a-e).\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eThe gene expression is a highly spatiotemporal event influenced by both external and internal factors. The profiling of gene expression patterns is among the most sought-after preliminary experiments in both animal and plant sciences [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Majorly done by qRT-PCR, there are chances of error due to pipetting and RNA quality affecting the imprecise estimation of gene expression level [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. These probable errors are taken care of by using a good reference gene for the normalization of the data generated in the qRT-PCR experiments [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Specifically, the expression of the reference gene needs to be consistent across the spatiotemporal conditions and the different development stages. Therefore, while doing a gene expression analysis, selecting the most trustworthy reference genes has become a crucial research endeavor [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. Limited efforts have been made to identify the suitable genes in legume crops, including \u003cem\u003eLupin angustifolius\u003c/em\u003e, \u003cem\u003eVigna angularis\u003c/em\u003e, \u003cem\u003eLens culinaris\u003c/em\u003e, \u003cem\u003eVicia faba\u003c/em\u003e, and \u003cem\u003eCicer arietinum\u003c/em\u003e [\u003cspan additionalcitationids=\"CR30 CR31 CR32\" citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eSo far only a handful of studies involving expression profiling have been conducted in \u003cem\u003eV. mungo\u003c/em\u003e for different biotic and abiotic stresses [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. The current study holds the potential to facilitate future research aiming to accurately explore the expression profiling of any selected gene across the developmental stages. To ensure the preciseness of the study, due care has been taken right from RNA isolation to qRT-PCR. Further, the majority of the candidate genes selected for the present study were novel to \u003cem\u003eV. mungo\u003c/em\u003e and have not been explored to the best of our knowledge.\u003c/p\u003e \u003cp\u003eA clear consensus on which HKGs to be used for data normalization in \u003cem\u003eV. mungo\u003c/em\u003e is lacking. Therefore, comparing the expression Ct values obtained from profiling multiples samples using different statistical algorithms will pave the way forward. Hence, the current study was designed to profile a set of 14 HKGs in 17 different tissues representing different stages right from seed germination to seed maturity. The genes used in the current study represent diverse metabolic and developmental pathways, including floral development, cellular cytoskeleton, RNA and protein processing, and protein trafficking. Diversity in the metabolic pathways and developmental stages used in the current study makes the study more reliable and of better utility for the \u003cem\u003eV. mungo\u003c/em\u003e researcher community. Further, five different algorithms, including geNorm [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e], NormFinder [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e], BestKeeper [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e], the comparative ΔCt method [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e], and RefFinder were utilized in the current study [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eInterestingly, two genes, RPS34 and RHA, were found to be the two most stable HKGs using the different statistical algorithms (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), with the exception of the BestKeeper algorithm. Likewise, RBSC3, EfTu and UFO were found to be least stable across the different developmental stages in the current study (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The ribosomal proteins have been reported to be major HKGs in all domains of life [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e]. Since RPS34 is an integral part of protein biosynthesis machinery, its stable expression could be explained. Surprisingly, any of the ribosomal proteins, although expressing stably, have not been explored much as a suitable HKGs for qRT-PCR normalization. The current study advocates to exploring the ribosomal proteins for gene normalization studies in other species. Likewise, RNA helicases are enzymes that are crucial for RNA metabolism. They unwind RNA secondary structures and/or displace RNA-binding proteins and are essential for various cellular processes involving RNA, including pre-mRNA splicing, stability and decay [\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. Since RNA helicase is ubiquitous and performs several functions vital for cell growth and development, their listing in the most stable HKGs is not surprising. Despite modest changes in the ranking of other genes, likely due to variation in the algorithms used, the findings clearly indicated that RPS34 and RHA consistently were most stable. Although ACT2 outperforms the RPS34 gene under all four abiotic stress conditions, RPS34 was still among the top stable HKGs.\u003c/p\u003e \u003cp\u003eKundu and his coworkers put an effort to identify the best suitable HKGs in \u003cem\u003eV. mungo\u003c/em\u003e under biotic (YMV) and abiotic stress (drought) but with only 7 HKGs in two stress conditions [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Comparably, our study not only uses more HKGs but also represents an array of developmental stages (17 tissues) along with the 4 stress conditions (salt, drought, aluminium and heat). The candidate HKGs identified in the current study were further demonstrated to be stably expressed and normalize the Ct value of the 5 test genes under independent experiment conditions. The HKGs identified in the current study i.e. RPS34, RHA, and ACT2 genes hold a greater potential to be widely used as HKGs in any gene expression experiment. Since this is the most comprehensive report on the identification of HKGs in Vigna species, the HKGs identified in the current study might be useful to the other Vigna species.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eThe current study explores the stability pattern of 14 HKGs representing diverse metabolic and development pathways in 17 different developmental stages from germination to seed maturity in \u003cem\u003eVigna mungo\u003c/em\u003e. The Ct values obtained were implemented to get a gene-wise ranking based on geNorm, NormFinder, BestKeeper, comparative ΔCt method and the RefFinder algorithm. Based on the comprehensive ranking, RPS34 and RHA were found to be most suitable for qRT-PCR data normalization in different developmental stages, while ACT2 and RPS34 were suitable under abiotic stress conditions. Based on the results obtained, the current study advocates using a combination of HKGs (RPS34 and RHA for developmental tissues and ACT2 and RPS34 for stress-treated tissues) for normalization of qRT-PCR data. The reference genes suggested here will ensure reliability and uniformity in qRT-PCR findings. Further, the expression of any gene is spatiotemporal and is influenced by both internal and external factors. And the genotypic background also plays an important role in it. The key findings of the current study need to be reconfirmed in independent background, including germplasm lines, to get more comprehensive insight. Parallelly, these genes further may also be tested and used in other related Vigna species such as \u003cem\u003eV. radiata\u003c/em\u003e, \u003cem\u003eV. unguiculata\u003c/em\u003e and \u003cem\u003eV. angularis\u003c/em\u003e, which share closer ancestry to \u003cem\u003eV. mungo\u003c/em\u003e.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eKK has conceptualized the whole study. PS, PK, SG, AAC, VK, KT performed the experiment and collected the plant tissues. PS, and AD performed the qRT-PCR. MR and KK performed statistical analysis. SKJ, KK, AD, MR, GPD wrote the draft manuscript. MR, GPD and KK edited the manuscript to prepare final draft. MR and GPD supervised and administered the whole study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe acknowledge the support provided by the Director, ICAR-IIPR, Kanpur. The research work carried out in this manuscript is a part of Institute funded project.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConflict of Interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThere is no conflict of interest, either financial or personal, with other people or organizations that could inappropriately influence or perceived to influence this work.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data that support the findings of this study are available in the supporting information of this article\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eSingh, N. K. et al. The first draft of the pigeonpea genome sequence. \u003cem\u003eJournal of plant biochemistry and biotechnology. \u003c/em\u003e\u003cstrong\u003e21\u003c/strong\u003e:98-112 (2012). https://doi.org/10.1007/s13562-011-0088-8\u003c/li\u003e\n\u003cli\u003eVarshney, R. K. et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. \u003cem\u003eNature biotechnology. \u003c/em\u003e31(\u003cstrong\u003e3\u003c/strong\u003e):240-246 (2013). https://doi.org/10.1038/nbt.2491 \u003c/li\u003e\n\u003cli\u003eJunaid, A. et al. Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo. \u003cem\u003eLife Science Alliance. \u003c/em\u003e7(\u003cstrong\u003e1\u003c/strong\u003e), (2024). https://doi.org/10.26508/lsa.202302074\u003c/li\u003e\n\u003cli\u003eGaikwad, K. et al. The chromosome-scale genome assembly of cluster bean provides molecular insight into edible gum (galactomannan) biosynthesis family genes. \u003cem\u003eScientific Reports.\u003c/em\u003e 13(\u003cstrong\u003e1\u003c/strong\u003e):9941 (2023). https://doi.org/10.1038/s41598-023-33762-3\u003c/li\u003e\n\u003cli\u003eAmeline-Torregrosa, C. et al. Transcriptomic approaches to unravel plant\u0026ndash;pathogen interactions in legumes. \u003cem\u003eEuphytica.\u003c/em\u003e\u003cstrong\u003e147\u003c/strong\u003e:25-36 (2006). https://doi.org/10.1007/s10681-006-6767-1 \u003c/li\u003e\n\u003cli\u003eGarg, R \u0026amp; Jain, M. Transcriptome analyses in legumes: A resource for functional genomics. \u003cem\u003eThe Plant Genome. \u003c/em\u003e6(\u003cstrong\u003e3\u003c/strong\u003e) (2013). https://doi.org/10.3835/plantgenome2013.04.0011\u003c/li\u003e\n\u003cli\u003eAbdelrahman, M. et al. Legume genetic resources and transcriptome dynamics under abiotic stress conditions. \u003cem\u003ePlant, cell \u0026amp; environment.\u003c/em\u003e 41(\u003cstrong\u003e9\u003c/strong\u003e):1972-1983 (2018). https://doi.org/10.1111/pce.13123\u003c/li\u003e\n\u003cli\u003eNicot, N. et al. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. \u003cem\u003eJournal of Experimental Botany.\u003c/em\u003e 56(\u003cstrong\u003e421\u003c/strong\u003e):2907-2914 (2005). https://doi.org/10.1093/jxb/eri285\u003c/li\u003e\n\u003cli\u003eUdvardi, M. K. Eleven golden rules of quantitative RT-PCR. \u003cem\u003eThe Plant Cell. \u003c/em\u003e20(\u003cstrong\u003e7\u003c/strong\u003e):1736-1737 (2008). https://doi.org/10.1105/tpc.108.061143\u003c/li\u003e\n\u003cli\u003eBustin, S. A. et al. MIQE precis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. \u003cem\u003eBMC Molecular Biology.\u003c/em\u003e\u003cstrong\u003e11\u003c/strong\u003e:1-5 (2010). https://doi.org/10.1186/1471-2199-11-74\u003c/li\u003e\n\u003cli\u003eManu, B. et al. Genomic designing for abiotic stress tolerance in mungbean and urdbean. In Genomic Designing for Abiotic Stress Resistant Pulse Crops (pp. 271-343). Cham: Springer International Publishing (2022). https://doi.org/10.1007/978-3-030-91039-6_8\u003c/li\u003e\n\u003cli\u003eAndersen, C. L. et al. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. \u003cem\u003eCancer research. \u003c/em\u003e64(\u003cstrong\u003e15\u003c/strong\u003e):5245-5250 (2004). https://doi.org/10.1158/0008-5472.CAN-04-0496\u003c/li\u003e\n\u003cli\u003eVandesompele, J. et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. \u003cem\u003eGenome Biology \u003c/em\u003e\u003cstrong\u003e3\u003c/strong\u003e:1-12 (2002). https://doi.org/10.1186/gb-2002-3-7-research0034\u003c/li\u003e\n\u003cli\u003ePfaffl, M. W. et al. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper\u0026ndash;Excel-based tool using pair-wise correlations. \u003cem\u003eBiotechnology letters. \u003c/em\u003e\u003cstrong\u003e26\u003c/strong\u003e:509-515 (2004). https://doi.org/10.1023/b:bile.0000019559.84305.47\u003c/li\u003e\n\u003cli\u003eSilver, N. et al. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. \u003cem\u003eBMC Molecular Biology.\u003c/em\u003e\u003cstrong\u003e7\u003c/strong\u003e:1-9 (2006). https://doi.org/10.1186/1471-2199-7-33\u003c/li\u003e\n\u003cli\u003eXie, F. et al. RefFinder: a web-based tool for comprehensively analyzing and identifying reference genes. \u003cem\u003eFunctional \u0026amp; Integrative Genomics.\u003c/em\u003e 23(\u003cstrong\u003e2\u003c/strong\u003e):125 (2023). https://doi.org/10.1007/s10142-023-01055-7\u003c/li\u003e\n\u003cli\u003eDantas, B. F. et al. Germination, initial growth and cotyledon protein content of bean cultivars under salinity stress. \u003cem\u003eRevista Brasileira de Sementes.\u003c/em\u003e\u003cstrong\u003e29\u003c/strong\u003e:106-110 (2007). https://doi.org/10.1590/S0101-31222007000200014\u003c/li\u003e\n\u003cli\u003eRoso, R. et al. Effects of Aluminum and pH on Germination of Echium plantagineum L. \u003cem\u003eJournal of Agricultural Studies\u003c/em\u003e. 8(\u003cstrong\u003e1\u003c/strong\u003e):313-327 (2020). https://doi.org/10.5296/jas.v8i1.16179\u003c/li\u003e\n\u003cli\u003eIslam, M. R. Evaluation of Polyethylene Glycol (PEG-6000) Induced Drought Stress Tolerant Mungbean Genotypes by using Correlation, Principal Component, Hierarchical Clustering and Multi-Trait Genotype-Ideotypes Distance Index Analysis. \u003cem\u003eCurrent Applied Science And Technology.\u003c/em\u003e e0260871-e0260871 (2023). https://doi.org/10.55003/cast.2024.260871\u003c/li\u003e\n\u003cli\u003eManasa, L. S. et al. Cold tolerance mechanisms in mungbean (Vigna radiata L.) genotypes during germination. \u003cem\u003eAgriculture.\u003c/em\u003e 13(\u003cstrong\u003e2\u003c/strong\u003e):315 (2023). https://doi.org/10.3390/agriculture13020315\u003c/li\u003e\n\u003cli\u003eKumar, K. et al. Identification and characterization of NHX gene family for their role under salt stress in Vigna mungo. \u003cem\u003ePhysiologia Plantarum.\u003c/em\u003e 176(\u003cstrong\u003e5\u003c/strong\u003e):e14563 (2024). https://doi.org/10.1111/ppl.14563\u003c/li\u003e\n\u003cli\u003eMarchler-Bauer, A. et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. \u003cem\u003eNucleic Acids Research.\u003c/em\u003e 45(\u003cstrong\u003eD1\u003c/strong\u003e), D200-D203 (2017). https://doi.org/10.1093/nar/gkw1129\u003c/li\u003e\n\u003cli\u003eKumar, S. et al. Gene expression and biochemical profiling of contrasting Vigna mungo genotypes against Mungbean Yellow Mosaic India Virus (MYMIV). \u003cem\u003eJournal of Food Legumes.\u003c/em\u003e 35(\u003cstrong\u003e2\u003c/strong\u003e):107-116 (2022).\u003c/li\u003e\n\u003cli\u003eKumar, K. et al. Single trait versus principal component based association analysis for flowering related traits in pigeonpea. \u003cem\u003eScientific Reports.\u003c/em\u003e 12(\u003cstrong\u003e1\u003c/strong\u003e):10453 (2022). https://doi.org/10.1038/s41598-022-14568-1\u003c/li\u003e\n\u003cli\u003eVijayakumar S \u0026amp; Sakuntala M. Validation of reference gene stability for normalization of RT-qPCR in Phytophthora capsici Leonian during its interaction with Piper nigrum L. \u003cem\u003eScientific Reports.\u003c/em\u003e 14(1):7331 (2024). https://doi.org/10.1038/s41598-024-58139-y\u003c/li\u003e\n\u003cli\u003eNolan, T. Hands, R. E. \u0026amp; Bustin, S. A. Quantification of mRNA using real-time RT-PCR. \u003cem\u003eNature Protocols.\u003c/em\u003e 1(\u003cstrong\u003e3\u003c/strong\u003e):1559-1582 (2006). https://doi.org/10.1038/nprot.2006.236\u003c/li\u003e\n\u003cli\u003eHu, R. et al. Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR. \u003cem\u003eBMC Molecular Biology.\u003c/em\u003e 10:1-12 (2009). https://doi.org/10.1186/1471-2199-10-93\u003c/li\u003e\n\u003cli\u003ede Vega-Bartol, J. J. et al. Normalizing gene expression by quantitative PCR during somatic embryogenesis in two representative conifer species: Pinus pinaster and Picea abies. \u003cem\u003ePlant Cell Reports. \u003c/em\u003e32:715-729 (2013). https://doi.org/10.1007/s00299-013-1407-4\u003c/li\u003e\n\u003cli\u003eTaylor, C. M. et al. Identifying stable reference genes for qRT-PCR normalisation in gene expression studies of narrow-leafed lupin (\u003cem\u003eLupinus angustifolius\u003c/em\u003e L.). \u003cem\u003ePLoS One.\u003c/em\u003e 11(2):e0148300 (2016). https://doi.org/10.1371/journal.pone.0148300\u003c/li\u003e\n\u003cli\u003eChi, C. et al. Selection and validation of reference genes for gene expression analysis in Vigna angularis using quantitative real-time RT-PCR. \u003cem\u003ePloS one.\u003c/em\u003e 11(12):e0168479 (2016). https://doi.org/10.1371/journal.pone.0168479\u003c/li\u003e\n\u003cli\u003eSinha, R. Sharma, T. R., Singh, A. K. Validation of reference genes for qRT-PCR data normalisation in lentil (Lens culinaris) under leaf developmental stages and abiotic stresses. \u003cem\u003ePhysiology and Molecular Biology of Plants.\u003c/em\u003e 25(\u003cstrong\u003e1\u003c/strong\u003e):123-134 (2019). https://doi.org/10.1007/s12298-018-0609-1\u003c/li\u003e\n\u003cli\u003eGutierrez, N. et al. Assessment of candidate reference genes for expression studies in Vicia faba L. by real-time quantitative PCR. \u003cem\u003eMolecular Breeding.\u003c/em\u003e\u003cstrong\u003e28\u003c/strong\u003e:13-24 (2011). https://doi.org/10.1007/s11032-010-9456-7\u003c/li\u003e\n\u003cli\u003eReddy, D. S. et al. Identification and validation of reference genes and their impact on normalized gene expression studies across cultivated and wild Cicer species. \u003cem\u003ePloS one.\u003c/em\u003e 11(2):e0148451 (2016). https://doi.org/10.1371/journal.pone.0148451 \u003c/li\u003e\n\u003cli\u003eKundu, A. et al. Complex molecular mechanisms underlying MYMIV-resistance in Vigna mungo revealed by comparative transcriptome profiling. \u003cem\u003eScientific Reports.\u003c/em\u003e 9(\u003cstrong\u003e1\u003c/strong\u003e):8858 (2019). https://doi.org/10.1038/s41598-019-45383-w\u003c/li\u003e\n\u003cli\u003eKumar J \u0026amp; Kumar V. Ribosome proteins\u0026mdash;Their balanced production. In: Emerging Concepts in Ribosome Structure, Biogenesis, and Function (pp. 47-87). Academic Press (2021). https://doi.org/10.1016/B978-0-12-816364-1.00003-2\u003c/li\u003e\n\u003cli\u003eBourgeois, C. F., Mortreux, F. \u0026amp; Auboeuf, D. The multiple functions of RNA helicases as drivers and regulators of gene expression. \u003cem\u003eNature Reviews Molecular Cell Biology\u003c/em\u003e\u003cem\u003e.\u003c/em\u003e 17(\u003cstrong\u003e7\u003c/strong\u003e):426-438 (2016). https://doi.org/10.1038/nrm.2016.50\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003e\u003cstrong\u003eTable 1:\u003c/strong\u003e \u003cstrong\u003eList of genes used in the current study and their function\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"99%\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eS.no.\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eArabidopsis gene\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eGene name\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eVigna mungo ortholog used in the current study\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e\u003cstrong\u003ePercent homology with Arabidopsis homolog\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eFunction\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003col\u003e\n \u003cli\u003e\u0026nbsp;\u003c/li\u003e\n \u003c/ol\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT1G30950.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eUnusual Floral Organ (UFO)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ23601.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e66.743\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes Unusual floral organ genes which is required for the proper identity of the floral meristem and is constitutively expressed in the floral tissues\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003col start=\"2\"\u003e\n \u003cli\u003e\u0026nbsp;\u003c/li\u003e\n \u003c/ol\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT5G62690.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eTUBULIN BETA CHAIN 2 (TUB2)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVY95928.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e95.778\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes tubulin beta-2/beta-3 chain, which are components of microtubules in the cytoskeleton and play essential roles in various cellular processes\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003col start=\"3\"\u003e\n \u003cli\u003e\u0026nbsp;\u003c/li\u003e\n \u003c/ol\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT3G18780.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eACTIN 2 (ACT2)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVY99582.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e92.042\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes an actin that is constitutively expressed in vegetative structures\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003col start=\"4\"\u003e\n \u003cli\u003e\u0026nbsp;\u003c/li\u003e\n \u003c/ol\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT1G58050.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eRNA helicase (RHA)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ12136.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e47.337\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes RNA helicase that unwind RNA secondary structures and/or displace RNA-binding proteins in an ATP-dependent manner, that is essential for various cellular processes involving RNA\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003col start=\"5\"\u003e\n \u003cli\u003e\u0026nbsp;\u003c/li\u003e\n \u003c/ol\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT1G13340.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ09401.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e51.507\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eISTL6 gene is crucial for the proper trafficking and degradation of membrane proteins within the cell. Proper functioning of this pathway is vital for maintaining cellular homeostasis and responding to environmental stressors\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e6.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT2G36130.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eCyclophilin (CyP)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ05807.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e87.805\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eCyclophilins are a family of highly conserved proteins with peptidyl-prolyl cis-trans isomerase (PPIase) activity, which catalyzes the isomerization of peptide bonds at proline residues that is crucial for protein folding, and function\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e7.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT2G20560.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ01865.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e79.769\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eDNAJ heat shock family protein are molecular chaperones involved in various cellular processes and play essential roles in protein folding, stabilization, degradation, and trafficking\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e8.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT5G38410.3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eRubisco small subunit 3b (RBCS3B)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ13955.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e67.935\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes a member of the Rubisco small subunit (RBCS) multi-gene family\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e9.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT5G58990.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003e28S ribosomal S34 protein (RPS34)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ02079.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e84.286\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003e(RPS34) is a component of the small subunit of the ribosome in eukaryotes that are responsible for protein synthesis, and their subunits (small and large) work together to translate mRNA into polypeptides\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e10.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT4G27320.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eUniversal Stress Protein 21 (USP21)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ04305.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e70.892\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eProtein that contains a universal stress protein domain. They are essential for plant stress tolerance as well as cellular homeostasis\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e11.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT5G51170.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eU6 snRNA phosphodiesterase like protein (USB1)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ06716.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e65.351\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes U6 snRNA phosphodiesterase-like protein, that are involved in RNA processing, specifically within the spliceosome machinery, where they contribute to pre-mRNA splicing\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e12.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT1G07920.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eElongation Factor-TU (EfTu)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVY95914.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e96.197\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eGTP binding Elongation factor Tu family protein\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e13.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT4G05320.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eUbiquitin 10 (UBQ10)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVZ24872.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e99.738\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes Ubiquitin 10 that is a member of the ubiquitin protein family, which plays a central role in post-translational modification of proteins through ubiquitination\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 4px;\"\u003e\n \u003cp\u003e14.\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eAT4G27960.2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 12px;\"\u003e\n \u003cp\u003eUbiquitin Conjugating Enzyme 9 (UBQE9)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003eWVY93904.1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 10px;\"\u003e\n \u003cp\u003e94.667\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 51px;\"\u003e\n \u003cp\u003eEncodes UBIQUITIN CONJUGATING ENZYME 9 is a key enzyme in the ubiquitination pathway, primarily involved in SUMOylation\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003eTable 2: Gene stability ranking of the 14 HKGs used in the current study through different algorithms using 17 developmental tissues\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"100%\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"15\" valign=\"top\" style=\"width: 100px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eRanking Order (Better--Good--Average)\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003eMethod\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003e1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003e6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003e8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003e12\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003e13\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003e14\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003eDelta CT\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRPS34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eRHA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eCyP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eACT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQC9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQ10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUSB1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eTUB2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUSP21\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eEfTu\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRBCS3B\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eUFO\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003eBestKeeper\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eUSB1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eRHA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eACT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUFO\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eCyP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eRPS34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQC9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUBQ10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eTUB2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUSP21\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRBCS3B\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eEfTu\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003eNormFinder\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRPS34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eACT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eCyP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eRHA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQC9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQ10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUSB1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eTUB2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUSP21\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eEfTu\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRBCS3B\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eUFO\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003egeNorm\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRPS34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eRHA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eCyP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eACT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUBQC9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUBQ10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUSB1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eTUB2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUSP21\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eEfTu\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRBCS3B\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eUFO\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 9px;\"\u003e\n \u003cp\u003eComprehensive ranking\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRPS34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eRHA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eACT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eDNAJ\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eCyP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eUSB1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUBQC9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 6px;\"\u003e\n \u003cp\u003eISTL6\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 8px;\"\u003e\n \u003cp\u003eUBQ10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eTUB2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUFO\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 5px;\"\u003e\n \u003cp\u003eUSP21\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eEfTu\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 7px;\"\u003e\n \u003cp\u003eRBCS3B\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Housekeeping gene, Normalization, qRT-PCR, urdbean","lastPublishedDoi":"10.21203/rs.3.rs-5819031/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5819031/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe qRT-PCR has emerged as the most sought-after local gene expression profiling technique. The absolute quantification of the mRNA is highly affected by the handling artifacts. In this regard, the internal control genes serve a vital role in negating the handling errors. So far, the most stable housekeeping genes to be used for the normalization of qRT-PCR data are lacking in \u003cem\u003eV. mungo\u003c/em\u003e, commonly known as urdbean or blackgram. Therefore, in the current study we intend to explore the stability of the UFO, TUB2, ACT2, RHA, ISTL6, CyP, DNAJ, RBCS3B, RPS34, USP21, USB1, EfTu, UBQ10, and UBQE9 genes of \u003cem\u003eV. mungo\u003c/em\u003e in 17 different developmental stages and 4 different abiotic stresses. The well-acknowledged algorithms such as geNorm, NormFinder, BestKeeper, the comparative ΔCt method and comprehensive RefFinder were used to assess the expression stability. According to the comprehensive ranking, RPS34 and RHA were identified as the most appropriate genes for qRT-PCR data normalization throughout all developmental stages, whereas ACT2 and RPS34 were optimal under abiotic stress conditions. Further, the suitability of these three candidate housekeeping genes for normalization has been demonstrated in an independent experiment. Based on the findings obtained, the present research advocates the use of a combination of HKGs (RPS34 and RHA for developing tissues and ACT2 and RPS34 for stress-treated tissues) for the normalization of qRT-PCR data. The proposed reference genes will guarantee the reliability and consistency of qRT-PCR results.\u003c/p\u003e","manuscriptTitle":"Identification of best housekeeping gene for normalization in qRT-PCR data under different development stages and stress conditions in Vigna mungo","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-04-07 06:59:43","doi":"10.21203/rs.3.rs-5819031/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-05-19T07:06:51+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-05-09T12:51:54+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-05-01T18:24:36+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"55601850229176891707027745026188832803","date":"2025-04-18T04:42:31+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"313327421580686674079026752602968553877","date":"2025-04-16T18:16:44+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-04-16T15:49:09+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-04-10T10:42:20+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2025-03-29T09:43:22+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"7e68063f-814a-47c9-85fc-9bd38f8477eb","owner":[],"postedDate":"April 7th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":46392665,"name":"Biological sciences/Biological techniques"},{"id":46392666,"name":"Biological sciences/Biotechnology"},{"id":46392667,"name":"Biological sciences/Plant sciences"}],"tags":[],"updatedAt":"2025-09-29T15:58:55+00:00","versionOfRecord":{"articleIdentity":"rs-5819031","link":"https://doi.org/10.1038/s41598-025-10295-5","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2025-09-26 15:56:52","publishedOnDateReadable":"September 26th, 2025"},"versionCreatedAt":"2025-04-07 06:59:43","video":"","vorDoi":"10.1038/s41598-025-10295-5","vorDoiUrl":"https://doi.org/10.1038/s41598-025-10295-5","workflowStages":[]},"version":"v1","identity":"rs-5819031","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5819031","identity":"rs-5819031","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}