Anthocyanins, which are colored pigments, have long been used as food and pharmaceutical ingredients due to their potential health benefits, but the intermediate signals through which environmental or developmental cues regulate anthocyanin biosynthesis remains poorly understood. Fleshy fruits have become a good system for studying the regulation of anthocyanin biosynthesis, and exploring the mechanism underlying pigment metabolism is valuable for controlling fruit ripening. The present study revealed that ABA accumulated during Lycium fruit ripening, and this accumulation was positively correlated with the anthocyanin contents and the LbNCED1 transcript levels. The application of exogenous ABA and of the ABA biosynthesis inhibitor fluridon increased and decreased the content of anthocyanins in Lycium fruit, respectively. This is the first report to show that ABA promotes the accumulation of anthocyanins in Lycium fruits. The variations in the anthocyanin content were consistent with the variations in the expression of the genes encoding the MYB-bHLH-WD40 transcription factor complex or anthocyanin biosynthesis-related enzymes. Virus-induced LbNCED1 gene silencing significantly slowed fruit coloration and decreased both anthocyanin and ABA accumulation during Lycium fruit ripening. An qRT-PCR analysis showed that LbNCED1 gene silencing clearly reduced the transcript levels of both structural and regulatory genes in the flavonoid biosynthetic pathway. Based on the results, a model of ABA-mediated development-dependent anthocyanin biosynthesis and fruit coloration during Lycium fruit maturation was proposed. In this model, the developmental cues transcriptionally activates LbNCED1 and thus enhances accumulation of the phytohormone ABA, and the accumulated ABA stimulates transcription of the MYB-bHLH-WD40 transcription factor complex to upregulate the expression of structural genes in the flavonoid biosynthetic pathway and thereby promoting anthocyanin production and fruit coloration. Our results provide a valuable strategy that could be used in practice to regulate the ripening and quality of fresh fruit in medicinal and edible plants by modifying the phytohormone ABA.

This study aimed at assessing the climatic factors influencing the wolfberry fruit morphology, and the composition of its nutritious metabolites. The cultivar Ningqi1, widely grown in Northwest China was collected from three typical ecological growing counties with contrasting climatic conditions: Ningxia Zhongning (NF), Xinjiang Jinghe (XF) and Qinghai Nomuhong (QF). During the ripening period, 45 fruits from different plantations at each location were sampled. A total of 393 metabolites were detected in all samples through the widely targeted metabolomics approach and grouped into 19 known classes. Fruits from QF were the biggest followed by those from XF and NF. The altitude, relative humidity and light intensity had negative and strong correlations with most of the metabolites, suggesting that growing wolfberry in very high altitudes and under high light intensity is detrimental for the fruit nutritional quality. Soil moisture content is highly and negatively correlated with vitamins, organic acids and carbohydrates while moderately and positively correlated with other classes of metabolites. In contrast, air and soil temperatures exhibited positive correlation with majority of the metabolites. Overall, our results suggest high soil and air temperatures, low altitude and light intensity and moderate soil moisture, as the suitable conditions to produce Lycium fruits with high content of nutritious metabolites.

The R2R3-MYB is a large gene family involved in various plant functions, including carotenoid biosynthesis. However, this gene family lacks a comprehensive analysis in wolfberry (Lycium barbarum L.) and other Solanaceae species. The recent sequencing of the wolfberry genome provides an opportunity for investigating the organization and evolutionary characteristics of R2R3-MYB genes in wolfberry and other Solanaceae species. A total of 610 R2R3-MYB genes were identified in five Solanaceae species, including 137 in wolfberry. The LbaR2R3-MYB genes were grouped into 31 subgroups based on phylogenetic analysis, conserved gene structures, and motif composition. Five groups only of Solanaceae R2R3-MYB genes were functionally divergent during evolution. Dispersed and whole duplication events are critical for expanding the R2R3-MYB gene family. There were 287 orthologous gene pairs between wolfberry and the other four selected Solanaceae species. RNA-seq analysis identified the expression level of LbaR2R3-MYB differential gene expression (DEGs) and carotenoid biosynthesis genes (CBGs) in fruit development stages. The highly expressed LbaR2R3-MYB genes are co-expressed with CBGs during fruit development. A quantitative Real-Time (qRT)-PCR verified seven selected candidate genes. Thus, Lba11g0183 and Lba02g01219 are candidate genes regulating carotenoid biosynthesis in wolfberry. This study elucidates the evolution and function of R2R3-MYB genes in wolfberry and the four Solanaceae species.

Wolfberry (Lycium Linn. 2n = 24) fruit, Gouqizi, is a perennial shrub, traditional food and medicinal plant resource in China. Leaf and fruit related characteristics are economically important traits that are the focus for genetic improvement, but few studies into the molecular genetics of this crop have been reported to date. Here, an F₁ population (302 individuals) derived from a cross between “NO.1 Ningqi” (Lycium barbarum L.) and “Chinese gouqi” (Lycium chinese Mill.) was constructed. We recorded fruit weight, longitude, diameter and index along with leaf length, width and index for three consecutive years from 2015 to 2017. Based on this population and these phenotypic data, we constructed the first high-density genetic map of Lycium using specific length amplified fragment sequencing (SLAF-seq) and analyzed quantitative trait loci (QTLs). The map contains 6733 single nucleotide polymorphisms and 12 linkage groups (LG) with a total map distance of 1702.45 cM and an average map distance of 0.253 cM. A total of 55 QTLs were mapped for more than 2 years, of which 18 stable QTLs for fruit index on LG 11, spanning an interval of 73.492–90.945 cM, were detected. qFI11-15 for fruit index was an impressive QTL with logarithm of odds (LOD) and proportion of variance explained (PEV) values reaching 11.07 and 19.7%, respectively. The QTLs on LG 11 were gathered tightly, having an average interval of less than 1 cM per QTL, suggesting that there might be a cluster region controlling fruit index. Remarkably, qLI10-2 and qLI11-2 for leaf index were detectable for 3 years. These results give novel insight into the genetic control of leaf and fruit related traits in Lycium and provide robust support for undertaking further positional cloning studies and implementing marker-assisted selection in seedlings.