The Penium margaritaceum Genome: Hallmarks of the Origins of Land Plants.

The evolutionary features and molecular innovations that enabled plants to first colonize land are not well understood. Here, insights are provided through our report of the genome sequence of the unicellular alga Penium margaritaceum, a member of the Zygnematophyceae, the sister lineage to land plants. The genome has a high proportion of repeat sequences that are associated with massive segmental gene duplications, likely facilitating neofunctionalization. Compared with representatives of earlier diverging algal lineages, P. margaritaceum has expanded repertoires of gene families, signaling networks, and adaptive responses that highlight the evolutionary trajectory toward terrestrialization. These encompass a broad range of physiological processes and protective cellular features, such as flavonoid compounds and large families of modifying enzymes involved in cell wall biosynthesis, assembly, and remodeling. Transcriptome profiling further elucidated adaptations, responses, and selective pressures associated with the semi-terrestrial ecosystems of P. margaritaceum, where a simple body plan would be an advantage.

Abscisic acid biosynthesis is necessary for full auxin effects on hypocotyl elongation.

In concert with other phytohormones, auxin regulates plant growth and development. However, how auxin and other phytohormones coordinately regulate distinct processes is not fully understood. In this work, we uncover an auxin-abscisic acid (ABA) interaction module in Arabidopsis that is specific to coordinating activities of these hormones in the hypocotyl. From our forward genetics screen, we determine that ABA biosynthesis is required for the full effects of auxin on hypocotyl elongation. Our data also suggest that ABA biosynthesis is not required for the inhibitory effects of auxin treatment on root elongation. Our transcriptome analysis identified distinct auxin-responsive genes in root and shoot tissues, which is consistent with differential regulation of growth in these tissues. Further, our data suggest that many gene targets repressed upon auxin treatment require an intact ABA pathway for full repression. Our results support a model in which auxin stimulates ABA biosynthesis to fully regulate hypocotyl elongation.

HY5 functions as a systemic signal by integrating BRC1-dependent hormone signaling in tomato bud outgrowth.

Light plays an important role in determining plant architecture, which greatly influences crop yield. However, the precise mechanisms by which light signaling regulates bud outgrowth remain to be identified. Here, we show that light regulates bud outgrowth via both HY5 and brassinosteroid (BR)-dependent pathways in tomato. Inactivation of the red-light photoreceptor PHYB, or deficiencies in PHYB or the blue-light photoreceptor CRY1a, inhibits bud outgrowth and leads to decreased accumulation of HY5 protein and increased transcript level of BRANCHED1 (BRC1), a central integrator of branching signals. HY5, functioning as a mobile systemic signal from leaves, promotes bud outgrowth by directly suppressing BRC1 transcript and activating the transcript of BR biosynthesis genes within the lateral buds in tomato. Furthermore, BRC1 prevents the accumulation of cytokinin (CK) and gibberellin (GA) by directly inhibiting the transcript of CK synthesis gene LOG4, while increasing the transcript levels of CK and GA degradation genes (CKX7, GA2ox4, and GA2ox5), leading to an arrest of bud outgrowth. Moreover, bud outgrowth occurs predominantly in the day due to the suppression of BRC1 transcript by HY5. These findings demonstrate that light-inducible HY5 acts as a systemic signaling factor in fine-tuning the bud outgrowth of tomato.

Structural insights into rice KAI2 receptor provide functional implications for perception and signal transduction.

KAI2 receptors, classified as plant α/ß hydrolase enzymes, are capable of perceiving smoke-derived butenolide signals and endogenous yet unidentified KAI2-ligands (KLs). While the number of functional KAI2 receptors varies among land plant species, rice has only one KAI2 gene. Rice, a significant crop and representative of grasses, relies on KAI2-mediated Arbuscular mycorrhiza (AM) symbioses to flourish in traditionally arid and nutrient-poor environments. This study presents the first crystal structure of an active rice (Oryza sativa, Os) KAI2 hydrolase receptor. Our structural and biochemical analyses uncover grass-unique pocket residues influencing ligand sensitivity and hydrolytic activity. Through structure-guided analysis, we identify a specific residue whose mutation enables the increase or decrease of ligand perception, catalytic activity, and signal transduction. Furthermore, we investigate OsKAI2-mediated signaling by examining its ability to form a complex with its binding partner, the F-box protein DWARF3 (D3) ubiquitin ligase and subsequent degradation of the target substrate OsSMAX1, demonstrating the significant role of hydrophobic interactions in the OsKAI2-D3 interface. This study provides new insights into the diverse and pivotal roles of the OsKAI2 signaling pathway in the plant kingdom, particularly in grasses.

Fidelity in plant hormone modifications catalyzed by Arabidopsis GH3 acyl acid amido synthetases.

GRETCHEN HAGEN 3 (GH3) acyl acid amido synthetases conjugate amino acids to acyl acid hormones to either activate or inactivate the hormone molecule. The largest subgroup of GH3 proteins modify the growth-promoting hormone auxin (indole-3-acetic acid; IAA) with the second largest class activating the defense hormone jasmonic acid (JA). The two-step reaction mechanism of GH3 proteins provides a potential proofreading mechanism to ensure fidelity of hormone modification. Examining pyrophosphate release in the first-half reaction of Arabidopsis GH3 proteins that modify IAA (AtGH3.2/YDK2, AtGH3.5/WES1, AtGH3.17/VAS2), JA (AtGH3.11/JAR1), and other acyl acids (AtGH3.7, AtGH3.12/PBS3) indicates that acyl acid-AMP intermediates are hydrolyzed into acyl acid and AMP in the absence of the amino acid, a typical feature of pre-transfer editing mechanisms. Single-turnover kinetic analysis of AtGH3.2/YDK2 and AtGH3.5/WES1 shows that non-cognate acyl acid-adenylate intermediates are more rapidly hydrolyzed than the cognate IAA-adenylate. In contrast, AtGH3.11/JAR1 only adenylates JA, not IAA. While some of the auxin-conjugating GH3 proteins in Arabidopsis (i.e., AtGH3.5/WES1) accept multiple acyl acid substrates, others, like AtGH3.2/YDK2, are specific for IAA; however, both these proteins share similar active site residues. Biochemical analysis of chimeric variants of AtGH3.2/YDK2 and AtGH3.5/WES1 indicates that the C-terminal domain contributes to selection of cognate acyl acid substrates. These findings suggest that the hydrolysis of non-cognate acyl acid-adenylate intermediates, or proofreading, proceeds via a slowed structural switch that provides a checkpoint for fidelity before the full reaction proceeds.

Ratiometric gibberellin biosensors for the analysis of signaling dynamics and metabolism in plant protoplasts.

Gibberellins (GAs) are major regulators of developmental and growth processes in plants. Using the degradation-based signaling mechanism of GAs, we have built transcriptional regulator (DELLA)-based, genetically encoded ratiometric biosensors as proxies for hormone quantification at high temporal resolution and sensitivity that allow dynamic, rapid and simple analysis in a plant cell system, i.e. Arabidopsis protoplasts. These ratiometric biosensors incorporate a DELLA protein as a degradation target fused to a firefly luciferase connected via a 2A peptide to a renilla luciferase as a co-expressed normalization element. We have implemented these biosensors for all five Arabidopsis DELLA proteins, GA-INSENSITIVE, GAI; REPRESSOR-of-ga1-3, RGA; RGA-like1, RGL1; RGL2 and RGL3, by applying a modular design. The sensors are highly sensitive (in the low pm range), specific and dynamic. As a proof of concept, we have tested the applicability in three domains: the study of substrate specificity and activity of putative GA-oxidases, the characterization of GA transporters, and the use as a discrimination platform coupled to a GA agonists' chemical screening. This work demonstrates the development of a genetically encoded quantitative biosensor complementary to existing tools that allow the visualization of GA in planta.

Genetic Networks in Plant Vascular Development.

Understanding the development of vascular tissues in plants is crucial because the evolution of vasculature enabled plants to thrive on land. Various systems and approaches have been used to advance our knowledge about the genetic regulation of vasculature development, from the scale of single genes to networks. In this review, we provide a perspective on the major approaches used in studying plant vascular development, and we cover the mechanisms and genetic networks underlying vascular tissue specification, patterning, and differentiation.

Developmental responses of roots to limited phosphate availability: research progress and application in cereals.

Phosphorus (P), an essential macronutrient, is crucial for plant growth and development. However, available inorganic phosphate (Pi) is often scarce in soil, and its limited mobility exacerbates P deficiency in plants. Plants have developed complex mechanisms to adapt to Pi-limited soils. The root, the primary interface of the plant with soil, plays an essential role in plant adaptation to Pi-limited soil environments. Root system architecture significantly influences Pi acquisition via the dynamic modulation of primary root and/or crown root length, lateral root proliferation and length, root hair development, and root growth angle in response to Pi availability. This review focuses on the physiological, anatomical, and molecular mechanisms underpinning changes in root development in response to Pi starvation in cereals, mainly focusing on the model monocot plant rice (Oryza sativa). We also review recent efforts to modify root architecture to enhance P uptake efficiency in crops and propose future research directions aimed at the genetic improvement of Pi uptake and use efficiency in crops based on root system architecture.

Systems Phytohormone Responses to Mitochondrial Proteotoxic Stress.

Mitochondrial function is controlled by two separate genomes. This feature makes mitochondria prone to proteotoxic stress when a stoichiometric imbalance occurs in the protein complexes that perform oxidative phosphorylation, which consist of both nuclear- and mitochondrial-encoded proteins. Such a proteotoxic stress is known to induce the mitochondrial unfolded protein response (UPRmt) in animals. It is unknown whether UPRmt occurs in plants. Here, we induced a mitonuclear protein imbalance in Arabidopsis through chemical or genetic interference. Mitochondrial proteotoxic stress activated a plant-specific UPRmt and impaired plant growth and development. The plant UPRmt pathway is triggered by a transient oxidative burst, activating MAPK and hormonal (involving ethylene and auxin) signaling, which are all geared to repair proteostasis. This also establishes phytohormones as bona fide plant mitokines. Our data ascertain that mitochondrial protein quality control pathways, such as the UPRmt, are conserved in plants and that hormone signaling is an essential mediator that regulates mitochondrial proteostasis.

AUXIN RESPONSE FACTOR protein accumulation and function.

Auxin is a key regulator of plant developmental processes. Its effects on transcription are mediated by the AUXIN RESPONSE FACTOR (ARF) family of transcription factors. ARFs tightly control specific auxin responses necessary for proper plant growth and development. Recent research has revealed that regulated ARF protein accumulation and ARF nucleo-cytoplasmic partitioning can determine auxin transcriptional outputs. In this review, we explore these recent findings and consider the potential for regulated ARF accumulation in driving auxin responses in plants.

Ligand diversity contributes to the full activation of the jasmonate pathway in Marchantia polymorpha.

In plants, jasmonate signaling regulates a wide range of processes from growth and development to defense responses and thermotolerance. Jasmonates, such as jasmonic acid (JA), (+)-7-iso-jasmonoyl-l-isoleucine (JA-Ile), 12-oxo-10,15(Z)-phytodienoic acid (OPDA), and dinor-12-oxo-10,15(Z)-phytodienoic acid (dn-OPDA), are derived from C18 (18 Carbon atoms) and C16 polyunsaturated fatty acids (PUFAs), which are found ubiquitously in the plant kingdom. Bryophytes are also rich in C20 and C22 long-chain polyunsaturated fatty acids (LCPUFAs), which are found only at low levels in some vascular plants but are abundant in organisms of other kingdoms, including animals. The existence of bioactive jasmonates derived from LCPUFAs is currently unknown. Here, we describe the identification of an OPDA-like molecule derived from a C20 fatty acid (FA) in the liverwort Marchantia polymorpha (Mp), which we term (5Z,8Z)-10-(4-oxo-5-((Z)-pent-2-en-1-yl)cyclopent-2-en-1-yl)deca-5,8-dienoic acid (C20-OPDA). This molecule accumulates upon wounding and, when applied exogenously, can activate known Coronatine Insensitive 1 (COI1) -dependent and -independent jasmonate responses. Furthermore, we identify a dn-OPDA-like molecule (Δ4-dn-OPDA) deriving from C20-OPDA and demonstrate it to be a ligand of the jasmonate coreceptor (MpCOI1-Mp Jasmonate-Zinc finger inflorescence meristem domain [MpJAZ]) in Marchantia. By analyzing mutants impaired in the production of LCPUFAs, we elucidate the major biosynthetic pathway of C20-OPDA and Δ4-dn-OPDA. Moreover, using a double mutant compromised in the production of both Δ4-dn-OPDA and dn-OPDA, we demonstrate the additive nature of these molecules in the activation of jasmonate responses. Taken together, our data identify a ligand of MpCOI1 and demonstrate LCPUFAs as a source of bioactive jasmonates that are essential to the immune response of M. polymorpha.

Integrated morphological, physiological and transcriptomic analyses reveal the responses of Toona sinensis seedlings to low-nitrogen stress.

Nitrogen is one of the most essential elements for plant growth and development. In this study, the growth, physiology, and transcriptome of Toona sinensis (A. Juss) Roem seedlings were compared between low-nitrogen (LN) and normal-nitrogen (NN) conditions. These results indicate that LN stress adversely influences T. sinensis seedling growth. The activities of key enzymes related to nitrogen assimilation and phytohormone contents were altered by LN stress. A total of 2828 differentially expressed genes (DEGs) in roots and 1547 in leaves were identified between the LN and NN treatments. A differential enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that nitrogen and sugar metabolism, flavonoid biosynthesis, plant hormone signal transduction, and ABC transporters, were strongly affected by LN stress. In summary, this research provides information for further understanding the response of T. sinensis to LN stress.

RAV1 family members function as transcriptional regulators and play a positive role in plant disease resistance.

Phytopathogens pose a severe threat to agriculture and strengthening the plant defense response is an important strategy for disease control. Here, we report that AtRAV1, an AP2 and B3 domain-containing transcription factor, is required for basal plant defense in Arabidopsis thaliana. The atrav1 mutant lines demonstrate hyper-susceptibility against fungal pathogens (Rhizoctonia solani and Botrytis cinerea), whereas AtRAV1 overexpressing lines exhibit disease resistance against them. Enhanced expression of various defense genes and activation of mitogen-activated protein kinases (AtMPK3 and AtMPK6) are observed in the R. solani infected overexpressing lines, but not in the atrav1 mutant plants. An in vitro phosphorylation assay suggests AtRAV1 to be a novel phosphorylation target of AtMPK3. Bimolecular fluorescence complementation and yeast two-hybrid assays support physical interactions between AtRAV1 and AtMPK3. Overexpression of the native as well as phospho-mimic but not the phospho-defective variant of AtRAV1 imparts disease resistance in the atrav1 mutant A. thaliana lines. On the other hand, overexpression of AtRAV1 fails to impart disease resistance in the atmpk3 mutant. These analyses emphasize that AtMPK3-mediated phosphorylation of AtRAV1 is important for the elaboration of the defense response in A. thaliana. Considering that RAV1 homologs are conserved in diverse plant species, we propose that they can be gainfully deployed to impart disease resistance in agriculturally important crop plants. Indeed, overexpression of SlRAV1 (a member of the RAV1 family) imparts disease tolerance against not only fungal (R. solani and B. cinerea), but also against bacterial (Ralstonia solanacearum) pathogens in tomato, whereas silencing of the gene enhances disease susceptibility.

Genome-wide identification, characterization and expression pattern analysis of TIFY family members in Artemisia argyi.

BACKGROUND: Plant-specific TIFY proteins play crucial roles in regulating plant growth, development, and various stress responses. However, there is no information available about this family in Artemisia argyi, a well-known traditional medicinal plant with great economic value. RESULTS: A total of 34 AaTIFY genes were identified, including 4 TIFY, 22 JAZ, 5 PPD, and 3 ZML genes. Structural, motif scanning, and phylogenetic relationships analysis of these genes revealed that members within the same group or subgroup exhibit similar exon-intron structures and conserved motif compositions. The TIFY genes were unevenly distributed across the 15 chromosomes. Tandem duplication events and segmental duplication events have been identified in the TIFY family in A. argyi. These events have played a crucial role in the gene multiplication and compression of different subfamilies within the TIFY family. Promoter analysis revealed that most AaTIFY genes contain multiple cis-elements associated with stress response, phytohormone signal transduction, and plant growth and development. Expression analysis of roots and leaves using RNA-seq data revealed that certain AaTIFY genes showed tissue-specific expression patterns, and some AaTIFY genes, such as AaTIFY19/29, were found to be involved in regulating salt and saline-alkali stresses. In addition, RT-qPCR analysis showed that TIFY genes, especially AaTIFY19/23/27/29, respond to a variety of hormonal treatments, such as MeJA, ABA, SA, and IAA. This suggested that TIFY genes in A. argyi regulate plant growth and respond to different stresses by following different hormone signaling pathways. CONCLUSION: Taken together, our study conducted a comprehensive identification and analysis of the TIFY gene family in A. argyi. These findings suggested that TIFY might play an important role in plant development and stress responses, which laid a valuable foundation for further understanding the function of TIFY genes in multiple stress responses and phytohormone crosstalk in A. argyi.

Genome-wide identification of B-box zinc finger (BBX) gene family in Medicago sativa and their roles in abiotic stress responses.

BACKGROUND: B-box (BBX) family is a class of zinc finger transcription factors (TFs) that play essential roles in regulating plant growth, development, as well as abiotic stress. However, no systematic analysis of BBX genes has yet been conducted in alfalfa (Medica go sativa L.), and their functions have not been elucidated up to now. RESULTS: In this study, 28 MsBBX genes were identified from the alfalfa genome, which were clustered into 4 subfamilies according to an evolutionary tree of BBX proteins. Exon-intron structure and conserved motif analysis reflected the evolutionary conservation of MsBBXs in alfalfa. Collinearity analysis showed that segmental duplication promoted the expansion of the MsBBX family. Analysis of cis-regulatory elements suggested that the MsBBX genes possessed many growth/development-, light-, phytohormone-, and abiotic stress-related elements. MsBBX genes were differentially expressed in leaves, flowers, pre-elongated stems, elongated stems, roots and nodules, and most MsBBXs were remarkably induced by drought, salt and various plant growth regulators (ABA, JA, and SA). Further functional verification demonstrated that overexpressing of the MsBBX11 gene clearly promoted salt tolerance in transgenic Arabidopsis by regulating growth and physiological processes of seedlings. CONCLUSIONS: This research provides insights into further functional research and regulatory mechanisms of MsBBX family genes under abiotic stress of alfalfa.

The Diterpene Isopimaric Acid Modulates the Phytohormone Pathway to Promote Oryza sativa L. Rice Seedling Growth.

Many plant secondary metabolites are active and important in the regulation of plant growth. Certain plant-derived diterpenes are known to promote plant growth, but the pathways by which this promotion occurs are still unknown. Activity screening revealed that the plant-derived diterpene isopimaric acid exhibits growth-promoting activity in rice (Oryza sativa L.) seedlings. Furthermore, 25 µg/mL of isopimaric acid promoted the growth of 15 self-incompatible associated populations from different rice lineages to different extents. Quantitative analyses revealed a significant decrease in the concentration of the defense-related phytohormone abscisic acid (ABA) following treatment with isopimaric acid. Correlation analysis of the phytohormone concentrations with growth characteristics revealed that the length of seedling shoots was significantly negatively correlated with concentrations of 3-indole-butyric acid (IBA). Moreover, the total root weight was not only negatively correlated with ABA concentrations but also negatively correlated with concentrations of isopentenyl adenine (iP). These data suggest that isopimaric acid is able to influence the phytohormone pathway to balance energy allocation between growth and defense in rice seedlings and also alter the correlation between the concentrations of phytohormones and traits such as shoot and root length and weight. We provide a theoretical basis for the development and utilization of isopimaric acid as a plant growth regulator for rice.

Genome-wide identification and expression analysis of GRAS gene family in Eucalyptus grandis.

BACKGROUND: The GRAS gene family is a class of plant-specific transcription factors with important roles in many biological processes, such as signal transduction, disease resistance and stress tolerance, plant growth and development. So far, no information available describes the functions of the GRAS genes in Eucalyptus grandis. RESULTS: A total of 82 GRAS genes were identified with amino acid lengths ranging from 267 to 817 aa, and most EgrGRAS genes had one exon. Members of the GRAS gene family of Eucalyptus grandis are divided into 9 subfamilies with different protein structures, while members of the same subfamily have similar gene structures and conserved motifs. Moreover, these EgrGRAS genes expanded primarily due to segmental duplication. In addition, cis-acting element analysis showed that this family of genes was involved involved in the signal transduction of various plant hormones, growth and development, and stress response. The qRT-PCR data indicated that 18 EgrGRAS genes significantly responded to hormonal and abiotic stresses. Among them, the expression of EgrGRAS13, EgrGRAS68 and EgrGRAS55 genes was significantly up-regulated during the treatment period, and it was hypothesised that members of the EgrGRAS family play an important role in stress tolerance. CONCLUSIONS: In this study, the phylogenetic relationship, conserved domains, cis-elements and expression patterns of GRAS gene family of Eucalyptus grandis were analyzed, which filled the gap in the identification of GRAS gene family of Eucalyptus grandis and laid the foundation for analyzing the function of EgrGRAS gene in hormone and stress response.

RhMED15a-like, a subunit of the Mediator complex, is involved in the drought stress response in Rosa hybrida.

BACKGROUND: Rose (Rosa hybrida) is a globally recognized ornamental plant whose growth and distribution are strongly limited by drought stress. The role of Mediator, a multiprotein complex crucial for RNA polymerase II-driven transcription, has been elucidated in drought stress responses in plants. However, its physiological function and regulatory mechanism in horticultural crop species remain elusive. RESULTS: In this study, we identified a Tail module subunit of Mediator, RhMED15a-like, in rose. Drought stress, as well as treatment with methyl jasmonate (MeJA) and abscisic acid (ABA), significantly suppressed the transcript level of RhMED15a-like. Overexpressing RhMED15a-like markedly bolstered the osmotic stress tolerance of Arabidopsis, as evidenced by increased germination rate, root length, and fresh weight. In contrast, the silencing of RhMED15a-like through virus induced gene silencing in rose resulted in elevated malondialdehyde accumulation, exacerbated leaf wilting, reduced survival rate, and downregulated expression of drought-responsive genes during drought stress. Additionally, using RNA-seq, we identified 972 differentially expressed genes (DEGs) between tobacco rattle virus (TRV)-RhMED15a-like plants and TRV controls. Gene Ontology (GO) analysis revealed that some DEGs were predominantly associated with terms related to the oxidative stress response, such as 'response to reactive oxygen species' and 'peroxisome'. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment highlighted pathways related to 'plant hormone signal transduction', in which the majority of DEGs in the jasmonate (JA) and ABA signalling pathways were induced in TRV-RhMED15a-like plants. CONCLUSION: Our findings underscore the pivotal role of the Mediator subunit RhMED15a-like in the ability of rose to withstand drought stress, probably by controlling the transcript levels of drought-responsive genes and signalling pathway elements of stress-related hormones, providing a solid foundation for future research into the molecular mechanisms underlying drought tolerance in rose.

New insights into azelaic acid-induced resistance against Alternaria Solani in tomato plants.

BACKGROUND: The effect of azelaic acid (Aza) on the response of tomato plants to Alternaria solani was investigated in this study. After being treated with Aza, tomato plants were infected with A. solani, and their antioxidant, biochemical, and molecular responses were analyzed. RESULTS: The results demonstrated that H2O2 and MDA accumulation increased in control plants after pathogen infection. Aza-treated plants exhibited a remarkable rise in peroxidase (POD) and catalase (CAT) activities during the initial stages of A. solani infection. Gene expression analysis revealed that both Aza treatment and pathogen infection altered the expression patterns of the SlNPR1, SlERF2, SlPR1, and SlPDF1.2 genes. The expression of SlPDF1.2, a marker gene for the jasmonic acid/ethylene (JA/ET) signaling pathway, showed a remarkable increase of 4.2-fold upon pathogen infection. In contrast, for the SlNPR1, a key gene in salicylic acid (SA) pathway, this increased expression was recorded with a delay at 96 hpi. Also, the phytohormone analysis showed significantly increased SA accumulation in plant tissues with disease development. It was also revealed that tissue accumulation of JA in Aza-treated plants was increased following pathogen infection, while it was not increased in plants without pathogen inoculation. CONCLUSION: The results suggest that the resistance induced by Aza is mainly a result of modulations in both SA and JA pathways following complex antioxidant and molecular defense responses in tomato plants during A. solani infection. These findings provide novel information regarding inducing mechanisms of azelaic acid which would add to the current body of knowledge of SAR induction in plants as result of Aza application.

Cold-stress induced metabolomic and transcriptomic changes in leaves of three mango varieties with different cold tolerance.

BACKGROUND: Mango (Mangifera indica L.) is grown in Hainan, Guangdong, Yunnan, Sichuan, and Fujian provinces and Guanxi autonomous region of China. However, trees growing in these areas suffer severe cold stress during winter, which affects the yield. To this regard, data on global metabolome and transcriptome profiles of leaves are limited. Here, we used combined metabolome and transcriptome analyses of leaves of three mango cultivars with different cold stress tolerance, i.e. Jinhuang (J)-tolerant, Tainung (T) and Guiremang No. 82 (G)-susceptible, after 24 (LF), 48 (MF) and 72 (HF) hours of cold. RESULTS: A total of 1,323 metabolites belonging to 12 compound classes were detected. Of these, amino acids and derivatives, nucleotides and derivatives, and lipids accumulated in higher quantities after cold stress exposure in the three cultivars. Notably, Jinhuang leaves showed increasing accumulation trends of flavonoids, terpenoids, lignans and coumarins, and alkaloids with exposure time. Among the phytohormones, jasmonic acid and abscisic acid levels decreased, while N6-isopentenyladenine increased with cold stress time. Transcriptome analysis led to the identification of 22,526 differentially expressed genes. Many genes enriched in photosynthesis, antenna proteins, flavonoid, terpenoid (di- and sesquiterpenoids) and alkaloid biosynthesis pathways were upregulated in Jihuang leaves. Moreover, expression changes related to phytohormones, MAPK (including calcium and H2O2), and the ICE-CBF-COR signalling cascade indicate involvement of these pathways in cold stress responses. CONCLUSION: Cold stress tolerance in mango leaves is associated with regulation of primary and secondary metabolite biosynthesis pathways. Jasmonic acid, abscisic acid, and cytokinins are potential regulators of cold stress responses in mango leaves.