Zanthoxylum-specific whole genome duplication and recent activity of transposable elements in the highly repetitive paleotetraploid Z. bungeanum genome

Citrus genomes: past, present and future

Over the past decade, genome sequencing and assembly approaches have been greatly improved, resulting in the assembly of many genomes for citrus, including wild, domesticated, and citrus-related genomes. Improvements in technologies have led to assembled genomes with higher completeness, contiguity, quality, and accuracy that have greatly facilitated annotation and analysis. This review summarizes the evolution of the sequencing, assembly, and annotation technologies leading to citrus genomes over the past 11 years, a comprehensive evaluation of their quality, contiguity, and completeness, and the major findings and applications. Of the 50 genomes now available, 35 have been assembled to chromosome level and 15 to draft level, and 14 were haplotype-resolved assemblies. To date there have been four pangenome-wide studies for citrus. The very recent genomes assembled with long-read sequencing have achieved >99% and >98% assembly and annotation completeness (BUSCO), respectively. However, some early genomes are not of the same high quality as more recently sequenced genomes and would benefit from re-sequencing. A more comprehensive pangenome based upon a larger set of species and genotypes assembled at the haplotype level would allow genomics to deliver the maximum benefits for citrus improvement and research.

Identification and Characterization of SQUAMOSA Promoter Binding Protein-like Transcription Factor Family Members in Zanthoxylum bungeanum and Their Expression Profiles in Response to Abiotic Stresses

Plant-specific transcription factors known as SQUAMOSA promoter binding protein-like (SPL) genes are essential for development, growth, and abiotic stress responses. While the SPL gene family has been extensively studied in various plant species, a systematic characterization in Zanthoxylum bungeanum (Zb) is lacking. This study used transcriptomic and bioinformatics data to conduct a thorough genomic identification and expression investigation of the ZbSPL gene family. Eight subfamilies including 73 ZbSPL members were identified, most of which are predicted to be localized in the nucleus. Ka/Ks ratio analysis indicates that most ZbSPL genes have undergone purifying selection. According to evolutionary research, segmental duplication is a major factor in the amplification of the ZbSPL gene family. Gene structures, conserved motifs, and domains were found to be highly conserved among paralogs. Cis-element research revealed that ZbSPLs may be implicated in hormone and abiotic stress responses. Codon usage pattern analysis showed that the ZbSPL gene family was more inclined to A/T base endings; the higher the A/T content, the stronger the preference of the codons; and the use pattern was mainly affected by natural selection. Additionally, 36 ZbSPLs were found to be potential targets of miR156. RNA-seq demonstrated that SPL genes in Zb are differentially expressed in response to distinct abiotic stressors. Two ZbSPL genes (ZbSPL10 and ZbSPL17) were implicated in the response to salt stress, while four ZbSPL genes (ZbSPL06, ZbSPL43, ZbSPL60, and ZbSPL61) showed response to drought stress, based on a qRT-PCR investigation of the ZbSPL genes under various abiotic stress conditions. This study will help us gain a deeper understanding of the functions of ZbSPLs and lay a genetic foundation for future breeding of high-quality, highly abiotic resistant varieties of Z. bungeanum.

Metabolome and transcriptome association study reveals biosynthesis of specialized benzylisoquinoline alkaloids in Phellodendron amurense

Objective

Benzylisoquinoline alkaloids (BIAs) have pharmacological functions and clinical use. BIAs are mainly distributed in plant species across the order Ranunculales and the genus Phellodendron from Sapindales. The BIA biosynthesis has been intensively investigated in Ranunculales species. However, the accumulation mechanism of BIAs in Phellodendron is largely unknown. The aim of this study is to unravel the biosynthetic pathways of BIAs in Phellodendron amurens.

Methods

The transcriptome and metabolome data from 18 different tissues of P. amurense were meticulously sequenced and subsequently subjected to a thorough analysis. Weighted gene co-expression network analysis (WGCNA), a powerful systems biology approach that facilitates the construction and subsequent analysis of co-expression networks, was utilized to identify candidate genes involved in BIAs biosynthesis. Following this, recombinant plasmids containing candidate genes were expressed in Escherichia coli, a widely used prokaryotic expression system. The purpose of this genetic engineering endeavor was to express the candidate genes within the bacteria, thereby enabling the assessment of the resultant enzyme activity.

Results

The synonymous substitutions per synonymous site for paralogs indicated that at least one whole genome duplication event has occurred. The potential BIA biosynthetic pathway of P. amurense was proposed, and two PR10/Bet v1 members, 14 CYP450s, and 33 methyltransferases were selected as related to BIA biosynthesis. One PR10/Bet v1 was identified as norcoclaurine synthase, which could catalyze dopamine and 4-hydroxyphenylacetaldehyde into (S)-norcoclaurine.

Conclusion

Our studies provide important insights into the biosynthesis and evolution of BIAs in non-Ranunculales species.

The genomes of Australian wild limes

Australian wild limes occur in highly diverse range of environments and are a unique genetic resource within the genus Citrus. Here we compare the haplotype-resolved genome assemblies of six Australian native limes, including four new assemblies generated using PacBio HiFi and Hi-C sequencing data. The size of the genomes was between 315 and 391 Mb with contig N50s from 29.5 to 35 Mb. Gene completeness of the assemblies was estimated to be from 98.4 to 99.3% and the annotations from 97.7 to 98.9% based upon BUSCO, confirming the high contiguity and completeness of the assembled genomes. High collinearity was observed among the genomes and the two haplotype assemblies for each species. Gene duplication and evolutionary analysis demonstrated that the Australian citrus have undergone only one ancient whole-genome triplication event during evolution. The highest number of species-specific and expanded gene families were found in C. glauca and they were primarily enriched in purine, thiamine metabolism, amino acids and aromatic amino acids metabolism which might help C. glauca to mitigate drought, salinity, and pathogen attacks in the drier environments in which this species is found. Unique genes related to terpene biosynthesis, glutathione metabolism, and toll-like receptors in C. australasica, and starch and sucrose metabolism genes in both C. australis and C. australasica might be important candidate genes for HLB tolerance in these species. Expanded gene families were not lineage specific, however, a greater number of genes related to plant-pathogen interactions, predominantly disease resistant protein, was found in C. australasica and C. australis.

Genome-wide analysis of AP2/ERF transcription factors that regulate fruit development of Chinese prickly ash

BACKGROUND

AP2/ERF is a large family of plant transcription factor proteins that play essential roles in signal transduction, plant growth and development, and responses to various stresses. The AP2/ERF family has been identified and verified by functional analysis in various plants, but so far there has been no comprehensive study of these factors in Chinese prickly ash. Phylogenetic, motif, and functional analyses combined with transcriptome analysis of Chinese prickly ash fruits at different developmental stages (30, 60, and 90 days after anthesis) were conducted in this study.

RESULTS

The analysis identified 146 ZbAP2/ERF genes that could be classified into 15 subgroups. The motif analysis revealed the presence of different motifs or elements in each group that may explain the functional differences between the groups. ZbERF13.2, ZbRAP2-12, and ZbERF2.1 showed high levels of expression in the early stages of fruit development. ZbRAP2-4, and ZbERF3.1 were significantly expressed at the fruit coloring stage (R2 and G2). ZbERF16 were significantly expressed at fruit ripening and expression level increased as the fruit continued to develop. Relative gene expression levels of 6 representative ZbAP2/ERFs assessed by RT-qPCR agreed with transcriptome analysis results.

CONCLUSIONS

These genes identified by screening can be used as candidate genes that affect fruit development. The results of the analysis can help guide future genetic improvement of Chinese prickly ash and enrich our understanding of AP2/ERF transcription factors and their regulatory functions in plants.

Genome-Wide Identification and Characterization of the RWP-RK Proteins in Zanthoxylum armatum

Apomixis is a common reproductive characteristic of Zanthoxylum plants, and RWP-RKs are plant-specific transcription factors known to regulate embryonic development. However, the genome-wide analysis and function prediction of RWP-RK family genes in Z. armatum are unclear. In this study, 36 ZaRWP-RK transcription factors were identified in the genome of Z. armatum, among which 15 genes belonged to the RKD subfamily and 21 belonged to the NLP subfamily. Duplication events of ZaRWP-RK genes were mainly segmental duplication, and synteny analysis revealed a close phylogenetic relationship between Z. armatum and Arabidopsis. The analysis of cis-elements indicated that ZaRWP-RK genes may be involved in the regulation of the embryonic development of Z. armatum by responding to plant hormones such as abscisic acid, auxin, and gibberellin. Results of a real-time PCR showed that the expression levels of most ZaRWP-RK genes were significantly increased from flowers to young fruits. Protein-protein interaction network analysis further revealed the potential roles of the ZaRWP-RK proteins in apomixis. Collectively, this study is expected to improve our understanding of ZaRWP-RK transcription factors and provide a theoretical basis for future investigations into the ZaRWP-RK genes and their regulatory mechanisms in the apomixis process of Z. armatum.

Genome sequencing and comparative genomics reveal insights into pathogenicity and evolution of Fusarium zanthoxyli, the causal agent of stem canker in prickly ash

BACKGROUND

Fusarium zanthoxyli is a destructive pathogen causing stem canker in prickly ash, an ecologically and economically important forest tree. However, the genome lack of F. zanthoxyli has hindered research on its interaction with prickly ash and the development of precise control strategies for stem canker.

RESULTS

In this study, we sequenced and annotated a relatively high-quality genome of F. zanthoxyli with a size of 43.39 Mb, encoding 11,316 putative genes. Pathogenicity-related factors are predicted, comprising 495 CAZymes, 217 effectors, 156 CYP450s, and 202 enzymes associated with secondary metabolism. Besides, a comparative genomics analysis revealed Fusarium and Colletotrichum diverged from a shared ancestor approximately 141.1 ~ 88.4 million years ago (MYA). Additionally, a phylogenomic investigation of 12 different phytopathogens within Fusarium indicated that F. zanthoxyli originated approximately 34.6 ~ 26.9 MYA, and events of gene expansion and contraction within them were also unveiled. Finally, utilizing conserved domain prediction, the results revealed that among the 59 unique genes, the most enriched domains were PnbA and ULP1. Among the 783 expanded genes, the most enriched domains were PKc_like kinases and those belonging to the APH_ChoK_Like family.

CONCLUSION

This study sheds light on the genetic basis of F. zanthoxyli's pathogenicity and evolution which provides valuable information for future research on its molecular interactions with prickly ash and the development of effective strategies to combat stem canker.

Genomic survey and expression analysis of cellulose synthase superfamily and COBRA-like gene family in Zanthoxylum bungeanum stipule thorns

The Cellulose Synthase gene (CS) superfamily and COBRA-like (COBL) gene family are essential for synthesizing cellulose and hemicellulose, which play a crucial role in cell wall biosynthesis and the hardening of plant tissues. Our study identified 126 ZbCS and 31 ZbCOBL genes from the Zanthoxylum bungeanum (Zb) genome. Phylogenetic analysis and conservative domain analysis unfolded that ZbCS and ZbCOBL genes were divided into seven and two subfamilies, respectively. Gene duplication data suggested that more than 75% of these genes had tandem and fragment duplications. Codon usage patterns analysis indicated that the ZbCS and ZbCOBL genes prefer ending with A/T base, with weak codon preference. Furthermore, seven key ZbCS and five key ZbCOBL genes were identified based on the content of cellulose and hemicellulose and the expression characteristics of ZbCS and ZbCOBL genes in various stages of stipule thorns. Altogether, these results improve the understanding of CS and COBL genes and provide valuable reference data for cultivating Zb with soft thorns.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01432-x.

Study on SH-SY5Y autophagy inhibition and apoptosis induced by methanol extract of Zanthoxylum armatum DC. based on mTOR signal pathway

Background

Zanthoxylum armatum DC. (ZADC) is a novel food raw material resource, offering both edible and medicinal properties. Recent research has unveiled the toxic nature of ZADC, particularly its close association with the nervous system. In a prior study, we observed that administering methanol extract of Zanthoxylum armatum DC. (MZADC) to rats via gavage at a dose of 1.038 g/kg resulted in various neurotoxicity symptoms, including excessive salivation, reduced mobility, unsteady gait, muscle twitching, and altered respiratory rates.

Materials and methods

We conducted cell-based research to assess the safety of ZADC and elucidate its potential toxic mechanism. In addition, we used experimental methods such as Cell Counting Kit-8, Western blot, and Flow cytometry to detect cytotoxicity in SH-SY5Y cells after intervention with MZADC.

Results

Following exposure of SY-SY5Y cells with MZADC, a substantial decline in cell viability was observed, accompanied by a concentration-dependent increase in intracellular reactive oxygen species (ROS) levels. Additionally, MZADC induced cellular oxidative stress, leading to elevated malonic dialdehyde (MDA) and superoxide dismutase (SOD) concentrations while decreasing glutathione (GSH) levels. Furthermore, MZADC induced apoptosis at varying doses (20, 40, and 60 μg/mL), and this effect was associated with increased Caspase-3, Bax expressions, and reduced Bcl2 and Bcl2/Bax expressions. In addition, the investigation revealed that MZADC induced autophagy inhibition in SH-SY5Y cells by activating the mTOR signaling pathway, resulting in a decrease in LC3II/LCI and Beclin-1, while increasing p-mTOR/mTOR, p62.

Conclusion

Consequently, this study suggests that MZADC triggers the mTOR pathway through oxidative stress in SH-SY5Y cells, ultimately leading to apoptosis. Understanding the toxicity mechanisms associated with ZADC can offer a valuable theoretical and experimental basis for its development and utilization.

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.

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.

Karyotype and Phylogenetic Relationship Analysis of Five Varieties and Cultivars of Zanthoxylum armatum Based on Oligo-FISH

Green prickly ash (Zanthoxylum armatum) has edible and medicinal value and is an economically significant plant in many countries. Z. armatum has many cultivars and varieties with similar phenotypes that are difficult to distinguish via traditional methods. In this study, we utilized oligo-FISH to distinguish five varieties and cultivars of Z. armatum on the basis of three oligonucleotide probes of 5S rDNA, (AG₃T₃)₃, and (GAA)₆. Karyotype analysis of the five varieties and cultivars of Z. armatum showed that the Z. armatum 'Tengjiao' karyotype formula was 2n = 2x = 98m with karyotype type 1C and an arm ratio of 4.3237, including two pairs of 5S rDNA signals and five pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'Youkangtengjiao' was 2n = 2x = 128m + 8sm with karyotype type 2B and an arm ratio of 3.5336, including three pairs of 5S rDNA signals and 17 pairs of (GAA)₆ signals. The karyotype formula of Z. armatum var. novemfolius was 2n = 2x = 134m + 2sm with karyotype type 1C and an arm ratio of 5.5224, including two pairs of 5S rDNA signals and eight pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'YT-03' was 2n = 2x = 2M + 128m + 4sm + 2st with karyotype type 2C and an arm ratio of 4.1829, including three pairs of 5S rDNA signals and nine pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'YT-06' was 2n = 2x = 126m + 10sm with cytotype 2B and an arm ratio of 3.3011, including three pairs of 5S rDNA signals and two pairs of (GAA)₆ signals. The five varieties and cultivars of Z. armatum had (AG₃T₃)₃ signals on all chromosomes. The chromosomal symmetry of Z. armatum 'Tengjiao' was high, whereas the chromosomal symmetry of Z. armatum 'YT-03' was low, with the karyotypes of the five materials showing a trend toward polyploid evolution. The phylogenetic relationship between Z. armatum 'Tengjiao' and Z. armatum var. novemfolius was the closest, while that between Z. armatum 'YT-03' and Z. armatum 'YT-06' was closer than with Z. armatum 'Youkangtengjiao' according to oligo-FISH. The results provided a karyotype profile and a physical map that contributes to the distinction of varieties and cultivars of Z. armatum and provides strategies for distinguishing other cultivated species.

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.

Functional validation of ZbFAD2 and ZbFAD3 in the alkylamide biosynthesis pathway from Zanthoxylum bungeanum Maxim

The spicy taste and medicinal properties of Zanthoxylum bungeanum are imparted by several alkylamides. Although most studies have focused on their isolation and identification, few have reported their biosynthesis pathways. Among the differentially expressed genes (DEGs) reported in the numerous varieties of Z. bungeanum, some might contribute to alkylamide biosynthesis. However, they are not yet functionally validated. The present study explored the function of two genes, ZbFAD2 and ZbFAD3, in the alkylamide biosynthesis pathway, and their stable and transient expression in Arabidopsis thaliana and Nicotiana benthamiana were also analyzed. As compared with the wild-type (WT), the fatty acid content analysis indicated that ZbFAD2-A. thaliana transgenic seeds had lower oleic acid and higher linoleic acid contents, while the ZbFAD3-A. thaliana transgenic seeds showed lower linoleic acid and higher α-linolenic acid levels. Moreover, hydroxy-α-sanshool, a major alkylamide, was considerably higher in the ZbFAD2-N. benthamiana transgenic plants (0.2167 ± 0.0026 mg/g) than in the WT (0.0875 ± 0.0049 mg/g), while it was lower in the ZbFAD3-N. benthamiana transgenic plants (0.0535 ± 0.0037 mg/g). These results suggest that both ZbFAD2 and ZbFAD3 are vital alkylamide biosynthesis enzymes in Z. bungeanum. Our study not only helps to scale up the alkylamide production, but also establishes the role of the uncharacterized genes.

Plant Hormone Response to Low-Temperature Stress in Cold-Tolerant and Cold-Sensitive Varieties of Zanthoxylum bungeanum Maxim

Plant growth and survival in nature, its growth process, will be affected by various factors from the environment, among which temperature has a greater impact. In recent years, extreme weather has frequently appeared, and the growth of crops has been increasingly affected by the environment. As an important flavoring and Chinese herbal medicine crop, Zanthoxylum bungeanum is also facing the harm of low-temperature stress. Plant hormones play a vital role in the response of plants to low temperatures. In this study, ultra-performance liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to determine the hormone components of cold-tolerant and cold-sensitive varieties of Z. bungeanum. Combined with chemometric analysis and weighted gene co-expression network analysis (WGCNA), the hormone component differences and hormone response strategies of Z. bungeanum under low-temperature stress were comprehensively studied. The results showed that 45 hormones were detected in Z. bungeanum. Among them, there were 7 kinds of components with high content and were detected in both two varieties. At the late stage of low-temperature stress, the contents of abscisic acid (ABA) and ABA-glucosyl ester (ABA-GE) in Fuguhuajiao (FG) were significantly increased, and the latter served as the storage of the former to supplement the active ABA. Orthogonal partial least squares discriminant analysis (OPLS-DA) found that indole-3-carboxylic acid (ICA), indole-3-carboxaldehyde (ICAld), meta-Topolin riboside (mTR), cis-Zeatin-O-glucoside riboside (cZROG), and N6-isopentenyladenosine (IPR) in FG were the upregulated important difference components, and IPR and 2-methylthio-cis-zeatin riboside (2MeScZR) in Fengxiandahongpao (FX) were the upregulated important difference components. There were common crossing points and independent response pathways in response to low temperature in two varieties. WGCNA analysis found that the main hormone components were associated with multiple metabolic pathways including carbon, fatty acid, amino acid, and sugar metabolism, indicating that hormone regulation plays an important role in the response of Z. bungeanum to low temperature. This study clarified the hormone response mechanism of Z. bungeanum under low-temperature stress and provided a reference and basis for further improving the cold resistance of Z. bungeanum and cultivating new varieties.

Genome-Wide Identification of the NAC Gene Family in Zanthoxylum bungeanum and Their Transcriptional Responses to Drought Stress

NAC (NAM, ATAF1/2, and CUC2) transcription factors (TFs) are one of the largest plant-specific TF families and play a pivotal role in adaptation to abiotic stresses. The genome-wide analysis of NAC TFs is still absent in Zanthoxylum bungeanum. Here, 109 ZbNAC proteins were identified from the Z. bungeanum genome and were classified into four groups with Arabidopsis NAC proteins. The 109 ZbNAC genes were unevenly distributed on 46 chromosomes and included 4 tandem duplication events and 17 segmental duplication events. Synteny analysis of six species pairs revealed the closely phylogenetic relationship between Z. bungeanum and C. sinensis. Twenty-four types of cis-elements were identified in the ZbNAC promoters and were classified into three types: abiotic stress, plant growth and development, and response to phytohormones. Co-expression network analysis of the ZbNACs revealed 10 hub genes, and their expression levels were validated by real-time quantitative polymerase chain reaction (qRT-PCR). Finally, ZbNAC007, ZbNAC018, ZbNAC047, ZbNAC072, and ZbNAC079 were considered the pivotal NAC genes for drought tolerance in Z. bungeanum. This study represented the first genome-wide analysis of the NAC family in Z. bungeanum, improving our understanding of NAC proteins and providing useful information for molecular breeding of Z. bungeanum.

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.