Green prickly ash (Zanthoxylum armatum) and red prickly ash (Zanthoxylum bungeanum) fruit have unique flavor and aroma characteristics that affect consumers’ purchasing preferences. However, differences in aroma components and relevant biosynthesis genes have not been systematically investigated in green and red prickly ash. Here, through the analysis of differentially expressed genes (DEGs), differentially abundant metabolites, and terpenoid biosynthetic pathways, we characterize the different aroma components of green and red prickly ash fruits and identify key genes in the terpenoid biosynthetic pathway. Gas chromatography-mass spectrometry (GC-MS) was used to identify 41 terpenoids from green prickly ash and 61 terpenoids from red prickly ash. Piperitone was the most abundant terpenoid in green prickly ash fruit, whereas limonene was most abundant in red prickly ash. Intergroup correlation analysis and redundancy analysis showed that HDS2, MVK2, and MVD are key genes for terpenoid synthesis in green prickly ash, whereas FDPS2 and FDPS3 play an important role in the terpenoid synthesis of red prickly ash. In summary, differences in the composition and content of terpenoids are the main factors that cause differences in the aromas of green and red prickly ash, and these differences reflect contrasting expression patterns of terpenoid synthesis genes.

Zanthoxylum bungeanum is an important spice and medicinal plant that is unique for its accumulation of abundant secondary metabolites, which create a characteristic aroma and tingling sensation in the mouth. Owing to the high proportion of repetitive sequences, high heterozygosity, and increased chromosome number of Z. bungeanum, the assembly of its chromosomal pseudomolecules is extremely challenging. Here, we present a genome sequence for Z. bungeanum, with a dramatically expanded size of 4.23 Gb, assembled into 68 chromosomes. This genome is approximately tenfold larger than that of its close relative Citrus sinensis. After the divergence of Zanthoxylum and Citrus, the lineage-specific whole-genome duplication event η-WGD approximately 26.8 million years ago (MYA) and the recent transposable element (TE) burst ~6.41 MYA account for the substantial genome expansion in Z. bungeanum. The independent Zanthoxylum-specific WGD event was followed by numerous fusion/fission events that shaped the genomic architecture. Integrative genomic and transcriptomic analyses suggested that prominent species-specific gene family expansions and changes in gene expression have shaped the biosynthesis of sanshools, terpenoids, and anthocyanins, which contribute to the special flavor and appearance of Z. bungeanum. In summary, the reference genome provides a valuable model for studying the impact of WGDs with recent TE activity on gene gain and loss and genome reconstruction and provides resources to accelerate Zanthoxylum improvement.

Zanthoxylum bungeanum Maxim. leaves (ZBLs) are rich in flavonoids and have become popular in nutrition, foods and medicine. However, the flavonoid components in ZBLs and the mechanism of flavonoid biosynthesis under drought stress have received little attention. Here, we performed an integrative analysis of the metabolome and transcriptome of ZBLs from HJ (Z. bungeanum cv. “Hanjiao”) and FJ (Z. bungeanum cv. “Fengjiao”) at four drought stages. A total of 231 individual flavonoids divided into nine classes were identified and flavones and flavonols were considered the most abundant flavonoid components in ZBLs. The total flavonoid content of ZBLs was higher in FJ; it increased in FJ under drought stress but decreased in HJ. Nine-quadrant analysis identified five and eight differentially abundant flavonoids in FJ and HJ leaves, respectively, under drought stress. Weighted gene correlation network analysis (WGCNA) identified nine structural genes and eight transcription factor genes involved in the regulation of flavonoid biosynthesis. Moreover, qRT-PCR results verified the accuracy of the transcriptome data and the reliability of the candidate genes. Taken together, our results reveal the flavonoid components of ZBLs and document changes in flavonoid metabolism under drought stress, providing valuable information for nutrition value and food utilization of ZBLs.

Chinese pepper, mainly includingZanthoxylum bungeanumandZanthoxylum armatum, is an economically important crop popular in Asian countries due to its unique taste characteristics and potential medical uses. Numerous cultivars of Chinese pepper have been developed in China through long-term domestication. To better understand the population structure, demographic history, and speciation of Chinese pepper, we performed a comprehensive analysis at a genome-wide level by analyzing 38,395 genomic SNPs that were identified in 112 cultivated and wild accessions using a high-throughput genome-wide genotyping-by-sequencing (GBS) approach. Our analysis provides genetic evidence of multiple splitting events occurring between and within species, resulting in at least four clades inZ. bungeanumand two clades inZ. armatum. Despite no evidence of recent admixture between species, we detected substantial gene flow within species. Estimates of demographic dynamics and species distribution modeling suggest that climatic oscillations during the Pleistocene (including the Penultimate Glaciation and the Last Glacial Maximum) and recent domestication events together shaped the demography and evolution of Chinese pepper. Our analyses also suggest that southeastern Gansu province is the most likely origin ofZ. bungeanumin China. These findings provide comprehensive insights into genetic diversity, population structure, demography, and adaptation inZanthoxylum.

Real-time reverse transcription quantitative PCR has become a common method for studying gene expression, however, the optimal selection of stable reference genes is a prerequisite for obtaining accurate quantification of transcript abundance. Suitable reference genes for RT-qPCR have not yet been identified for Chinese prickly ash (Zanthoxylum bungeanum Maxim.). Chinese prickly ash is the source of an important food seasoning in China. In recent years, Chinese prickly ash has also been developed as a medicinal plant. The expression stabilities of ten genes (18S, 28S, EF, UBA, UBQ, TIF, NTB, TUA, RPS, and TIF5A) were evaluated in roots, stems, leaves, flowers and fruits at five developmental stages and also under stress from cold, drought, and salt. To do this we used three different statistical algorithms: geNorm, NormFinder and BestKeeper. Among the genes investigated, UBA and UBQ were found to be most stable for the different cultivars and different tissues examined, UBQ and TIF for fruit developmental stage. Meanwhile, EF and TUA were most stable under cold treatment, EF and UBQ under drought treatment and NTB and RPS under salt treatment. UBA and UBQ for all samples evaluated were most stably expressed, but 18S, TUA and RPS were found to be generally unreliable as reference genes. Our results provide a basis for the future selection of reference genes for biological research with Chinese prickly ash, under a variety of conditions.

Zanthoxylum bungeanum is one of the most important medicinal and edible homologous plants because of its potential health benefits and unique flavors. The chemical components in compositions and contents vary with plant genotype variations and various environmental stress conditions. Fatty acids participate in various important metabolic pathways in organisms to resist biotic and abiotic stresses. To determine the variations in metabolic profiling and genotypes, the fatty acid profiling and key differential genes under low temperature stress in two Z. bungeanum varieties, cold-tolerant (FG) and sensitive (FX), were investigated. Twelve main fatty acids were found in two Z. bungeanum varieties under cold stress. Results showed that the contents of total fatty acids and unsaturated fatty acids in FG were higher than those in FX, which made FG more resistant to low temperature. Based on the result of orthogonal partial least squares discriminant analysis, palmitic acid, isostearic acid, linolenic acid and eicosenoic acid were the important differential fatty acids in FG under cold stress, while isomyristic acid, palmitic acid, isostearic acid, stearic acid, oleic acid, linolenic acid and eicosenoic acid were the important differential fatty acids in FX. Furthermore, fatty acid synthesis pathway genes fatty acyl-ACP thioesterase A (FATA), Delta (8)-fatty-acid desaturase 2 (SLD2), protein ECERIFERUM 3 (CER3), fatty acid desaturase 3 (FAD3) and fatty acid desaturase 5 (FAD5) played key roles in FG, and SLD2, FAD5, 3-oxoacyl-[acyl-carrier-protein] synthase I (KAS I), fatty acyl-ACP thioesterase B (FATB) and acetyl-CoA carboxylase (ACC) were the key genes responding to low temperature in FX. The variation and strategies of fatty acids in two varieties of Z. bungeanum were revealed at the metabolic and molecular level. This work provides a reference for the study of chemical components in plant stress resistance.