Integrated Transcriptome and Metabolome Analysis Revealed That Flavonoid Biosynthesis May Dominate the Resistance of Zanthoxylum bungeanum against Stem Canker

Stem canker of Zanthoxylum bungeanum is a devastating disease that seriously affects the plantation and industrial development of Z. bungeanum due to a lack of effective control measures. The objective of this study was to screen out resistant Z. bungeanum varieties and further explore their resistance mechanisms against stem canker. Results showed that the most resistant and susceptible varieties were, respectively, Doujiao (DJ) and Fengxian Dahongpao (FD). Combining transcriptomic and metabolomic analyses, we found that the genes and metabolites associated with the phenylpropanoid metabolism, especially flavonoid biosynthesis, were highly significantly enriched in DJ following pathogen infection compared with that in FD, which indicated that the flavonoid metabolism may positively dominate the resistance of Z. bungeanum. This finding was further confirmed by quantitative real-time polymerase chain reaction analysis, through which higher expression levels of core genes involved in flavonoid metabolism in resistant variety were observed. Moreover, by analyzing the differences in the flavonoid content in the stems of resistant and susceptible varieties and the antifungal activities of flavonoids extracted from Z. bungeanum stems, the conclusion that flavonoid metabolism positively regulates the resistance of Z. bungeanum was further supported. Our results not only aid in better understanding the resistance mechanisms of Z. bungeanum against stem canker but also promote the breeding and utilization of resistant varieties.

Expression Stabilities of Ten Candidate Reference Genes for RT-qPCR in Zanthoxylum bungeanum Maxim

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.

Comparative Transcriptome Analysis and Expression of Genes Reveal the Biosynthesis and Accumulation Patterns of Key Flavonoids in Different Varieties of Zanthoxylum bungeanum Leaves

Zanthoxylum bungeanum (Rutaceae), a popular food flavoring and traditional Chinese medicine ingredient, is an important cash crop. Its leaves are rich in flavonoids with multiple bioactivities. However, the transcriptional sequencing has not been investigated, and the molecular basis for the flavonoid biosynthesis remains unclear in this plant. This paper, the key flavonoids (epicatechin, rutin, hyperoside, trifolin, quercitrin, and afzelin) contents were determined in the leaves of 10 Z. bungeanum varieties from a common garden. Results show the leaves of Z. bungeanum mainly contained hyperoside (11.410-21.721 mg/g) and quercitrin (9.401-18.016 mg/g). The total content of these key components was the highest in Fengxian Dahongpao (66.012 mg/g) and the lowest in Fugu (32.223 mg/g). Three varieties (Hancheng stingless, Fugu, and Fengxian Dahongpao) with significant differences in the total content of key flavonoids were selected for transcriptome analysis to obtain flavonoid biosynthesis-related genes. In total, 83 522 unigenes were obtained, 40 668 (48.69%) unigenes were annotated, and 6656 differentially expressed genes (DEGs) were identified. Comparison of the other two varieties, Fugu had many differentially expressed genes indicating the particularity of its variety. Flavonoid-related DEGs of 22 structural genes, including three PALs, one CYP73A, three 4CLs, six CHSs, one CHI, one F3H, one DFR, two ANSs, one ANR, one FLS, and two CYP75B1s, as well as nine MYBs were obtained. These structural genes had different expression patterns in different Z. bungeanum varieties. It is worth noting that the genes expressing the flavonoid 3'5' hydroxylase are absent in Z. bungeanum. Furthermore, quantitative real-time PCR experiment showed consistent results in transcriptome analysis. The RNA-Seq data set of this study sheds lights on the molecular mechanism of flavonoid biosynthesis in Z. bungeanum, provides valuable information for the metabolic regulation of flavonoids, and may serve as a guide for future breeding programs.

miRNAs and their target genes regulate the antioxidant system of Zanthoxylum bungeanum under drought stress

Plants can accumulate a large amount of reactive oxygen species under adverse conditions such as drought and high temperature, which seriously affect the normal growth and development of plants. The antioxidant system can scavenge the reactive oxygen species produced under drought conditions and so mitigate oxidative damage. However, the regulation patterns of many miRNAs under drought stress are still unclear. The content of antioxidant enzymes and the expression patterns of miRNAs and their target genes related to antioxidant systems were studied under drought stress in Zanthoxylum bungeanum. The results indicate that under drought stress, POD, CAT, APX, proline, MDA and related genes all show positive responses to drought, while SOD and its genes showed a negative response. It is indicated that in the antioxidant process of Z. bungeanum, POD, CAT, and APX play a major role, and SOD plays a supporting role. In addition, GUS histochemical and RT-qPCR experimental results show that the expression levels of miRNAs and their target genes are basically negatively correlated, indicating that miRNAs can inhibit the expression of related genes and are also important regulators in the antioxidant system of Z. bungeanum. According to the expression patterns of antioxidant enzymes, miRNA and its target genes under drought stress, combined with previous research results, a model of plant antioxidant mechanism was constructed to provide a reference for further understanding of plant antioxidant mechanism.

Low Temperature Affects Fatty Acids Profiling and Key Synthesis Genes Expression Patterns in Zanthoxylum bungeanum Maxim

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.

The complex genome and adaptive evolution of polyploid Chinese pepper (Zanthoxylum armatum and Zanthoxylum bungeanum)

Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.

Authentication of Zanthoxylum Species Based on Integrated Analysis of Complete Chloroplast Genome Sequences and Metabolite Profiles

We performed chloroplast genome sequencing and comparative analysis of two Rutaceae species, Zanthoxylum schinifolium (Korean pepper tree) and Z. piperitum (Japanese pepper tree), which are medicinal and culinary crops in Asia. We identified more than 837 single nucleotide polymorphisms and 103 insertions/deletions (InDels) based on a comparison of the two chloroplast genomes and developed seven DNA markers derived from five tandem repeats and two InDel variations that discriminated between Korean Zanthoxylum species. Metabolite profile analysis pointed to three metabolic groups, one with Korean Z. piperitum samples, one with Korean Z. schinifolium samples, and the last containing all the tested Chinese Zanthoxylum species samples, which are considered to be Z. bungeanum based on our results. Two markers were capable of distinguishing among these three groups. The chloroplast genome sequences identified in this study represent a valuable genomics resource for exploring diversity in Rutaceae, and the molecular markers will be useful for authenticating dried Zanthoxylum berries in the marketplace.

Transcriptome mining of genes in Zanthoxylum armatum revealed ZaMYB86 as a negative regulator of prickly development

Zanthoxylum armatum DC. is an important economic tree species. Prickle is a type of trichome with special morphology, and there are a lot of prickles on the leaves of Z. armatum, which seriously restricts the development of Z. armatum industry. In this study, the leaves of Z. armatum cv. Zhuye (ZY) and its budding variety 'Rongchangwuci' (WC) (A less prickly mutant variety) at different developmental stages were used as materials, and the transcriptome sequencing data were analyzed. A total of 96,931 differentially expressed genes (DEGs) were identified among the samples, among which 1560 were candidate DEGs that might be involved in hormone metabolism. The contents of JA, auxin and CK phytohormones in ZY leaves were significantly higher than those in WC leaves. Combined with weighted gene co-expression network analysis, eight genes (MYC, IAA, ARF, CRE/AHK, PP2C, ARR-A, AOS and LOX) were identified, including 25 transcripts, which might affect the metabolism of the three hormones and indirectly participate in the formation of prickles. Combining with the proteins successfully reported in other plants to regulate trichome formation, ZaMYB86, a transcription factor of R2R3 MYB family, was identified through local Blast and phylogenetic tree analysis, which might regulate prickle formation of Z. armatum. Overexpression of ZaMYB86 in mutant A. thaliana resulted in the reduction of trichomes in A. thaliana leaves, which further verified that ZaMYB86 was involved in the formation of pickles. Yeast two-hybrid results showed that ZaMYB86 interacted with ZaMYB5. Furthermore, ZaMYB5 was highly homologous to AtMYB5, a transcription factor that regulated trichomes development, in MYB DNA binding domain. Taken together, these results indicated that ZaMYB86 and ZaMYB5 act together to regulate the formation of prickles in Z. armatum. Our findings provided a new perspective for revealing the molecular mechanism of prickly formation.

Regulation mechanisms of flavonoids biosynthesis of Hancheng Dahongpao peels (Zanthoxylum bungeanum Maxim) at different development stages by integrated metabolomics and transcriptomics analysis

BACKGROUND

Flavonoids have strong free radical scavenging and antioxidant capacity. The high abundance of flavonoids in Chinese prickly ash peels have many benefits to human health. In this study, 'Hancheng Dahongpao', a main cultivar, was taken as materials to investigate the flavonoids biosynthesis mechanism of Zanthoxylum bungeanum Maxim at three key development stages by integration of metabolomics and transcriptomics analysis.

RESULTS

A total of 19 differentially accumulated metabolites were identified, the key flavonoids compounds were kaempferol, quercetin and their glycoside derivatives, and two major anthocyanins (peonidin O-hexoside and peonidin 3-O-glucoside). 5 gene networks/modules including 15 important candidate genes were identified, which was highly correlated with flavonoids. Among these genes, ZM-163828 and ZM-184209 were strongly correlated with kaempferol and quercetin, and ZM-125833 and ZM-97481 were controlled the anthocyanins biosynthesis. Moreover, it was shown that MYB-ZM1, MYB-ZM3, MYB-ZM5, MYB-ZM6 and MYB-ZM7 coordinately controlled flavonoids accumulation through regulating the structural genes.

CONCLUSIONS

Generally, this study systematically revealed the flavonoids metabolic pathways and candidate genes involved in flavonoids biosynthesis and laid a foundation for the potential targets for the breeding of new valuable Chinese prickly ash cultivars.

Genetic structure and putative selective sweep in the pioneer tree, Zanthoxylum ailanthoides

Zanthoxylum ailanthoides Siebold & Zucc. is one of the most frequently encountered pioneer trees in Japanese warm-temperate evergreen oak forests. Our previous study in one region of Japan suggested high levels of population differentiation and putative natural selection acting on one of the nuclear loci analyzed. Here, we extend our analysis to study the genetic structure of 10 populations of Z. ailanthoides across Japan using 9 simple sequence repeat (SSR) loci for a better understanding of its genetic structure. First, the southernmost population (Kagoshima) in the samples was found to have the highest genetic diversity, suggesting there was a glacial refugium at or near the location of the population. Second, relatively strong genetic differentiation was found among populations, and there was a positive correlation between genetic distances and geographic distances (Mantel test; P < 0.001). Based on this information, we analyzed nucleotide variation at the putatively selected locus homologous to the gene encoding the ADP-glucose pyrophosphorylase large subunit (agpL). Despite the strong genetic differentiation among populations suggested by the SSR loci, the agpL locus was monomorphic in almost all populations analyzed. The results of this study strongly supported the possibility of a selective sweep at or near the agpL locus.

Target enrichment improves phylogenetic resolution in the genus Zanthoxylum (Rutaceae) and indicates both incomplete lineage sorting and hybridization events

BACKGROUND AND AIMS

Zanthoxylum is the only pantropical genus within Rutaceae, with a few species native to temperate eastern Asia and North America. Efforts using Sanger sequencing failed to resolve the backbone phylogeny of Zanthoxylum. In this study, we employed target-enrichment high-throughput sequencing to improve resolution. Gene trees were examined for concordance and sectional classifications of Zanthoxylum were evaluated. Off-target reads were investigated to identify putative single-copy markers for bait refinement, and low-copy markers for evidence of putative hybridization events.

METHODS

A custom bait set targeting 354 genes, with a median of 321 bp, was designed for Zanthoxylum and applied to 44 Zanthoxylum species and one Tetradium species as the outgroup. Illumina reads were processed via the HybPhyloMaker pipeline. Phylogenetic inferences were conducted using coalescent and maximum likelihood methods based on concatenated datasets. Concordance was assessed using quartet sampling. Additional phylogenetic analyses were performed on putative single and low-copy genes extracted from off-target reads.

KEY RESULTS

Four major clades are supported within Zanthoxylum: the African clade, the Z. asiaticum clade, the Asian-Pacific-Australian clade and the American-eastern Asian clade. While overall support has improved, regions of conflict are similar to those previously observed. Gene tree discordances indicate a hybridization event in the ancestor of the Hawaiian lineage, and incomplete lineage sorting in the American backbone. Off-target putative single-copy genes largely confirm on-target results, and putative low-copy genes provide additional evidence for hybridization in the Hawaiian lineage. Only two of the five sections of Zanthoxylum are resolved as monophyletic.

CONCLUSIONS

Target enrichment is suitable for assessing phylogenetic relationships in Zanthoxylum. Our phylogenetic analyses reveal that current sectional classifications need revision. Quartet tree concordance indicates several instances of reticulate evolution. Off-target reads are proven useful to identify additional phylogenetically informative regions for bait refinement or gene tree based approaches.

SMRT and Illumina RNA sequencing reveal the complexity of terpenoid biosynthesis in Zanthoxylum armatum

Sichuan pepper (Zanthoxylum armatum DC) is a popular spice and is often prescribed in traditional Chinese medicine to treat vomiting, diarrhea, ascariasis and eczema, among other conditions. Volatile oils from Z. armatum leaves contain active ingredients, with terpenoids being one of the main components. In the present study, the combination of sequencing data of Z. armatum from PacBio single molecule real time (SMRT) and Illumina RNA sequencing (RNA-Seq) platforms facilitated an understanding of the gene regulatory network of terpenoid biosynthesis in pepper leaves. The leaves of three developmental stages from two Z. armatum cultivars, 'Rongchangwuci' (WC) and 'Zhuye' (ZY), were selected as test materials to construct sequencing libraries. A total of 143,122 predictions of unique coding sequences, 105,465 simple sequence repeats, 20,145 transcription factors and 4719 long non-coding RNAs (lncRNAs) were identified, and 142,829 transcripts were successfully annotated. The occurrence of alternative splicing events was verified by reverse transcription PCR, and quantitative real-time PCR was used to confirm the expression pattern of six randomly selected lncRNAs. A total of 96,931 differentially expressed genes were filtered from different samples. According to functional annotation, a total of 560 candidate genes were involved in terpenoid synthesis, of which 526 were differentially expressed genes (DEGs). To identify the key genes involved in terpenoid biosynthesis, the module genes in different samples, including structural and transcription factors genes, were analyzed using the weighted gene co-expression network method, and the co-expression network of genes was constructed. Thirty-one terpenoids were identified by gas chromatography-mass spectrometry. The correlation between 18 compounds with significantly different contents and genes with high connectivity in the module was jointly analyzed in both cultivars, yielding 12 candidate DEGs presumably involved in the regulation of terpenoid biosynthesis. Our findings showed that full-length transcriptome SMRT and Illumina RNA-Seq can play an important role in studying organisms without reference genomes and elucidating the gene regulation of a biosynthetic pathway.

Development of conventional PCR and real-time PCR assays to discriminate the origins of Chinese pepper oil and herbal materials from Zanthoxylum

BACKGROUND

To ensure the safety, quality and therapeutic efficacy of processed foods and herbal medicines, it is important to identify and discriminate economically motivated adulterants. Zanthoxylum schinifolium is sold at a higher price than other Zanthoxylum species and is frequently adulterated with closely related Zanthoxylum species because of its high demand as a Korean food ingredient and medicinal material in markets. In addition, the pericarps of three Zanthoxylum species (Z. schinifolium, Z. bungeanum and Z. piperitum) are defined as herbal medicine Zanthoxyli Pericarpium in Korean pharmacopoeias, but not Z. piperitum in Chinese pharmacopoeias. Further confusion arises in the morphological similarity between Z. armatum (adulterant) and Z. bungeanum. Therefore, the aim of this study was to develop a sequence characterized amplified region (SCAR) marker for discrimination of four Zanthoxylum species.

RESULTS

With the goal of developing rapid and reliable tools for genetic discrimination of authentic Zanthoxyli Pericarpium, we designed species-specific SCAR markers, based on ITS2 sequences, that generate amplicons of less than 200 bp. Using these markers, we established both conventional and real-time PCR assay methods capable of differentiating samples at the species level. We validated the ability of SCAR markers to authenticate edible oil and herbal medicine, and confirmed that some herbal medicines contaminated with Z. armatum are being distributed as Zanthoxyli Pericarpium in Korean and Chinese markets.

CONCLUSIONS

The SCAR markers and PCR methods described represent powerful tools for protecting against adulteration and ensuring standardization of processed foods and herbal medicine. © 2018 The Authors. Journal of the Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Insights into the catalytic and regulatory mechanisms of dihydroflavonol 4-reductase, a key enzyme of anthocyanin synthesis in Zanthoxylum bungeanum

Accumulation of anthocyanins largely determines the fruit color, and dihydroflavonol 4-reductase (DFR) is a key enzyme involved in the formation of anthocyanins. However, the catalytic and regulatory mechanisms of DFR are unclear. In this study, the gene encoding DFR from Zanthoxylum bungeanum Maxim. was cloned and ZbDFR was analyzed in detail. The ZbDFR accepted dihydrokaempferol, dihydroquercetin and dihydromyricetin as substrates. Flavonols such as myricetin, quercetin and kaempferol significantly inhibited the activity of ZbDFR, while quercitrin and isoquercitrin slightly increased the activity. Quercetin was a competitive inhibitor at low concentrations, and it had a combined effect of competitive and noncompetitive inhibition at high concentrations, which was consistent with ZbDFR having two inhibitor binding sites. In addition, the content of different types of flavonoids in Z. bungeanum peel at green, semi-red and red stage was analyzed, and the in vivo results could be explained by the regulation of ZbDFR activity in vitro. Site-directed mutagenesis combined with enzyme activity experiments showed that Ser128, Tyr163, Phe164 and Lys167 are the key catalytic amino acid residues. The Ser128, Tyr163 and Lys167 were crucial for the hydrogen transfer reaction, and mutation of these amino acids resulted in the loss of all or most of the activity. Phe164 was found to be important for the regulation of ZbDFR by flavonols. Accordingly, ZbDFR is a node at which flavonoids regulate the synthesis of anthocyanins and proanthocyanins.

Monoterpene synthases contribute to the volatile production in tana (Zanthoxylum ailanthoides) through indigenous cultivation practices

Tana (Zanthoxylum ailanthoides), a perennial deciduous species in the Rutaceae family, possesses leaves with a unique fragrance that indigenous peoples incorporate into their traditional cuisine. In Kalibuan, the cultivated tana trees were pruned repeatedly to maintain a shorter height, which led to the growth of new leaves that were spicier and pricklier. Tana leaves contain a range of volatile terpenoids, and the pungent aroma may arise from the presence of monoterpenoids. To gain insight into the biosynthetic pathway, five candidate monoterpene synthase genes were cloned and characterized using a purified recombinant protein assay. The main product of Za_mTPS1, Za_mTPS2, and Za_mTPS5 is sabinene, geraniol, and (E)-β-ocimene, respectively. The main product of Za_mTPS3 and Za_mTPS4 is linalool. Real-time PCR analysis revealed that Za_mTPS1 and Za_mTPS5 are expressed at higher levels in prickly leaves of cultivated tana, suggesting that they may contribute to the distinctive aroma of this plant.

Physiological and transcriptomic analyses provide insight into thermotolerance in desert plant Zygophyllum xanthoxylum

BACKGROUND

Heat stress has adverse effects on the growth and reproduction of plants. Zygophyllum xanthoxylum, a typical xerophyte, is a dominant species in the desert where summer temperatures are around 40 °C. However, the mechanism underlying the thermotolerance of Z. xanthoxylum remained unclear.

RESULTS

Here, we characterized the acclimation of Z. xanthoxylum to heat using a combination of physiological measurements and transcriptional profiles under treatments at 40 °C and 45 °C, respectively. Strikingly, moderate high temperature (40 °C) led to an increase in photosynthetic capacity and superior plant performance, whereas severe high temperature (45 °C) was accompanied by reduced photosynthetic capacity and inhibited growth. Transcriptome profiling indicated that the differentially expressed genes (DEGs) were related to transcription factor activity, protein folding and photosynthesis under heat conditions. Furthermore, numerous genes encoding heat transcription shock factors (HSFs) and heat shock proteins (HSPs) were significantly up-regulated under heat treatments, which were correlated with thermotolerance of Z. xanthoxylum. Interestingly, the up-regulation of PSI and PSII genes and the down-regulation of chlorophyll catabolism genes likely contribute to improving plant performance of Z. xanthoxylum under moderate high temperature.

CONCLUSIONS

We identified key genes associated with of thermotolerance and growth in Z. xanthoxylum, which provide significant insights into the regulatory mechanisms of thermotolerance and growth regulation in Z. xanthoxylum under high temperature conditions.

Small RNA sequencing provides candidate miRNA-target pairs for revealing the mechanism of apomixis in Zanthoxylum bungeanum

BACKGROUND

Apomixis is a form of asexual reproduction that produces offspring without the need for combining male and female gametes, and the offspring have the same genetic makeup as the mother. Therefore, apomixis technology has great application potential in plant breeding. To identify the apomixis types and critical period, embryonic development at different flower development stages of Zanthoxylum bungeanum was observed by cytology.

RESULTS

The results show that the S3 stage is the critical period of apomixis, during which the nucellar cells develop into an adventitious primordial embryo. Cytological observations showed that the type of apomixis in Z. bungeanum is sporophytic apomixis. Furthermore, miRNA sequencing, miRNA-target gene interaction, dual luciferase reporter assay, and RT-qPCR verification were used to reveal the dynamic regulation of miRNA-target pairs in Z. bungeanum apomixis. The miRNA sequencing identified 96 mature miRNAs, of which 40 were known and 56 were novel. Additionally, 29 differentially expressed miRNAs were screened according to the miRNAs expression levels at the different developmental stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses showed that the target genes of the differentially expressed miRNAs were mainly enriched in plant hormone signal transduction, RNA biosynthetic process, and response to hormone pathways.

CONCLUSIONS

During the critical period of apomictic embryonic development, miR172c significantly reduces the expression levels of TOE3 and APETALA 2 (AP2) genes, thereby upregulating the expression of the AGAMOUS gene. A molecular regulation model of miRNA-target pairs was constructed based on their interactions and expression patterns to further understand the role of miRNA-target pairs in apomixis. Our data suggest that miR172c may regulates AGAMOUS expression by inhibiting TOE3 in the critical period of apomixis.

A high-quality chromosome-level genome assembly of the traditional Chinese medicinal herb Zanthoxylum nitidum

Dried roots of Zanthoxylum nitidum (2n=2x=70, family Rutaceae) was referred as "Liang-Mian-Zhen" in Chinese, acting as a valuable species due to its notable pharmacological activities. Herein, we combined PacBio HiFi data together with Hi-C mapping technology to construct a chromosome-scale reference genome assembly for Z. nitidum. The assembly reached a length of 2.24 Gb, successfully anchoring 99.31% of sequences onto 35 pseudo-chromosomes. Among these, 26 chromosomes achieved telomere-to-telomere assembly, and 11 chromosomes were gap-free. The contigs N50 and scaffolds N50 reached 57.61 Mb and 78.00 Mb, respectively. Transposable elements comprised 81.44% of the Z. nitidum genome, and over 78% of them were long-terminal repeat retrotransposon elements. Furthermore, 32,737 protein-coding genes were identified and 99.38% of all were functionally annotated. The completeness of the genome assembly and final gene sets reached 97.83% and 96.47% based on Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. Taken together, our results provided a high-quality chromosome-level assembly of Z. nitidum genome and will be a valuable resource that will facilitate breeding varieties with higher alkaloids content.

Transcriptomics integrated with metabolomics to characterize key pigment compounds and genes related to anthocyanin biosynthesis in Zanthoxylum bungeanum peel

Zanthoxylum bungeanum is an important condiment with high economic value and its peel color is one of the main quality indexes. However, the key pigment compounds and related genes are still unclear affecting the quality control of the plants. In this study, the contents of four types of pigments were measured in Z. bungeanum and flavonoids were identified as the most important pigments. Based on the targeted flavonoid metabolomics of Z. bungeanum peels, 14 key pigment compounds were screened out from 152 flavonoids, among which cyanidin-3-O-rutinoside and cyanidin-3-O-glucoside were the most critical compounds for peel color. They were further verified to be present in nine varieties of Z. bungeanum by HPLC fingerprints. The 14 compounds were all associated with flavonoid and anthocyanin biosynthesis pathways and the 39 differentially expressed genes related to these pathways were annotated and screened based on transcriptomics. The genes ZbDFR, ZbANS, and ZbUFGT were identified as three key genes for anthocyanin synthesis in Z. bungeanum peels. Further qRT-PCR results confirmed the reliability of transcriptomics and the accuracy of gene screening. Subsequent protein induced expression demonstrated that ZbANS and ZbUFGT were expressed after 12 h induced by IPTG while ZbDFR was expressed after 15 h. Further transient and stable transformation analysis confirmed that both anthocyanin content and the expression of ZbDFR were significantly increased in overexpression Z. bungeanum leaves and Nicotiana benthamiana. The functional effect of stable transformation of ZbDFR was more significant than that of transient transformation with a 7.67-fold/1.49-fold difference in total anthocyanin content and a 42.37-fold/12.32-fold difference in the expression of ZbDFR. This study provides new insights into the chemical composition and the molecular mechanisms of Z. bungeanum peel color and lays an effective foundation for the color quality control, multi-purpose utilization of Z. bungeanum and the creation of new germplasm.

Genome-wide identification of HD-ZIP gene family and screening of genes related to prickle development in Zanthoxylum armatum

Zanthoxylum armatum is an important cash crop for medicinal and food purposes in Asia. However, its stems and leaves are covered with a large number of prickles, which cause many problems in the production process. The homeodomain leucine zipper (HD-ZIP) gene family is a class of transcription factors unique to plants that play an important role in biological processes such as morphogenesis, signal transduction, and secondary metabolite synthesis. However, little is known about HD-ZIP gene information that may be involved in prickle development of Z. armatum. Here, we identified 76 ZaHDZ genes from the Z. armatum genome and classified them into four subfamilies (I-IV) based on phylogenetic analysis, a classification further supported by gene structure and conserved motif analysis. Seventy-six ZaHDZ genes were unevenly distributed on chromosomes. Evolutionary analysis revealed that the expansion of ZaHDZ genes mainly were due to whole-genome duplication (WGD) or segmental duplication, and they experienced strong purifying selection pressure in the process of evolution. A total of 47 cis-elements were identified in the promoter region of ZaHDZ genes. Quantitative real-time polymerase chain reaction analysis was performed on subfamily IV ZaHDZ gene expression levels in five tissues and under four hormone treatments. Finally, ZaHDZ16 was predicted to be the candidate gene most likely to be involved in prickle development of Z. armatum. These results contribute to a better understanding of the characteristics of HD-ZIP gene family and lay a foundation for further study on the function of genes related to prickle development of Z. armatum.

SSR molecular marker developments and genetic diversity analysis of Zanthoxylum nitidum (Roxb.) DC

Zanthoxylum nitidum (Roxb.) is a commonly used traditional Chinese medicine. However, the collection and protection of wild germplasm resources of Z. nitidum are still insufficient, and there is limited research on its genetic diversity and fingerprint. In the present study, 15 simple sequence repeat (SSR) markers were developed by genotyping based on multiplexed shotgun sequencing. The genetic diversity of 51 populations (142 individuals) of Z. nitidum was evaluated using these 15 SSRs. A total of 245 alleles (Na) were detected, with an average value of 16.333, and the average polymorphism information content was 0.756. The genetic distance among 51 populations was 0.164~1.000, with an average of 0.659. Analysis of molecular variance showed low genetic differentiation (40%) and high genetic differentiation (60%) between populations and individuals, respectively. The genetic differentiation coefficient (Fst) of the population was 0.338, indicating that 66.2% of the genetic variation occurred within the population, and the gene flow (Nm) was 0.636, demonstrating that the gene exchange between populations was low. Clustering analysis revealed that the genetic similarity coefficient was 0.30, dividing the 51 populations into 4 groups of 2, 17, 3, and 29 populations. There was no specific relationship between geographical location differences and genetic distance. The genetic diversity level of Z. nitidum is relatively high, and our results provide a theoretical basis for the rapid identification of Z. nitidum germplasm resources and variety selection.

A NAC transcription factor ZaNAC93 confers floral initiation, fruit development, and prickle formation in Zanthoxylum armatum

Zanthoxylum armatum is a dioecious prickly plant which developed apomictic reproduction. The increases in male flowers and prickle density in female plants lead to low yield and picking efficiency. However, little is known concerning the mechanisms of floral development and prickle formation. NAC is a well-known transcription factor that participates in multiple aspects of plant growth and development. Herein, we characterize the functions and regulatory mechanisms of candidate NACs controlling both traits in Z. armatum. A total of 159 ZaNACs were identified, and 16 of these were male-biased, represented by the NAP subfamily members ZaNAC93 and ZaNAC34, orthologs of AtNAC025 and AtNARS1/NAC2 respectively. Overexpression of ZaNAC93 in tomato led to modifications in flower and fruit development, including earlier flowering, increased numbers of lateral shoots and flowers, accelerated plant senescence, and reduced size and weight of fruits and seeds. In addition, the trichome density in leaves and inflorescences was dramatically reduced in ZaNAC93-OX lines. Overexpression of ZaNAC93 resulted in the up-/downregulation of genes associated with GA, ABA and JA signaling pathways, such as GAI, PYL and JAZ, as well as several TFs, including bZIP2, AGL11, FBP24 and MYB52. Yeast two-hybrid analysis revealed that ZaNAC93 protein could interact with AP1, GAI, bZIP2 and AGL11 in Z. armatum, which might contribute to floral induction, fruit growth, and trichome initiation. This work provides new insights into the molecular mechanisms of ZaNAC93 in reproductive development and prickle formation in Z. armatum.

Integrated Transcriptome and Metabolome Analysis Reveals the Molecular Mechanism of Rust Resistance in Resistant (Youkang) and Susceptive (Tengjiao) Zanthoxylum armatum Cultivars

Chinese pepper rust is a live parasitic fungal disease caused by Coleosporium zanthoxyli, which seriously affects the cultivation and industrial development of Z. armatum. Cultivating and planting resistant cultivars is considered the most economical and environmentally friendly strategy to control this disease. Therefore, the mining of excellent genes for rust resistance and the analysis of the mechanism of rust resistance are the key strategies to achieve the targeted breeding of rust resistance. However, there is no relevant report on pepper rust resistance at present. The aim of the present study was to further explore the resistance mechanism of pepper by screening the rust-resistant germplasm resources in the early stage. Combined with the analysis of plant pathology, transcriptomics, and metabolomics, we found that compared with susceptible cultivar TJ, resistant cultivar YK had 2752 differentially expressed genes (DEGs, 1253 up-, and 1499 downregulated) and 321 differentially accumulated metabolites (DAMs, 133 up- and 188 down-accumulated) after pathogen infection. And the genes and metabolites related to phenylpropanoid metabolism were highly enriched in resistant varieties, which indicated that phenylpropanoid metabolism might mediate the resistance of Z. armatum. This finding was further confirmed by a real-time quantitative polymerase chain reaction analysis, which revealed that the expression levels of core genes involved in phenylpropane metabolism in disease-resistant varieties were high. In addition, the difference in flavonoid and MeJA contents in the leaves between resistant and susceptible varieties further supported the conclusion that the flavonoid pathway and methyl jasmonate may be involved in the formation of Chinese pepper resistance. Our research results not only help to better understand the resistance mechanism of Z. armatum rust but also contribute to the breeding and utilization of resistant varieties.

The UV-B-Induced Transcription Factor HY5 Regulated Anthocyanin Biosynthesis in Zanthoxylum bungeanum

Pericarp color is an important economic characteristic of Zanthoxylum bungeanum. Anthocyanins are the main reason for the pericarp's red appearance in Z. bungeanum. In this study, through the combined analysis of the metabolome and transcriptome, HY5, whose expression is highly correlated to changes in the anthocyanin content, was screened and identified. Under natural ripening conditions, the Z. bungeanum fruit gradually changed in color from green to red, while bagging resulted in the fruit maintaining its green color. After unbagging, the fruit gradually turned red, and the ZbHY5 expression and anthocyanin content increased. In addition, the leaves changed from green to red after exposure to UV-B radiation, and the ZbHY5 expression and anthocyanin content increased. The transient overexpression of ZbHY5 deepened the redness of the Z. bungeanum leaves and promoted the expression of ZbHY5 and ZbMYB113 as well as anthocyanin accumulation. Bimolecular fluorescence complementation (BIFC) showed that there was an interaction between ZbHY5 and ZbMYB113. These results revealed that under UV-B irradiation, ZbHY5 might regulate the expression levels of the structural genes related to anthocyanin biosynthesis through combination with ZbMYB113, thereby affecting anthocyanin accumulation. This finding provides useful insights for further studies focusing on UV-B-induced anthocyanin accumulation in Z. bungeanum.

The ERF transcription factor ZbERF3 promotes ethylene-induced anthocyanin biosynthesis in Zanthoxylum bungeanum

Ethylene regulates fruit ripening, and in Zanthoxylum bungeanum, fruit color deepened with increasing of ethylene during fruit ripening. However, the molecular mechanism of this physiological process was still unclear. In this study, through the combined analysis of transcriptome and metabolome, it was found that ethylene release was consistent with anthocyanin synthesis, and ethylene response factors (ERFs) were significantly related to anthocyanin biosynthesis during the fruit ripening of Z. bungeanum. Ethylene treatment significantly induced fruit coloration and promoted anthocyanin synthesis and the expression of ZbERF3. Furthermore, Yeast one-hybrid assays and Luciferase reporter assays demonstrated that ZbERF3 directly bound to the promoter of ZbMYB17 and transcriptionally activated its expression. What's more, it was demonstrated that ZbMYB17 directly bound to the promoter of ZbANS, promoting anthocyanin biosynthesis. Overall, this study revealed the mechanism of ERF and MYB synergistically regulating the coloration of Z. bungeanum fruit.