WRKY genes family study reveals tissue-specific and stress-responsive TFs in wild potato species

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

Identification of WRKY gene family in Dioscorea opposita Thunb. reveals that DoWRKY71 enhanced the tolerance to cold and ABA stress

WRKY transcription factors constitute one of the largest plant-specific gene families, regulating various aspects of plant growth, development, physiological processes, and responses to abiotic stresses. This study aimed to comprehensively analyze the WRKY gene family of yam (Dioscorea opposita Thunb.), to understand their expression patterns during the growth and development process and their response to different treatments of yam and analyze the function of DoWRKY71 in detail. A total of 25 DoWRKY genes were identified from the transcriptome of yam, which were divided into six clades (I, IIa, IIc, IId, IIe, III) based on phylogenetic analysis. The analysis of conserved motifs revealed 10 motifs, varying in length from 16 to 50 amino acids. Based on real-time quantitative PCR (qRT-PCR) analysis, DoWRKY genes were expressed at different stages of growth and development and responded differentially to various abiotic stresses. The expression level of DoWRKY71 genes was up-regulated in the early stage and then down-regulated in tuber enlargement. This gene showed responsiveness to cold and abiotic stresses, such as abscisic acid (ABA) and methyl jasmonate (MeJA). Therefore, further study was conducted on this gene. Subcellular localization analysis revealed that the DoWRKY71 protein was localized in the nucleus. Moreover, the overexpression of DoWRKY71 enhanced the cold tolerance of transgenic tobacco and promoted ABA mediated stomatal closure. This study presents the first systematic analysis of the WRKY gene family in yam, offering new insights for studying WRKY transcription factors in yam. The functional study of DoWRKY71 lays theoretical foundation for further exploring the regulatory function of the DoWRKY71 gene in the growth and development related signaling pathway of yam.

Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

Comparative Transcriptome Analysis of High- and Low-Growth Genotypes of Eucalyptus urophylla in Response to Long-Term Nitrogen Deficiency

Nutrients play important roles in the growth and development of most plant species. However, in perennial trees, the function of nutrients in different genotypes is poorly understood. Three different nutrient levels (low, sufficient, and high nutrient levels) were applied to two contrasting Eucalyptus urophylla cultivars (a high-growth cultivar ZQUA44 and a low-growth cultivar ZQUB15), and growth and expression levels were analyzed. Although the growth traits of both genotypes under nutrient starvation treatment were much lower than under abundant nutrients, tree height, crown width, and biomass of different ZQUA44 tissues were much higher than those of ZQUB15 at all three nutrient levels. Differentially expressed genes (DEGs) clustered into six subclusters based on their expression patterns, and functional annotation showed that the DEGs involved in glutathione metabolism and flavonoid biosynthesis may be responsible for nutrient starvation across different genotypes, while the DEGs involved in carotenoid biosynthesis and starch and sucrose metabolism may have a range of functions in different genotypes. The DEGs encoding the MYB-related family may be responsible for nutrient deficiency in all genotypes, while B3 may have different functions in different genotypes. Our results demonstrate that different genotypes may form different pathways to coordinate plant survival when they face abiotic stresses.

WRKY transcription factors: evolution, regulation, and functional diversity in plants

The recent advancements in sequencing technologies and informatic tools promoted a paradigm shift to decipher the hidden biological mysteries and transformed the biological issues into digital data to express both qualitative and quantitative forms. The transcriptomic approach, in particular, has added new dimensions to the versatile essence of plant genomics through the large and deep transcripts generated in the process. This has enabled the mining of super families from the sequenced plants, both model and non-model, understanding their ancestry, diversity, and evolution. The elucidation of the crystal structure of the WRKY proteins and recent advancement in computational prediction through homology modeling and molecular dynamic simulation has provided an insight into the DNA-protein complex formation, stability, and interaction, thereby giving a new dimension in understanding the WRKY regulation. The present review summarizes the functional aspects of the high volume of sequence data of WRKY transcription factors studied from different species, till date. The review focuses on the dynamics of structural classification and lineage in light of the recent information. Additionally, a comparative analysis approach was incorporated to understand the functions of the identified WRKY transcription factors subjected to abiotic (heat, cold, salinity, senescence, dark, wounding, UV, and carbon starvation) stresses as revealed through various sets of studies on different plant species. The review will be instrumental in understanding the events of evolution and the importance of WRKY TFs under the threat of climate change, considering the new scientific evidences to propose a fresh perspective.

Genome-Wide Identification of Sweet Orange WRKY Transcription Factors and Analysis of Their Expression in Response to Infection by Penicillium digitatum

WRKY transcription factors (TFs) play a vital role in plant stress signal transduction and regulate the expression of various stress resistance genes. Sweet orange (Citrus sinensis) accounts for a large proportion of the world's citrus industry, which has high economic value, while Penicillium digitatum is a prime pathogenic causing postharvest rot of oranges. There are few reports on how CsWRKY TFs play their regulatory roles after P. digitatum infects the fruit. In this study, we performed genome-wide identification, classification, phylogenetic and conserved domain analysis of CsWRKY TFs, visualized the structure and chromosomal localization of the encoded genes, explored the expression pattern of each CsWRKY gene under P. digitatum stress by transcriptome data, and made the functional prediction of the related genes. This study provided insight into the characteristics of 47 CsWRKY TFs, which were divided into three subfamilies and eight subgroups. TFs coding genes were unevenly distributed on nine chromosomes. The visualized results of the intron-exon structure and domain are closely related to phylogeny, and widely distributed cis-regulatory elements on each gene played a global regulatory role in gene expression. The expansion of the CSWRKY TFs family was probably facilitated by twenty-one pairs of duplicated genes, and the results of Ka/Ks calculations indicated that this gene family was primarily subjected to purifying selection during evolution. Our transcriptome data showed that 95.7% of WRKY genes were involved in the transcriptional regulation of sweet orange in response to P. digitatum infection. We obtained 15 differentially expressed genes and used the reported function of AtWRKY genes as references. They may be involved in defense against P. digitatum and other pathogens, closely related to the stress responses during plant growth and development. Two interesting genes, CsWRKY2 and CsWRKY14, were expressed more than 60 times and could be used as excellent candidate genes in sweet orange genetic improvement. This study offers a theoretical basis for the response of CSWRKY TFs to P. digitatum infection and provides a vital reference for molecular breeding.

Transcriptome-wide identification of WRKY transcription factors and their expression profiles under different stress in Cynanchum thesioides

Cynanchum thesioides (Freyn) K. Schum. is an important economic and medicinal plant widely distributed in northern China. WRKY transcription factors (TFs) play important roles in plant growth, development and regulating responses. However, there is no report on the WRKY genes in Cynanchum thesioides. A total of 19 WRKY transcriptome sequences with complete ORFs were identified as WRKY transcriptome sequences by searching for WRKYs in RNA sequencing data. Then, the WRKY genes were classified by phylogenetic and conserved motif analysis of the WRKY family in Cynanchum thesioides and Arabidopsis thaliana. qRT-PCR was used to determine the expression patterns of 19 CtWRKY genes in different tissues and seedlings of Cynanchum thesioides under plant hormone (ABA and ETH) and abiotic stresses (cold and salt). The results showed that 19 CtWRKY genes could be divided into groups I-III according to their structure and phylogenetic characteristics, and group II could be divided into five subgroups. The prediction of CtWRKY gene protein interactions indicates that CtWRKY is involved in many biological processes. In addition, the CtWRKY gene was differentially expressed in different tissues and positively responded to abiotic stress and phytohormone treatment, among which CtWRKY9, CtWRKY18, and CtWRKY19 were significantly induced under various stresses. This study is the first to identify the WRKY gene family in Cynanchum thesioides, and the systematic analysis lays a foundation for further identification of the function of WRKY genes in Cynanchum thesioides.

Comprehensive Genomic Survey, Structural Classification, and Expression Analysis of WRKY Transcription Factor Family in Rhododendron simsii

(1) Rhododendron is one of the top ten traditional flowers in China, with both high ornamental and economic values. However, with the change of the environment, Rhododendron suffers from various biological stresses. The WRKY transcription factor is a member of the most crucial transcription factor families, which plays an essential regulatory role in a variety of physiological processes and developmental stresses. (2) In this study, 57 RsWRKYs were identified using genome data and found to be randomly distributed on 13 chromosomes. Based on gene structure and phylogenetic relationships, 57 proteins were divided into three groups: I, II, and III. Multiple alignments of RsWRKYs with Arabidopsis thaliana homologous genes revealed that WRKY domains in different groups had different conserved sites. RsWRKYs have a highly conserved domain, WRKYGQK, with three variants, WRKYGKK, WRKYGEK, and WRKYGRK. Furthermore, cis-acting elements analysis revealed that all of the RsWRKYs had stress and plant hormone cis-elements, with figures varying by group. Finally, the expression patterns of nine WRKY genes treated with gibberellin acid (GA), methyl jasmonate (MeJA), heat, and drought in Rhododendron were also measured using quantitative real-time PCR (qRT-PCR). The results showed that the expression levels of the majority of RsWRKY genes changed in response to multiple phytohormones and abiotic stressors. (3) This current study establishes a theoretical basis for future studies on the response of RsWRKY transcription factors to various hormone and abiotic stresses as well as a significant foundation for the breeding of new stress-tolerant Rhododendron varieties.

Uncovering anthocyanin diversity in potato landraces (Solanum tuberosum L. Phureja) using RNA-seq

Potato (Solanum tuberosum L.) is the third largest source of antioxidants in the human diet, after maize and tomato. Potato landraces have particularly diverse contents of antioxidant compounds such as anthocyanins. We used this diversity to study the evolutionary and genetic basis of anthocyanin pigmentation. Specifically, we analyzed the transcriptomes and anthocyanin content of tubers from 37 landraces with different colorations. We conducted analyses of differential expression between potatoes with different colorations and used weighted correlation network analysis to identify genes whose expression is correlated to anthocyanin content across landraces. A very significant fraction of the genes identified in these two analyses had annotations related to the flavonoid-anthocyanin biosynthetic pathway, including 18 enzymes and 5 transcription factors. Importantly, the causal genes at the D, P and R loci governing anthocyanin accumulation in potato cultivars also showed correlations to anthocyanin production in the landraces studied here. Furthermore, we found that 60% of the genes identified in our study were located within anthocyanin QTLs. Finally, we identified new candidate enzymes and transcription factors that could have driven the diversification of anthocyanins. Our results indicate that many anthocyanins biosynthetic genes were manipulated in ancestral potato breeding and can be used in future breeding programs.

WRKY transcription factors: a promising way to deal with arsenic stress in rice

Arsenic (As) is a global carcinogenic contaminant, and is one of the significant environmental constraints that limits the development and yield of crop plants. It is always tagged along with rice as rice takes up As and tends to accumulate it in grains. This amassment makes a way for As to get into the food chain that leads to unforeseen human health risks. Being viewed as parallel with toxicity, As in rice is an important global risk that calls for an urgent solution. WRKY Transcription Factors (TFs) seems to be promising in this area. The classical and substantial progress in the molecular mechanism of WRKY TFs, strengthened the understanding of innovative solutions for dealing with As in rice. Here, we review the potential of WRKY TFs under As stressed rice as a genetic solution and also provide insights into As and rice. Further, we develop an understanding of WRKY TF gene family and its regulation in rice. To date, studies on the role of WRKY TFs under As stressed rice are lacking. This area needs to be explored more so that this gene family can be utilized as an effective genetic tool that can break the As cycle to develop low or As free rice cultivar.

Transcriptome and metabolome profiling in different stages of infestation of Eucalyptus urophylla clones by Ralstonia solanacearum

Eucalyptus urophylla is an economically important tree species that widely planted in tropical and sub-tropical areas around the world, which suffers significant losses due to Ralstonia solanacearum. However, little is known about the molecular mechanism of pathogen-response of Eucalyptus. We collected the vascular tissues of a E. urophylla clone infected by R. solanacearum in the laboratory, and combined transcriptome and metabolome analysis to investigate the defense responses of Eucalyptus. A total of 11 flavonoids that differentially accumulated at the first stage or a later stage after infection. The phenylpropanoid of p-coumaraldehyde, the two alkaloids trigonelline and DL-ephedrine, two types of traditional Chinese medicine with patchouli alcohol and 3-dihydrocadambine, and the amino acid phenylalanine were differentially accumulated after infection, which could be biomarkers indicating a response to R. solanacearum. Differentially expressed genes involved in plant hormone signal transduction, phenylpropanoids, flavonoids, mitogen-activated protein kinase (MAPK) signaling, and amino acid metabolism were activated at the first stage of infection or a later stage, indicating that they may participate in the defense against infection. This study is expected to deliver several insights into the molecular mechanism in response to pathogens in E. urophylla, and the findings have far-reaching implications in the control of E. urophylla pathogens.

Whole-genome resequencing reveals genomic footprints of Italian sweet and hot pepper heirlooms giving insight into genes underlying key agronomic and qualitative traits

Background

Pepper is a major crop species of the Solanaceae family, largely appreciated for its high nutritional and healthy contribution to human diets. In the Mediterranean basin, the favorable pedoclimatic conditions enhanced the selection of several diversified landraces cultivated pepper (Capsicum annuum), for whom Italy can be considered a main pole of diversification. Hence, a survey of traditional C. annuum genetic resources is essential for deep understanding of such diversity and for applications in genomics assisted breeding. Here, we report whole-genome resequencing analyses of two sweet and two pungent genotypes highly diffused in South Italy and representative of the variability for shape, colour and nutritional properties.

Results

The four genomes were reconstructed at a chromosomal scale using a reference-guided approach, based on a dataset of 2.6 billion paired-end reads, corresponding to 20× genome coverage and a mapping rate above 99% for a final genomes size of approximately 3 Gb. After five iterations of variant calling, a total of 29,258,818 single nucleotide polymorphisms (SNPs) and 1,879,112 InDels, were identified. Substantial differences were observed among the four genomes based on geographical origin, with chromosomes 9 and 11 showing more polymorphisms in the accessions with higher fruit weight and absence of pungency. Among the identified variants, a small private indel (T - > TA) shared between sweet and big fruits accessions induces a frameshift with the generation of a new stop codon in a gene annotated as extensin, whereas two private SNPs within hot types were identified in 1-aminocyclopropane-1-carboxylate oxidase (ACO), a key gene involved in fruit ripening. The estimation of repetitive elements highlights a preponderant presence of Long Terminal Repeats (LTRs), the majority of which belonged to Gypsy superfamily. By comparing the four genomes with publicly available references including ‘CM334’ and Zunla-1 highlight the presence of 49,475 shared gene families.

Conclusions

The new genomic sequences aim to enrich the whole genome information of pepper local varieties, providing a valuable tool for precision gene mapping, marker discovery, comparative studies. Such knowledge widens the frontiers to understand the selection history of Italian pepper landraces toward the recognition of specificity local agri-food products marks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01039-9.

Approaches Involved in the Vegetable Crops Salt Stress Tolerance Improvement: Present Status and Way Ahead

Salt stress is one of the most important abiotic stresses as it persists throughout the plant life cycle. The productivity of crops is prominently affected by soil salinization due to faulty agricultural practices, increasing human activities, and natural processes. Approximately 10% of the total land area (950 Mha) and 50% of the total irrigated area (230 Mha) in the world are under salt stress. As a consequence, an annual loss of 12 billion US$ is estimated because of reduction in agriculture production inflicted by salt stress. The severity of salt stress will increase in the upcoming years with the increasing world population, and hence the forced use of poor-quality soil and irrigation water. Unfortunately, majority of the vegetable crops, such as bean, carrot, celery, eggplant, lettuce, muskmelon, okra, pea, pepper, potato, spinach, and tomato, have very low salinity threshold (ECt, which ranged from 1 to 2.5 dS m-1 in saturated soil). These crops used almost every part of the world and lakes' novel salt tolerance gene within their gene pool. Salt stress severely affects the yield and quality of these crops. To resolve this issue, novel genes governing salt tolerance under extreme salt stress were identified and transferred to the vegetable crops. The vegetable improvement for salt tolerance will require not only the yield influencing trait but also target those characters or traits that directly influence the salt stress to the crop developmental stage. Genetic engineering and grafting is the potential tool which can improve salt tolerance in vegetable crop regardless of species barriers. In the present review, an updated detail of the various physio-biochemical and molecular aspects involved in salt stress have been explored.

Pathogenesis-related protein 1 (PR-1) genes in soybean: Genome-wide identification, structural analysis and expression profiling under multiple biotic and abiotic stresses

Plant pathogenesis-related (PR) proteins are a large group of proteins, classified in 17 families, that are induced by pathological conditions. Here, we characterized the soybean PR-1 (GmPR-1) gene repertoire at the sequence, structural and expression levels. We found 24 GmPR-1 genes, clustered in two phylogenetic groups. GmPR-1 genes are under strong purifying selection, particularly those that emerged by tandem duplications. GmPR-1 promoter regions are abundant in cis-regulatory elements associated with major stress-related transcription factor families, namely WRKY, ERF, HD-Zip, C2H2, NAC, and GATA. We observed that 23 GmPR-1 genes are induced by stress conditions or exclusively expressed upon stress. We explored 1972 transcriptome samples, including 26 stress conditions, revealing that most GmPR-1 genes are differentially expressed in a plethora of biotic and abiotic stresses. Our findings highlight stress-responsive GmPR-1 genes with potential biotechnological applications, such as the development of transgenic lines with increased resistance to biotic and abiotic stresses.

WRKY Transcription Factor Response to High-Temperature Stress

Plant growth and development are closely related to the environment, and high-temperature stress is an important environmental factor that affects these processes. WRKY transcription factors (TFs) play important roles in plant responses to high-temperature stress. WRKY TFs can bind to the W-box cis-acting elements of target gene promoters, thereby regulating the expression of multiple types of target genes and participating in multiple signaling pathways in plants. A number of studies have shown the important biological functions and working mechanisms of WRKY TFs in plant responses to high temperature. However, there are few reviews that summarize the research progress on this topic. To fully understand the role of WRKY TFs in the response to high temperature, this paper reviews the structure and regulatory mechanism of WRKY TFs, as well as the related signaling pathways that regulate plant growth under high-temperature stress, which have been described in recent years, and this paper provides references for the further exploration of the molecular mechanisms underlying plant tolerance to high temperature.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Progressive Genomic Approaches to Explore Drought- and Salt-Induced Oxidative Stress Responses in Plants under Changing Climate

Drought and salinity are the major environmental abiotic stresses that negatively impact crop development and yield. To improve yields under abiotic stress conditions, drought- and salinity-tolerant crops are key to support world crop production and mitigate the demand of the growing world population. Nevertheless, plant responses to abiotic stresses are highly complex and controlled by networks of genetic and ecological factors that are the main targets of crop breeding programs. Several genomics strategies are employed to improve crop productivity under abiotic stress conditions, but traditional techniques are not sufficient to prevent stress-related losses in productivity. Within the last decade, modern genomics studies have advanced our capabilities of improving crop genetics, especially those traits relevant to abiotic stress management. This review provided updated and comprehensive knowledge concerning all possible combinations of advanced genomics tools and the gene regulatory network of reactive oxygen species homeostasis for the appropriate planning of future breeding programs, which will assist sustainable crop production under salinity and drought conditions.

The whole-genome and expression profile analysis of WRKY and RGAs in Dactylis glomerata showed that DG6C02319.1 and DgWRKYs may cooperate in the immunity against rust

Orchardgrass (Dactylis glomerata) is one of the top four perennial forages worldwide and, despite its large economic advantages, often threatened by various environmental stresses. WRKY transcription factors (TFs) can regulate a variety of plant processes, widely participate in plant responses to biotic and abiotic stresses, and are one of the largest gene families in plants. WRKYs can usually bind W-box elements specifically. In this study, we identified a total of 93 DgWRKY genes and 281 RGAs, including 65, 169 and 47 nucleotide-binding site-leucine-rich repeats (NBS-LRRs), leucine-rich repeats receptor-like protein kinases (LRR-RLKs), and leucine-rich repeats receptor-like proteins (LRR-RLPs), respectively. Through analyzing the expression of DgWRKY genes in orchardgrass under different environmental stresses, it was found that many DgWRKY genes were differentially expressed under heat, drought, submergence, and rust stress. In particular, it was found that the greatest number of genes were differentially expressed under rust infection. Consistently, GO and KEGG enrichment analysis of all genes showed that 78 DgWRKY TFs were identified in the plant–pathogen interaction pathway, with 59 of them differentially expressed. Through cis-acting element prediction, 154 RGAs were found to contain W-box elements. Among them, DG6C02319.1 (a member of the LRR-RLK family) was identified as likely to interact with 14 DGWRKYs. Moreover, their expression levels in susceptible plants after rust inoculation were first up-regulated and then down-regulated, while those in the resistant plants were always up-regulated. In general, DgWRKYs responded to both biotic stress and abiotic stress. DgWRKYs and RGAs may synergistically respond to the response of orchardgrass to rust. This study provides meaningful insight into the molecular mechanisms of WRKY proteins in orchardgrass.

Pathway-based analysis of anthocyanin diversity in diploid potato

Anthocyanin biosynthesis is one of the most studied pathways in plants due to the important ecological role played by these compounds and the potential health benefits of anthocyanin consumption. Given the interest in identifying new genetic factors underlying anthocyanin content we studied a diverse collection of diploid potatoes by combining a genome-wide association study and pathway-based analyses. By using an expanded SNP dataset, we identified candidate genes that had not been associated with anthocyanin variation in potatoes, namely a Myb transcription factor, a Leucoanthocyanidin dioxygenase gene and a vacuolar membrane protein. Importantly, a genomic region in chromosome 10 harbored the SNPs with strongest associations with anthocyanin content in GWAS. Some of these SNPs were associated with multiple anthocyanin compounds and therefore could underline the existence of pleiotropic genes or anthocyanin biosynthetic clusters. We identified multiple anthocyanin homologs in this genomic region, including four transcription factors and five enzymes that could be governing anthocyanin variation. For instance, a SNP linked to the phenylalanine ammonia-lyase gene, encoding the first enzyme in the phenylpropanoid biosynthetic pathway, was associated with all of the five anthocyanins measured. Finally, we combined a pathway analysis and GWAS of other agronomic traits to identify pathways related to anthocyanin biosynthesis in potatoes. We found that methionine metabolism and the production of sugars and hydroxycinnamic acids are genetically correlated to anthocyanin biosynthesis. The results contribute to the understanding of anthocyanins regulation in potatoes and can be used in future breeding programs focused on nutraceutical food.

The chromosome-level reference genome assembly for Panax notoginseng and insights into ginsenoside biosynthesis

Panax notoginseng, a perennial herb of the genus Panax in the family Araliaceae, has played an important role in clinical treatment in China for thousands of years because of its extensive pharmacological effects. Here, we report a high-quality reference genome of P. notoginseng, with a genome size up to 2.66 Gb and a contig N50 of 1.12 Mb, produced with third-generation PacBio sequencing technology. This is the first chromosome-level genome assembly for the genus Panax. Through genome evolution analysis, we explored phylogenetic and whole-genome duplication events and examined their impact on saponin biosynthesis. We performed a detailed transcriptional analysis of P. notoginseng and explored gene-level mechanisms that regulate the formation of characteristic tubercles. Next, we studied the biosynthesis and regulation of saponins at temporal and spatial levels. We combined multi-omics data to identify genes that encode key enzymes in the P. notoginseng terpenoid biosynthetic pathway. Finally, we identified five glycosyltransferase genes whose products catalyzed the formation of different ginsenosides in P. notoginseng. The genetic information obtained in this study provides a resource for further exploration of the growth characteristics, cultivation, breeding, and saponin biosynthesis of P. notoginseng.