MPK3 and MPK6 control salicylic acid signaling by up-regulating NLR receptors during pattern- and effector-triggered immunity

Genome-wide identification and analysis of abiotic stress responsiveness of the mitogen-activated protein kinase gene family in Medicago sativa L

BACKGROUND

The mitogen-activated protein kinase (MAPK) cascade is crucial cell signal transduction mechanism that plays an important role in plant growth and development, metabolism, and stress responses. The MAPK cascade includes three protein kinases, MAPK, MAPKK, and MAPKKK. The three protein kinases mediate signaling to downstream response molecules by sequential phosphorylation. The MAPK gene family has been identified and analyzed in many plants, however it has not been investigated in alfalfa.

RESULTS

In this study, Medicago sativa MAPK genes (referred to as MsMAPKs) were identified in the tetraploid alfalfa genome. Eighty MsMAPKs were divided into four groups, with eight in group A, 21 in group B, 21 in group C and 30 in group D. Analysis of the basic structures of the MsMAPKs revealed presence of a conserved TXY motif. Groups A, B and C contained a TEY motif, while group D contained a TDY motif. RNA-seq analysis revealed tissue-specificity of two MsMAPKs and tissue-wide expression of 35 MsMAPKs. Further analysis identified MsMAPK members responsive to drought, salt, and cold stress conditions. Two MsMAPKs (MsMAPK70 and MsMAPK75) responds to salt and cold stresses; two MsMAPKs (MsMAPK60 and MsMAPK73) responds to cold and drought stresses; four MsMAPKs (MsMAPK1, MsMAPK33, MsMAPK64 and MsMAPK71) responds to salt and drought stresses; and two MsMAPKs (MsMAPK5 and MsMAPK7) responded to all three stresses.

CONCLUSION

This study comprehensively identified and analysed the alfalfa MAPK gene family. Candidate genes related to abiotic stresses were screened by analysing the RNA-seq data. The results provide key information for further analysis of alfalfa MAPK gene functions and improvement of stress tolerance.

A signaling cascade mediating fruit trait development via phosphorylation-modulated nuclear accumulation of JAZ repressor

It is generally accepted that jasmonate-ZIM domain (JAZ) repressors act to mediate jasmonate (JA) signaling via CORONATINE-INSENSITIVE1 (COI1)-mediated degradation. Here, we report a cryptic signaling cascade where a JAZ repressor, FvJAZ12, mediates multiple signaling inputs via phosphorylation-modulated subcellular translocation rather than the COI1-mediated degradation mechanism in strawberry (Fragaria vesca). FvJAZ12 acts to regulate flavor metabolism and defense response, and was found to be the target of FvMPK6, a mitogen-activated protein kinase that is capable of responding to multiple signal stimuli. FvMPK6 phosphorylates FvJAZ12 at the amino acid residues S179 and T183 adjacent to the PY residues, thereby attenuating its nuclear accumulation and relieving its repression for FvMYC2, which acts to control the expression of lipoxygenase 3 (FvLOX3), an important gene involved in JA biosynthesis and a diverse array of cellular metabolisms. Our data reveal a previously unreported mechanism for JA signaling and decipher a signaling cascade that links multiple signaling inputs with fruit trait development.

Integrated multi-omics and genetic analyses reveal molecular determinants underlying Arabidopsis snap33 mutant phenotype

The secretory pathway is essential for plant immunity, delivering diverse antimicrobial molecules into the extracellular space. Arabidopsis thaliana soluble N-ethylmaleimide-sensitive-factor attachment protein receptor SNAP33 is a key actor of this process. The snap33 mutant displays dwarfism and necrotic lesions, however the molecular determinants of its macroscopic phenotypes remain elusive. Here, we isolated several new snap33 mutants that exhibited constitutive cell death and H2O2 accumulation, further defining snap33 as an autoimmune mutant. We then carried out quantitative transcriptomic and proteomic analyses showing that numerous defense transcripts and proteins were up-regulated in the snap33 mutant, among which genes/proteins involved in defense hormone, pattern-triggered immunity, and nucleotide-binding domain leucine-rich-repeat receptor signaling. qRT-PCR analyses and hormone dosages supported these results. Furthermore, genetic analyses elucidated the diverse contributions of the main defense hormones and some nucleotide-binding domain leucine-rich-repeat receptor signaling actors in the establishment of the snap33 phenotype, emphasizing the preponderant role of salicylic acid over other defense phytohormones. Moreover, the accumulation of pattern-triggered immunity and nucleotide-binding domain leucine-rich-repeat receptor signaling proteins in the snap33 mutant was confirmed by immunoblotting analyses and further shown to be salicylic acid-dependent. Collectively, this study unveiled molecular determinants underlying the Arabidopsis snap33 mutant phenotype and brought new insights into autoimmunity signaling.

Pathogen perception and signaling in plant immunity

Plant diseases are a constant and serious threat to agriculture and ecological biodiversity. Plants possess a sophisticated innate immunity system capable of detecting and responding to pathogen infection to prevent disease. Our understanding of this system has grown enormously over the past century. Early genetic descriptions of plant disease resistance and pathogen virulence were embodied in the gene-for-gene hypothesis, while physiological studies identified pathogen-derived elicitors that could trigger defense responses in plant cells and tissues. Molecular studies of these phenomena have now coalesced into an integrated model of plant immunity involving cell surface and intracellular detection of specific pathogen-derived molecules and proteins culminating in the induction of various cellular responses. Extracellular and intracellular receptors engage distinct signaling processes but converge on many similar outputs with substantial evidence now for integration of these pathways into interdependent networks controlling disease outcomes. Many of the molecular details of pathogen recognition and signaling processes are now known, providing opportunities for bioengineering to enhance plant protection from disease. Here we provide an overview of the current understanding of the main principles of plant immunity, with an emphasis on the key scientific milestones leading to these insights.

A new spider mite elicitor triggers plant defence and promotes resistance to herbivores

Herbivore-associated elicitors (HAEs) are active molecules produced by herbivorous insects. Recognition of HAEs by plants induces defence that resist herbivore attacks. We previously demonstrated that the tomato red spider mite Tetranychus evansi triggered defence in Nicotiana benthamiana. However, our knowledge of HAEs from T. evansi remains limited. Here, we characterize a novel HAE, Te16, from T. evansi and dissect its function in mite-plant interactions. We investigate the effects of Te16 on spider mites and plants by heterologous expression, virus-induced gene silencing assay, and RNA interference. Te16 induces cell death, reactive oxygen species (ROS) accumulation, callose deposition, and jasmonate (JA)-related responses in N. benthamiana leaves. Te16-mediated cell death requires a calcium signalling pathway, cytoplasmic localization, the plant co-receptor BAK1, and the signalling components SGT1 and HSP90. The active region of Te16-induced cell death is located at amino acids 114-293. Moreover, silencing Te16 gene in T. evansi reduces spider mite survival and hatchability, but expressing Te16 in N. benthamiana leaves enhances plant resistance to herbivores. Finally, Te16 gene is specific to Tetranychidae species and is highly conserved in activating plant immunity. Our findings reveal a novel salivary protein produced by spider mites that elicits plant defence and resistance to insects, providing valuable clues for pest management.

Identification of genetic loci and candidate genes underlying freezing tolerance in wheat seedlings

KEY MESSAGE

Ten stable loci for freezing tolerance (FT) in wheat were detected by genome-wide association analysis. The putative candidate gene TaRPM1-7BL underlying the major locus QFT.ahau-7B.2 was identified and validated. Frost damage restricts wheat growth, development, and geographical distribution. However, the genetic mechanism of freezing tolerance (FT) remains unclear. Here, we evaluated FT phenotypes of 245 wheat varieties and lines, and genotyped them using a Wheat 90 K array. The association analysis showed that ten stable loci were significantly associated with FT (P < 1 × 10-4), and explained 6.45-26.33% of the phenotypic variation. In particular, the major locus QFT.ahau-7B.2 was consistently related to all nine sets of FT phenotypic data. Based on five cleaved amplified polymorphic sequence (CAPS) markers closely linked to QFT.ahau-7B.2, we narrowed down the target region to the 570.67-571.16 Mb interval (0.49 Mb) on chromosome 7B, in which four candidate genes were annotated. Of these, only TaRPM1-7BL exhibited consistent differential expression after low temperature treatment between freezing-tolerant and freezing-sensitive varieties. The results of cloning and whole-exome capture sequencing indicated that there were two main haplotypes for TaRPM1-7BL, including freezing-tolerant Hap1 and freezing-sensitive Hap2. Based on the representative SNP (+1956, A/G), leading to an amino acid change in the NBS domain, a CAPS marker (CAPS-TaRPM1-7BL) was developed and validated in 431 wheat varieties (including the above 245 materials) and 318 F2 lines derived from the cross of 'Annong 9267' (freezing-tolerant) × 'Yumai 9' (freezing-sensitive). Subsequently, the TaRPM1-7BL gene was silenced in 'Yumai 9' by virus-induced gene silencing (VIGS), and these silenced wheat seedlings exhibited enhanced FT phenotypes, suggesting that TaRPM1-7BL negatively regulates FT. These findings are valuable for understanding the complex genetic basis of FT in wheat.

Rice NLR protein XinN1, induced by a pattern recognition receptor XA21, confers enhanced resistance to bacterial blight

KEY MESSAGE

Rice CC-type NLR XinN1, specifically induced by a PRR XA21, activates defense pathways against Xoo. Plants have evolved two layers of immune systems regulated by two different types of immune receptors, cell surface located pattern recognition receptors (PRRs) and intracellular nucleotide-binding domain leucine-rich repeat-containing receptors (NLRs). Plant PRRs recognize conserved molecular patterns from diverse pathogens, resulting in pattern-triggered immunity (PTI), whereas NLRs are activated by effectors secreted by pathogens into plant cells, inducing effector-triggered immunity (ETI). Rice PRR, XA21, recognizes a tyrosine-sulfated RaxX peptide (required for activation of XA21-mediated immunity X) as a molecular pattern secreted by Xanthomonas oryzae pv. oryzae (Xoo). Here, we identified a rice NLR gene, XinN1, that is specifically induced during the XA21-mediated immune response against Xoo. Transgenic rice plants overexpressing XinN1 displayed increased resistance to infection by Xoo with reduced lesion length and bacterial growth. Overexpression of autoactive mutant of XinN1 (XinN1D543V) also displayed increased resistance to Xoo, accompanied with severe growth retardation and cell death. In rice protoplast system, overexpression of XinN1 or XinN1D543V significantly elevated reactive oxygen species (ROS) production and cytosolic-free calcium (Ca2+) accumulations. In addition, XinN1 overexpression additionally elevated the ROS burst caused by the interaction between XA21 and RaxX-sY and induced the transcription of PTI signaling components, including somatic embryogenesis receptor kinases (OsSERKs) and receptor-like cytoplasmic kinases (OsRLCKs). Our results suggest that XinN1 induced by the PRR XA21 activates defense pathways and provides enhanced resistance to Xoo in rice.

Editing Metabolism, Sex, and Microbiome: How Can We Help Poplar Resist Pathogens?

Poplar (Populus) is a genus of woody plants of great economic value. Due to the growing economic importance of poplar, there is a need to ensure its stable growth by increasing its resistance to pathogens. Genetic engineering can create organisms with improved traits faster than traditional methods, and with the development of CRISPR/Cas-based genome editing systems, scientists have a new highly effective tool for creating valuable genotypes. In this review, we summarize the latest research data on poplar diseases, the biology of their pathogens and how these plants resist pathogens. In the final section, we propose to plant male or mixed poplar populations; consider the genes of the MLO group, transcription factors of the WRKY and MYB families and defensive proteins BbChit1, LJAMP2, MsrA2 and PtDef as the most promising targets for genetic engineering; and also pay attention to the possibility of microbiome engineering.

Integrated proteomic analysis reveals interactions between phosphorylation and ubiquitination in rose response to Botrytis infection

As two of the most abundant post-translational modifications, phosphorylation and ubiquitination play a significant role in modulating plant-pathogen interactions and increasing evidence indicates their crosstalk in plant immunity. Rose (Rosa sp.) is one of the most important ornamental plants and can be seriously infected by Botrytis cinerea. Here, integrated proteomics analysis was performed to detect global proteome, phosphorylation, and ubiquitination changes in rose upon B. cinerea infection and investigate the possible phosphorylation and ubiquitination crosstalk. A total of 6165 proteins, 11 774 phosphorylation and 10 582 ubiquitination sites, and 77 phosphorylation and 13 ubiquitination motifs were identified. Botrytis cinerea infection resulted in 169 up-regulated and 122 down-regulated proteins, 291 up-regulated and 404 down-regulated phosphorylation sites, and 250 up-regulated and 634 down-regulated ubiquitination sites. There were 12 up-regulated PR10 proteins and half of them also showed reduced ubiquitination. A lot of kinases probably involved in plant pattern-triggered immunity signaling were up-regulated phosphoproteins. Noticeably, numerous kinases and ubiquitination-related proteins also showed a significant change in ubiquitination and phosphorylation, respectively. A cross-comparison of phosphoproteome and ubiquitylome indicated that both of two post-translational modifications of 104 proteins were dynamically regulated, and many putative pattern-triggered immunity signaling components in the plant plasma membrane were co-regulated. Moreover, five selected proteins, including four PR10 proteins and a plasma membrane aquaporin, were proven to be involved in rose resistance to B. cinerea. Our study provides insights into the molecular mechanisms underlying rose resistance to B. cinerea and also increases the database of phosphorylation and ubiquitination sites in plants.

Cell specialization and coordination in Arabidopsis leaves upon pathogenic attack revealed by scRNA-seq

Plant defense responses involve several biological processes that allow plants to fight against pathogenic attacks. How these different processes are orchestrated within organs and depend on specific cell types is poorly known. Here, using single-cell RNA sequencing (scRNA-seq) technology on three independent biological replicates, we identified several cell populations representing the core transcriptional responses of wild-type Arabidopsis leaves inoculated with the bacterial pathogen Pseudomonas syringae DC3000. Among these populations, we retrieved major cell types of the leaves (mesophyll, guard, epidermal, companion, and vascular S cells) with which we could associate characteristic transcriptional reprogramming and regulators, thereby specifying different cell-type responses to the pathogen. Further analyses of transcriptional dynamics, on the basis of inference of cell trajectories, indicated that the different cell types, in addition to their characteristic defense responses, can also share similar modules of gene reprogramming, uncovering a ubiquitous antagonism between immune and susceptible processes. Moreover, it appears that the defense responses of vascular S cells, epidermal cells, and mesophyll cells can evolve along two separate paths, one converging toward an identical cell fate, characterized mostly by lignification and detoxification functions. As this divergence does not correspond to the differentiation between immune and susceptible cells, we speculate that this might reflect the discrimination between cell-autonomous and non-cell-autonomous responses. Altogether our data provide an upgraded framework to describe, explore, and explain the specialization and the coordination of plant cell responses upon pathogenic challenge.

NLR signaling in plants: from resistosomes to second messengers

Nucleotide binding and leucine-rich repeat-containing receptors (NLRs) have a critical role in plant immunity through direct or indirect recognition of pathogen effectors. Recent studies have demonstrated that such recognition induces formation of large protein complexes called resistosomes to mediate NLR immune signaling. Some NLR resistosomes activate Ca2+ influx by acting as Ca2+-permeable channels, whereas others function as active NADases to catalyze the production of nucleotide-derived second messengers. In this review we summarize these studies on pathogen effector-induced assembly of NLR resistosomes and resistosome-mediated production of the second messengers of Ca2+ and nucleotide derivatives. We also discuss downstream events and regulation of resistosome signaling.

The responses of poplars to fungal pathogens: A review of the defensive pathway

Long-lived tree species need to cope with changing environments and pathogens during their lifetime. Fungal diseases cause damage to trees growth and forest nurseries. As model system for woody plants, poplars are also hosts of a large variety of fungus. The defense strategies to fungus are generally associated with the type of fungus, therefore, the defense strategies of poplar against necrotrophic and biotrophic fungus are different. Poplars initiate constitutive defenses and induced defenses based on recognition of the fungus, hormone signaling network cascades, activation of defense-related genes and transcription factors and production of phytochemicals. The means of sensing fungus invasion in poplars are similar with herbs, both of which are mediated by receptor proteins and resistance (R) proteins, leading to pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), but poplars have evolved some unique defense mechanisms compared with Arabidopsis due to their longevity. In this paper, current researches on poplar defensive responses to necrotrophic and biotrophic fungus, which mainly include the physiological and genetic aspects, and the role of noncoding RNA (ncRNA) in fungal resistance are reviewed. This review also provides strategies to enhance poplar disease resistance and some new insights into future research directions.

Comprehensive analysis of MAPK gene family in Populus trichocarpa and physiological characterization of PtMAPK3-1 in response to MeJA induction

Mitogen-activated protein kinases (MAPKs) play important roles in plant growth and development, as well as hormone and stress responses by signaling to eukaryotic cells, through MAPK cascade, the presence of various cues; thereby, regulating various responses. The MAPK cascade consists mainly of three gene families, MAPK, MAPKK, and MAPKKK, which activate downstream signaling pathways through sequential phosphorylation. Although the MAPK cascade gene family has been reported in several species, there is a lack of comprehensive analysis in poplar. We identified 21 MAPK genes, 11 MAPKK genes, and 104 MAPKKK genes in Populus trichocarpa. The phylogenetic classification was supported by conservative motif, gene structure and motif analysis. Whole genome duplication has an important role in the expansion of MAPK cascade genes. Analysis of promoter cis-elements and expression profiles indicates that MAPK cascade genes have important roles in plant growth and development, abiotic and biotic stresses, and phytohormone response. Expression profiling revealed a significant upregulation of PtMAPK3-1 expression in response to drought, salt and disease stresses. Poplar transiently overexpressing PtMAPK3-1 and treated with methyl jasmonic acid (MeJA) had higher catalase and peroxidase levels than non-overexpressing poplar. This work represents the first complete inventory of the MAPK cascade in P. trichocarpa, which reveals that PtMAPK3-1 is induced by the MeJA hormone and participates in the MeJA-induced enhancement of the antioxidant enzyme system.

RING-Type E3 Ubiquitin Ligases AtRDUF1 and AtRDUF2 Positively Regulate the Expression of PR1 Gene and Pattern-Triggered Immunity

The importance of E3 ubiquitin ligases from different families for plant immune signaling has been confirmed. Plant RING-type E3 ubiquitin ligases are members of the E3 ligase superfamily and have been shown to play positive or negative roles during the regulation of various steps of plant immunity. Here, we present Arabidopsis RING-type E3 ubiquitin ligases AtRDUF1 and AtRDUF2 which act as positive regulators of flg22- and SA-mediated defense signaling. Expression of AtRDUF1 and AtRDUF2 is induced by pathogen-associated molecular patterns (PAMPs) and pathogens. The atrduf1 and atrduf2 mutants displayed weakened responses when triggered by PAMPs. Immune responses, including oxidative burst, mitogen-activated protein kinase (MAPK) activity, and transcriptional activation of marker genes, were attenuated in the atrduf1 and atrduf2 mutants. The suppressed activation of PTI responses also resulted in enhanced susceptibility to bacterial pathogens. Interestingly, atrduf1 and atrduf2 mutants showed defects in SA-mediated or pathogen-mediated PR1 expression; however, avirulent Pseudomonas syringae pv. tomato DC3000-induced cell death was unaffected. Our findings suggest that AtRDUF1 and AtRDUF2 are not just PTI-positive regulators but are also involved in SA-mediated PR1 gene expression, which is important for resistance to P. syringae.

Molecular mechanisms of resistance to Myzus persicae conferred by the peach Rm2 gene: A multi-omics view

The transcriptomic and metabolomic responses of peach to Myzus persicae infestation were studied in Rubira, an accession carrying the major resistance gene Rm2 causing antixenosis, and GF305, a susceptible accession. Transcriptome and metabolome showed both a massive reconfiguration in Rubira 48 hours after infestation while GF305 displayed very limited changes. The Rubira immune system was massively stimulated, with simultaneous activation of genes encoding cell surface receptors involved in pattern-triggered immunity and cytoplasmic NLRs (nucleotide-binding domain, leucine-rich repeat containing proteins) involved in effector-triggered immunity. Hypersensitive reaction featured by necrotic lesions surrounding stylet punctures was supported by the induction of cell death stimulating NLRs/helpers couples, as well as the activation of H2O2-generating metabolic pathways: photorespiratory glyoxylate synthesis and activation of the futile P5C/proline cycle. The triggering of systemic acquired resistance was suggested by the activation of pipecolate pathway and accumulation of this defense hormone together with salicylate. Important reduction in carbon, nitrogen and sulphur metabolic pools and the repression of many genes related to cell division and growth, consistent with reduced apices elongation, suggested a decline in the nutritional value of apices. Finally, the accumulation of caffeic acid conjugates pointed toward their contribution as deterrent and/or toxic compounds in the mechanisms of resistance.

Connecting the dots between cell surface- and intracellular-triggered immune pathways in plants

Plants can detect microbial molecules via surface-localized pattern-recognition receptors (PRRs) and intracellular immune receptors from the nucleotide-binding, leucine-rich repeat receptor (NLR) family. The corresponding pattern-triggered (PTI) and effector-triggered (ETI) immunity were long considered separate pathways, although they converge on largely similar cellular responses, such as calcium influx and overlapping gene reprogramming. A number of studies recently uncovered genetic and molecular interconnections between PTI and ETI, highlighting the complexity of the plant immune network. Notably, PRR- and NLR-mediated immune responses require and potentiate each other to reach an optimal immune output. How PTI and ETI connect to confer robust immunity in different plant species, including crops will be an exciting future research area.

The TGA Transcription Factors from Clade II Negatively Regulate the Salicylic Acid Accumulation in Arabidopsis

Salicylic acid (SA) is a hormone that modulates plant defenses by inducing changes in gene expression. The mechanisms that control SA accumulation are essential for understanding the defensive process. TGA transcription factors from clade II in Arabidopsis, which include the proteins TGA2, TGA5, and TGA6, are known to be key positive mediators for the transcription of genes such as PR-1 that are induced by SA application. However, unexpectedly, stress conditions that induce SA accumulation, such as infection with the avirulent pathogen P. syringae DC3000/AvrRPM1 and UV-C irradiation, result in enhanced PR-1 induction in plants lacking the clade II TGAs (tga256 plants). Increased PR-1 induction was accompanied by enhanced isochorismate synthase-dependent SA production as well as the upregulation of several genes involved in the hormone’s accumulation. In response to avirulent P. syringae, PR-1 was previously shown to be controlled by both SA-dependent and -independent pathways. Therefore, the enhanced induction of PR-1 (and other defense genes) and accumulation of SA in the tga256 mutant plants is consistent with the clade II TGA factors providing negative feedback regulation of the SA-dependent and/or -independent pathways. Together, our results indicate that the TGA transcription factors from clade II negatively control SA accumulation under stress conditions that induce the hormone production. Our study describes a mechanism involving old actors playing new roles in regulating SA homeostasis under stress.

Concerted actions of PRR- and NLR-mediated immunity

Plants utilise cell-surface immune receptors (functioning as pattern recognition receptors, PRRs) and intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) to detect pathogens. Perception of pathogens by these receptors activates immune signalling and resistance to infections. PRR- and NLR-mediated immunity have primarily been considered parallel processes contributing to disease resistance. Recent studies suggest that these two pathways are interdependent and converge at multiple nodes. This review summarises and provides a perspective on these convergent points.

Distinct Responses to Pathogenic and Symbionic Microorganisms: The Role of Plant Immunity

Plants must balance both beneficial (symbiotic) and pathogenic challenges from microorganisms, the former benefitting the plant and agriculture and the latter causing disease and economic harm. Plant innate immunity describes a highly conserved set of defense mechanisms that play pivotal roles in sensing immunogenic signals associated with both symbiotic and pathogenic microbes and subsequent downstream activation of signaling effector networks that protect the plant. An intriguing question is how the innate immune system distinguishes “friends” from “foes”. Here, we summarize recent advances in our understanding of the role and spectrum of innate immunity in recognizing and responding to different microbes. In addition, we also review some of the strategies used by microbes to manipulate plant signaling pathways and thus evade immunity, with emphasis on the use of effector proteins and micro-RNAs (miRNAs). Furthermore, we discuss potential questions that need addressing to advance the field of plant–microbe interactions.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses

INTRODUCTION

Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown.

OBJECTIVES

We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species.

METHODS

Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues.

RESULTS

Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors.

CONCLUSION

We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.