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.

Physiological and transcriptome analyses reveal the photosynthetic response to drought stress in drought-sensitive (Fengjiao) and drought-tolerant (Hanjiao) Zanthoxylum bungeanum cultivars.

As an important economical plant, Zanthoxylum bungeanum is widely cultivated in arid and semi-arid areas. The studies associated with photosynthesis under drought stress were widely carried out, but not yet in Z. bungeanum. Here, the photosynthesis of two Z. bungeanum cultivars (FJ, Z. bungeanum cv. "Fengjiao"; HJ, Z. bungeanum cv. "Hanjiao") was analyzed under drought stress using physiological indicators and transcriptome data. Drought decreased stomatal aperture and stomatal conductance (Gsw), reduced transpiration rate (E) and sub-stomatal CO2 concentration (Ci), and lowered chlorophyll and carotenoid content, which reduced the net photosynthetic rate (Pn) of Z. bungeanum. The higher photosynthetic rate in HJ stemmed from its higher chlorophyll content, larger stomatal aperture and Gsw, and higher Ci. Weighted gene co-expression network analysis (WGCNA) identified several ABA signal transduction genes (PYL4, PYL9, and PYR1), LCH-encoding genes (LHCB4.3), and chlorophyll metabolism genes (CRD1, PORA, and CHLH). Additionally, seven transcription factor genes were identified as important factors regulating photosynthesis under drought conditions. In general, a photosynthetic response model under drought stress was built firstly in Z. bungeanum, and the key genes involved in photosynthesis under drought stress were identified. Therefore, the results in our research provide important information for photosynthesis under drought and provided key clues for future molecular breeding in Z. bungeanum.

Widely targeted metabolomic profiling combined with transcriptome analysis provides new insights into amino acid biosynthesis in green and red pepper fruits.

The type and content of amino acids in pepper are important indicators to reflect its nutritional value, largely affecting the purchasing behavior of consumers. Understanding the biosynthesis of amino acids in pepper fruit is beneficial to the development of pepper functional food. Widely targeted metabolomics, transcriptome analysis, correlation analysis, weighted gene co-expression network analysis (WGCNA), and canonical correlation analysis (CCA) were used to evaluate the quality characteristics of green and red pepper amino acids. The results showed 78 kinds of amino acids and their derivatives in the fruit of pepper. The essential amino acids were comprehensive and abundant. Especially, the contents of lysine and tryptophan were high. However, significant differences were found in the ratio of essential amino acids to total amino acids in green and red pepper. The ratio of essential amino acids to total amino acids in red pepper was up to 28.88%, while that in green pepper was up to 17.69%. WGCNA and CCA analyses were performed in combination with amino acid metabolism profiling and transcriptome analysis to further identify the main contributors to amino acid synthesis in green and red pepper. The results showed PK and PFK were the genes in the backbone of the amino acid biosynthesis pathway, which had a direct impact on the synthesis of various amino acids, and were the main genes for amino acid synthesis in pepper fruit. In this study, the amino acid biosynthesis rules for two kinds of pepper were analyzed by amino acid metabolism profiling and transcriptome analysis, which provided the basis for the development of amino acid nutritional supplements and pepper functional food.

UPLC-MS/MS Profile Combined With RNA-Seq Reveals the Amino Acid Metabolism in Zanthoxylum bungeanum Leaves Under Drought Stress.

Zanthoxylum bungeanum leaves have a unique taste and incomparable nutritional value and hence are popular as a food item and traditional medicine in China. However, the studies on the metabolites in Z. bungeanum leaves are quite limited, especially for amino acids. Therefore, this study explored the amino acid component in Z. bungeanum leaves and also the accumulation mechanism under drought stress in two Z. bungeanum cultivars using the widely targeted metabolome combined with transcriptome analysis. A total of 56 amino acids and their derivatives were identified in Z. bungeanum leaves, including eight essential amino acids. The total amino acid content with most individual amino acids increased under progressive drought stress. More differentially accumulated amino acids (DAAs) and differentially expressed genes (DEGs) were found in FJ (Z. bungeanum cv. 'Fengjiao') than in HJ (Z. bungeanum cv. 'Hanjiao'). The orthogonal projections to latent structures discriminant analysis identified nine and seven indicator DAAs in FJ and HJ leaves, respectively. The weighted gene co-expression network analysis (WGCNA) showed that the green module was significantly correlated with most indicator DAAs and revealed the important role of FBA3, DELTA-OAT, PROC, and 15 transcription factor genes in regulating the amino acid synthesis. Furthermore, the correlation analysis and redundancy analysis (RDA) identified four candidate synthesis genes (ASNS, AK, ASPS, and PK) in amino acid biosynthesis pathway. This study provided useful information for the development of Z. bungeanum leaves in food and nutrition industry and also laid the foundations for future molecular breeding.

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.

Integrated LC-MS/MS and Transcriptome Sequencing Analysis Reveals the Mechanism of Color Formation During Prickly Ash Fruit Ripening.

Prickly ash peel is one of the eight major condiments in China and is widely used in cooking because of its unique fragrance and numbing taste. The color of prickly ash fruit is the most intuitive quality that affects consumer choice. However, the main components and key biosynthetic genes responsible for prickly ash fruit color have not yet been determined. To better understand the biosynthetic mechanisms and accumulation of prickly ash fruit color components, we performed an integrated transcriptomic and metabolomic analysis of red and green prickly ash fruit at different growth periods. The transcriptome analysis identified 17,269 differentially expressed genes (DEGs) between fruit of red and green prickly ash: 7,236 upregulated in green fruit and 10,033 downregulated. Liquid chromatography tandem mass spectrometry (LC-MS/MS) identified 214 flavonoids of 10 types. Flavonoids and flavonols are the main flavonoids in prickly ash, and the total flavonoid content of red prickly ash is higher than that of green prickly ash. Comprehensive analysis showed that the main colored metabolites that differed between green and red prickly ash were cyanidin-3-O-galactoside and cyanidin-3-O-glucoside, and differences in the contents of these metabolites were due mainly to differences in the expression of ANS and UFGT. Our results provide insight into the mechanisms underlying color differences in red and green prickly ash and will be useful for improving the quality of prickly ash fruit.

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.

ZbAGL11, a class D MADS-box transcription factor of Zanthoxylum bungeanum, is involved in sporophytic apomixis.

Apomixis is a reproductive model that bypasses sexual reproduction, so it does not require the combination of paternal and maternal gametes but instead results in the production of offspring directly from maternal tissues. This reproductive mode results in the same genetic material in the mother and the offspring and has significant applications in agricultural breeding. Molecular and cytological methods were used to identify the reproductive type of Zanthoxylum bungeanum (ZB). Fluorescence detection of the amplified products of 12 pairs of polymorphic SSR primers showed consistent fluorescence signals for mother and offspring, indicating that no trait separation occurred during reproduction. In addition, the cytological observation results showed differentiation of ZB embryos (2n) from nucellar cells (2n) to form indefinite embryonic primordia and then form adventitious embryos (2n), indicating that the apomictic type of ZB is sporophytic apomixis. The MADS-box transcription factor ZbAGL11 was highly expressed during the critical period of nucellar embryo development in ZB. Unpollinated ZbAGL11-OE Arabidopsis produced fertile offspring and exhibited an apomictic phenotype. The overexpression of ZbAGL11 increased the callus induction rate of ZB tissue. In addition, the results of the yeast two-hybrid experiment showed that ZbAGL11 could interact with the ZbCYP450 and ZbCAD11 proteins. Our results demonstrate that ZbAGL11 can cause developmental disorders of Arabidopsis flower organs and result in apomixis-like phenotypes.

Transcriptome and Metabolome Dynamics Explain Aroma Differences between Green and Red Prickly Ash Fruit.

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.

Pollination promotes ABA synthesis but not sexual reproduction in the apomictic species Zanthoxylum bungeanum Maxim.

Apomixis is a form of reproduction that does not involve the fertilization of female gametes by male gametes but instead involves the production of offspring directly from the female parent. The offspring of apomixis are genetically identical to the female parent and inherit its traits. Therefore, apomixis has great potential for application to agricultural genetic breeding. However, it remains unclear whether apomictic species require pollination, and the impacts of pollination on such species are poorly understood. We investigated the effects of pollination on the apomictic species Zanthoxylum bungeanum Maxim. by analyzing its fertilization process, assembling its transcriptome, and measuring hormone concentrations, fruit setting rate and gene expression levels. Transcriptome sequencing of pollinated and unpollinated fruits resulted in a total of 69,131 PacBio reads. Of these, 7102 genes were up-regulated and 6491 genes were down-regulated. Analysis of the differentially expressed genes (DEGs) and construction of a weighted gene co-expression network showed that many DEGs were involved in plant hormone signal transduction, suggesting that hormonal signaling during development differs between pollinated and unpollinated fruit. The germination rate of Z. bungeanum pollen in vitro was only 11%, and pollen could not germinate in the embryo sac to complete fertilization. Although pollination did not enable Z. bungeanum to complete the sexual reproduction process, it significantly increased abscisic acid (ABA) concentration and fruit setting rate. Spraying 100 µg l-1 ABA also significantly increased the fruit setting rate. Therefore, ABA appears to be a key factor in the regulation of fruit setting in apomictic Z. bungeanum. Based on these results, we suggest that some male plants be cultivated in Z. bungeanum plantations or exogenous ABA be sprayed to increase the likelihood of pollination and thereby increase the fruit setting rate.

Integrated Analysis of Metabolome and Transcriptome Data for Uncovering Flavonoid Components of Zanthoxylum bungeanum Maxim. Leaves Under Drought Stress.

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.

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.

Effects of location, climate, soil conditions and plant species on levels of potentially toxic elements in Chinese Prickly Ash pericarps from the main cultivation regions in China.

The presence of potentially toxic mineral elements is often reported in soil, fruits and vegetables. The concentrations of these are influenced by location, climate and soil conditions, and plant species. This study reports levels of aluminium (Al), arsenic (As), cadmium (Cd), lead (Pb), manganese (Mn) and nickel (Ni) in the soils where Chinese Prickly Ash - Zanthoxylum bungeanum (ZB) and Z. armatum (ZA) are grown. Soil and pericarp samples were taken from 72 ZB and ZA plantations in 12 provinces. The soil pollution index (PI) was generally low for As (0.04-0.35), Pb (0.01-0.74) and Ni (0.15-1.06). The PI was higher for Cd (0.06-6.99) and was assessed at three categorical levels: 'low' (in 26% of soils), 'moderate' (in 67% of soils) and high (in 7% of soils). The majority (94%) of soils were slightly polluted with Cd with the integrated pollution index in the range 0.17-2.02. More than 60% of pericarp samples posed a high non-cancer health risk based on the hazard index, while 43% of pericarp samples posed a considerable cancer risk based on the incremental lifetime cancer risk. The mean temperature and the mean relative humidity contributed most significantly (28% and 10%, respectively) to variation in the concentrations of potentially toxic elements in the pericarps. This study provides a comprehensive report on the concentrations of potentially toxic elements in the soils and associated pericarps, and the effects of location, climate, and soil on the levels of potentially toxic elements in the soil and pericarps, including the effects of tree species.

The steps from sexual reproduction to apomixis.

MAIN CONCLUSION: In this paper, an interaction model of apomixis-related genes was constructed to analyze the emergence of apomictic types. It is speculated that apomixis technology will be first implemented in gramineous plants. Apomixis (asexual seed formation) is a phenomenon in which a plant bypasses the most fundamental aspects of sexual reproduction-meiosis and fertilization-to form a viable seed. Plants can form seeds without fertilization, and the seed genotype is consistent with the female parent. The development of apomictic technology would be revolutionary for agriculture and for food production as it would reduce costs and breeding times and also avoid many complications typical of sexual reproduction (e.g. incompatibility barriers) and of vegetative propagation (e.g. viral transfer). The application of apomictic reproductive technology has the potential to revolutionize crop breeding. This article reviews recent advances in apomixis in cytology and molecular biology. The general idea of identifying apomixis was proposed and the process of the emergence of non-fusion types was discussed. To better understand the apomixis mechanism, an apomixis regulatory model was established. At the same time, the realization of apomixis technology is proposed, which provides reference for the research and application of apomixis.

Patterns of Drought Response of 38 WRKY Transcription Factors of Zanthoxylum bungeanum Maxim.

The WRKY family of transcription factors (TFs) includes a number of transcription-specific groupings that play important roles in plant growth and development and in plant responses to various stresses. To screen for WRKY transcription factors associated with drought stress in Zanthoxylum bungeanum, a total of 38 ZbWRKY were identified and these were then classified and identified with Arabidopsis WRKY. Using bioinformatics analyses based on the structural characteristics of the conservative domain, 38 WRKY transcription factors were identified and categorized into three groups: Groups I, II, and III. Of these, Group II can be divided into four subgroups: subgroups IIb, IIc, IId, and IIe. No ZbWRKY members of subgroup IIa were found in the sequencing data. In addition, 38 ZbWRKY were identified by real-time PCR to determine the behavior of this family of genes under drought stress. Twelve ZbWRKY transcription factors were found to be significantly upregulated under drought stress and these were identified by relative quantification. As predicted by the STRING website, the results show that the WRKYs are involved in four signaling pathways-the jasmonic acid (JA), the salicylic acid (SA), the mitogen-activated protein kinase (MAPK), and the ethylene signaling pathways. ZbWRKY33 is the most intense transcription factor in response to drought stress. We predict that WRKY33 binds directly to the ethylene synthesis precursor gene ACS6, to promote ethylene synthesis. Ethylene then binds to the ethylene activator release signal to activate a series of downstream genes for cold stress and osmotic responses. The roles of ZbWRKY transcription factors in drought stress rely on a regulatory network center on the JA signaling pathway.

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.