Wolfberry genomes and the evolution of Lycium (Solanaceae)

Metabolite-based genome-wide association studies enable the dissection of the genetic bases of flavonoids, betaine and spermidine in wolfberry (Lycium)

Wolfberry is a plant with medicinal and food values. However, its bioactive ingredients and the corresponding genetic bases have not been determined. Here, we de novo generated a chromosome-level genome assembly for wolfberry, yielding a genome sequence of ~1.77 Gb with contig N50 of 50.55 Mb and 39 224 predicted gene models. A variation map, using 307 re-sequenced accessions, was called based on this genome assembly. Furthermore, the fruit metabolome of these accessions was profiled using 563 annotated metabolites, which separated Lycium barbarum L. and non-L. barbarum L. The flavonoids, coumarins, alkaloids and nicotinic acid contents were higher in the former than in the latter. A metabolite-based genome-wide association study mapped 156 164 significant single nucleotide polymorphisms corresponding to 340 metabolites. This included 19 219 unique lead single nucleotide polymorphisms in 1517 significant association loci, of which three metabolites, flavonoids, betaine and spermidine, were highlighted. Two candidate genes, LbUGT (evm.TU.chr07.2692) and LbCHS (evm.TU.chr07.2738), with non-synonymous mutations, were associated with the flavonoids content. LbCHS is a structural gene that interacts with a nearby MYB transcription factor (evm.TU.chr07.2726) both in L. barbarum and L. ruthenicum. Thus, these three genes might be involved in the biosynthesis/metabolism of flavonoids. LbSSADH (evm.TU.chr09.627) was identified as possibly participating in betaine biosynthesis/metabolism. Four lycibarbarspermidines (E-G and O) were identified, and only the lycibarbarspermidines O content was higher in L. barbarum varieties than in non-L. barbarum varieties. The evm.TU.chr07.2680 gene associated with lycibarbarspermidines O was annotated as an acetyl-CoA-benzylalcohol acetyltransferase, suggesting that it is a candidate gene for spermidine biosynthesis. These results provide novel insights into the specific metabolite profile of non-L. barbarum L. and the genetic bases of flavonoids, betaine and spermidine biosynthesis/metabolism.

Mapping quantitative trait loci associated with self-(in)compatibility in goji berries (Lycium barbarum)

BACKGROUND

Goji (Lycium barbarum L.) is a perennial deciduous shrub widely distributed in arid and semiarid regions of Northwest China. It is highly valued for its medicinal and functional properties. Most goji varieties are naturally self-incompatible, posing challenges in breeding and cultivation. Self-incompatibility is a complex genetic trait, with ongoing debates regarding the number of self-incompatible loci. To date, no genetic mappings has been conducted for S loci or other loci related to self-incompatibility in goji.

RESULTS

We used genome resequencing to create a high-resolution map for detecting de novo single-nucleotide polymorphisms (SNP) in goji. We focused on 229 F1 individuals from self-compatible '13-19' and self-incompatible 'new 9' varieties. Subsequently, we conducted a quantitative trait locus (QTL) analysis on traits associated with self-compatibility in goji berries. The genetic map consisted of 249,327 SNPs distributed across 12 linkage groups (LGs), spanning a total distance of 1243.74 cM, with an average interval of 0.002 cM. Phenotypic data related to self-incompatibility, such as average fruit weight, fruit rate, compatibility index, and comparable compatibility index after self-pollination and geitonogamy, were collected for the years 2021-2022, as well as for an extra year representing the mean data from 2021 to 2022 (2021/22). A total of 43 significant QTL, corresponding to multiple traits were identified, accounting for more than 11% of the observed phenotypic variation. Notably, a specific QTL on chromosome 2 consistently appeared across different years, irrespective of the relationship between self-pollination and geitonogamy. Within the localization interval, 1180 genes were annotated, including Lba02g01102 (annotated as an S-RNase gene), which showed pistil-specific expression. Cloning of S-RNase genes revealed that the parents had two different S-RNase alleles, namely S1S11 and S2S8. S-genotype identification of the F1 population indicated segregation of the four S-alleles from the parents in the offspring, with the type of S-RNase gene significantly associated with self-compatibility.

CONCLUSIONS

In summary, our study provides valuable insights into the genetic mechanism underlying self-compatibility in goji berries. This highlights the importance of further positional cloning investigations and emphasizes the importance of integration of marker-assisted selection in goji breeding programs.

Comparative metabolomics combined with genome sequencing provides insights into novel wolfberry-specific metabolites and their formation mechanisms

Wolfberry (Lycium, of the family Solanaceae) has special nutritional benefits due to its valuable metabolites. Here, 16 wolfberry-specific metabolites were identified by comparing the metabolome of wolfberry with those of six species, including maize, rice, wheat, soybean, tomato and grape. The copy numbers of the riboflavin and phenyllactate degradation genes riboflavin kinase (RFK) and phenyllactate UDP-glycosyltransferase (UGT1) were lower in wolfberry than in other species, while the copy number of the phenyllactate synthesis gene hydroxyphenyl-pyruvate reductase (HPPR) was higher in wolfberry, suggesting that the copy number variation of these genes among species may be the main reason for the specific accumulation of riboflavin and phenyllactate in wolfberry. Moreover, the metabolome-based neighbor-joining tree revealed distinct clustering of monocots and dicots, suggesting that metabolites could reflect the evolutionary relationship among those species. Taken together, we identified 16 specific metabolites in wolfberry and provided new insight into the accumulation mechanism of species-specific metabolites at the genomic level.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum)

The characterization of the PYL/RCAR ABA receptors in a great deal of plant species has dramatically advanced the study of ABA functions involved in key physiological processes. However, the genes in this family are still unclear in Lycium (Goji) plants, one of the well-known economically, medicinally, and ecologically valuable fruit crops. In the present work, 12 homologs of Arabidopsis PYL/RCAR ABA receptors were first identified and characterized from Lycium (L.) barbarum (LbPYLs). The quantitative real-time PCR (qRT-PCR) analysis showed that these genes had clear tissue-specific expression patterns, and most of them were transcribed in the root with the largest amount. Among the three subfamilies, while the Group I and Group III members were down-regulated by extraneous ABA, the Group II members were up-regulated. At 42 °C, most transcripts showed a rapid and violent up-regulation response to higher temperature, especially members of Group II. One of the genes in the Group II members, LbPYL10, was further functionally validated by virus-induced gene silencing (VIGS) technology. LbPYL10 positively regulates heat stress tolerance in L. barbarum by alleviating chlorophyll degradation, thus maintaining chlorophyll stability. Integrating the endogenous ABA level increase following heat stress, it may be concluded that LbPYL-mediated ABA signaling plays a vital role in the thermotolerance of L. barbarum plants. Our results highlight the strong potential of LbPYL genes in breeding genetically modified L. barbarum crops that acclimate to climate change.

Wolfberry genome database: integrated genomic datasets for studying molecular biology

Wolfberry, also known as goji berry or Lycium barbarum, is a highly valued fruit with significant health benefits and nutritional value. For more efficient and comprehensive usage of published L. barbarum genomic data, we established the Wolfberry database. The utility of the Wolfberry Genome Database (WGDB) is highlighted through the Genome browser, which enables the user to explore the L. barbarum genome, browse specific chromosomes, and access gene sequences. Gene annotation features provide comprehensive information about gene functions, locations, expression profiles, pathway involvement, protein domains, and regulatory transcription factors. The transcriptome feature allows the user to explore gene expression patterns using transcripts per kilobase million (TPM) and fragments per kilobase per million mapped reads (FPKM) metrics. The Metabolism pathway page provides insights into metabolic pathways and the involvement of the selected genes. In addition to the database content, we also introduce six analysis tools developed for the WGDB. These tools offer functionalities for gene function prediction, nucleotide and amino acid BLAST analysis, protein domain analysis, GO annotation, and gene expression pattern analysis. The WGDB is freely accessible at https://cosbi7.ee.ncku.edu.tw/Wolfberry/. Overall, WGDB serves as a valuable resource for researchers interested in the genomics and transcriptomics of L. barbarum. Its user-friendly web interface and comprehensive data facilitate the exploration of gene functions, regulatory mechanisms, and metabolic pathways, ultimately contributing to a deeper understanding of wolfberry and its potential applications in agronomy and nutrition.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Genomic Analysis of Romanian Lycium Genotypes: Exploring BODYGUARD Genes for Stress Resistance Breeding

Goji berries, long valued in Traditional Chinese Medicine and Asian cuisine for their wide range of medicinal benefits, are now considered a 'superfruit' and functional food worldwide. Because of growing demand, Europe and North America are increasing their goji berry production, using goji berry varieties that are not originally from these regions. European breeding programs are focusing on producing Lycium varieties adapted to local conditions and market demands. By 2023, seven varieties of goji berries were successfully registered in Romania, developed using germplasm that originated from sources outside the country. A broader project focused on goji berry breeding was initiated in 2014 at USAMV Bucharest. In the present research, five cultivated and three wild L. barbarum genotypes were compared to analyse genetic variation at the whole genome level. In addition, a case study presents the differences in the genomic coding sequences of BODYGUARD (BDG) 3 and 4 genes from chromosomes 4, 8, and 9, which are involved in cuticle-related resistance. All three BDG genes show distinctive differences between the cultivated and wild-type genotypes at the SNP level. In the BDG 4 gene located on chromosome 8, 69% of SNPs differentiate the wild from the cultivated genotypes, while in BDG 3 on chromosome 4, 64% of SNPs could tell the difference between the wild and cultivated goji berry. The research also uncovered significant SNP and InDel differences between cultivated and wild genotypes, in the entire genome, providing crucial insights for goji berry breeders to support the development of goji berry cultivation in Romania.

Integrated metabolome and transcriptome analysis of differences in quality of ripe Lycium barbarum L. fruits harvested at different periods

BACKGROUND

Wolfberry is well-known for its high nutritional value and medicinal benefits. Due to the continuous ripening nature of Goji berries and the fact that they can be commercially harvested within a few weeks, their phytochemical composition may change during the harvesting process at different periods.

RESULT

The involved molecular mechanisms of difference in fruit quality of ripe Lycium barbarum L. harvested at four different periods were investigated by transcriptomic and metabolomics analyses for the first time. According to the results we obtained, it was found that the appearance quality of L. barbarum fruits picked at the beginning of the harvesting season was superior, while the accumulation of sugar substances in L. barbarum fruits picked at the end of the harvesting season was better. At the same time the vitamin C and carotenoids content of wolfberry fruits picked during the summer harvesting season were richer. Ascorbic acid, succinic acid, glutamic acid, and phenolic acids have significant changes in transcription and metabolism levels. Through the network metabolic map, we found that ascorbic acid, glutamic acid, glutamine and related enzyme genes were differentially accumulated and expressed in wolfberry fruits at different harvesting periods. Nevertheless, these metabolites played important roles in the ascorbate-glutathione recycling system. Ascorbic acid, phenolic substances and the ascorbate-glutathione recycling system have antioxidant effects, which makes the L. barbarum fruits harvested in the summer more in line with market demand and health care concepts.

CONCLUSION

This study laid the foundation for understanding the molecular regulatory mechanisms of quality differences of ripe wolfberry fruits harvested at different periods, and provides a theoretical basis for enhancing the quality of L. barbarum fruits.

Effect of long-term deficit irrigation on tomato and goji berry quality: from fruit composition to in vitro bioaccessibility of carotenoids

Drought is a persistent challenge for horticulture, affecting various aspects of fruit development and ultimately fruit quality, but the effect on nutritional value has been under-investigated. Here, fruit quality was studied on six tomato genotypes and one goji cultivar under deficit irrigation (DI), from fruit composition to in vitro bioaccessibility of carotenoids. For both species, DI concentrated most health-related metabolites in fresh fruit. On a dry mass basis, DI increased total phenolic and sugar concentration, but had a negative or insignificant impact on fruit ascorbic acid, organic acid, and alcohol-insoluble matter contents. DI also reduced total carotenoids content in tomato (-18.7% on average), especially β-carotene (-32%), but not in goji berry DW (+15.5% and +19.6%, respectively). DI reduced the overall in vitro bioaccessibility of carotenoids to varying degrees depending on the compound and plant species. Consequently, mixed micelles produced by digestion of fruits subjected to DI contained either the same or lesser quantities of carotenoids, even though fresh fruits could contain similar or higher quantities. Thus, DI effects on fruit composition were species and genotype dependent, but an increase in the metabolite concentration did not necessarily translate into greater bioaccessibility potentially due to interactions with the fruit matrix.

Multiple independent losses of the biosynthetic pathway for two tropane alkaloids in the Solanaceae family

Hyoscyamine and scopolamine (HS), two valuable tropane alkaloids of significant medicinal importance, are found in multiple distantly related lineages within the Solanaceae family. Here we sequence the genomes of three representative species that produce HS from these lineages, and one species that does not produce HS. Our analysis reveals a shared biosynthetic pathway responsible for HS production in the three HS-producing species. We observe a high level of gene collinearity related to HS synthesis across the family in both types of species. By introducing gain-of-function and loss-of-function mutations at key sites, we confirm the reduced/lost or re-activated functions of critical genes involved in HS synthesis in both types of species, respectively. These findings indicate independent and repeated losses of the HS biosynthesis pathway since its origin in the ancestral lineage. Our results hold promise for potential future applications in the artificial engineering of HS biosynthesis in Solanaceae crops.

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (Solanum lycopersicum)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato.

Approaches to increase the validity of gene family identification using manual homology search tools

Identifying homologs is an important process in the analysis of genetic patterns underlying traits and evolutionary relationships among species. Analysis of gene families is often used to form and support hypotheses on genetic patterns such as gene presence, absence, or functional divergence which underlie traits examined in functional studies. These analyses often require precise identification of all members in a targeted gene family. Manual pipelines where homology search and orthology assignment tools are used separately are the most common approach for identifying small gene families where accurate identification of all members is important. The ability to curate sequences between steps in manual pipelines allows for simple and precise identification of all possible gene family members. However, the validity of such manual pipeline analyses is often decreased by inappropriate approaches to homology searches including too relaxed or stringent statistical thresholds, inappropriate query sequences, homology classification based on sequence similarity alone, and low-quality proteome or genome sequences. In this article, we propose several approaches to mitigate these issues and allow for precise identification of gene family members and support for hypotheses linking genetic patterns to functional traits.

Transcriptional deciphering of the metabolic pathways associated with the bioactive ingredients of wolfberry species with different quality characteristics

BACKGROUND

Wolfberry is rich in carotenoids, flavonoids, vitamins, alkaloids, betaines and other bioactive ingredients. For over 2,000 years, wolfberry has been used in China as a medicinal and edible plant resource. Nevertheless, the content of bioactive ingredients varies by cultivars, resulting in uneven quality across wolfberry cultivars and species. To date, research has revealed little about the underlying molecular mechanism of the metabolism of flavonoids, carotenoids, and other bioactive ingredients in wolfberry.

RESULTS

In this context, the transcriptomes of the Lycium barbarum L. cultivar 'Ningqi No. 1' and Lycium chinense Miller were compared during the fruit maturity stage using the Illumina NovaSeq 6000 sequencing platform, and subsequently, the changes of the gene expression profiles in two types of wolfberries were analysed. In total, 256,228,924 clean reads were obtained, and 8817 differentially expressed genes (DEGs) were identified, then assembled by Basic Local Alignment Search Tool (BLAST) similarity searches and annotated using Gene Ontology (GO), Clusters of Orthologous Groups of proteins (KOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). By combining these transcriptome data with data from the PubMed database, 36 DEGs related to the metabolism of bioactive ingredients and implicated in the metabolic pathway of carotenoids, flavonoids, terpenoids, alkaloids, vitamins, etc., were identified. In addition, among the 9 differentially expressed transcription factors, LbAPL, LbPHL11 and LbKAN4 have raised concerns. The protein physicochemical properties, structure prediction and phylogenetic analysis indicated that LbAPL and LbPHL11 may be good candidate genes involved in regulating the flavonoid metabolism pathway in wolfberry.

CONCLUSIONS

This study provides preliminary evidence for the differences in bioactive ingredient content at the transcription level among different wolfberry species, as well as a research and theoretical basis for the screening, cloning and functional analysis of key genes involved in the metabolism of bioactive ingredients in wolfberry.

Comparative Analysis of the GATA Transcription Factors in Five Solanaceae Species and Their Responses to Salt Stress in Wolfberry (Lycium barbarum L.)

GATA proteins are a class of zinc-finger DNA-binding proteins that participate in diverse regulatory processes in plants, including the development processes and responses to environmental stresses. However, a comprehensive analysis of the GATA gene family has not been performed in a wolfberry (Lycium barbarum L.) or other Solanaceae species. There are 156 GATA genes identified in five Solanaceae species (Lycium barbarum L., Solanum lycopersicum L., Capsicum annuum L., Solanum tuberosum L., and Solanum melongena L.) in this study. Based on their phylogeny, they can be categorized into four subfamilies (I-IV). Noticeably, synteny analysis revealed that dispersed- and whole-genome duplication contributed to the expansion of the GATA gene family. Purifying selection was a major force driving the evolution of GATA genes. Moreover, the predicted cis-elements revealed the potential roles of wolfberry GATA genes in phytohormone, development, and stress responses. Furthermore, the RNA-seq analysis identified 31 LbaGATA genes with different transcript profiling under salt stress. Nine candidate genes were then selected for further verification using quantitative real-time PCR. The results revealed that four candidate LbaGATA genes (LbaGATA8, LbaGATA19, LbaGATA20, and LbaGATA24) are potentially involved in salt-stress responses. In conclusion, this study contributes significantly to our understanding of the evolution and function of GATA genes among the Solanaceae species, including wolfberry.

Complete chloroplast genome assembly and phylogenetic analysis of blackcurrant (Ribes nigrum), red and white currant (Ribes rubrum), and gooseberry (Ribes uva-crispa) provide new insights into the phylogeny of Grossulariaceae

Background

Blackcurrant (Ribes nigrum), red currant (R. rubrum), white currant (R. rubrum), and gooseberry (R. uva-crispa) belong to Grossulariaceae and are popular small-berry crops worldwide. The lack of genomic data has severely limited their systematic classification and molecular breeding.

Methods

The complete chloroplast (cp) genomes of these four taxa were assembled for the first time using MGI-DNBSEQ reads, and their genome structures, repeat elements and protein-coding genes were annotated. By genomic comparison of the present four and previous released five Ribes cp genomes, the genomic variations were identified. By phylogenetic analysis based on maximum-likelihood and Bayesian methods, the phylogeny of Grossulariaceae and the infrageneric relationships of the Ribes were revealed.

Results

The four cp genomes have lengths ranging from 157,450 to 157,802 bp and 131 shared genes. A total of 3,322 SNPs and 485 Indels were identified from the nine released Ribes cp genomes. Red currant and white currant have 100% identical cp genomes partially supporting the hypothesis that white currant (R. rubrum) is a fruit color variant of red currant (R. rubrum). The most polymorphic genic and intergenic region is ycf1 and trnT-psbD, respectively. The phylogenetic analysis demonstrated the monophyly of Grossulariaceae in Saxifragales and the paraphyletic relationship between Saxifragaceae and Grossulariaceae. Notably, the Grossularia subgenus is well nested within the Ribes subgenus and shows a paraphyletic relationship with the co-ancestor of Calobotrya and Coreosma sections, which challenges the dichotomous subclassification of the Ribes genus based on morphology (subgenus Ribes and subgenus Grossularia). These data, results, and insights lay a foundation for the phylogenetic research and breeding of Ribes species.

Changes in m6A RNA methylation are associated with male sterility in wolfberry

BACKGROUND

N6-methyladenosine (m6A) modification is the most abundant type of RNA modification in eukaryotic cells, playing pivotal roles in multiple plant growth and development processes. Yet the potential role of m6A in conferring the trait of male sterility in plants remains unknown.

RESULTS

In this study, we performed RNA-sequencing (RNA-Seq) and m6A-sequencing (m6A-Seq) of RNAs obtained from the anther tissue of two wolfberry lines: 'Ningqi No.1' (LB1) and its natural male sterile mutant 'Ningqi No.5' (LB5). Based on the newly assembled transcriptome, we established transcriptome-wide m6A maps for LB1 and LB5 at the single nucleus pollen stage. We found that the gene XLOC_021201, a homolog of m6A eraser-related gene ALKBH10 in Arabidopsis thaliana, was significantly differentially expressed between LB1 and LB5. We also identified 1642 and 563 m6A-modified genes with hypermethylated and hypomethylated patterns, respectively, in LB1 compared with LB5. We found the hypermethylated genes significantly enriched in biological processes related to energy metabolism and lipid metabolism, while hypomethylation genes were mainly linked to cell cycle process, gametophyte development, and reproductive process. Among these 2205 differentially m6A methylated genes, 13.74% (303 of 2205) were differentially expressed in LB1 vis-à-vis LB5.

CONCLUSIONS

This study constructs the first m6A transcriptome map of wolfberry and establishes an association between m6A and the trait of male sterility in wolfberry.

Genome-Wide Analysis of MBF1 Family Genes in Five Solanaceous Plants and Functional Analysis of SlER24 in Salt Stress

Multiprotein bridging factor 1 (MBF1) is an ancient family of transcription coactivators that play a crucial role in the response of plants to abiotic stress. In this study, we analyzed the genomic data of five Solanaceae plants and identified a total of 21 MBF1 genes. The expansion of MBF1a and MBF1b subfamilies was attributed to whole-genome duplication (WGD), and the expansion of the MBF1c subfamily occurred through transposed duplication (TRD). Collinearity analysis within Solanaceae species revealed collinearity between members of the MBF1a and MBF1b subfamilies, whereas the MBF1c subfamily showed relative independence. The gene expression of SlER24 was induced by sodium chloride (NaCl), polyethylene glycol (PEG), ABA (abscisic acid), and ethrel treatments, with the highest expression observed under NaCl treatment. The overexpression of SlER24 significantly enhanced the salt tolerance of tomato, and the functional deficiency of SlER24 decreased the tolerance of tomato to salt stress. SlER24 enhanced antioxidant enzyme activity to reduce the accumulation of reactive oxygen species (ROS) and alleviated plasma membrane damage under salt stress. SlER24 upregulated the expression levels of salt stress-related genes to enhance salt tolerance in tomato. In conclusion, this study provides basic information for the study of the MBF1 family of Solanaceae under abiotic stress, as well as a reference for the study of other plants.

Transcriptomics and Metabolomics Reveal the Critical Genes of Carotenoid Biosynthesis and Color Formation of Goji (Lycium barbarum L.) Fruit Ripening

Carotenoids in goji (Lycium barbarum L.) have excellent health benefits, but the underlying mechanism of carotenoid synthesis and color formation in goji fruit ripening is still unclear. The present study uses transcriptomics and metabolomics to investigate carotenoid biosynthesis and color formation differences in N1 (red fruit) and N1Y (yellow fruit) at three stages of ripening. Twenty-seven carotenoids were identified in N1 and N1Y fruits during the M1, M2, and M3 periods, with the M2 and M3 periods being critical for the difference in carotenoid and color between N1 and N1Y fruit. Weighted gene co-expression network analysis (WGCNA), gene trend analysis, and correlation analysis suggest that PSY1 and ZDS16 may be important players in the synthesis of carotenoids during goji fruit ripening. Meanwhile, 63 transcription factors (TFs) were identified related to goji fruit carotenoid biosynthesis. Among them, four TFs (CMB1-1, WRKY22-1, WRKY22-3, and RAP2-13-like) may have potential regulatory relationships with PSY1 and ZDS16. This work sheds light on the molecular network of carotenoid synthesis and explains the differences in carotenoid accumulation in different colored goji fruits.

Front-face excitation-emission matrix fluorescence spectroscopy combined with interpretable deep learning for the rapid identification of the storage year of Ningxia wolfberry

Ningxia wolfberry stored for many years may be disguised as fresh wolfberry by unscrupulous traders and sold for huge profits. In this work, the front-face excitation-emission matrix (FF-EEM) fluorescence spectroscopy coupled with interpretable deep learning was proposed to identify the storage year of Ningxia wolfberry in a lossless, fast and accurate way. Alternating trilinear decomposition (ATLD) algorithm was used to decompose the three-way data array obtained by Ningxia wolfberry samples, extracting the chemically meaningful information. Meanwhile, a convolutional neural network (CNN) model for the identification of the storage year of Ningxia wolfberry, called EEMnet, was proposed. The model successfully classified wolfberry samples from different storage years by extracting the subtle feature differences of the spectra, and the correct classification rate of the training set, test set and prediction set was more than 98%. In addition, a series of interpretability analyses were implemented to break the "black box" of the deep learning model. These results indicated that the method based on FF-EEM fluorescence spectroscopy combined with EEMnet could quickly and accurately identify the year of Ningxia wolfberry in a green way, providing a new idea for the identification of the storage years of Chinese medicinal materials.

Chromosome-level genome assembly and annotation of the prickly nightshade Solanum rostratum Dunal

The prickly nightshade Solanum rostratum, an annual malignant weed, is native to North America and has globally invaded 34 countries, causing serious threats to ecosystems, agriculture, animal husbandry, and human health. In this study, we constructed a chromosome-level genome assembly and annotation of S. rostratum. The contig-level genome was initially assembled in 898.42 Mb with a contig N50 of 62.00 Mb from PacBio high-fidelity reads. With Hi-C sequencing data scaffolding, 96.80% of the initially assembled sequences were anchored and orientated onto 12 pseudo-chromosomes, generating a genome of 869.69 Mb with a contig N50 of 72.15 Mb. We identified 649.92 Mb (72.26%) of repetitive sequences and 3,588 non-coding RNAs in the genome. A total of 29,694 protein-coding genes were predicted, with 28,154 (94.81%) functionally annotated genes. We found 99.5% and 91.3% complete embryophyta_odb10 genes in the pseudo-chromosomes genome and predicted gene datasets by BUSCO assessment. The present genomic resource provides essential information for subsequent research on the mechanisms of environmental adaptation of S. rostratum and host shift in Colorado potato beetles.

Integrated Metabolome and Transcriptome during Fruit Development Reveal Metabolic Differences and Molecular Basis between Lycium barbarum and Lycium ruthenicum

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.

Identification of Terpenoid Compounds and Toxicity Assays of Essential Oil Microcapsules from Artemisia stechmanniana

Plant essential oils, as biological pesticides, have been reviewed from several perspectives and play a key role in chemical ecology. However, plant essential oils show rapid degradation and vulnerability during actual usage. In this study, we conducted a detailed analysis of the compounds present in the essential oils of A. stechmanniana using gas chromatography-mass spectrometry (GC-MS). The results showed seventeen terpenoid compounds in the A. stechmanniana oil, with four major terpenoid compounds, i.e., eucalyptol (15.84%), (+)-2-Bornanone (16.92%), 1-(1,2,3-Trimethyl-cyclopent-2-enyl)-ethanone (25.63%), and (-)-Spathulenol (16.38%), in addition to an amount of the other terpenoid compounds (25.26%). Indoor toxicity assays were used to evaluate the insecticidal activity of Artemisia stechmanniana essential oil against Aphis gossypii, Frankliniella occidentalis, and Bactericera gobica in Lycium barbarum. The LC₅₀/LD₅₀ values of A. stechmanniana essential oils against A. gossypii, F. occidentalis, and B. gobica were 5.39 mg/mL, 0.34 mg/L, and 1.40 μg/insect, respectively, all of which were highly efficient compared with azadirachtin essential oil. Interestingly, A. stechmanniana essential oil embedded in β-cyclodextrin (microencapsule) remained for only 21 days, whereas pure essential oils remained for only 5 days. A field efficacy assay with the A. stechmanniana microencapsule (AM) and doses at three concentrations was conducted in Lycium barbarum, revealing that the insecticidal activities of AM showed high efficiency, maintained a significant control efficacy at all concentrations tested, and remained for 21 days. Our study identified terpenoid compounds from untapped Artemisia plants and designed a novel method against pests using a new biopesticide on L. barbarum.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Genome-Wide Identification and Expression Analysis of the Aquaporin Gene Family in Lycium barbarum during Fruit Ripening and Seedling Response to Heat Stress

Plant−water relations mediated by aquaporins (AQPs) play vital roles in both key plant growth processes and responses to environmental challenges. As a well-known medicinal and edible plant, the harsh natural growth habitat endows Lycium plants with ideal materials for stress biology research. However, the details of their molecular switch for water transport remain unclear. In the present work, we first identified and characterized AQP family genes from Lycium (L.) barbarum at the genome scale and conducted systemic bioinformatics and expression analyses. The results showed that there were 38 Lycium barbarum AQPs (LbAQPs) in L. barbarum, which were classified into four subfamilies, including 17 LbPIP, 9 LbTIP, 10 LbNIP, and 2 LbXIP. Their encoded genes were unevenly distributed on all 12 chromosomes, except chromosome 10. Three of these genes encoded truncated proteins and three genes underwent clear gene duplication events. Cis-acting element analysis indicated that the expression of LbAQPs may be mainly regulated by biotic/abiotic stress, phytohormones and light. The qRT-PCR assay indicated that this family of genes presented a clear tissue-specific expression pattern, in which most of the genes had maximal transcript levels in roots, stems, and leaves, while there were relatively lower levels in flowers and fruits. Most of the LbAQP genes were downregulated during L. barbarum fruit ripening and presented a negative correlation with the fruit relative water content (RWC). Most of their transcripts presented a quick and sharp upregulation response to heat stress following exposure of the 2-month-old seedlings to a 42 °C temperature for 0, 1, 3, 12, or 24 h. Our results proposed that LbAQPs were involved in L. barbarum key development events and abiotic stress responses, which may lay a foundation for further studying the molecular mechanism of the water relationship of Lycium plants, especially in harsh environments.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

Genome-Wide Identification and Analysis of the BBX Gene Family and Its Role in Carotenoid Biosynthesis in Wolfberry (Lycium barbarum L.)

The B-box proteins (BBXs) are a family of zinc-finger transcription factors with one/two B-Box domain(s) and play important roles in plant growth and development as well as stress responses. Wolfberry (Lycium barbarum L.) is an important traditional medicinal and food supplement in China, and its genome has recently been released. However, comprehensive studies of BBX genes in Lycium species are lacking. In this study, 28 LbaBBX genes were identified and classified into five clades by a phylogeny analysis with BBX proteins from Arabidopsis thaliana and the LbaBBXs have similar protein motifs and gene structures. Promoter cis-regulatory element prediction revealed that LbaBBXs might be highly responsive to light, phytohormone, and stress conditions. A synteny analysis indicated that 23, 20, 8, and 5 LbaBBX genes were orthologous to Solanum lycopersicum, Solanum melongena, Capsicum annuum, and Arabidopsis thaliana, respectively. The gene pairs encoding LbaBBX proteins evolved under strong purifying selection. In addition, the carotenoid content and expression patterns of selected LbaBBX genes were analyzed. LbaBBX2 and LbaBBX4 might play key roles in the regulation of zeaxanthin and antheraxanthin biosynthesis. Overall, this study improves our understanding of LbaBBX gene family characteristics and identifies genes involved in the regulation of carotenoid biosynthesis in wolfberry.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Five Solanaceae Species and Identification of Members Regulating Carotenoid Biosynthesis in Wolfberry

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.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium species.

The branch-thorn occurrence of Lycium ruthenicum is associated with leaf DNA hypermethylation in response to soil water content

BACKGROUND

Lycium ruthenicum is an eco-economic shrub which can exist in two forms, thorny and thornless under varying soil moisture conditions. The aim of this study was to determine if the two forms of L. ruthenicum were influenced by soil water content (SWC) and to test the three-way link among SWC, occurrence of branch-thorn and DNA methylation modification.

METHODS AND RESULTS

Here, pot experiment was carried out to reveal the influence of SWC on the occurrence of branch-thorn and then paraffin sections, scanning electron microscope and methylation-sensitive amplification polymorphism(MSAP) analysis were used to determine the three-way link among SWC, branch-thorn occurrence and DNA methylation. The results showed that (a) soil drought promoted the development of thorn primordium into branch-thorn and (b) branch-thorn covered axillary bud to protect it against drought and other stresses; (c) the branch-thorn occurrence response to drought was correlated with hypermethylation of CCGG sites and (d) thorny and thornless plants of a clone were distinguished successfully based on the MSAP profiles of their leaves.

CONCLUSIONS

Branch-thorns of the L. ruthenicum clone, which occurred in response to drought, covered axillary buds to protect them against drought and other stresses; thorn primordium of the clone did not develop into branch-thorn under the adequate soil moisture condition. The occurrence and absence of the branch-thorns were correlated with the hyper- and hypo-methylation, respectively. We proposed that the branch-thorn plasticity might be an adjustment strategy for the environment, which seems to support the theory of "Use in, waste out".