Selection and Validation of Reference Genes for Gene Expression Studies by RT-PCR in Dalbergia odorifera

Dietary restriction reveals sex-specific expression of the mTOR pathway genes in Japanese quails

Limited resources affect an organism's physiology through the conserved metabolic pathway, the mechanistic target of rapamycin (mTOR). Males and females often react differently to nutritional limitation, but whether it leads to differential mTOR pathway expression remains unknown. Recently, we found that dietary restriction (DR) induced significant changes in the expression of mTOR pathway genes in female Japanese quails (Coturnix japonica). We simultaneously exposed 32 male and female Japanese quails to either 20%, 30%, 40% restriction or ad libitum feeding for 14 days and determined the expression of six key genes of the mTOR pathway in the liver to investigate sex differences in the expression patterns. We found that DR significantly reduced body mass, albeit the effect was milder in males compared to females. We observed sex-specific liver gene expression. DR downregulated mTOR expression more in females than in males. Under moderate DR, ATG9A and RPS6K1 expressions were increased more in males than in females. Like females, body mass in males was correlated positively with mTOR and IGF1, but negatively with ATG9A and RS6K1 expressions. Our findings highlight that sexes may cope with nutritional deficits differently and emphasise the importance of considering sexual differences in studies of dietary restriction.

Changes in the physiological activity of parenchyma cells in Dalbergia odorifera xylem and its relationship with heartwood formation

BACKGROUND

The formation of a tree's heartwood gives the wood properties such as natural decay resistance and aesthetic color, and often directly determines the value of wood products. Regulating the quantity and quality of heartwood is of great importance to the use of wood. However, the mechanism of heartwood formation has been poorly understood.

RESULTS

Using Dalbergia odorifera as the study species, the number of starch grains, the morphology of the nuclei, the changes in the content of water and secondary metabolites were observed continuously in the radial direction of the xylem. The results show that from the outer toward inner sapwood, the starch grains are abundant, the length to diameter ratio of the nuclei is decreasing, and the morphology changes from elongated elliptical and then to round. In the outer transition zone, the starch grains begin to decrease abruptly and the nuclei shrink at a slower rate, with a radial width of approximately 2 mm. In the inner transition zone, the heartwood color begins to appear, the starch grains disappear and a few nuclei with reduced fluorescence are present, with a radial width of approximately 1 mm. Heartwood formation after complete disappearance of the nuclei. The moisture content of the heartwood is higher than that of the sapwood, and the inner transition zone is where the content rises. The secondary metabolites of the heartwood begin to accumulate in large quantities in the inner transition zone.

CONCLUSION

Based on the physiological changes of parenchyma cells in the xylem, the radial width of the transition zone of Dalbergia odorifera is clearly defined as approximately 3 mm. Both the water and secondary metabolite abrupt changes occur at the final stage of programmed cell death, and neither is a direct cause of programmed cell death in parenchyma cells.

Function of the R2R3-MYB Transcription Factors in Dalbergia odorifera and Their Relationship with Heartwood Formation

R2R3-MYB transcription factors (TFs) form one of the most important TF families involved in regulating various physiological functions in plants. The heartwood of Dalbergia odorifera is a kind of high-grade mahogany and valuable herbal medicine with wide application. However, the role of R2R3-MYB genes in the growth and development of D. odorifera, especially their relevance to heartwood formation, has not been revealed. A total of 126 R2R3-MYBs were screened from the D. odorifera genome and named DodMYB1-126 based on their location on 10 chromosomes. The collinearity results showed that purification selection was the main driving force for the evolution of the R2R3-MYB TFs family, and whole genome/fragment replication event was the main form for expanding the R2R3-MYB family, generating a divergence of gene structure and function. Comparative phylogenetic analysis classified the R2R3-MYB TFs into 33 subfamilies. S3-7,10,12-13,21 and N4-7 were extensively involved in the metabolic process; S9,13,16-19,24-25 and N1-3,8 were associated with the growth and development of D. odorifera. Based on the differential transcriptional expression levels of R2R3-MYBs in different tissues, DodMYB32, DodMYB55, and DodMYB89 were tentatively screened for involvement in the regulatory process of heartwood. Further studies have shown that the DodMYB89, localized in the nucleus, has transcriptional activation activity and is involved in regulating the biosynthesis of the secondary metabolites of heartwood by activating the promoters of the structural genes DodI2'H and DodCOMT. This study aimed to comprehensively analyze the functions of the R2R3-MYB TFs and screen for candidate genes that might be involved in heartwood formation of D. odorifera.

Selection and verification of reliable internal reference genes in stem development of herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) has emerged in the cut flower market due to its beautiful appearance. The bending flower stems caused by a lack of mechanical strength is the main problem restricting the development of the cut P. lactiflora industry. So it is of great worth to reveal the basis of stem development changes in P. lactiflora to improve its cut flower quality. Quantitative research on gene expression characteristics can provide clues for understanding their biological functions, and the screening of relatively stable expression genes is a prerequisite for the quantitative study of gene expression characteristics. Thus, it is necessary to find appropriate genes during stem development so as to analyze the qRT‒PCR results. In this study, 10 genes were screened, and these genes expressed stably in stems of different stem strengths at three different developmental stages. Then, their expressions were evaluated by RefFinder, BestKeeper, NormFinder, and GeNorm programs. The results demonstrated that γ-tubulin (γ-TUB) was the most suitable gene, followed by α-tubulin (α-TUB) and β-D-glucosidase (β-GUS), whereas histone H3 (His) was the least suitable gene. Besides, the temporal and spatial expression characteristics of PlCOMT1, the key gene concerned with the synthesis of cell wall fillers in P. lactiflora, were also used to evaluate the suitability of genes. Consequently, γ-TUB and α-TUB are the two best combinations during stem development, and their combination can be used for the stem development of P. lactiflora. These findings will provide a reference for the selection of genes related to stem development and the study of molecular mechanisms related to stem development in P. lactiflora.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01325-5.

Genome-wide identification, characterization, and genetic diversity of CCR gene family in Dalbergia odorifera

INTRODUCTION

Lignin is a complex aromatic polymer plays major biological roles in maintaining the structure of plants and in defending them against biotic and abiotic stresses. Cinnamoyl-CoA reductase (CCR) is the first enzyme in the lignin-specific biosynthetic pathway, catalyzing the conversion of hydroxycinnamoyl-CoA into hydroxy cinnamaldehyde. Dalbergia odorifera T. Chen is a rare rosewood species for furniture, crafts and medicine. However, the CCR family genes in D. odorifera have not been identified, and their function in lignin biosynthesis remain uncertain.

METHODS AND RESULTS

Here, a total of 24 genes, with their complete domains were identified. Detailed sequence characterization and multiple sequence alignment revealed that the DoCCR protein sequences were relatively conserved. They were divided into three subfamilies and were unevenly distributed on 10 chromosomes. Phylogenetic analysis showed that seven DoCCRs were grouped together with functionally characterized CCRs of dicotyledons involved in developmental lignification. Synteny analysis showed that segmental and tandem duplications were crucial in the expansion of CCR family in D. odorifera, and purifying selection emerged as the main force driving these genes evolution. Cis-acting elements in the putative promoter regions of DoCCRs were mainly associated with stress, light, hormones, and growth/development. Further, analysis of expression profiles from the RNA-seq data showed distinct expression patterns of DoCCRs among different tissues and organs, as well as in response to stem wounding. Additionally, 74 simple sequence repeats (SSRs) were identified within 19 DoCCRs, located in the intron or untranslated regions (UTRs), and mononucleotide predominated. A pair of primers with high polymorphism and good interspecific generality was successfully developed from these SSRs, and 7 alleles were amplified in 105 wild D. odorifera trees from 17 areas covering its whole native distribution.

DISCUSSION

Overall, this study provides a basis for further functional dissection of CCR gene families, as well as breeding improvement for wood properties and stress resistance in D. odorifera.

Iso-Seq and RNA-Seq analyses uncover the molecular response of Dalbergia odorifera T. Chen to low temperature

In previous studies, we found young plants propagated from acclimated D. odorifera T. Chen trees displayed increased cold tolerance. To understand the molecular and physiological mechanisms underlying this phenomenon, we prepared acclimated and non-acclimated young D. odorifera T. Chen trees for physiological, RNA-Seq, and Iso-Seq analyses. First, both the acclimated and non-acclimated young trees were grown in 10 ℃ condition, then treated at -3 ℃ for one day, and finally recovered at 10 ℃. The measurement of physiological parameters, including superoxide dismutase, peroxidase, malondialdehyde, thiobarbituric acid, and soluble sugar, showed that the physiological change of acclimated plants is smoother than non-acclimated plants. The RNA-Seq analyses pointed out that cold acclimation already fixed the different gene expression patterns of D. odorifera T. Chen trees. The hormone-related, secondary metabolic, and signal transduction related biological genes tend to show different expressions between the acclimated and non-acclimated D. odorifera T. Chen trees. Moreover, the change of gene expression for some biological processes, such as alpha-Linolenic acid metabolism and its response to hydrogen peroxide, seems to occur earlier in non-acclimated than acclimated plants. The ISO-Seq analyses pointed out that alternative splicing (AS) of some genes was also found, and these AS events were predicted to play important roles in regulating different expression patterns between non-acclimated and acclimated plants. Therefore, according to the performance of this study, we are able to provide some novel understanding of cold tolerance enhanced by cold acclimation in perennial trees.

Integrated mRNA and Small RNA Sequencing Reveals microRNAs Associated With Xylem Development in Dalbergia odorifera

Dalbergia odorifera is a rare and precious rosewood specie, whose wood is a very high-quality material for valuable furniture and carving crafts. However, limited information is available about the process of wood formation in D. odorifera. To determine genes that might be closely associated with the xylem differentiation process, we analyzed the differentially expressed genes (DEGs) and microRNAs (miRNAs) from specific xylem tissues of D. odorifera by RNA sequencing (RNA-seq) and small RNA sequencing (small RNA-seq). In total, we obtained 134,221,955 clean reads from RNA-seq and 90,940,761 clean reads from small RNA-seq. By comparing the transition zone (Dotz) and sapwood (Dosw) samples, a total of 395 DEGs were identified. Further analysis revealed that DEGs encoded for WRKY transcription factors (eight genes), lignin synthesis (PER47, COMT, CCR2), cell wall composition (UXS2), gibberellin synthesis (KAO2, GA20OX1), jasmonic acid synthesis (OPR2, CYP74A), and synthesis of flavonoids (PAL2) and terpenoids (CYP71A1). Subsequently, a preliminary analysis by small RNA-seq showed that the expressions of 14 miRNAs (such as miR168a-5p, miR167f-5p, miR167h-5p, miR167e, miR390a, miR156g, novel_52, and novel_9) were significantly different between Dotz and Dosw. Further analysis revealed that the target genes of these differentially expressed miRNAs were enriched in the GO terms "amino acid binding," "cellulase activity," and "DNA beta-glucosyltransferase activity". Further, KEGG pathway annotation showed significant enrichment in "fatty acid elongation" and "biosynthesis of unsaturated fatty acids". These processes might be participating in the xylem differentiation of D. odorifera. Next, expression correlation analysis showed that nine differentially expressed miRNAs were significantly negatively associated with 21 target genes, which encoded for proteins such as pyrH, SPL6, SPL12, GCS1, and ARF8. Overall, this is the first study on miRNAs and their potential functions in the xylem development of D. odorifera, which provides a stepping stone for a detailed functional investigation of D. odorifera miRNAs.

Establishment of reference (housekeeping) genes via quantitative real-time PCR for investigation of the genomic basis of abiotic stress resistance in Psammochloa villosa (Poaceae)

Psammochloa villosa is a desert plant growing in Northwest China with considerable resistance to abiotic stress, including drought, cold, and salt. To facilitate future studies of stress resistance in Psammochloa villosa, we sought to establish a suite of reference (or housekeeping) genes for utilization within future gene expression studies. Specifically, we selected nine candidate genes based on prior studies and new transcriptomic data for P. villosa, and we evaluated their expression stability in three different tissues of P. villosa under different treatments simulating abiotic stress conditions using four different bioinformatics assessments. Our results showed that TIP41 (TIP41-like family protein) was the most stable reference gene in drought- and salt-stressed leaves and salt-stressed stems, ELF-1α (elongation factor 1-α) was the most stable in cold-stressed leaves and drought- and salt-stressed roots, ACT (actin) was the most stable in drought-stressed stems, TUA (α-tubulin) was the most stable in cold-stressed stems, and 18S rRNA (18S ribosomal RNA) was the most stable in cold-stressed roots. Additionally, we tested the utility of these candidate reference genes to detect the expression pattern of P5CS (Δ¹-pyrroline-5-carboxylate synthetase), which is a drought-related gene. This study is the first report on selecting and validating reference genes of P. villosa under various stress conditions and will benefit future investigations of the genomic mechanisms of stress resistance in this ecologically important species.

Selection and optimization of reference genes for RT-qPCR normalization: A case study in Solanum lycopersicum exposed to UV-B

Quantitative RT- PCR is one of the most common methods to study gene expression in response to stress. Therefore, it is crucial to have suitable reference genes (RGs) for result normalization. Although several reports describe UV-B-modulated gene expression in Solanum lycopersicum, there are no suitable RGs identified until now. The aim of this work was to evaluate the suitability of seven traditional genes: actin (ACT), tubulin (TUB), ubiquitin (UBI), glyceraldehyde- 3 phosphate dehydrogenase (GAPDH), elongation factor 1α (EF1α), phosphatase 2A catalytic subunit (PP2A) and GAGA binding transcriptional activator (GAGA); and two non-traditional genes: thioredoxin h1 (TRX h1) and UV-B RESISTANCE LOCUS 8 (UVR8), as candidate RGs for their potential use as reliable internal controls in leaves, stems and roots of tomato seedlings exposed to acute and chronic UV-B. The stability of these genes expression was evaluated using five statistical algorithms: geNorm, NormFinder, BestKeeper, Delta Ct and ANOVA. Considering the comprehensive stability ranking, we recommend ACT+TUB as the best pair of RGs for leaves, PP2A+GAPDH+TRX h1 for stems and TUB+UVR8 for roots. The reliability of the selected RGs for each tissue was verified amplifying tomato chalcone synthase 1 (CHS1) and cyclobutane pyrimidine dimer (CPD) photolyase (PHR1-LIKE). Under UV-B treatment, CHS1 was upregulated in leaves, stems and roots whereas PHR1-LIKE was only upregulated in leaves and stems. This interpretation differs when the most and least stable RGs are chosen. This is the first report regarding suitable RGs selection for accurate normalization of gene expression in tomato seedlings exposed to UV-B irradiation.

Molecular Mechanism Underlying Mechanical Wounding-Induced Flavonoid Accumulation in Dalbergia odorifera T. Chen, an Endangered Tree That Produces Chinese Rosewood

Dalbergia odorifera, a critically endangered tree species, produces heartwood containing a vast variety of flavonoids. This heartwood, also known as Chinese rosewood, has high economic and medicinal value, but its formation takes several decades. In this study, we showed that discolored wood induced by pruning displays similar color, structure, and flavonoids content to those of natural heartwood, suggesting that wounding is an efficient method for inducing flavonoid production in D. odorifera. Transcriptome analysis was performed to investigate the mechanism underlying wounding-induced flavonoids production in D. odorifera heartwood. Wounding upregulated the expression of 90 unigenes, which covered 19 gene families of the phenylpropanoid and flavonoid pathways, including PAL, C4H, 4CL, CHS, CHI, 6DCS, F3’5’H, F3H, FMO, GT, PMAT, CHOMT, IFS, HI4’OMT, HID, IOMT, I2’H, IFR, and I3’H. Furthermore, 47 upregulated unigenes were mapped to the biosynthesis pathways for five signal molecules (ET, JA, ABA, ROS, and SA). Exogenous application of these signal molecules resulted in the accumulation of flavonoids in cell suspensions of D. odorifera, supporting their role in wounding-induced flavonoid production. Insights from this study will help develop new methods for rapidly inducing the formation of heartwood with enhanced medicinal value.