Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The Complete Mitochondrial Genomes of Penthe kochi (Coleoptera: Tetratomidae) with Its Phylogenetic Implications
Previous Article in Journal
Improvement of Quality and Disease Resistance for a Heavy-Panicle Hybrid Restorer Line, R600, in Rice (Oryza sativa L.) by Gene Pyramiding Breeding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of the Optimal Quantitative RT-PCR Reference Gene for Paper Mulberry (Broussonetia papyrifera)

Zhejiang Key Laboratory of Forest Genetics and Breeding, Zhejiang Academy of Forestry, Hangzhou 310023, China
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2024, 46(10), 10779-10794; https://doi.org/10.3390/cimb46100640
Submission received: 13 August 2024 / Revised: 10 September 2024 / Accepted: 24 September 2024 / Published: 26 September 2024
(This article belongs to the Section Molecular Plant Sciences)

Abstract

:
Paper Mulberry (Broussonetia papyrifera) possesses medicinal, economic, and ecological significance and is extensively used for feed production, papermaking, and ecological restoration due to its ease of propagation, rapid growth rate, and strong stress resistance. The recent completion of the sequencing of the Paper Mulberry genome has prompted further research into the genetic breeding and molecular biology of this important species. A highly stable reference gene is essential to enhance the quantitative analysis of functional genes in Paper Mulberry; however, none has been identified. Accordingly, in this study, the leaves, stems, roots, petioles, young fruits, and mature fruits of Paper Mulberry plants were selected as experimental materials, and nine candidate reference genes, namely, α-TUB1, α-TUB2, β-TUB, H2A, ACT, DnaJ, UBQ, CDC2, and TIP41, were identified by RT-qPCR. Their stability was assessed using the geNorm, Normfinder, Delta Ct, BestKeeper, and RefFinder algorithms, identifying ACT and UBQ as showing the greatest stability. The expression of BpMYB090, which regulates the production of trichomes, was examined in the leaves of plants of the wild type (which have more trichomes) and mutant (which have fewer trichomes) at various developmental stages to validate the results of this study. As a result, their identification addresses a critical gap in the field of Paper Mulberry research, providing a solid foundation for future research that will concentrate on the characterization of pertinent functional genes in this economically valuable species.

1. Introduction

The deciduous tree or shrub, Paper Mulberry (Broussonetia papyrifera), belongs to the Moraceae family and is distributed over various regions of China, East Asia, Southeast Asia, and northern South Asia. It is considered a pioneering plant and a significant representative of local tree species [1,2]. Due to their rapid reproduction, growth, and remarkable resistance to stressful conditions, Paper Mulberry trees are extensively used in environmental restoration, paper manufacturing, and animal feed production. The trees contain high levels of protein (between 18 and 24% of the dry weight), similar to those of alfalfa hay and greater than those of grass or leguminous roughage. It is generally regarded as a woody feed with substantial development opportunities [3]. The bark of Paper Mulberry possesses elongated fibers with beneficial morphological characteristics, higher cellulose content, and lower lignin levels. It may be easily pulped, making it suitable for weaving or papermaking [4,5]. The tree is also a valuable pioneer in establishing or restoring plant communities. As it can reproduce both asexually through root sprouting and via the production of seeds, the rapid growth of Paper Mulberry roots and trees in damaged ecosystems is important as a means of preventing soil erosion and improving the quality of the soil microenvironment [4,5]. Paper Mulberry is extensively dispersed, highly versatile, possesses a small genome, and displays various phenotypic characteristics. Therefore, it is an exemplary model organism for a functional study of genes in woody plants. The sequenced genome of the Paper Mulberry has recently been published, allowing the application of biotechnological strategies for the genetic breeding and molecular study of this important woody plant species.
Quantifying gene expression is necessary for research centered on gene functions. Among the techniques established for such investigations, RT-qPCR is the most commonly employed due to its superior specificity, reproducibility, and sensitivity, facilitating quantification with a high level of precision [6,7]. The accuracy of RT-qPCR is influenced by factors such as the quality of DNA and RNA, amplification efficiency, and the original sample size. As a result, reference genes are essential to ensure precision and reduce the possibility of bias in interpreting gene expression patterns [8,9]. Optimal internal reference genes exhibit stable intracellular expression without interference by external factors. Housekeeping genes are commonly used, including genes encoding actin (ACT), TIP41-like protein (TIP41), β-Tubulin (β-TUB), ubiquitin (UBQ), and 18S ribosomal RNA (18S rRNA) [10,11,12,13,14]. It is essential that reference genes are present across several developmental stages and tissues and that they remain unaffected by environmental influences and unidentified stressors [15]. Therefore, it is essential to identify genes that meet these criteria to ensure the most accurate assessment of gene expression.
Reference gene expression has been extensively characterized in species including Arabidopsis thaliana [16], Oryza sativa [17], Zingiber officinale [18], Liquidambar formosana [19], Sorghum bicolor [20], Citrullus lanatus [21], Cucumis melo [22], Paeonia suffruticosa Andr. [23], Populus [24], Eucalyptus [25], and Dendrobium huoshanense [26]. Results from these studies suggest that, even within a given species, the most appropriate reference genes vary among tissue types or stressor conditions. In Paper Mulberry, BpDOUB and BpNADH were recognized as suitable reference genes in roots and leaves under Cd stress; however, no reference genes showing stable expression across several tissues have been discovered [27]. The instability in gene expression may lead to substantial inaccuracies in data related to different experimental settings, affecting the precision of target gene detection initiatives. A swift surge of studies has emerged regarding the functionality of Paper Mulberry genes following the completion of its entire genome sequencing [28]. Therefore, it is necessary to identify appropriate reference genes expressed across many tissues in Paper Mulberry to ensure the reliability of expression data.
Six different Paper Mulberry tissue types (leaves, stems, roots, petioles, young fruits, and mature fruits) were selected as experimental materials. Using genomic and transcriptomic data, ACT, UBQ, α-TUB1, α-TUB2, β-TUB, TIP41, CDC2, H2A, and DnaJ were initially identified as reference genes through the cloning of homologous genes from the Paper Mulberry genome. The expression of these genes in the six tissue types was analyzed using RT-qPCR, and stability was assessed using the GeNorm, NormFinder, Delta Ct, and BestKeeper algorithms, followed by the geometric mean and the ReFinder network program to determine the gene combination with the greatest stability. The expression levels of the two genes demonstrating the highest and lowest stability were later employed to evaluate the expression of the BpMYB090 gene, which regulates trichome development, in leaves at different developmental stages and in mutant plants with reduced trichome density as an additional measure of stability. The discovered genes collectively provide a crucial molecular foundation for precise RT-qPCR-based assessments of Paper Mulberry gene expression across many tissues and organs, thus underpinning functional research and transgenic breeding initiatives.

2. Materials and Methods

2.1. Plant Materials

Specimens of Paper Mulberry were cultivated at the Zhejiang Academy of Forestry. In August 2023, three trees of comparable height demonstrating rapid development were selected for sample collection, during which leaves, petioles, immature fruits, mature fruits, stems, and roots were systematically collected from each specimen (Figure 1). Three replicate samples per tree were harvested for each tissue type. After collection, these samples were transferred into centrifuge tubes, frozen in liquid nitrogen, and stored at −80 °C.

2.2. Extraction of RNA and cDNA Synthesis

Paper Mulberry tissue samples were fully ground, after which an RNAprep Pure Polysaccharide Polyphenol Plant Total RNA Extraction Kit (TIANGEN, Nanjing, China) was used for RNA isolation. RNA integrity was analyzed via 1% agarose gel electrophoresis (AGE), and a NanoDrop™ One instrument (Thermo Fisher, Waltham, MA, USA) was used to assess RNA quality and quantity. Following this, cDNA was synthesized for each sample employing the One-Step gDNA Removal and cDNA Synthesis SuperMix kit (Trans, Nanjing, China), using a reaction volume of 20 μL and 1 μg of total RNA. RNA from individual samples was used for cDNA synthesis and subsequently stored at −20 °C.

2.3. Reference Gene Selection and Primer Design

In total, nine potential reference genes were chosen (ACT, UBQ, α-TUB1, α-TUB2, β-TUB, TIP41, CDC2, H2A and DnaJ), while BpMYB090 was selected as a representative target gene (Table 1). Full-length and mRNA sequences for these reference genes were determined using the Paper Mulberry genome annotation file, and Primer Premier 5.0 was used for primer design with the following constraints: primer length, 18–24 bp; G.C. content, 45–55%; melting temperature, ~60 °C; and product size, 100–250 bp. The specificities of the primers were evaluated by BLAST comparisons with the Paper Mulberry genome, following which the final primers were synthesized by Shenggong Biotechnology Co., Ltd (Shanghai, China). An equal amount of the cDNA samples were combined and subjected to a five-fold dilution with ddH2O. Amplified PCR products were imaged following agarose gel electrophoresis, with only primers exhibiting sufficient specificity, products of the expected size, and an absence of primer dimers being selected for further use. Then, RT-qPCR was used to test the specificity of these primers. Primers must show a single melt curve peak, free of any negative control peaks or signs of impurities, to be considered appropriate for use in future experiments. These outcomes guided the selection of final primers for additional use.

2.4. Reference Gene RT-qPCR Analyses

Following the guidelines, RT-qPCR assays were conducted using a QIAquant 96 2plex equipment (QIAGEN, Hilden, Germany) and the PowerUpTM SYBRTM Green Master Mix (Thermo Fisher, Waltham, MA, USA). Each reaction contained 2 μL of diluted cDNA templates, 0.2 μL of each of the forward and reverse primers (10 μM), 10 μL of PowerUp™ SYBR™ Green Master Mix, and ddH2O to a final 20 μL volume. Then, RT-qPCR analyses were performed as follows: 95 °C for 3 min, 40 cycles of 95 °C for 15 s, 60 °C for 15 s, and 72 °C for 30 s. A default melt curve analysis was conducted from 55–95 °C with 0.3 °C increments for 60 s, recording the corresponding fluorescent signals. Sample analyses were repeated in triplicate.

2.5. Construction of Standard Curves for Reference Gene Primers

To enable standard curve construction, a five-fold cDNA dilution series was prepared for each sample (dilutions: 1, 1/5, 1/25, 1/125, 1/625), with ddH2O serving as a negative control. As previously mentioned, RT-qPCR was used to determine the Ct values for each gene, using the QIAquant 96 2plex system. A standard curve was generated on the graph by placing log-transformed dilutions on the horizontal axis and average Ct values on the vertical axis. This allowed for the computation of the correlation coefficient (R2) and slope (K). Amplification efficiency (E) values for each of these reference genes were calculated as E = [5(1/−K) − 1] × 100%, enabling further analyses of whether the specific primers, templates, and amplification program were suitable for RT-qPCR analyses using the following standards: R2 > 0.99, 90% ≤ E ≤ 110%. Analyses were repeated in triplicate to ensure reliability.

2.6. Data Processing

geNorm evaluates the stability values (M) for reference genes, and those with the highest expression stability can be chosen for normalization [9]. The M value is used to compare reference genes that show good stability; a cut-off value of 1.5 is used, meaning that M values less than 1.5 indicate that a reference gene is suitable, and smaller M values imply higher reference gene stability. This algorithm can also determine optimal reference gene combination numbers based on measures of the Vn/n+1 value. When paired Vn/n+1 < 0.15, the optimal reference gene number is n, whereas the n + 1st gene should be introduced if Vn/n+1 > 0.15 [29].
The model-based NormFinder approach allows for the calculation of expression stability (S) values according to the 2−ΔCt method [30], allowing for appropriate reference gene selection, with smaller S values indicating greater stability [31].
The ΔCt algorithm assesses gene pairs’ relative expression levels to identify suitably stable reference genes in a given sample. Consistency between ΔCt values for two genes across different samples indicates expression stability. In contrast, fluctuating ΔCt values indicate unstable expression for one or both genes. Furthermore, 3 to 5 additional genes may be integrated into the study to provide further insights. As a result, genes are categorized to identify the most suitable reference gene candidates [32].
Gene stability can be evaluated using BestKeeper, which analyses repeated paired correlations and regression based on the SDs and CVs for paired Ct values. Smaller values in these studies correspond to more stability [33].
RefFinder is an algorithm that facilitates the thorough evaluation and ranking of findings obtained from the abovementioned four techniques, allocating suitable weights to the genes and calculating the respective geometric mean values. The genes are then allocated scores, facilitating the selection of the most optimal genes [34].

2.7. Stability Validation

BpMYB090 was selected as the target to confirm the potential genes’ stability for RT-qPCR applications. The expression of BpMYB090 in Paper Mulberry was assessed using the two most stable genes (ACT and UBQ), whereas the two genes (H2A and CDC2) showing the least stability served as controls. The relative expression of BpMYB090 was determined by the 2−ΔΔCt method [35].

3. Results

3.1. Assessment of RNA Quality and Primer Specificities

The RNA yields in samples derived from the six target Paper Mulberry tissue types (Figure 1) were examined using 1% AGE and a Thermo Fisher NanoDropTM One device. AGE analyses revealed clear RNA bands in the prepared tissue extracts, with the 28S rRNA band exhibiting a brightness level ~1.5–2.0-fold greater than that of the 18S rRNA band. These total RNA samples had OD260/280 values between 1.80 and 2.10, which is an excellent indication of high-quality RNA that can be used for further analysis. BLAST alignments with the Paper Mulberry genome were used to identify the genes (Table 1). The cDNA was then amplified using gene-specific primers. The 1% AGE results in every case showed that the target genes had been successfully amplified, as seen by the formation of clear, distinct bands with sizes between 117 and 196 bp, consistent with the length of these amplicons (Figure S1). A reference gene standard curve was also successfully generated, revealing a reasonably minimal disparity in the linear correlation coefficients (R2) between the candidate and target genes, exceeding 0.99. Primer amplification efficiency (E) ranged from 98.92% to 109.48% (Table 1) and was thus within the acceptable range (90–120%). The RT-qPCR melting curve analyses for the genes revealed single peaks devoid of primer dimer activity (Figure S2). Excellent overlap between the repeat sample amplification curves and the highly specific reference gene primer amplification was observed (Figure S3). The results demonstrated that the RT-qPCR primers developed for the nine genes were appropriate. The product specificity and efficiency values were adequate for further evaluation of reference gene stability.

3.2. Expression of Candidate Reference Genes

Ct values show an inverse correlation with the expression levels of specific genes in designated samples. Significant variations in the Ct values of the nine genes were observed among samples (Figure 2), with average values spanning from 18.00 to 28.93 (Figure 2A), falling within the Ct value range requisite for RT-qPCR analysis. The average Ct values for the individual genes were above 20, with H2A exhibiting the highest expression (average Ct: 20.79) and TIP41 exhibiting the lowest (average Ct: 25.11). Ct distributions in all test samples are shown in Figure 2B. Reduced Ct value ranges signify enhanced stability for the expression of that gene. The thorough investigations indicated that the UBQ gene demonstrated the least variation in expression among the reference gene candidates, ranging from 19.43 to 23.95. At the same time, the Ct values for α-TUB1 varied most substantially, ranging from 17.99 to 26.91. The data indicate that these candidate genes have predominantly high expression levels with minimal changes, making them appropriate as reference gene candidates. However, gene levels demonstrated irregular fluctuations across several tissues, necessitating a thorough stability evaluation employing specific approaches for reference gene identification.

3.3. Stability Analyses

The expression levels of candidate genes in various tissues were assessed employing the BestKeeper, geNorm, NormFinder, and ΔCt algorithms. The total gene stability in the samples was evaluated using RefFinder.

3.3.1. BestKeeper Analyses

The BestKeeper program facilitates direct analysis of Ct values derived from RT-qPCR, incorporating stability evaluations based on paired standard deviations (SD) and coefficients of variation (CV) for each value. Larger correlation coefficients and smaller CV ± SD values signify enhanced stability. A reference gene with an SD over 1 is deemed unstable. In these BestKeeper analyses, the SD values for α-TUB1, β-TUB, CDC2, and H2A were all found to exceed 1 because they were considered unstable reference genes. UBQ and ACT exhibited the highest level of stability, and the rank order of the genes in terms of CV and SD values was UBQ > ACT > α-TUB2 > DnaJ > TIP41 > H2A > CDC2 > β-TUB > α-TUB1 (Table 2).

3.3.2. geNorm Analyses

The geNorm program was used to evaluate the genes’ stable expression values (M). A negative association between M and stability is present when M is less than 1.5, indicating that higher M values correspond to reduced gene stability levels. Genes with an M value less than 0.5 demonstrate promise as reference genes. The M values were all less than 1.5. It is important to note that ACT, UBQ, and α-TUB1 exhibited M values of less than 0.5 in all tissue types, indicating high stability (Figure 3). On the other hand, CDC2 exhibited the least stable expression. The number of genes with M < 0.5 was only three in leaves, but it increased to 8 in stems, underscoring the extent to which gene stability depends upon tissue type (Figure 3).
The geNorm version 2003 software can be applied to find the optimal numbers of references using the values of paired variation (Vn/n+1) for different tissues when Vn/n+1 < 0.15, n represents the ideal number of references without needing n+1 references. When analyzing Paper Mulberry leaves, V2/3 = 0.16 > 0.15, while when a third gene was introduced, V3/4 = 0.13 < 0.15, demonstrating the need for accurate normalization of these genes, with α-TUB1, UBQ, and ACT showing the greatest stability in the leaves. In contrast, V2/3 = 0.12 < 0.15 in stems, consistent with the necessity of only two reference genes. TIP41 and H2A exhibit the highest stability. The optimal value of the two genes was indicated by the V2/3 values in roots, immature fruits, mature fruits, and petioles, which were all less than 0.15. These comprehensive analyses of reference genes across Paper Mulberry tissue types revealed ACT and UBQ to be the most stable, with the stability ranking being ACT = UBQ > α-TUB1 > β-TUB > TIP41 > DnaJ > α-TUB2 > H2A > CDC2 (Figure 4).

3.3.3. NormFinder Analysis

NormFinder computes stability (S) values using the 2−ΔΔCt method, thereby identifying the optimal reference gene with the highest stability, with higher S values indicating greater instability. The nine candidate genes showed variations in instability. The most stable components of the plant were the leaves, roots, and mature fruits, which ACT exhibited. UBQ was the most stable component of the plant in the petioles and stems. Overall, ACT demonstrated the greatest stability. The stability ranking was ACT > UBQ > DnaJ > β-TUB >TIP41 > H2A > α-TUB2 > α-TUB1 > CDC2 (Table 3).

3.3.4. ΔCt Analysis

Table 4 lists the nine potential genes based on their stably expressed levels in various Paper Mulberry tissues, as determined by the ΔCt technique. ACT was the most stable in leaves, early fruits, mature fruits, and all samples. DnaJ and UBQ were the most stable in stems, roots, and petioles, and β-TUB was the second most stable in the three samples. For both individual and overall sample analysis, CDC2 was the least stable. In ΔCt analyses, the ranking was ACT > UBQ > DnaJ > β-TUB > TIP41 > α-TUB2 > α-TUB1 > H2A > CDC2. Thus, ACT and UBQ performed the best as references.

3.3.5. RefFinder Analysis

Gene stability rankings varied in proportion to the differences in implementation between the four techniques used above. Therefore, the geometric means of the algorithm outputs were determined using the RefFinder tool. This showed that the best combination was that of ACT and UBQ for analyzing stem, leaf, and petiole samples, whereas ACT and DnaJ were the best combination in young and mature fruits, and ACT and α-TUB1 were the best in root samples. When assessing all available samples, ACT and UBQ exhibited the highest levels of stability (Table 5). Figure 5 shows that, in at least three algorithms, both genes exhibited notable consistency among the top five rankings. Furthermore, it was discovered that the candidates with the lowest stability were H2A and CDC2.

3.4. Verification of Reference Gene Stability

BpMYB090 was selected as a representative target to confirm the use of the suggested reference genes. The homologous gene of BpMYB090 has been shown to have functional conservation and control the formation of trichomes in various species [36,37,38,39,40]. Its expression patterns were verified in leaves at different developmental stages in wild-type (containing more trichomes) (Figure 6A–G) and mutant (with fewer trichomes) (Figure 6H–N) Paper Mulberry, including topmost, unfolded, and mature leaves. Relative BpMYB090 levels were normalized against ACT and UBQ, individually or in combination, and the results revealed consistent expression trends. BpMYB090 was highly expressed in the upper leaves early in the trichome growth process. The trichome density gradually decreases as the leaves unfold, and BpMYB090 expression declines. No expression of the BpMYB090 gene was detected in mature leaves.
BpMYB090 is involved in the development of leaf trichomes, as evidenced by the very low expression observed in the uppermost leaves of the mutant with few trichomes and the absence of expression in unfolded and mature leaves. However, when H2A and CDC2 were used for normalization, the results were inconsistent and exhibited marked deviations (Figure 6O). As a result, these findings validated the stability of ACT and UBQ genes in RT-qPCR analysis, providing a foundation for further research on Paper Mulberry genes.

4. Discussion

Paper Mulberry wood and bark are raw materials in papermaking, while the leaves represent a protein-rich feed source [41,42]. Paper Mulberry seeds, stems, leaves, roots, and fruits also contain high flavonoid concentrations and are often employed in medicinal applications [43,44,45]. It is an ideal tree species for tailings repair and ecological restoration initiatives. Inadequate crude protein presents a significant challenge in animal husbandry, whereas ecological restoration is a worldwide concern [46,47]. Paper Mulberry-related research efforts have thus shifted to molecular biology-based analyses aimed at understanding the unique physiological characteristics of this important woody plant species. Stably expressed reference genes must be identified to facilitate data standardization and ensure the validity of quantitative analyses in this setting. There are currently few reference genes available in Paper Mulberry and other woody plant species, and no reference genes that exhibit consistent expression under various conditions have been found [48,49]. Thus, it is vital for future research that systematic efforts be used to identify appropriate reference genes.
This study used RT-qPCR to assess the stable expression of nine potential genes in six different types of Paper Mulberry tissue using the BestKeeper, geNorm, NormFinder, and ΔCt algorithms. These analyses revealed slight variations in the algorithms’ assessments of stability. While ACT showed the greatest stability in all four instances, BestKeeper revealed that UBQ was relatively stably expressed. Even greater variations were observed among these four algorithms when assessing individual tissue types. For example, BestKeeper found that UBQ showed the greatest stability in leaf samples, while geNorm instead identified α-TUB1, followed by UBQ. Even with these variations, however, a large degree of overlap was observed, and these findings were in line with those from prior reports identifying reference genes in Populus euphratica [13], Salix suchowensis [50], Eucalyptus camaldulensis [51], and other plant species. GeNorm and NormFinder showed comparable rankings for these candidate genes across all samples when comparing the overall expression stability rankings for these genes. RefFinder was used to compute geometric means from the results of the four algorithms and is applied for the thorough evaluation of stability in reference genes [52,53]. Based on the abovementioned analyses, ACT and UBQ were ultimately the best combination when analyzing Paper Mulberry samples.
Software-based analyses alone cannot ensure the stable expression of reference genes, highlighting the necessity to validate this stability across various tissues by examining target genes of interest. Due to their elevated crude protein levels, paper mulberry leaves are extensively utilized as woody feed [54]. The elevated trichome density leads to poor palatability and reduced digestion and absorption when employed as animal feed, resulting in low feed consumption efficiency. Therefore, the research concentrates on the selective breeding of Paper Mulberry, which lacks trichomes. This study identified the BpMYB090 gene, which regulates the growth of leaf trichomes, as the target gene for assessing reference gene stability. It was found that, when ACT and UBQ were used for standardization, BpMYB090 levels were consistent, with high expression in the topmost leaves of the wild-type Paper Mulberry that contain high trichome densities. The trichomes matured, the density dropped, and the expression of BpMYB090 reduced as the leaves unfolded. In the highest leaves of the Paper Mulberry mutant, which has very few trichomes, the expression of BpMYB090 is minimal. Future research on the gene may find that its expression pattern, which gradually reduced as the trichome density dropped, was similar to the mutant phenotype of BpMYB090. In contrast, using H2A and CDC2, the two least stably expressed reference gene candidates, resulted in aberrant BpMYB090 expression patterns. Thus, these findings add to an existing body of evidence supporting the validity of employing the ACT and UBQ reference gene combination in analyzing various Paper Mulberry tissues.
ACT and UBQ were found in this investigation to be reference genes for various Paper Mulberry tissues. This is the first report published on reference genes for stable expression in different Paper Mulberry tissues. These data establish an evidence-based platform for future quantitative research on these functional genes and support the development of Paper Mulberry agronomic characteristics through advanced scientific and technological approaches. In addition, the data also provide a reference framework for selecting reference genes for Paper Mulberry-related species such as Broussonetia kaempferi var. australis, Broussonetia kazinoki Sieb.et Zucc., and Broussonetia monoica Hance. Although homologous reference genes typically exhibit conservation in closely related species, differences in expression stability between species are also evident. As a result, a thorough study incorporating several reference gene types is necessary for the precise assessment of target gene expression and the quantitative analysis outcomes across diverse species.

5. Conclusions

In the present study, nine putative reference genes were expressed in Paper Mulberry using RT-qPCR, and their stability was assessed using five distinct statistical algorithms. We conclude that ACT and UBQ are the most suitable RT-qPCR reference genes for various tissues of Paper Mulberry, and BpMYB090 was used as a target gene to validate the stability of the selected reference genes. To our knowledge, the current research is the first to discover suitable reference genes for normalizing gene expression analyses via RT-qPCR in several tissues of Paper Mulberry. This establishes the basis for gene expression and associated molecular biology research on Paper Mulberry.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb46100640/s1. Figure S1: PCR products shown on 1% agarose gels; Figure S2: Single peaks in melting curves for 10 reference genes; Figure S3: Amplification plots of the ten candidate reference genes.

Author Contributions

Conceptualization, S.Y.; methodology, F.Z.; software, L.X. and F.W.; validation, F.Z., and L.X.; formal analysis, F.Z. and L.X.; investigation, F.Z. and C.S.; resources, F.Z. and C.S.; data curation, F.W. and L.X.; writing–original draft, F.Z.; writing–review and editing, F.Z., S.Y.; visualization, L.X. and C.S.; supervision, S.Y.; project administration, F.Z. and S.Y.; funding acquisition, L.X. and S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Zhejiang Special Project for Provincial Research Institutes (grant number 2024F1065-4), Project of Development and Research on Unordinary Ornamental Trees (grant number, 2018HX014).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Morgan, E.C.; Overholt, W.A. Wildland Weeds: Paper Mulberry, Broussonetia papyrifera; University of Florida, IFAS: Gainesville, FL, USA, 2004; Extension: 2. [Google Scholar]
  2. Chang, C.S.; Liu, H.L.; Moncada, X.; Chung, K.F. A holistic picture of Austronesian migrations revealed by phylogeography of Pacific Paper Mulberry. Proc. Natl. Acad. Sci. USA 2015, 112, 13537–13542. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, M.; Mu, L.; Wang, Z.L.; Liu, J.Y.; Liu, T.L.; Wanapat, M.; Huang, B.Z. Assessment of mulberry leaf as a potential feed supplement for animal feeding in P.R. China. Asian-Australas. J. Anim. Sci. 2019, 32, 1145–1152. [Google Scholar] [CrossRef] [PubMed]
  4. Li, G.Y.; Hu, N.; Ding, D.X.; Zheng, J.F.; Liu, Y.L.; Wang, Y.D.; Nie, X.Q. Screening of plant species for phytoremediation of uranium, thorium, barium, nickel, strontium and lead contaminated soils from a uranium mill tailings repository in South China. Bull. Environ. Contam. Toxicol. 2011, 86, 646–652. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, J.; Liu, J.; Peng, X.; Ni, Z.; Wang, G.; Shen, S. Applied hybrid Paper Mulberry in ecological virescence of the coastal saline. Tianjin Agr. Sci. 2014, 20, 95–101. [Google Scholar]
  6. Bustin, S.; Mueller, R. Real-time reverse transcription PCR (RT-qPCR) and its potential use in clinical diagnosis. Clin. Sci. 2005, 109, 365–379. [Google Scholar] [CrossRef]
  7. Gibson, U.; Heid, C.; Williams, P. A novel method for real time quantitative RT-PCR. Genome Res. 1996, 6, 995–1001. [Google Scholar] [CrossRef]
  8. Dash, P.K.; Rai, R.; Pradhan, S.K.; Shivaraj, S.M.; Deshmukh, R.; Sreevathsa, R.; Singh, N.K. Drought and Oxidative Stress in Flax (Linum usitatissimum L.) Entails Harnessing Non-Canonical Reference Gene for Precise Quantification of qRT-PCR Gene Expression. Antioxidants 2023, 12, 950. [Google Scholar] [CrossRef]
  9. Wang, M.; Wang, Z.W.; Wei, S.S.; Xie, J.; Huang, J.Q.; Li, D.D.; Hu, W.B.; Li, H.L.; Tang, H. Identification of RT-qPCR reference genes suitable for gene function studies in the pitaya canker disease pathogen Neoscytalidium dimidiatum. Sci. Rep. 2022, 12, 22357. [Google Scholar] [CrossRef]
  10. Fiume, E.; Fletcher, J.C. Regulation of Arabidopsis embryo and endosperm development by the polypeptide signaling molecule CLE8. Plant Cell 2012, 24, 1000–1012. [Google Scholar] [CrossRef]
  11. Albuquerque, G.M.R.; Fonseca, F.C.; Boiteux, L.S.; Borges, R.C.; Miller, R.N.; Lopes, C.A.; Souza, E.B.; Fonseca, M.E.N. Stability analysis of reference genes for RT-qPCR assays involving compatible and incompatible Ralstonia solanacearum-tomato ‘Hawaii 7996’ interactions. Sci. Rep. 2021, 11, 18719. [Google Scholar] [CrossRef]
  12. Ciesielska, A.; Stączek, P. Selection and validation of reference genes for qPCR in the human dermatophyte Trichophyton rubrum exposed to different carbon sources which promote adhesion-inducing conditions. Mycoses 2021, 64, 300–308. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, H.L.; Li, L.; Tang, S.; Yuan, C.; Tian, Q.; Su, Y.; Li, H.G.; Zhao, L.; Yin, W.L.; Zhao, R.; et al. Evaluation of Appropriate Reference Genes for Reverse Transcription-Quantitative PCR Studies in Different Tissues of a Desert Poplar via Comparision of Different Algorithms. Int. J. Mol. Sci. 2015, 16, 20468–20491. [Google Scholar] [CrossRef] [PubMed]
  14. Sakuraba, Y.; Jeong, J.; Kang, M.Y.; Kim, J.; Paek, N.C.; Choi, G. Phytochrome-interacting transcription factors PIF4 and PIF5 induce leaf senescence in Arabidopsis. Nat. Commun. 2014, 5, 4636. [Google Scholar] [CrossRef] [PubMed]
  15. Gutierrez, L.; Mauriat, M.; Guénin, S.; Pelloux, J.; Lefebvre, J.F.; Louvet, R.; Rusterucci, C.; Moritz, T.; Guerineau, F.; Bellini, C.; et al. The lack of a systematic validation of reference genes: A serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnol. J. 2008, 6, 609–618. [Google Scholar] [CrossRef] [PubMed]
  16. Remans, T.; Smeets, K.; Opdenakker, K.; Mathijsen, D.; Vangronsveld, J.; Cuypers, A. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 2008, 227, 1343–1349. [Google Scholar] [CrossRef]
  17. Jain, M.; Nijhawan, A.; Tyagi, A.K.; Khurana, J.P. Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem. Biophys. Res. Commun. 2006, 345, 646–651. [Google Scholar] [CrossRef]
  18. Lv, Y.; Li, Y.Y.; Liu, X.H.; Xu, K. Identification of Ginger (Zingiber officinale Roscoe) Reference Genes for Gene Expression Analysis. Front. Genet. 2020, 11, 586098. [Google Scholar] [CrossRef]
  19. Zhou, F.W.; Xu, L.; Shi, C.G.; Yang, S.Z.; Chen, Y.H. Selection and Validation of Reliable Reference Genes for Liquidambar formosana Leaves with Different Leaf Colors. Curr. Issues Mol. Biol. 2024, 46, 9449–9462. [Google Scholar] [CrossRef]
  20. Reddy, P.S.; Reddy, D.S.; Sivasakthi, K.; Bhatnagar-Mathur, P.; Vadez, V.; Sharma, K.K. Evaluation of sorghum Sorghum bicolor (L.) reference genes in various tissues and under abiotic stress conditions for quantitative real-time PCR data normalization. Front. Plant Sci. 2016, 7, 529. [Google Scholar]
  21. Kong, Q.S.; Yuan, J.X.; Gao, L.Y.; Zhao, S.; Jiang, W.; Huang, Y.; Bie, Z.L. Identification of suitable reference genes for gene expression normalization in RT-qPCR analysis in watermelon. PLoS ONE 2014, 9, e90612. [Google Scholar]
  22. Kong, Q.S.; Yuan, J.X.; Niu, P.H.; Xie, J.J.; Jiang, W.; Huang, Y.; Bie, Z.L. Screening suitable reference genes for normalization in reverse transcription quantitative real-time PCR analysis in melon. PLoS ONE 2014, 9, e87197. [Google Scholar] [CrossRef] [PubMed]
  23. Li, J.; Han, J.G.; Hu, Y.H.; Yang, J. Selection of reference Genes for quantitative real-time PCR during flower development in Tree Peony (Paeonia suffruticosa Andr.). Front. Plant Sci. 2016, 7, 516. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, H.L.; Chen, J.H.; Tian, Q.Q.; Wang, S.; Xia, X.L.; Yin, W.L. Identification and validation of reference genes for Populus euphratica gene expression analysis during abiotic stresses by quantitative real-time PCR. Physiol. Plant. 2014, 152, 529–545. [Google Scholar] [CrossRef] [PubMed]
  25. Oliveira, L.A.; Breton, M.C.; Bastolla, F.M.; Camargo, S.D.S.; Margis, R.; Frazzon, J.; Pasquali, G. Reference genes for the normalization of gene expression in Eucalyptus species. Plant Cell Physiol. 2012, 53, 405–422. [Google Scholar] [CrossRef] [PubMed]
  26. Yi, S.Y.; Lu, H.B.; Tian, C.J.; Xu, T.; Song, C.; Wang, W.; Wei, P.P.; Gu, F.L.; Liu, D.; Cai, Y.P.; et al. Selection of Suitable Reference Genes for Gene Expression Normalization Studies in Dendrobium huoshanense. Genes 2022, 13, 1486. [Google Scholar] [CrossRef]
  27. Li, H.Y.; Chen, M.D.; Wang, Z.B.; Hao, Z.Y.; Liu, L.C.; Zhao, X.P.; Ni, J.W. Screening and validation of qRT PCR reference genes in Broussonetia papyrifera under cadmium stress. For. Sci. Res. 2023, 36, 129–138. [Google Scholar]
  28. Peng, X.J.; Liu, H.; Chen, P.L.; Tang, F.; Hu, Y.M.; Wang, F.F.; Pi, Z.; Zhao, M.L.; Chen, N.Z.; Chen, H.; et al. A Chromosome-Scale Genome Assembly of Paper Mulberry (Broussonetia papyrifera) Provides New Insights into Its Forage and Papermaking Usage. Mol. Plant 2019, 12, 661–677. [Google Scholar] [CrossRef]
  29. Hellemans, J.; Mortier, G.; De, P.A.; Speleman, F.; Vandesompele, J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8, R19. [Google Scholar] [CrossRef]
  30. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. 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 Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef]
  31. Li, L.; Li, N.; Fang, H.L.; Qi, X.W.; Zhou, Y.F. Selection and Validation of Reference Genes for Normalisation of Gene Expression in Glehnia littoralis. Sci. Rep. 2020, 10, 7374. [Google Scholar] [CrossRef]
  32. Silver, N.; Best, S.; Jiang, J.; Thein, S.L. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 2006, 7, 33. [Google Scholar] [CrossRef] [PubMed]
  33. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper-Excel-based tool using pair-wise correlations. Biotechnol. Lett. 2004, 26, 509–515. [Google Scholar] [CrossRef] [PubMed]
  34. Xie, F.L.; Xiao, P.; Chen, D.L.; Xu, L.; Zhang, B.H. miRDeepFinder: A miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 2012, 80, 75–84. [Google Scholar] [CrossRef] [PubMed]
  35. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−∆∆ct method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  36. Brockington, S.F.; Ruben, A.F.; Landis, J.B.; Alcorn, K.; Walker, R.H.; Thomas, M.M.; Hileman, L.C.; Glover, B.J. Evolutionary Analysis of the MIXTA Gene Family Highlights Potential Targets for the Study of Cellular Differentiation. Mol. Biol. Evol. 2013, 30, 526–540. [Google Scholar] [CrossRef]
  37. Zhou, F.W.; Wu, H.T.; Chen, Y.N.; Wang, M.X.; Tuskan, G.A.; Yin, T.M. Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells. Int. J. Biol. Macromol. 2023, 242 Pt 2, 124743. [Google Scholar] [CrossRef]
  38. Wu, H.T.; Tian, Y.; Wan, Q.; Fang, L.; Guan, X.Y.; Chen, J.D.; Hu, Y.; Ye, W.X.; Zhang, H.; Guo, W.Z.; et al. Genetics and evolution of MIXTA genes regulating cotton lint fiber development. New Phytol. 2018, 217, 883–895. [Google Scholar] [CrossRef]
  39. Yan, T.X.; Li, L.; Xie, L.H.; Chen, M.H.; Shen, Q.; Pan, Q.F.; Fu, X.Q.; Shi, P.; Tang, Y.L.; Huang, H.Y.; et al. A novel HD-ZIP IV/MIXTA complex promotes glandular trichome initiation and cuticle development in Artemisia annua. New Phytol. 2018, 218, 567–578. [Google Scholar] [CrossRef]
  40. Yang, S.; Cai, Y.L.; Liu, X.W.; Dong, M.M.; Zhang, Y.Q.; Chen, S.Y.; Zhang, W.B.; Li, Y.J.; Tang, M.; Zhai, X.L.; et al. CsMYB6-CsTRY module regulates fruit trichome initiation in cucumber. J. Exp. Bot. 2018, 69, 1887–1902. [Google Scholar] [CrossRef]
  41. Si, B.W.; Xu, W.C.; Guo, J.P.; Cui, K.; Wang, S.Y.; Tu, Y.; Diao, Q.Y. Effect of Broussonetia papyrifera L. (Paper Mulberry) silage on growth performance, biochemical indexes of serum and fatty acid composition in the longissimus dorsi muscle in dorper×thin-tailed han crossbred sheep. Chinese J. Anim. Vet. Sci. 2019, 50, 1424–1432. [Google Scholar]
  42. Du, Z.; Sun, L.; Chen, C.; Lin, J.; Yang, F.; Cai, Y. Exploring microbial community structure and metabolic gene clusters during silage fermentation of Paper Mulberry, a high-protein woody plant. Anim. Feed. Sci. Technol. 2020, 275, 114766. [Google Scholar] [CrossRef]
  43. Ko, H.H.; Chang, W.L.; Lu, T.M. Antityrosinase and antioxidant effects of ent-kaurane diterpenes from leaves of Broussonetia papyrifera. J. Nat. Prod. 2008, 71, 1930–1933. [Google Scholar] [CrossRef] [PubMed]
  44. Mei, R.Q.; Wang, Y.H.; Du, G.H.; Liu, G.M.; Zhang, L.; Cheng, Y.X. Antioxidant lignans from the fruits of Broussonetia papyrifera. J. Nat. Prod. 2009, 72, 621–625. [Google Scholar] [CrossRef]
  45. Lee, J.M.; Choi, S.S.; Park, M.H.; Jang, H.; Lee, Y.H.; Khim, K.W.; Oh, S.R.; Park, J.; Ryu, H.W.; Choi, J.H. Broussonetia papyrifera Root Bark Extract Exhibits Anti-inflammatory Effects on Adipose Tissue and Improves Insulin Sensitivity Potentially via AMPK Activation. Nutrients 2020, 12, 773. [Google Scholar] [CrossRef]
  46. Ye, B.; Li, J.; Xu, L.J.; Liu, H.; Yang, M.J. Metabolomic Effects of the Dietary Inclusion of Hermetia illucens Larva Meal in Tilapia. Metabolites 2022, 12, 286. [Google Scholar] [CrossRef] [PubMed]
  47. Salazar-Fillippo, A.A.; Teunkens, B.; Leirs, H.; Frouz, J.; Diggelen, R.; Miko, L. Quantitative assessment of the dispersal of soil-dwelling oribatid mites via rodents in restored heathlands. Ecol. Evol. 2022, 12, e9653. [Google Scholar] [CrossRef]
  48. Narsai, R.; Ivanova, A.; Ng, S.; Whelan, J. Defining reference genes in Oryza sativa using organ, development, biotic and abiotic transcriptome datasets. BMC Plant Biol. 2010, 10, 56. [Google Scholar] [CrossRef]
  49. Xu, Z.; Dai, J.Y.; Su, W.J.; Wu, H.X.; Shah, K.; Xing, L.B.; Ma, J.J.; Zhang, D.; Zhao, C.P. Selection and Validation of Reliable Reference Genes for Gene Expression Studies in Different Genotypes and TRV-Infected Fruits of Peach (Prunus persica L. Batsch) during Ripening. Genes 2022, 13, 160. [Google Scholar] [CrossRef]
  50. Zhou, F.W.; Chen, Y.N.; Wu, H.T.; Yin, T.M. Selection of Reliable Reference Genes for Gene Expression Analysis of Female and Male Flowers of Salix suchowensis. Plants 2022, 11, 647. [Google Scholar] [CrossRef]
  51. Unnikrishnan, B.; Shankaranarayan, D.G.; Sharma, N. Selection of reference gene in Eucalyptus camaldulensis for real-time RT-qPCR. BMC Proc. 2011, 5 (Suppl. S7), P125. [Google Scholar] [CrossRef]
  52. Ciesielska, A.; Oleksak, B.; Stączek, P. Reference genes for accurate evaluation of expression levels in Trichophyton interdigitale grown under different carbon sources, pH levels and phosphate levels. Sci. Rep. 2019, 9, 5566. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, K.; Fan, W.; Chen, D.; Jiang, L.; Li, Y.; Yao, Z.Q.; Yang, Y.; Qiu, D. Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp. Sci. Rep. 2020, 10, 22205. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, X.K.; Cao, X.X.; Liu, H.; Guo, L.N.; Lin, Y.L.; Liu, X.J.; Xiong, Y.; Ni, K.K.; Yang, F.Y. Effects of Lactic Acid Bacteria on Microbial Metabolic Functions of Paper Mulberry Silage: A BIOLOG ECO Microplates Approach. Front. Microbiol. 2021, 12, 689174. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Different Paper Mulberry tissues. (A) leaves, (B) young fruit, (C) mature fruit, (D) stems, (E) roots, and (F) petioles. Scale bar: 1 cm.
Figure 1. Different Paper Mulberry tissues. (A) leaves, (B) young fruit, (C) mature fruit, (D) stems, (E) roots, and (F) petioles. Scale bar: 1 cm.
Cimb 46 00640 g001
Figure 2. Ct values for the nine candidate reference genes in different sample types. (A) Average Ct values. (B) Distributions of Ct values in the test samples. Medians are represented by the horizontal lines, upper quartiles by the top of the boxes, and lower quartiles by the bottom of the boxes, while the upper lines represent the maximum Ct values and the outermost lines the minimum values.
Figure 2. Ct values for the nine candidate reference genes in different sample types. (A) Average Ct values. (B) Distributions of Ct values in the test samples. Medians are represented by the horizontal lines, upper quartiles by the top of the boxes, and lower quartiles by the bottom of the boxes, while the upper lines represent the maximum Ct values and the outermost lines the minimum values.
Cimb 46 00640 g002
Figure 3. The M values for the nine candidate reference genes were analyzed using the geNorm algorithm on leaves, stems, roots, young fruits, mature fruits, petioles, and overall samples.
Figure 3. The M values for the nine candidate reference genes were analyzed using the geNorm algorithm on leaves, stems, roots, young fruits, mature fruits, petioles, and overall samples.
Cimb 46 00640 g003
Figure 4. Optimal reference gene numbers, assessed by paired variation analyses (Vn/n+1). The Vn/n+1 values were determined using geNorm using a critical Vn/n+1 value of 0.15, where n represents the optimal number of references.
Figure 4. Optimal reference gene numbers, assessed by paired variation analyses (Vn/n+1). The Vn/n+1 values were determined using geNorm using a critical Vn/n+1 value of 0.15, where n represents the optimal number of references.
Cimb 46 00640 g004
Figure 5. Venn diagrams showing the five reference genes according to the BestKeeper, geNorm, NormFinder, and Delta Ct analyses. Blue, purple, green, and pink ovals, respectively, indicate the five genes with the greatest stability in samples from (A) leaves, (B) stems, (C) young fruits, (D) petioles, (E) roots, (F) mature fruits, and (G) all samples. Genes in overlapping regions of these diagrams were identified as being among the five most stable according to more than one of the utilized algorithms.
Figure 5. Venn diagrams showing the five reference genes according to the BestKeeper, geNorm, NormFinder, and Delta Ct analyses. Blue, purple, green, and pink ovals, respectively, indicate the five genes with the greatest stability in samples from (A) leaves, (B) stems, (C) young fruits, (D) petioles, (E) roots, (F) mature fruits, and (G) all samples. Genes in overlapping regions of these diagrams were identified as being among the five most stable according to more than one of the utilized algorithms.
Cimb 46 00640 g005
Figure 6. Patterns of BpMYB090 expression. (AF) The front and back of the topmost unfolded and mature leaves of wild-type Paper Mulberry; (G) Microscopic observation of trichomes in the topmost leaves of wild-type Paper Mulberry (Magnification: ×100.0); (HM) The front and back of the topmost, unfolded, and mature leaves of the mutant Paper Mulberry; (N) Microscopic observation of trichomes in the topmost leaves of the mutant Paper Mulberry (Magnification: ×100.0); (O) RT-qPCR was used to assess BpMYB090 expression in different Paper Mulberry tissues, normalizing gene expression using ACT, UBQ, or a combination of the two. Genes with the lowest stability, H2A, and CDC2, were also included as a comparison. Gene expression was analyzed via the. 2−ΔΔCt approach.
Figure 6. Patterns of BpMYB090 expression. (AF) The front and back of the topmost unfolded and mature leaves of wild-type Paper Mulberry; (G) Microscopic observation of trichomes in the topmost leaves of wild-type Paper Mulberry (Magnification: ×100.0); (HM) The front and back of the topmost, unfolded, and mature leaves of the mutant Paper Mulberry; (N) Microscopic observation of trichomes in the topmost leaves of the mutant Paper Mulberry (Magnification: ×100.0); (O) RT-qPCR was used to assess BpMYB090 expression in different Paper Mulberry tissues, normalizing gene expression using ACT, UBQ, or a combination of the two. Genes with the lowest stability, H2A, and CDC2, were also included as a comparison. Gene expression was analyzed via the. 2−ΔΔCt approach.
Cimb 46 00640 g006
Table 1. Candidate reference genes, their primer sequences, and amplification by RT-qPCR.
Table 1. Candidate reference genes, their primer sequences, and amplification by RT-qPCR.
GeneGene DescriptionGene IDPrimer Sequence F/R(5′-3′)Product Size (bp)Efficiency (%)R2
α-TUB1Alpha-tubulin1Bp02g1361.1GGTGGAGCCATACAACAGT17898.920.9951
GGCAGTGAGCGATGAGAT
α-TUB2Alpha-tubulin2Bp11g0006.1CACTGGACTACAAGGGTTT138102.450.9971
GGTGAAGGGTAGATGGTG
β-TUBβ-tubulinBp06g0215.1CGGATGATGCTCACCTTCTC117105.830.9947
CATACACTCGTCGGCGTTCT
TIP41Type 2A phosphatase activatorBp11g1240.1TGCTTATGAGACTAAGGGAC136107.770.9955
GCGGAATCAGTAGGGTAT
ACTActinBp01g3225.1AATGGTGAAGGCTGGGTT178101.280.9954
ACCGTGCTCAATGGGATA
CDC2Cyclin dependent kinase-putativeBp07g1967.1GGCATTGCTTACTGTCATTC175106.730.9913
GTGCTCTGTACCAAAGGGTC
H2AHistone H2ABp05g0572.1TGGCTGCTGAGGTGCTAGAGTT174107.610.9961
GGAGGAGGTTGTGGATGTTGG
DnaJChaperone protein DnaJ 49Bp05g0847.1AGAGGCGAACCACGAGACAT196109.480.9920
ACATTGAACCCACCAGACCC
UBQUbiquitin familyBp10g1600.1CCCTCGCCGACTACAACA189109.120.9955
TCAGCCTCTGGACCTTGC
Target gene
BpMYB090MYBBp12g0904.1TACCTGACGGCTTGGCTAC185103.210.9942
ATCCTCAATCCACCGCTCT
Table 2. BestKeeper calculations of SDs and CVs corresponding to Ct values for the specific genes.
Table 2. BestKeeper calculations of SDs and CVs corresponding to Ct values for the specific genes.
α-TUB1α-TUB2β-TUBDnaJACTTIP41H2AUBQCDC2
Leaf
SD1.511.511.120.740.640.790.960.571.74
CV6.222.535.523.23.43.584.282.725.15
Stem
SD20.820.950.7311.530.611.95
CV7.323.67.453.753.383.735.782.636.99
Young fruit
SD1.960.871.910.950.821.071.460.690.94
CV7.083.876.933.823.773.965.52.993.94
Petiole
SD2.340.771.070.750.770.831.540.61.31
CV8.473.484.533.083.573.125.712.554.54
Root
SD0.830.581.320.30.480.891.120.581.62
CV2.812.594.631.212.23.354.122.445.67
Mature fruit
SD1.320.7520.890.61.121.560.491.3
CV5.093.47.543.592.754.195.92.075.82
Total
SD2.010.731.930.80.70.981.430.641.78
CV7.263.267.053.253.213.665.392.746.27
Deeper green coloration is indicative of higher levels of stability.
Table 3. Determination of the S values of the 9 candidate genes using NormFinder.
Table 3. Determination of the S values of the 9 candidate genes using NormFinder.
Rank123456789
LeafGeneACTUBQβ-TUBDnaJTIP41H2Aα-TUB2α-TUB1CDC2
stability0.1190.1200.3170.3750.4180.4460.4930.5410.723
StemGeneUBQACTDnaJTIP41β-TUBα-TUB1α-TUB2H2ACDC2
stability0.0950.1780.2030.2260.2300.2810.3430.3470.562
Young fruitGeneDnaJUBQACTβ-TUBα-TUB2H2Aα-TUB1TIP41CDC2
stability0.0640.0640.1070.2160.3050.3330.4450.4500.776
PetioleGeneUBQDnaJACTβ-TUBα-TUB2TIP41H2Aα-TUB1CDC2
stability0.2170.2270.2390.2750.3100.3230.3500.4200.650
RootGeneACTUBQDnaJβ-TUBα-TUB2TIP41H2Aα-TUB1CDC2
stability0.2170.2270.2390.2750.3100.3230.3500.4200.650
Mature fruitGeneACTUBQDnaJα-TUB1β-TUBTIP41α-TUB2H2ACDC2
stability0.1100.1230.2600.3140.3340.3490.3640.4951.070
TotalGeneACTUBQDnaJβ-TUBTIP41H2Aα-TUB2α-TUB1CDC2
stability0.1330.2290.2870.2910.3040.3320.3370.3510.718
Table 4. ΔCt analyses of gene average STEDV for the 9 candidate genes.
Table 4. ΔCt analyses of gene average STEDV for the 9 candidate genes.
Rank123456789
LeafGeneDnaJβ-TUBUBQACTH2ATIP41α-TUB2α-TUB1CDC2
STEDV0.640.650.730.770.820.830.860.91.11
StemGeneACTβ-TUBH2AUBQTIP41α-TUB1α-TUB2DnaJCDC2
STEDV0.470.480.510.520.530.570.60.60.86
Young fruitGeneACTDnaJα-TUB2β-TUBTIP41H2AUBQα-TUB1CDC2
STEDV0.50.520.550.580.610.680.730.741.15
PetioleGeneUBQβ-TUBα-TUB2TIP41ACTα-TUB1H2ADnaJCDC2
STEDV0.590.590.610.620.640.640.710.751.01
RootGeneUBQDnaJα-TUB1ACTα-TUB2β-TUBH2ATIP41CDC2
STEDV0.550.550.570.610.630.720.750.871.07
Mature fruitGeneACTβ-TUBα-TUB1DnaJα-TUB2TIP41UBQH2ACDC2
STEDV0.60.580.650.650.650.680.680.951.56
TotalGeneACTUBQDnaJβ-TUBTIP41α-TUB2α-TUB1H2ACDC2
STEDV0.63 0.58 0.64 0.67 0.68 0.69 0.70 0.75 1.13
Table 5. Comprehensive RefFinder rankings of the stability of the candidate genes.
Table 5. Comprehensive RefFinder rankings of the stability of the candidate genes.
Method123456789
Ranking order in Leaf (Better–Good–Average)
BestKeeperUBQACTDnaJTIP41H2Aβ-TUBα-TUB1α-TUB2CDC2
geNormα-TUB1/
UBQ
ACTTIP41β-TUBDnaJα-TUB2H2ACDC2
NormFinderACTUBQβ-TUBDnaJTIP41H2Aα-TUB2α-TUB1CDC2
Delta CTDnaJβ-TUBUBQACTH2ATIP41α-TUB2α-TUB1CDC2
RefFinderUBQACTDnaJβ-TUBTIP41α-TUB1H2Aα-TUB2CDC2
Ranking order in Stem (Better–Good–Average)
BestKeeperUBQACTα-TUB2DnaJTIP41H2ACDC2β-TUBα-TUB1
geNormTIP41/
H2A
UBQACTβ-TUBα-TUB1α-TUB2DnaJCDC2
NormFinderUBQACTDnaJTIP41β-TUBα-TUB1α-TUB2H2ACDC2
Delta CTACTβ-TUBH2AUBQTIP41α-TUB1α-TUB2DnaJCDC2
RefFinderUBQACTTIP41DnaJβ-TUBH2Aα-TUB1α-TUB2CDC2
Ranking order in Young Fruit (Better–Good–Average)
BestKeeperUBQACTα-TUB2CDC2DnaJTIP41H2Aβ-TUBα-TUB1
geNormDnaJ/
TIP41
β-TUBα-TUB2ACTUBQα-TUB1H2ACDC2
NormFinderDnaJUBQACTβ-TUBα-TUB2H2Aα-TUB1TIP41CDC2
Delta CTACTDnaJα-TUB2β-TUBTIP41H2AUBQα-TUB1CDC2
RefFinderDnaJACTUBQβ-TUBα-TUB2TIP41H2Aα-TUB1CDC2
Ranking order in Petiole (Better–Good–Average)
BestKeeperUBQDnaJα-TUB2ACTTIP41β-TUBCDC2H2Aα-TUB1
geNormACT/
UBQ
α-TUB1β-TUBTIP41α-TUB2DnaJH2ACDC2
NormFinderUBQDnaJACTβ-TUBα-TUB2TIP41H2Aα-TUB1CDC2
Delta CTUBQβ-TUBα-TUB2TIP41ACTα-TUB1H2ADnaJCDC2
RefFinderUBQACTDnaJβ-TUBα-TUB2TIP41α-TUB1H2ACDC2
Ranking order in Root (Better–Good–Average)
BestKeeperα-TUB1DnaJACTα-TUB2UBQTIP41H2Aβ-TUBCDC2
geNormα-TUB1/
ACT
DnaJβ-TUBUBQα-TUB2H2ATIP41CDC2
NormFinderACTUBQDnaJβ-TUBα-TUB2TIP41H2Aα-TUB1CDC2
Delta CTUBQDnaJα-TUB1ACTα-TUB2β-TUBH2ATIP41CDC2
RefFinderα-TUB1ACTDnaJUBQβ-TUBα-TUB2H2ATIP41CDC2
Ranking order in Mature Fruit (Better–Good–Average)
BestKeeperUBQACTα-TUB2DnaJTIP41CDC2α-TUB1H2Aβ-TUB
geNormα-TUB1
/DnaJ
ACTTIP41UBQβ-TUBα-TUB2H2ACDC2
NormFinderACTUBQDnaJα-TUB1β-TUBTIP41α-TUB2H2ACDC2
Delta CTACTβ-TUBα-TUB1DnaJα-TUB2TIP41UBQH2ACDC2
RefFinderACTDnaJUBQα-TUB1TIP41β-TUBα-TUB2H2ACDC2
Ranking order in Total ( Better–Good–Average)
BestKeeperUBQACTα-TUB2DnaJTIP41H2ACDC2β-TUBα-TUB1
geNormACT/UBQ α-TUB1β-TUBTIP41DnaJα-TUB2H2ACDC2
NormFinderACTUBQDnaJβ-TUBTIP41H2Aα-TUB2α-TUB1CDC2
Delta CTACTUBQDnaJβ-TUBTIP41α-TUB2α-TUB1H2ACDC2
RefFinderACTUBQDnaJβ-TUBTIP41α-TUB2α-TUB1H2ACDC2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, F.; Xu, L.; Shi, C.; Wu, F.; Yang, S. Identification of the Optimal Quantitative RT-PCR Reference Gene for Paper Mulberry (Broussonetia papyrifera). Curr. Issues Mol. Biol. 2024, 46, 10779-10794. https://doi.org/10.3390/cimb46100640

AMA Style

Zhou F, Xu L, Shi C, Wu F, Yang S. Identification of the Optimal Quantitative RT-PCR Reference Gene for Paper Mulberry (Broussonetia papyrifera). Current Issues in Molecular Biology. 2024; 46(10):10779-10794. https://doi.org/10.3390/cimb46100640

Chicago/Turabian Style

Zhou, Fangwei, Liang Xu, Congguang Shi, Fengying Wu, and Shaozong Yang. 2024. "Identification of the Optimal Quantitative RT-PCR Reference Gene for Paper Mulberry (Broussonetia papyrifera)" Current Issues in Molecular Biology 46, no. 10: 10779-10794. https://doi.org/10.3390/cimb46100640

APA Style

Zhou, F., Xu, L., Shi, C., Wu, F., & Yang, S. (2024). Identification of the Optimal Quantitative RT-PCR Reference Gene for Paper Mulberry (Broussonetia papyrifera). Current Issues in Molecular Biology, 46(10), 10779-10794. https://doi.org/10.3390/cimb46100640

Article Metrics

Back to TopTop