Genomic analysis reveals the genetic diversity, population structure, evolutionary history and relationships of Chinese pepper

Analysis of population structure and genetic diversity of Camellia tachangensis in Guizhou based on SNP markers

BACKGROUND

Camellia tachangensis F. C. Zhang is a five-compartment species in the ovary of tea group plants, which represents the original germline of early differentiation of some tea group plants.

METHODS AND RESULTS

In this study, we analyzed single-nucleotide polymorphisms (SNPs) at the genome level, constructed a phylogenetic tree, analyzed the genetic diversity, and further investigated the population structure of 100 C. tachangensis accessions using the genotyping-by-sequencing (GBS) method. A total of 91,959 high-quality SNPs were obtained. Population structure analysis showed that the 100 C. tachangensis accessions clustered into three groups: YQ-1 (Village Group), YQ-2 (Forest Group) and YQ-3 (Transition Group), which was further consistent with the results of phylogenetic analysis and principal component analyses (PCA). In addition, a comparative analysis of the genetic diversity among the three populations (Forest, Village, and Transition Groups) detected the highest genetic diversity in the Transition Group and the highest differentiation between Forest and Village Groups.

CONCLUSIONS

C. tachangensis plants growing in the forest had different genetic backgrounds from those growing in villages. This study provides a basis for the effective protection and utilization of C. tachangensis populations and lays a foundation for future C. tachangensis breeding.

Leaf traits of prickly ash and its correlation with ecological and geographical factors of origin

The morphological, physiological, and biochemical characteristics of leaves result from the long-term adaptation of plants to their environment and are closely related to plant growth and development. In this study, 37 prickly ash germplasm resources from 18 production areas were utilized as the subjects of research. Logistic equations, principal component analysis, and cluster analysis were employed to comprehensively evaluate the leaf traits of prickly ash germplasm resources, with an analysis of their correlation with ecological and geographical factors in the production areas. The results showed that the leaf traits of prickly ash germplasms of different origins are substantially different and diverse. The coefficient of variation for the 14 leaf traits was greater than 10%. The coefficient of variation of the compound leaflet number was the highest among all the considered leaf traits, and the coefficient of variation of leaf thickness was the lowest, at 49.86% and 11.37%, respectively. The leaf traits of the prickly ash germplasm originating from Chongqing in Yongchuan, Chongqing in Rongchang, and Yunnan in Honghe ranked highest, whereas the leaf traits of the prickly ash germplasm from Henan in Jiaozuo, Gansu in Tianshui, and Shanxi in Yuncheng ranked lowest. The results of the correlation analysis showed that among the ecological and geographical factors of the origins, latitude had the strongest correlation with the leaf traits of the prickly ash germplasm. As latitude increased, the leaves of prickly ash gradually decreased in size, weight, and leaf shape index. The factor with the second strongest correlation was temperature. The leaves of the prickly ash germplasm originating from warmer climate areas were larger and heavier than those from areas with colder climates. Altitude and longitude did not significantly affect the leaf traits of the prickly ash germplasm, but at similar latitudes, the leaves of the prickly ash germplasm in high-altitude areas were smaller, and the leaves of the prickly ash germplasm in low-altitude areas were larger. These findings can provide valuable references for breeding and the sustainable utilization of new varieties of prickly ash resources.

Optimized Pepper Target SNP-Seq Applied in Population Structure and Genetic Diversity Analysis of 496 Pepper (Capsicum spp.) Lines

Peppers are a major vegetable crop worldwide. With the completion of additional genome assemblies, a multitude of single-nucleotide polymorphisms (SNPs) can be utilized for population structure and genetic diversity analysis. In this study, we used target SNP-sequencing as a new high-throughput sequencing technology, screening out 425 perfect SNPs for analyzing the genetic diversity and population structure among 496 pepper lines from five pepper species in China and abroad. The perfect SNP panel exhibited commendable discriminative ability, as indicated by the average values of polymorphism information content, observed heterozygosity, minor allele frequency, and genetic diversity, which were 0.346, 0.011, 0.371, and 0.449, respectively. Based on phylogenetic, population structure, and principal component analyses, 484 C. annuum lines were divided into four subpopulations according to the shape of fruit: blocky fruit, wide-horn fruit, narrow-horn fruit, and linear fruit. These subpopulations displayed clear clustering with minimal or no overlap. Moreover, F statistic (Fst) analysis revealed considerable distinctions among these subpopulations. Additionally, we established a set of 47 core SNPs that could effectively differentiate among all pepper lines. This core SNP set could precisely classify the C. annuum lines into four distinct fruit-shape groups. The blocky and narrow-horn fruit subpopulations displayed the lowest and highest genetic diversity, respectively. This study highlights the importance of fruit shape as a crucial trait in pepper breeding. Moreover, this work indicates the immense potential of optimized target SNP technology in the addition of foreground markers of important traits to improve molecular breeding efficiency, and demonstrates its broad application prospects in the genetic analysis and variety identification of peppers.

Six Express Sequence Tag-Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species

Zanthoxylum (Sichuan pepper), with its rich cultivars, has long been widely cultivated in China for its unique seasoning and medicinal uses, but most of its cultivars have similar morphological characteristics. Therefore, we hypothesized that the genetic diversity of Zanthoxylum cultivars is low because of their apomixis and long cultivation history. In this study, we aimed to investigate the genetic diversity of three Zanthoxylum species on the cultivar level based on express sequence tag-simple sequence repeat (EST-SSR) primers. In total, 121 samples of three Zanthoxylum species (Z. bungeanum, Z. armatum and Z. piperitum) were collected from different areas in China for genetic diversity analysis. A total of six specificity and polymorphism EST-SSR primers, which we selected from among 120 primers based on two transcriptomes (Z. bungeanum, Z. armatum) in our earlier study, were used to evaluate genetic diversity based on capillary electrophoresis technology. The results of our analysis using the unweighted pair group method with arithmetic mean (UPGMA) indicated that most of the samples are clustered in one clade in the UPGMA dendrogram, and the average genetic distance was 0.6409. Principal component analysis (PCA) showed that Z. piperitum may have a closer genetic relationship with Z. bungeanum than with Z. armatum. An analysis of molecular variation (AMOVA) showed that the genetic variation mainly stemmed from individuals within populations; the genetic differentiation coefficient (PhiPT) was 0.429, the gene flow (Nm) between populations was 0.333, and the differences among populations were not significant (p > 0.001). For the intraspecific populations of ZB, the percentage of genetic variation was 53% among populations and 47% within populations, with non-significant differences between populations (p > 0.001). The genetic differentiation coefficient (PhiT) was 0.529, and the gene flow (Nm) was 0.223. For the intraspecific populations of ZA, the results indicated that the percentage of genetic variation was 29% among populations and 71% within populations, with non-significant differences between populations (p > 0.001); the genetic differentiation coefficient (PhiPT) was 0.293, and the gene flow (Nm) was 0.223. Through genetic structure analysis (GSA), we predicted that these 121 samples belonged to two optimal subgroups, which means that all the samples probably originated from two gene pools. Above all, this indicated that the genetic diversity of the 121 Zanthoxylum samples was relatively low at both the species and cultivar levels, a finding which was consistent with our initial assumptions. This study provides a reference, with molecular-level data, for the further identification of Zanthoxylum species.

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.

α-Linolenic Acid Inhibits RANKL-Induced Osteoclastogenesis In Vitro and Prevents Inflammation In Vivo

Inflammation is an important risk factor for bone-destroying diseases. Our preliminary research found that Zanthoxylum bungeanum seed oil (ZBSO) is abundant in unsaturated fatty acids and could inhibit osteoclastogenesis in receptor activator of nuclear factor κB ligand (RANKL)-induced RAW264.7 cells. However, the key constituents in ZBSO in the prevention of osteoclastogenesis and its possible mechanism related to inflammation are still unclear. Therefore, in this study, oleic acid (OA), linoleic acid (LA), palmitoleic acid (PLA), and alpha-linolenic acid (ALA) in ZBSO, havingthe strongest effect on RANKL-induced osteoclastogenesis, were selected by a tartrate-resistant acid phosphatase (TRAP) staining method. Furthermore, the effects of the selected fatty acids on anti-inflammation and anti-osteoclastogenesis in vitro and in vivo were assessed using RT-qPCR. Among the four major unsaturated fatty acids we tested, ALA displayed the strongest inhibitory effect on osteoclastogenesis. The increased expression of free fatty acid receptor 4 (FFAR4) and β-arrestin2 (βarr2), as well as the decreased expression of nuclear factor-kappa B (NF-κB), tumor necrosis factor-α (TNF-α), nuclear factor of activated T-cells c1 (NFATc1), and tartrate-resistant acid phosphatase (TRAP) in RAW264.7 cells after ALA treatment were observed. Moreover, in ovariectomized osteoporotic rats with ALA preventive intervention, we found that the expression of TNF-α, interleukin-6 (IL-6), interleukin-1β (IL-1β), NFATc1, and TRAP were decreased, while with the ALA therapeutic intervention, downregulated expression of NF-κB, NFATc1, TRAP, and transforming growth factor beta-activated kinase 1 (TAK1) were noticed. These results indicate that ALA, as the major unsaturated fatty acid in ZBSO, could inhibit RANKL-induced osteoclastogenesis via the FFAR4/βarr2 signaling pathway and could prevent inflammation, suggesting that ZBSO may be a promising potential natural product of unsaturated fatty acids and a dietary supplement for the prevention of osteoclastogenesis and inflammatory diseases.

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.

Genome-wide characterization of the MBF1 gene family and its expression pattern in different tissues and stresses in Zanthoxylum armatum

Background

Multiprotein bridging factor 1 (MBF1) is a crucial transcriptional coactivator in animals, plants, and some microorganisms, that plays a necessary role in growth development and stress tolerance. Zanthoxylum armatum is an important perennial plant for the condiments and pharmaceutical industries, whereas the potential information in the genes related to stress resistance remains poorly understood in Z. armatum. 

Results

Herein, six representative species were selected for use in a genome-wide investigation of the MBF1 family, including Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Citrus sinensis, Ginkgo biloba, and Z. armatum. The results showed that the MBF1 genes could be divided into two groups: Group I contained the MBF1a and MBF1b subfamilies, and group II was independent of the MBF1c subfamily.. Most species have at least two different MBF1 genes, and MBF1c is usually an essential member. The three ZaMBF1 genes were respectively located on ZaChr26, ZaChr32, and ZaChr4 of Zanthoxylum chromosomes. The collinearity were occurred between three ZaMBF1 genes, and ZaMBF1c showed the collinearity between Z. armatum and both P. trichocarpa and C. sinensis. Moreover, many cis-elements associated with abiotic stress and phytohormone pathways were detected in the promoter regions of MBF1 of six representative species. The ERF binding sites were the most abundant targets in the sequences of the ZaMBF1 family, and some transcription factor sites related to floral differentiation were also identified in ZaMBF1c, such as MADS, LFY, Dof, and AP2. ZaMBF1a was observed to be very highly expressed in 25 different samples except in the seeds, and ZaMBF1c may be associated with the male and female floral initiation processes. In addition, expression in all the ZaMBF1 genes could be significantly induced by water-logging, cold stress, ethephon, methyl jasmonate, and salicylic acid treatments, especially in ZaMBF1c.

Conclusion

The present study carried out a comprehensive bioinformatic investigation related to the MBF1 family in six representative species, and the responsiveness of ZaMBF1 genes to various abiotic stresses and phytohormone inductions was also revealed. This work not only lays a solid foundation to uncover the biological roles of the ZaMBF1 family in Z. armatum, but also provides some broad references for conducting the MBF1 research in other plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08863-4.

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 CO₂ 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.

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.

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.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.

Apomixis: genetic basis and controlling genes

Apomixis is the phenomenon of clonal reproduction by seed. As apomixis can produce clonal progeny with exactly the same genotype as the maternal plant, it has an important application in genotype fixation and accelerating agricultural breeding strategies. The introduction of apomixis to major crops would bring many benefits to agriculture, including permanent fixation of superior genotypes and simplifying the procedures of hybrid seed production, as well as purification and rejuvenation of crops propagated vegetatively. Although apomixis naturally occurs in more than 400 plant species, it is rare among the major crops. Currently, with better understanding of apomixis, some achievements have been made in synthetic apomixis. However, due to prevailing limitations, there is still a long way to go to achieve large-scale application of apomixis to crop breeding. Here, we compare the developmental features of apomixis and sexual plant reproduction and review the recent identification of apomixis genes, transposons, epigenetic regulation, and genetic events leading to apomixis. We also summarize the possible strategies and potential genes for engineering apomixis into crop plants.

Volatile Oil Profile of Prickly Ash (Zanthoxylum) Pericarps from Different Locations in China

Volatile oils of prickly ash (Zanthoxylum) pericarps have various potential biological functions with considerable relevance to food, pharmacological, and industrial applications. The volatile profile of oils extracted from prickly ash pericarps obtained from 72 plantations in China was determined by gas chromatography and mass spectrometry. Several chemometric analyses were used to better understand the volatile oil profile differences among different pericarps and to determine the key factors that affected geographical variations in the main volatile constituents of oils. A total of 47 constituents were detected with D-limonene, alfa-myrcene, and linalool as the most abundant. The volatile profile of pericarp oils was significantly affected by prickly ash species and some environmental factors, and the key factors that affected volatile profile variations for different prickly ash species were diverse. Chemometric analyses based on the volatile oil profile could properly distinguish Z. armatum pericarps from other pericarps. This study provides comprehensive information on the volatile oil profile of pericarps from different prickly ash species and different plantations, and it can be beneficial to a system for evaluating of pericarp quality. Moreover, this study speculates on the key environmental factors that cause volatile oil variations for each species, and can help to obtain better prickly ash pericarp volatile oils by improving the cultivated environments.

Zanthoxylum Species: A Review of Traditional Uses, Phytochemistry and Pharmacology in Relation to Cancer, Infectious Diseases and Sickle Cell Anemia

The health benefits and toxicity of plant products are largely dependent on their secondary metabolite contents. These compounds are biosynthesized by plants as protection mechanisms against environmental factors and infectious agents. This review discusses the traditional uses, phytochemical constituents and health benefits of plant species in genus Zanthoxylum with a focus on cancer, microbial and parasitic infections, and sickle cell disease as reported in articles published from 1970 to 2021 in peer-reviewed journals and indexed in major scientific databases. Generally, Z. species are widely distributed in Asia, America and Africa, where they are used as food and for disease treatment. Several compounds belonging to alkaloids, flavonoids, terpenoids, and lignans, among others have been isolated from Z. species. This review discusses the biological activities reported for the plant species and their phytochemicals, including anticancer, antibacterial, antifungal, antiviral, anti-trypanosomal, antimalarial and anti-sickling properties. The safety profiles and suggestions for conservation of the Z. species were also discussed. Taken together, this review demonstrates that Z. species are rich in a wide range of bioactive phytochemicals with multiple health benefits, but more research is needed towards their practical application in the development of functional foods, nutraceuticals and lead compounds for new drugs.

Chromosome-level genome assembly of Sichuan pepper provides insights into apomixis, drought tolerance, and alkaloid biosynthesis

Sichuan pepper is a commonly used spice in Asian cuisine. Sanshools and wgx-50/gx-50 isolated from it have been shown to possess a wide spectrum of medicinal properties. Here we generated a chromosome-level genome assembly of one Sichuan pepper species Zanthoxylum armatum characterized by drought tolerance and apomixis. Analyses of functionally related genes suggested that increased gene copy number and expression level of drought-tolerant genes might play an important role in improving drought tolerance of Z. armatum. Moreover, a gene encoding an RWP-RK domain-containing protein was shown to contribute to apomixis in Z. armatum, which was further characterized by overexpression in Arabidopsis thaliana. Furthermore, based on gene homology searching and co-expression patterns of metabolite concentration and gene expressions, we identified a number of candidate genes involved in the biosynthesis of sanshools and wgx-50/gx-50. Taken together, our results yield valuable insights for understanding the evolution of apomixis, drought tolerance, and alkaloid biosynthesis in Z. armatum.

Alkylamide Profiling of Pericarps Coupled with Chemometric Analysis to Distinguish Prickly Ash Pericarps

Because of extensive cultivation areas, various cultivars, nonstandard naming notations, and morphology similarity among relative cultivars, adulteration and associated business fraud may happen in the marketplaces of prickly ash pericarps due to higher financial gain and high-frequency trading. This study presents variations in the chemical components and contents of different prickly ash species from different plantations. Alkylamide profiling of pericarps derived from Zanthoxylum armatum, Z. bungeanum, and some relative Zanthoxylum species from 72 plantations across China were tested using ultra-performance liquid chromatography. Then, several chemometrics were applied to classify the prickly ash pericarps to reveal potential indicators that distinguish prickly ash pericarps and to identify the key factors that affect pericarp alkylamide profiling. The dominating alkylamides in the prickly ash pericarps were Z. piperitum (ZP)-amide C (0–20.64 mg/g) and ZP-amide D (0–30.43 mg/g). Alkylamide profiling of prickly ash pericarps varied significantly across species and geographical variations. ZP-amide D in prickly ash pericarps was identified as a potential indicator to distinguish prickly ash species. Longitude and aluminum content in soils were identified as key factors that affected alkylamide profiling of prickly ash pericarps. This study provides a useful tool to classify prickly ash species based on pericarp alkylamide profiling and to determine the key influence factors on pericarp alkylamide variations.