Background: RT-qPCR is a preferred method for rapid and reliable quantification of gene expression studies. Appropriate application of RT-qPCR in such studies requires the use of reference gene(s) as an internal control to normalize mRNA levels between different samples for an exact comparison of gene expression level. However, recent studies have shown that no single reference gene is universal for all experiments. Thus, the identification of high quality reference gene(s) is of paramount importance for the interpretation of data generated by RT-qPCR. Only a few studies on reference genes have been done in plants and none in peach (Prunus persica L. Batsch). Therefore, the present study was conducted to identify suitable reference gene(s) for normalization of gene expression in peach. Results: In this work, eleven reference genes were investigated in different peach samples using RT-qPCR with SYBR green. These genes are: actin 2/7 (ACT), cyclophilin (CYP2), RNA polymerase II (RP II), phospholipase A2 (PLA2), ribosomal protein L13 (RPL13), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S ribosomal RNA (18S rRNA), tubblin beta (TUB), tubblin alpha (TUA), translation elongation factor 2 (TEF2) and ubiquitin 10 (UBQ10). All eleven reference genes displayed a wide range of C_q values in all samples, indicating that they expressed variably. The stability of these genes except for RPL13 was determined by three different descriptive statistics, geNorm, NormFinder and BestKeeper, which produced highly comparable results. Conclusion: Our study demonstrates that expression stability varied greatly between genes studied in peach. Based on the results from geNorm, NormFinder and BestKeeper analyses, for all the sample pools analyzed, TEF2, UBQ10 and RP II were found to be the most suitable reference genes with a very high statistical reliability, and TEF2 and RP II for the other sample series, while 18S rRNA, RPL13 and PLA2 were unsuitable as internal controls. GAPDH and ACT also performed poorly and were less stable in our analysis. To achieve accurate comparison of levels of gene expression, two or more reference genes must be used for data normalization. The combinations of TEF2/UBQ10/RP II and TEF2/RP II were suggested for use in all samples and subsets, respectively.

Chinese jujube (Ziziphus jujuba Mill.) fruit are much admired for their unique flavor and high nutritional value. This study tracks changes in skin color and antioxidant activity over six stages (S1 … S6) of fruit development in two cultivars of jujube, ‘Junzao’ and the color mutant ‘Tailihong’. The study records the changing levels of chlorophylls, carotenoids, anthocyanins, and phenolic compounds during fruit development. Contents of chlorophylls, β-carotenes and anthocyanins decreased throughout the developmental stages in two jujube cultivars, while lutein contents decreased at first and then increased to a maximum at S6. The levels of total phenolics, total flavonoids, total flavanols, total anthocyanins, procyanidin B1, procyanidin B2, procyanidin B3, (+)-catechin, quercetin, and ferulic acid are significantly higher in ‘Tailihong’ than in ‘Junzao’ before the onset ripening (before S3). However, after S3 the level differences of these components in the two cultivars are not significant. In both cultivars, antioxidant activity reduces gradually throughout fruit development. Our findings indicate how the skin color of jujube fruit during maturation is due to changes in the levels of flavonoids, carotenoids, and anthocyanins. The color changes are also associated with changes in antioxidant activity.

The MYB transcription factor (TF) superfamily, one of the largest gene superfamilies, regulates a variety of physiological processes in plants. Although many MYB superfamily genes have been identified in plants, the MYB TFs in Chinese jujube (Ziziphus jujuba Mill.) have not been fully identified and characterized. Additionally, the functions of these genes remain unclear. In total, we identified 171 MYB superfamily genes in jujube and divided them into five subfamilies containing 99 genes of the R2R3-MYB subfamily, 58 genes of the MYB-related subfamily, four genes of the R1R2R3-MYB subfamily, one gene of the 4R-MYB subfamily, and nine genes of the atypical MYB subfamily. The 99 R2R3-MYB genes of jujube were divided into 35 groups, C1–C35, and the 58 MYB-related genes were divided into the following groups: the R-R-type, CCA1-like, I-box-binding-like, TBP-like, CPC-like, and Chinese jujube-specific groups. ZjMYB genes in jujube were well supported by additional highly conserved motifs and exon/intron structures. Most R1 repeats of MYB-related proteins comprised the R2 repeat and had highly conserved EED and EEE residue groups in jujube. Three tandem duplicated gene pairs were found on 12 chromosomes in jujube. According to an expression analysis of 126 ZjMYB genes, MYB-related genes played important roles in jujube development and fruit-related biological processes. The total flavonoid content of jujube fruit decreased as ripening progressed. A total of 93 expressed genes were identified in the RNA-sequencing data from jujube fruit, and 56 ZjMYB members presented significant correlations with total flavonoid contents by correlation analysis. Five pairs of paralogous MYB genes within jujube were composed of nine jujube MYB genes. A total of 14 ZjMYB genes had the same homology to the MYB genes of Arabidopsis and peach, indicating that these 14 MYB genes and their orthologs probably existed before the ancestral divergence of the MYB superfamily. We used a synteny analysis of MYB genes between jujube and Arabidopsis to predict that the functions of the ZjMYBs involve flavonoid/phenylpropanoid metabolism, the light signaling pathway, auxin signal transduction, and responses to various abiotic stresses (cold, drought, and salt stresses). Additionally, we speculate that ZjMYB108 is an important TF involved in the flavonoid metabolic pathway. This comprehensive analysis of MYB superfamily genes in jujube lay a solid foundation for future comprehensive analyses of ZjMYB gene functions.

Wolfberry (Lycium barbarum) is a traditional cash crop in China and is well-known worldwide for its outstanding nutritional and medicinal value. Lycium ruthenicum is a close relative of Lycium barbarum but differs significantly in size, color, flavor and nutritional composition. To date, the metabolic differences between the fruits of the two wolfberry varieties and the genetic basis behind them are unclear. Here, we compared metabolome and transcriptome data of two kinds of wolfberry fruits at five stages of development. Metabolome results show that amino acids, vitamins and flavonoids had the same accumulation pattern in various developmental stages of fruit but that Lycium ruthenicum accumulated more metabolites than Lycium barbarum during the same developmental stage, including L-glutamate, L-proline, L-serine, abscisic acid (ABA), sucrose, thiamine, naringenin and quercetin. Based on the metabolite and gene networks, many key genes that may be involved in the flavonoid synthesis pathway in wolfberry were identified, including PAL, C4H, 4CL, CHS, CHI, F3H, F3’H and FLS. The expression of these genes was significantly higher in Lycium ruthenicum than in Lycium barbarum, indicating that the difference in the expression of these genes was the main reason for the variation in flavonoid accumulation between Lycium barbarum and Lycium ruthenicum. Taken together, our results reveal the genetic basis of the difference in metabolomics between Lycium barbarum and Lycium ruthenicum and provide new insights into the flavonoid synthesis of wolfberry.

Programmed cell death (PCD) in tapetum provides nutrients for pollen development. Once this process becomes abnormal, the pollen will be aborted, and sterile lines will be formed. Lycium barbarum L. is a well-known medicinal and edible plant, and male sterile lines play an important role in the cultivation of its new varieties by utilizing hybrid vigor. However, the male sterility mechanism of L. barbarum is poorly understood. In this study, the main characteristic changes during the pollen development of L. barbarum sterile line (Ningqi No. 5) and fertile line (Ningqi No. 1) were compared through anatomical observation, physiological detection and gene expression analysis. Anatomical observation showed that compared with that of the fertile line, the tapetum of the sterile line persisted during pollen development, the organelle morphology and number of tapetal cells changed remarkably, and the PCD was remarkably delayed. Membranous peroxidation caused by reactive oxygen species (ROS) in the sterile line occurred from the late tetrad to the pollen grain stage, and that in Ningqi No. 1 transpired in the tetrad stage. This difference in the dynamic changes of ROS affected the redox balance of microspore development. qPCR detection of DYT1 and MS1 genes regulating tapetum development showed that compared with those in the fertile line, the expression levels of both genes in the sterile line changed significantly from pollen mother cell stage to pollen grain stage. This finding may be associated with the start-up delay of tapetal PCD. All these results suggested that abnormal tapetal PCD is an important mechanism leading to male sterility in L. barbarum.

Pollen fertility is a critical factor in seed development and crop breeding. Extensive studies have explored the mechanisms of pollen fertility in model plants and economic crops. However, the mechanisms of pollen abortion in medicinal and edible plants, including Lycium barbarum, remain elusive. This study utilized transcriptome analysis to identify key genes and regulatory networks implicated in pollen fertility in L. barbarum. The results demonstrated differential expression of 12,185 genes (DEGs) between the sterile and fertile lines, encompassing 489 genes that exhibited variation across the five stages of pollen development. Additionally, GO and KEGG enrichment analyses indicated that the DEGs were predominantly associated with energy metabolism, carbohydrate metabolism, and notably, hydrolase activity. Co-expression network analysis unveiled two modules intimately associated with fertility, each comprising 908 and 756 hub genes, incorporating β-1,3-glucanase genes (Glu) and co-expressed transcription factors (TFs). Phylogenetic analysis implied that LbGlu1 was a potential candidate gene implicated in regulating pollen abortion in L. barbarum. This work advances a novel understanding of pollen abortion in L. barbarum and offers theoretical support for the utilization of sterility genes to enhance crop improvement.