Transcriptome and Co-expression Network Analyses Reveal Differential Gene Expression and Pathways in Response to Severe Drought Stress in Peanut (Arachis hypogaea L.)

Transcriptome sequencing of garlic reveals key genes related to the heat stress response

With global warming, heat stress has become an important factor that seriously affects crop yield and quality. Therefore, understanding plant responses to heat stress is important for agricultural practice, but the molecular mechanism of high-temperature tolerance in garlic remains unclear. In this study, 'Xusuan No. 6' was used as the experimental material. After heat stress for 0 (CK), 2 and 24 h, transcriptome sequencing was used to screen metabolic pathways and differentially expressed genes (DEGs) closely related to heat stress and was further verified by quantitative real-time polymerase chain reaction (qRT-PCR). A total of 86,110 unigenes obtained from the raw transcriptome sequencing data were spliced. After 2 h of heat treatment, the expression levels of 8898 genes increased, and 3829 genes were decreased in leaves. After 24 h, the expression levels of 7167 genes were upregulated, and 3176 genes were downregulated. Gene Ontology enrichment analysis showed that DEGs were mainly enriched in seven categories: cellular processes, metabolic processes, binging, catalytic activity, cellular anatomical entity and protein-containing complex response to stimulus. Kyoto Encyclopedia of Genes and Genomes pathway enrichment showed that DEGs are involved in protein processing in the endoplasmic reticulum, plant hormone signal transduction, phenylpropanoid biosynthesis, and photosynthetic antenna proteins. Six genes were selected and further verified by qRT-PCR. In this study, the full-length transcriptome of garlic was constructed, and the regulatory genes related to the heat resistance of garlic were studied. Taken together, these findings can provide a theoretical basis for the cloning of heat resistance genes in garlic and for the analysis of heat resistance mechanisms.

Transcriptome Analysis of Sesame (Sesamum indicum L.) Reveals the LncRNA and mRNA Regulatory Network Responding to Low Nitrogen Stress

Nitrogen is one of the important factors restricting the development of sesame planting and industry in China. Cultivating sesame varieties tolerant to low nitrogen is an effective way to solve the problem of crop nitrogen deficiency. To date, the mechanism of low nitrogen tolerance in sesame has not been elucidated at the transcriptional level. In this study, two sesame varieties Zhengzhi HL05 (ZZ, nitrogen efficient) and Burmese prolific (MD, nitrogen inefficient) in low nitrogen were used for RNA-sequencing. A total of 3964 DEGs (differentially expressed genes) and 221 DELs (differentially expressed lncRNAs) were identified in two sesame varieties at 3d and 9d after low nitrogen stress. Among them, 1227 genes related to low nitrogen tolerance are mainly located in amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and participate in the process of transporter activity and antioxidant activity. In addition, a total of 209 pairs of lncRNA-mRNA were detected, including 21 pairs of trans and 188 cis. WGCNA (weighted gene co-expression network analysis) analysis divided the obtained genes into 29 modules; phenotypic association analysis identified three low-nitrogen response modules; through lncRNA-mRNA co-expression network, a number of hub genes and cis/trans-regulatory factors were identified in response to low-nitrogen stress including GS1-2 (glutamine synthetase 1-2), PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase, CHS), CAB21 (chlorophyll a-b binding protein 21) and transcription factors MYB54, MYB88 and NAC75 and so on. As a trans regulator, lncRNA MSTRG.13854.1 affects the expression of some genes related to low nitrogen response by regulating the expression of MYB54, thus responding to low nitrogen stress. Our research is the first to provide a more comprehensive understanding of DEGs involved in the low nitrogen stress of sesame at the transcriptome level. These results may reveal insights into the molecular mechanisms of low nitrogen tolerance in sesame and provide diverse genetic resources involved in low nitrogen tolerance research.

The genus Arachis: an excellent resource for studies on differential gene expression for stress tolerance

Peanut Arachis hypogaea is a segmental allotetraploid in the section Arachis of the genus Arachis along with the Section Rhizomataceae. Section Arachis has several diploid species along with Arachis hypogaea and A. monticola. The section Rhizomataceae comprises polyploid species. Several species in the genus are highly tolerant to biotic and abiotic stresses and provide excellent sets of genotypes for studies on differential gene expression. Though there were several studies in this direction, more studies are needed to identify more and more gene combinations. Next generation RNA-seq based differential gene expression study is a powerful tool to identify the genes and regulatory pathways involved in stress tolerance. Transcriptomic and proteomic study of peanut plants under biotic stresses reveals a number of differentially expressed genes such as R genes (NBS-LRR, LRR-RLK, protein kinases, MAP kinases), pathogenesis related proteins (PR1, PR2, PR5, PR10) and defense related genes (defensin, F-box, glutathione S-transferase) that are the most consistently expressed genes throughout the studies reported so far. In most of the studies on biotic stress induction, the differentially expressed genes involved in the process with enriched pathways showed plant-pathogen interactions, phenylpropanoid biosynthesis, defense and signal transduction. Differential gene expression studies in response to abiotic stresses, reported the most commonly expressed genes are transcription factors (MYB, WRKY, NAC, bZIP, bHLH, AP2/ERF), LEA proteins, chitinase, aquaporins, F-box, cytochrome p450 and ROS scavenging enzymes. These differentially expressed genes are in enriched pathways of transcription regulation, starch and sucrose metabolism, signal transduction and biosynthesis of unsaturated fatty acids. These identified differentially expressed genes provide a better understanding of the resistance/tolerance mechanism, and the genes for manipulating biotic and abiotic stress tolerance in peanut and other crop plants. There are a number of differentially expressed genes during biotic and abiotic stresses were successfully characterized in peanut or model plants (tobacco or Arabidopsis) by genetic manipulation to develop stress tolerance plants, which have been detailed out in this review and more concerted studies are needed to identify more and more gene/gene combinations.

A 3-Ketoacyl-CoA Synthase 10 (KCS10) Homologue from Alfalfa Enhances Drought Tolerance by Regulating Cuticular Wax Biosynthesis

Cuticular wax, forming the first line of defense against adverse environmental stresses, comprises very long-chain fatty acids (VLCFAs) and their derivatives. 3-Ketoacyl-CoA synthase (KCS) is a rate-limiting enzyme for VLCFA biosynthesis. In this study, we isolated KCS10, a KCS gene from alfalfa, and analyzed the effect of gene expression on wax production and drought stress in transgenic plants. MsKCS10 overexpression increased compact platelet-like crystal deposition and promoted primary alcohol biosynthesis through acyl reduction pathways in alfalfa leaves. Overexpression of MsKCS10 induced the formation of coiled-rodlet-like crystals and increased n-alkane content through decarbonylation pathways in tobacco and tomato fruits. Overexpression of MsKCS10 enhanced drought tolerance by limiting nonstomatal water loss, improving photosynthesis, and maintaining osmotic potential under drought stress in transgenic tobacco. In summary, MsKCS10 plays an important role in wax biosynthesis, wax crystal morphology, and drought tolerance, although the mechanisms are different among the plant species. MsKCS10 can be targeted in future breeding programs to improve drought tolerance in plants.

Comprehensive Transcriptome and Metabolome Analyses Reveal Primary Molecular Regulation Pathways Involved in Peanut under Water and Nitrogen Co-Limitation

The yield and quality of peanut (Arachis hypogaea L.), an oil crop planted worldwide, are often limited by drought stress (DS) and nitrogen (N) deficiency. To investigate the molecular mechanism by which peanut counteracts DS and N deficiency, we conducted comprehensive transcriptomic and metabolomic analyses of peanut leaves. Herein, 829 known differentially accumulated metabolites, 324 differentially expressed transcription factors, and 5294 differentially expressed genes (DEGs) were identified under different water and N conditions. The transcriptome analysis demonstrated that drought-related DEGs were predominantly expressed in "glycolysis/gluconeogenesis" and "glycerolipid metabolism", while N-deficiency-related DEGs were mainly expressed in starch and sucrose metabolism, as well as in the biosynthesis of amino acid pathways. The biosynthesis, transport, and catabolism of secondary metabolites accounted for a large proportion of the 1317 DEGs present in water and N co-limitation. Metabolomic analysis showed that the metabolic accumulation of these pathways was significantly dependent on the stress conditions. Additionally, the roles of metabolites and genes in these pathways, such as the biosynthesis of amino acids and phenylpropanoid biosynthesis under different stress conditions, were discussed. The results demonstrated that different genes, metabolic pathways, and metabolites were related to DS and N deficiency. Thus, this study elucidates the metabolic pathways and functional genes that can be used for the improvement of peanut resistance to abiotic stress.

The basic helix-loop-helix transcription factor gene, OsbHLH38, plays a key role in controlling rice salt tolerance

The plant hormone abscisic acid (ABA) is crucial for plant seed germination and abiotic stress tolerance. However, the association between ABA sensitivity and plant abiotic stress tolerance remains largely unknown. In this study, 436 rice accessions were assessed for their sensitivity to ABA during seed germination. The considerable diversity in ABA sensitivity among rice germplasm accessions was primarily reflected by the differentiation between the Xian (indica) and Geng (japonica) subspecies and between the upland-Geng and lowland-Geng ecotypes. The upland-Geng accessions were most sensitive to ABA. Genome-wide association analyses identified four major quantitative trait loci containing 21 candidate genes associated with ABA sensitivity of which a basic helix-loop-helix transcription factor gene, OsbHLH38, was the most important for ABA sensitivity. Comprehensive functional analyses using knockout and overexpression transgenic lines revealed that OsbHLH38 expression was responsive to multiple abiotic stresses. Overexpression of OsbHLH38 increased seedling salt tolerance, while knockout of OsbHLH38 increased sensitivity to salt stress. A salt-responsive transcription factor, OsDREB2A, interacted with OsbHLH38 and was directly regulated by OsbHLH38. Moreover, OsbHLH38 affected rice abiotic stress tolerance by mediating the expression of a large set of transporter genes of phytohormones, transcription factor genes, and many downstream genes with diverse functions, including photosynthesis, redox homeostasis, and abiotic stress responsiveness. These results demonstrated that OsbHLH38 is a key regulator in plant abiotic stress tolerance.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Genome-wide identification of NF-Y gene family in maize (Zea mays L.) and the positive role of ZmNF-YC12 in drought resistance and recovery ability

Nuclear factor Y (NF-Y) genes play important roles in many biological processes, such as leaf growth, nitrogen nutrition, and drought resistance. However, the biological functions of these transcription factor family members have not been systematically analyzed in maize. In the present study, a total of 52 ZmNF-Y genes were identified and classified into three groups in the maize genome. An analysis of the evolutionary relationship, gene structure, and conserved motifs of these genes supports the evolutionary conservation of NF-Y family genes in maize. The tissue expression profiles based on RNA-seq data showed that all genes apart from ZmNF-Y16, ZmNF-YC15, and ZmNF-YC17 were expressed in different maize tissues. A weighted correlation network analysis was conducted and a gene co expression network method was used to analyze the transcriptome sequencing results; six core genes responding to drought and rewatering were identified. A real time fluorescence quantitative analysis showed that these six genes responded to high temperature, drought, high salt, and abscisic acid (ABA) treatments, and subsequent restoration to normal levels. ZmNF-YC12 was highly induced by drought and rewatering treatments. The ZmNF-YC12 protein was localized in the nucleus, and the Gal4-LexA/UAS system and a transactivation analysis demonstrated that ZmNF-YC12 in maize (Zea mays L.) is a transcriptional activator that regulates drought resistance and recovery ability. Silencing ZmNF-YC12 reduced net photosynthesis, chlorophyll content, antioxidant (superoxide dismutase, catalase, peroxidase and ascorbate peroxidase) system activation, and soluble protein and proline contents; it increased the malondialdehyde content, the relative water content, and the water loss rate, which weakened drought resistance and the recoverability of maize. These results provide insights into understanding the evolution of ZmNF-Y family genes in maize and their potential roles in genetic improvement. Our work provides a foundation for subsequent functional studies of the NF-Y gene family and provides deep insights into the role of the ZmNF-YC12 regulatory network in controlling drought resistance and the recoverability of maize.

Identification of drought-responsive phenolic compounds and their biosynthetic regulation under drought stress in Ligularia fischeri

Ligularia fischeri, a leafy edible plant found in damp shady regions, has been used as an herbal medicine and is also consumed as a horticultural crop. In this study, we investigated the physiological and transcriptomic changes, especially those involved in phenylpropanoid biosynthesis, induced by severe drought stress in L. fischeri plants. A distinguishing characteristic of L. fischeri is a color change from green to purple due to anthocyanin biosynthesis. We chromatographically isolated and identified two anthocyanins and two flavones upregulated by drought stress using liquid chromatography-mass spectrometry and nuclear magnetic resonance analyses in this plant for the first time. In contrast, all types of caffeoylquinic acids (CQAs) and flavonol contents were decreased under drought stress. Further, we performed RNA sequencing to examine the molecular changes in these phenolic compounds at the transcriptome level. In an overview of drought-inducible responses, we identified 2,105 hits for 516 distinct transcripts as drought-responsive genes. Moreover, differentially expressed genes (DEGs) associated with phenylpropanoid biosynthesis accounted for the greatest number of both up- and downregulated DEGs by Kyoto Encyclopedia of Genes and Genomes enrichment analysis. We identified 24 meaningful DEGs based on the regulation of phenylpropanoid biosynthetic genes. Potential drought-responsive genes included upregulated flavone synthase (LfFNS, TRINITY DN31661 c0 g1 i1) and anthocyanin 5-O-glucosyltransferase (LfA5GT1, TRINITY DN782 c0 g1 i1), which could contribute to the high levels of flavones and anthocyanins under drought stress in L. fischeri. In addition, the downregulated shikimate O-hydroxycinnamolytransferase (LfHCT, TRINITY DN31661 c0 g1 i1) and hydroxycinnamoyl-CoA quinate/shikimate transferase (LfHQT4, TRINITY DN15180 c0 g1 i1) genes led to a reduction in CQAs. Only one or two BLASTP hits for LfHCT were obtained for six different Asteraceae species. It is possible that the HCT gene plays a crucial role in CQAs biosynthesis in these species. These findings expand our knowledge of the response mechanisms to drought stress, particularly regarding the regulation of key phenylpropanoid biosynthetic genes in L. fischeri.

Comparative physiological, transcriptomic, and WGCNA analyses reveal the key genes and regulatory pathways associated with drought tolerance in Tartary buckwheat

Drought stress is one of the major abiotic stress factors that affect plant growth and crop productivity. Tartary buckwheat is a nutritionally balanced and flavonoid-rich pseudocereal crop and also has strong adaptability to different adverse environments including drought. However, little is known about its drought tolerance mechanism. In this study, we performed comparative physiological and transcriptomic analyses of two contrasting drought-resistant Tartary buckwheat genotypes under nature drought treatment in the reproductive stage. Under drought stress, the drought-tolerant genotype XZSN had significantly higher contents of relative water, proline, and soluble sugar, as well as lower relative electrolyte leakage in the leaves than the drought-susceptible LK3. A total of 5,058 (2,165 upregulated and 2,893 downregulated) and 5,182 (2,358 upregulated and 2,824 downregulated) potential drought-responsive genes were identified in XZSN and LK3 by transcriptome sequencing analysis, respectively. Among the potential drought-responsive genes of XZSN, 1,206 and 1,274 genes were identified to be potential positive and negative contributors for XZSN having higher drought resistance ability than LK3. Furthermore, 851 out of 1,206 positive drought-resistant genes were further identified to be the core drought-resistant genes of XZSN based on WGCNA analysis, and most of them were induced earlier and quicker by drought stress than those in LK3. Functional annotation of the 851 core drought-resistant genes found that a large number of stress-responsive genes were involved in TFs, abscisic acid (ABA) biosynthesis, signal transduction and response, non-ABA signal molecule biosynthesis, water holding, oxygen species scavenging, osmotic adjustment, cell damage prevention, and so on. Transcriptional regulatory network analyses identified the potential regulators of these drought-resistant functional genes and found that the HD-ZIP and MYB TFs might be the key downstream TFs of drought resistance in Tartary buckwheat. Taken together, these results indicated that the XZSN genotype was more drought-tolerant than the LK3 genotype as evidenced by triggering the rapid and dramatic transcriptional reprogramming of drought-resistant genes to reduce water loss, prevent cell damage, and so on. This research expands our current understanding of the drought tolerance mechanisms of Tartary buckwheat and provides important information for its further drought resistance research and variety breeding.

Exogenous application of acetic acid enhances drought tolerance by influencing the MAPK signaling pathway induced by ABA and JA in apple plants

The external application of acetic acid (AA) has been shown to improve drought survival in plants, such as Arabidopsis, rice, maize, wheat, rapeseed and cassava, and the application of AA also increased drought tolerance in perennial woody apple (Malus domestica) plants. An understanding of AA-induced drought tolerance in apple plants at the molecular level will contribute to the development of technology that can be used to enhance drought tolerance. In this study, the morphological, physiological and transcriptomic responses to drought stress were analyzed in apple plants after watering without AA (CK), watering with AA (AA), drought treatment (D) and drought treatment with AA (DA). The results suggested that the AA-treated apple plants had a higher tolerance to drought than water-treated plants. Higher levels of chlorophyll and carotenoids were found under the DA conditions than under D stress. The levels of abscisic acid (ABA), jasmonic acid (JA) and methyl jasmonate were increased in AA-treated apple plants. Transcriptomic profiling indicated the key biological pathways involved in metabolic processes, mitogen-activated protein kinase (MAPK) signaling, plant hormone signal transduction and the biosynthesis of secondary metabolites in response to different drought conditions. The 9-cis-epoxycarotenoid dioxygenase, (9S,13S)-cis-oxophytodienoic acid reductase, allene oxide synthase, allene oxide cyclase and lipoxygenase genes participate in the synthase of ABA and JA under drought and AA treatments. Collectively, the results showed that external application of AA enhanced drought tolerance in apple plants by influencing the ABA- and JA-induced MAPK signaling pathways. These data indicated that the application of AA in plants is beneficial for enhancing drought tolerance and decreasing growth inhibition in agricultural fields.

Arbuscular mycorrhizal fungi improve the growth and drought tolerance of Cinnamomum migao by enhancing physio-biochemical responses

Drought is the main limiting factor for plant growth in karst areas with a fragile ecological environment. Cinnamomum migao H.W. Li is an endemic medicinal woody plant present in the karst areas of southwestern China, and it is endangered due to poor drought tolerance. Arbuscular mycorrhizal fungi (AMF) are known to enhance the drought tolerance of plants. However, few studies have examined the contribution of AMF in improving the drought tolerance of C. migao seedlings. Therefore, we conducted a series of experiments to determine whether a single inoculation and coinoculation of AMF (Claroideoglomus lamellosum and Claroideoglomus etunicatum) enhanced the drought tolerance of C. migao. Furthermore, we compared the effects of single inoculation and coinoculation with different inoculum sizes (20, 40, 60, and 100 g; four replicates per treatment) on mycorrhizal colonization rate, plant growth, photosynthetic parameters, antioxidant enzyme activity, and malondialdehyde (MDA) and osmoregulatory substance contents. The results showed that compared with nonmycorrhizal plants, AMF colonization significantly improved plant growing status; net photosynthetic rate; superoxide dismutase, catalase, and peroxidase activities; and soluble sugar, soluble protein, and proline contents. Furthermore, AMF colonization increased relative water content and reduced MDA content in cells. These combined cumulative effects of AMF symbiosis ultimately enhanced the drought tolerance of seedlings and were closely related to the inoculum size. With an increase in inoculum size, the growth rate and drought tolerance of plants first increased and then decreased. The damage caused by drought stress could be reduced by inoculating 40-60 g of AMF, and the effect of coinoculation was significantly better than that of single inoculation at 60 g of AMF, while the effect was opposite at 40 g of AMF. Additionally, the interaction between AMF and inoculum sizes had a significant effect on drought tolerance. In conclusion, the inoculation of the AMF (Cl. lamellosum and Cl. etunicatum) improved photosynthesis, activated antioxidant enzymes, regulated cell osmotic state, and enhanced the drought tolerance of C. migao, enabling its growth in fragile ecological environments.

Genome-Wide Identification and Analysis of the NF-Y Transcription Factor Family Reveal Its Potential Roles in Salt Stress in Alfalfa (Medicago sativa L.)

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor that plays an important role in various biological processes in plants, such as flowering regulation, drought resistance, and salt stress. However, few in-depth studies investigated the alfalfa NF-Y gene family. In this study, in total, 60 MsNF-Y genes, including 9 MsNF-YAs, 26 MsNF-YBs, and 25 MsNF-YCs, were identified in the alfalfa genome. The genomic locations, gene structures, protein molecular weights, conserved domains, phylogenetic relationships, and gene expression patterns in different tissues and under different stresses (cold stress, drought stress, and salt stress) of these NF-Y genes were analyzed. The illustration of the conserved domains and specific domains of the different subfamilies of the MsNF-Y genes implicates the conservation and diversity of their functions in alfalfa growth, development, and stress resistance. The gene expression analysis showed that 48 MsNF-Y genes (7 MsNF-YAs, 22 MsNF-YBs, and 19 MsNF-YCs) were expressed in all tissues at different expression levels, indicating that these genes have tissue expression specificity and different biological functions. In total, seven, seven, six, and eight MsNF-Y genes responded to cold stress, the ABA treatment, drought stress, and salt stress in alfalfa, respectively. According to the WGCNA, molecular regulatory networks related to salt stress were constructed for MsNF-YB2, MsNF-YB5, MsNF-YB7, MsNF-YB15, MsNF-YC5, and MsNF-YC6. This study could provide valuable information for further elucidating the biological functions of MsNF-Ys and improving salt tolerance and other abiotic stress resistance in alfalfa.

Seed-Borne Erwinia persicina Affects the Growth and Physiology of Alfalfa (Medicago sativa L.)

Seed-borne Erwinia persicina can be transmitted globally via alfalfa (Medicago sativa L.) seed trade, but there is limited information about the impact of this plant-pathogenic bacterium on alfalfa plants. In this study, strain Cp2, isolated from alfalfa seeds, was confirmed by whole-genome sequencing to belong to E. persicina. Subsequently, the effects of Cp2 on alfalfa growth and physiology were evaluated by constructing a rhizosphere infection model. Strain Cp2 had a strong inhibitory effect on the elongation and growth of alfalfa roots, which was very unfavorable to these perennial plants. Furthermore, an increased number of leaf spots and yellowing symptoms were observed in plants of the Cp2 group from day 10 to day 21 and the strain Cp2 was re-isolated from these leaves. Correlation between growth and photosynthetic parameters was analyzed and the significant decreases in fresh weight and root and plant lengths in the Cp2 group were related to the marked reduction of chlorophyll b, carotenoid, transpiration rate, and stomatal conductance of leaves (r > 0.75). In addition, nine physiological indicators of root, stem, and leaf were measured in the plants 21 days after treatment with Cp2. The physiological response of root and leaf to Cp2 treatment was stronger than that of stem. The physiological indicators with the greatest response to Cp2 infection were further explored through principal component analysis, and superoxide dismutase, peroxidase, ascorbate peroxidase, and soluble protein showed the greatest changes in roots, stems, and leaves (P < 0.001). Among tissues, the commonality was the change of soluble protein. Therefore, soluble protein is speculated to be a physiological marker during alfalfa–E. persicina interactions. These findings indicate that once E. persicina spreads from alfalfa seeds to the rhizosphere, it can invade alfalfa roots and cause disease. This study demonstrates that this plant pathogenic bacterium may be a potential threat to new environment when they spread via seed trade and these “dangerous hitchhikers” warrant further attention, especially in the study of bacterial diseases in pasture-based production systems.

Genome-wide identification of MAPK gene family members in Fagopyrum tataricum and their expression during development and stress responses

BACKGROUND

Mitogen-activated protein kinases (MAPKs) plays essential roles in the development, hormone regulation and abiotic stress response of plants. Nevertheless, a comprehensive study on MAPK family members has thus far not been performed in Tartary buckwheat.

RESULTS

Here, we identified 16 FtMAPKs in the Fagopyrum tataricum genome. Phylogenetic analysis showed that the FtMAPK family members could be classified into Groups A, B, C and D, in which A, B and C members contain a Thr-Glu-Tyr (TEY) signature motif and Group D members contain a Thr-Asp-Tyr (TDY) signature motif. Promoter cis-acting elements showed that most ProFtMAPks contain light response elements, hormone response elements and abiotic stress response elements, and several ProFtMAPks have MYB-binding sites, which may be involved in the regulation of flavonoid biosynthesis-related enzyme gene expression. Synteny analysis indicated that FtMAPKs have a variety of biological functions. Protein interaction prediction suggested that MAPKs can interact with proteins involved in development and stress resistance. Correlation analysis further confirmed that most of the FtMAPK genes and transcription factors involved in the stress response have the same expression pattern. The transient transformation of FtMAPK1 significantly increased the antioxidant enzymes activity in Tartary buckwheat leaves. In addition, we also found that FtMAPK1 can respond to salt stress by up-regulating the transcription abundance of downstream genes.

CONCLUSIONS

A total of 16 MAPKs were identified in Tartary buckwheat, and the members of the MAPK family containing the TDY motif were found to have expanded. The same subfamily members have relatively conserved gene structures and similar protein motifs. Tissue-specific expression indicated that the expression of all FtMAPK genes varied widely in the roots, stems, leaves and flowers. Most FtMAPKs can regulate the expression of other transcription factors and participate in the abiotic stress response. Our findings comprehensively revealed the FtMAPK gene family and laid a theoretical foundation for the functional characterization of FtMAPKs.

Transcriptomic Analysis Reveals Regulatory Networks for Osmotic Water Stress and Rewatering Response in the Leaves of Ginkgo biloba

To elucidate the transcriptomic regulation mechanisms that underlie the response of Ginkgo biloba to dehydration and rehydration, we used ginkgo saplings exposed to osmotically driven water stress and subsequent rewatering. When compared with a control group, 137, 1453, 1148, and 679 genes were differentially expressed in ginkgo leaves responding to 2, 6, 12, and 24 h of water deficit, and 796 and 1530 genes were differentially expressed responding to 24 and 48 h of rewatering. Upregulated genes participated in the biosynthesis of abscisic acid, eliminating reactive oxygen species (ROS), and biosynthesis of flavonoids and bilobalide, and downregulated genes were involved in water transport and cell wall enlargement in water stress-treated ginkgo leaves. Under rehydration conditions, the genes associated with water transport and cell wall enlargement were upregulated, and the genes that participated in eliminating ROS and the biosynthesis of flavonoids and bilobalide were downregulated in the leaves of G. biloba. Furthermore, the weighted gene coexpression networks were established and correlated with distinct water stress and rewatering time-point samples. Hub genes that act as key players in the networks were identified. Overall, these results indicate that the gene coexpression networks play essential roles in the transcriptional reconfiguration of ginkgo leaves in response to water stress and rewatering.