Background: Blueberry (Vaccinium corymbosumL.) is an important species with a high content of flavonoids in fruits. As a perennial shrub, blueberry is characterized by shallow-rooted property and susceptible to drought stress. MYB transcription factor was reported to be widely involved in plant response to abiotic stresses, however, the role of MYB family in blueberry responding to drought stress remains elusive.Results: In this study, we conducted a comprehensive analysis of VcMYBs in blueberry based on the genome data under drought stress, including phylogenetic relationship, identification of differentially expressed genes (DEGs), expression profiling, conserved motifs, expression correlation and protein-protein interaction prediction, etc. The results showed that 229 non-redundant MYB sequences were identified in the blueberry genome, and divided into 23 subgroups. A total of 102 MYB DEGs with a significant response to drought stress were identified, of which 72 in leaves and 69 in roots, and 8 differential expression genes with a > 20-fold change in the level of expression. 17 DEGs had a higher expression correlation with other MYB members. The interaction partners of the key VcMYB proteins were predicted by STRING analysis and in combination with physiological and morphological observation. 10 keyVcMYBgenes such asVcMYB8, VcMYB102andVcMYB228were predicted to be probably involved in reactive oxygen species (ROS) pathway, and 7 keyVcMYBgenes (VcMYB41,VcMYB88andVcMYB100, etc..) probably participated in leaf regulation under drought treatment.Conclusions: Our studies provide a new understanding of the regulation mechanism ofVcMYBfamily in blueberry response to drought stress, and lay fundamental support for future studies on blueberry grown in regions with limited water supply for this crop.

Camellia (C.) oleifera Abel. is an evergreen small arbor with high economic value for producing edible oil that is well known for its high level of unsaturated fatty acids. The yield formation of tea oil extracted from fruit originates from the leaves, so leaf senescence, the final stage of leaf development, is an important agronomic trait affecting the production and quality of tea oil. However, the physiological characteristics and molecular mechanism underlying leaf senescence of C. oleifera are poorly understood. In this study, we performed physiological observation and de novo transcriptome assembly for annual leaves and biennial leaves of C. oleifera. The physiological assays showed that the content of chlorophyll (Chl), soluble protein, and antioxidant enzymes including superoxide dismutase, peroxide dismutase, and catalase in senescing leaves decreased significantly, while the proline and malondialdehyde concentration increased. By analyzing RNA-Seq data, we identified 4645 significantly differentially expressed unigenes (DEGs) in biennial leaves with most associated with flavonoid and phenylpropanoid biosynthesis and phenylalanine metabolism pathways. Among these DEGs, 77 senescence-associated genes (SAGs) including NOL, ATAF1, MDAR, and SAG12 were classified to be related to Chl degradation, plant hormone, and oxidation pathways. The further analysis of the 77 SAGs based on the Spearman correlation algorithm showed that there was a significant expression correlation between these SAGs, suggesting the potential connections between SAGs in jointly regulating leaf senescence. A total of 162 differentially expressed transcription factors (TFs) identified during leaf senescence were mostly distributed in MYB (myeloblastosis), ERF (Ethylene-responsive factor), WRKY, and NAC (NAM, ATAF1/2 and CUCU2) families. In addition, qRT-PCR analysis of 19 putative SAGs were in accordance with the RNA-Seq data, further confirming the reliability and accuracy of the RNA-Seq. Collectively, we provide the first report of the transcriptome analysis of C. oleifera leaves of two kinds of age and a basis for understanding the molecular mechanism of leaf senescence.