Integrated Physiological and Transcriptomic Analyses Responses to Altitude Stress in Oat (Avena sativa L.)

Oat is an annual gramineous forage grass with the remarkable ability to survive under various stressful environments. However, understanding the effects of high altitude stresses on oats is poor. Therefore, the physiological and the transcriptomic changes were analyzed at two sites with different altitudes, low (ca. 2,080 m) or high (ca. 2,918 m), respectively. Higher levels of antioxidant enzyme activity, reactive oxygen and major reductions in photosynthesis-related markers were suggested for oats at high altitudes. Furthermore, oat yields were severely suppressed at the high altitude. RNA-seq results showed that 11,639 differentially expressed genes were detected at both the low and the high altitudes in which 5,203 up-regulated and 6,436 down-regulated. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment tests were conducted and a group of major high altitude-responsive pigment metabolism genes, photosynthesis, hormone signaling, and cutin, suberine and wax biosynthesis were excavated. Using quantitative real-time polymerase chain response, we also confirmed expression levels of 20 DEGs (qRT-PCR). In summary, our study generated genome-wide transcript profile and may be useful for understanding the molecular mechanisms of Avena sativa L. in response to high altitude stress. These new findings contribute to our deeper relevant researches on high altitude stresses and further exploring new candidategenes for adapting plateau environment oat molecular breeding.

Introduction

Extreme environments provide natural laboratories for studies on the processes of speciation and adaptive evolution of organisms (Qiao et al., 2015). The Qinghai-Tibetan Plateau, the highest plateau in the world, plays an important role in climate changes in Asia and even the world (Chi et al., 2019). In particular, the Qinghai-Tibetan Plateau environment is known for its harsh conditions, characterized by severe coldness, intensive ultraviolet radiations, hypoxia, poor soils and low CO₂ pressure (Xin et al., 2011; Jia et al., 2017). Therefore, the survival in Qinghai-Tibetan Plateau is very challenging for most organisms. Nevertheless, many plant species can thrive in the cold and hypoxic conditions in high-alpine areas (Storz et al., 2010; Yang et al., 2015). Generally, in response to the bioclimatic conditions, plants can make corresponding molecular and physiological changes (Thiebaut et al., 2019). For example, the high-altitudinal gradient can restrict plant growth and reproduction via strong solar UV-B radiation, resulting in a reduction of photosynthetic rates by bleaching chlorophyll a (Chl a) and damaging the photosynthetic apparatus (Zhu and Yang, 2015). As such, some alpine plants are able to activate antioxidants such as APX, CAT, GR, proline and abscisic acid to confer plants with the tolerance to the alpine environments (Ch, 2005; Vinocur and Altman, 2005; Li et al., 2014).

In recent years, transcriptome sequencing has been proven to be an effective and efficient method for determining adaptive evolutions and differential gene expressions in high-altitude plants, e.g., Kobresia pygmaea, Potentilla saundersiana, Lamiophlomis rotate, and Lobelia. Some studies demonstrated that alpine plants had various morphological and physiological response strategies to adapt to high-elevation environments. For example, researchers had revealed that some genes were significantly involved in energy metabolism, hypoxia response under positive selection and rapid adaptation (Ma et al., 2015a,b; Zhao et al., 2016). However, up to now, few transcriptome-based investigations have been devoted to the molecular mechanisms of high-altitude adaptation in oats.

Oats have not only high yield, rich nutrition and good palatability, but also the characteristics of cold tolerance, drought resistance and strong adaptability. Oats can germinate when the temperature is 3∼4°C and seedlings can resist the low temperature of −3∼−4°C. Oats are the main forage species in the Qinghai-Tibet Plateau where they are high-quality forages for both winter and spring feeding and disaster-preserving livestock in local farming and pastoral areas. Planting oat grass is of great significance for solving short pasture supply in winter and spring grassland and developing herbivorous livestock and animal husbandry (Sharma et al., 2016). The oat has also been identified as an excellent species that can adapt to various environmental stresses including drought, salinity and pathogen attacks. Oats make morphological and physiological changes to adapt to environments under abiotic stress (Wu et al., 2017; Anwar et al., 2018). However, the adaptation mechanisms of oats in high altitude environments have been rarely studied.

In this study, the first step was designed to characterize the impact of altitudes on oat physiology and agriculture. It showed that high altitude influenced photosynthetic ability, antioxidant enzyme activity, stoma aperture and agricultural parameters of oats. The second step, RNA-seq was employed to study the changes at gene levels between low and high altitudes. According to extensive data analyses, many DEGs and metabolic pathways were identified and characterized which were involved in adaption of oats to high altitude stresses.

Results

High Altitude Environments Induced Physiological Changes in Oats

H₂O₂ and O₂^– are essential signals of reactive oxygen species (ROS). They were increased significantly when oats were exposed to high altitude stress (Figures 1A,B). Due to the overproduction of ROS in the cell, MDA is an important marker for lipid peroxidation. The MDA behaviors of high altitude oats were higher than those of low altitude oats (Figure 1C). To further investigate the antioxidant safety mechanism for oats, antioxidant enzyme functions have been investigated. The SOD, CAT, and GR behaviors of the enzymes improved dramatically in high altitude conditions relative to those at low altitude, except for APX (Figures 1D–G). However, the high altitude caused the aggregation of proline osmoregulation substances and soluble sugars in oats (Figures 1H,I).

FIGURE 1

Figure 1. Physiological changes of oats at low and high altitudes. (A) H₂O₂ activity, (B) O₂^– activity, (C) MDA activity, (D) SOD activity, (E) CAT activity, (F) APX activity, (G) GR activity, (H) proline content, (I) soluble sugars content, (J) chlorophyll content, (K) NPQ, (L) stomata density. Values represented the means (± SE) from three fully independent biological replicates. Significant differences from the low altitude were denoted by one or two asterisks corresponding to P < 0.05 and P < 0.01, respectively, by the Student’s t-test. (M) The photographs of stomata of oats at low and high altitudes.

The chlorophyll content and NPQ at the high altitude decreased and were significantly lower than those at the low altitude (Figures 1J,K). This result indicated that there was a large photosynthetic apparatus damage in high altitude conditions. Stomata control the carbon dioxide absorption and play a crucial role in photosynthesis. This result suggested that the stoma density had significantly decreased at high altitude (Figures 1L,M).

The Agronomic Traits of Oats Were Remarkably Affected by High Altitude

Oats displayed significant differences from two altitudes in terms of plant height, stem/leaf ratio, crude fat and total hay yield. Under high altitude conditions, the plant height of oats was 104.7 cm which was significantly shorter than 126.9 cm at the low altitude (Figure 2A). The stem/leaf ratio and crude fat at the low altitude were 1.18 and 2.7%, respectively, which were significantly different from 1.05 to 1.7% at the high altitude (Figures 2B,C). On the other hand, compared with the low altitude, the high altitude exhibited remarkable decreases in hay yield (Figure 2D). It was clear that oats altered their agronomic traits in two different regions.

FIGURE 2

Figure 2. Agronomic trait changes of oats at low and high altitudes. (A) plant height, (B) crude protein, (C) crude fat, (D) hay yield. Means were compared using one-way ANOVA. All parameters were shown as the mean ± SE. Significant differences from low altitude were denoted by one or two asterisks corresponding to P < 0.05 and P < 0.01, respectively by the Student’s t-test.

Transcriptome Profiling

To identify the transcriptomes and gene expression profiles of oats at low and high altitudes, six cDNA samples from oat leaves were prepared and sequenced using Illumina HiSeq 2000 platform. A total of 47 50 million raw reads were harvested from each cDNA library. After the low quality reads were removed, 39.49 Gb clean reads were obtained with an average of 6.66 Gb reads for each sample. The percentage of bases higher than Q30 in each sample was not less than 90.91% (Table 1).

TABLE 1

Table 1. Overview of the RNA-sequencing reads generated from each sample.

Using Trinity program, these clean reads were then assembled de novo, and a cluster deduplication was performed to get the final unigene for a subsequent analysis designated all-unigene. As a consequence, with a mean length of 950, a total of 274,639 All-unigenes were created. Meanwhile, with a GC content of 49.93%, a N50 value of 1,635 bp (Table 2).

TABLE 2

Table 2. Trinity statistics.

To predict and analyze the functions of the unigenes, we carried out functional annotations by using BLAST against multiple databases such as Nr, Nt, KEGG, Swissprot, KOG, Interpro, and GO. A total of 94,954 unigenes were successfully matched up in the databases mentioned above. Among them, 11,507 (12.12%), 38,504 (40.55%), 12,044 (12.68%), 11,801 (12.43%), 12,154 (12.80%), 15.199 (16.00%), 49,383 (52.00%) unigenes were found in Nr, Nt, KEGG, Swissprot, KOG, Interpro, and GO databases, respectively (Supplementary Table S1). Moreover, a species distribution map was drawn (Figure 3). The results indicated that oats showed the most similar to Aegilops tauschii subsp. tauschii, followed by Brachypodium distachyon and Hordeum vulgare subsp. vulgare.

FIGURE 3

Figure 3. Similarities of Oats to other species.

The Analyses of Differentially Expressed Genes (DEGs)

To obtain the differential expression genes response to high altitude, the Fragments Per Kilobase of transcript, FPKM method was used to analyze the expression abundance of unigenes. As a results, a total of 11,639 differentially expressed genes (|log 2 Fold Change > 2|) were differentially expressed between the low and the high altitudes (Supplementary Table S2 and Figure 4A). Among them, 5,203 DEGs were up-regulated and 6,436 DEGs were down-regulated (Figure 4B). Interestingly, the number of down-regulation genes were more than the up-regulation ones. It indicated that more genes were down-regulated in response to high altitude stresses.

FIGURE 4

Figure 4. Venn diagram and volcano plot of DEGs at low and high altitudes. (A) Venn diagram showed DEGs expressed at each of the two altitude treatments. (B) Volcano plot showed the numbers of DEGs identified in comparisons between pairs of libraries.

GO and KEGG Analyses of DEGs

To investigate the functions of these DEGs under high altitude stresses, we performed GO annotation enrichment analysis using the software package top GO (Supplementary Table S3 and Figure 5A). The results showed that 2,879 transcripts were characterized and annotated with the top GO terms for classified genes being “cellular process” (1,956 transcripts) for biological processes. Other categories with high numbers of genes included “metabolic process” (1,819 transcripts) and “biological regulation” (887 transcripts). These showed that substance synthesis and catabolism, cellular processes, and physiological regulation responses to external stresses were involved in responses of oats to high altitude environmental stresses. In addition, 3,308 transcripts were mapped to interior cell components in which “cell” (2,348 transcripts) and “cell part” (2,314 transcripts) were the most represented categories. The mechanism by which oats responding to high altitude was closely related to the structural composition of their cells, cell membranes and organelles. As for molecular functions, 4,011 transcripts were significantly annotated with the “binding” (2,533 transcripts) and “catalytic activity” (2,478 transcripts) pathways having the most DEGs, respectively. These indicated that the non-covalent bonding of molecules and the catalytic activity of enzymes in biochemical reactions played important roles in the physiological adaptability of oats to high altitude stresses.

FIGURE 5

Figure 5. Functional annotation of assembled transcriptome. (A) Gene ontology (GO) classifications of the annotated DEGs (up- and down-regulated). (B) Kyoto Encyclopedia of Genetics and Genomics (KEGG) database analyses of DEGs (up- and down-regulated) enriched in different biological pathways. The X-axis represented enriched pathways and Y-axis represented the total number of transcripts.

At the same time, the DEGs and the biological processes involved in altitude stresses were illuminated and all of the DEGs were queried in KEGG databases. In this study 4,831 transcripts were allocated to 30 pathways in KEGG databases (Supplementary Table S4 and Figure 5B) while “Biosynthesis of secondary metabolites” (1,060 transcripts), “plant-pathogen interaction” (541 transcripts) and “RNA transport” (504 transcripts) pathways were mostly enriched followed by “phenylpropanoid biosynthesis” (390 transcripts), “plant hormone signal transduction” (259 transcripts), “MAPK sigaling pathway” (253 transcripts), and “starch and sucrose metabolism” (215 transcripts). Interestingly most of the genes mapped in the first seven significantly enriched pathways had the trend in down-regulation. These results above indicated that various pathways were involved in responses of oats to different altitudes.

Analyses of Transcription Factors (TF) in the DEGs

TFs control the expression of numerous genes and play crucial functions in stress-induced signal transduction pathways. In this study, 125 TF genes were differentially expressed between two different altitudes (Supplementary Table S5). Among them, 53 bHLH transcripts, 36 AP2-EREBP transcripts, 11 C₂C₂-CO-like transcripts were the top three most TF families. The rest of the transcription factors were ABI3VP1, C₂C₂-Dof, C₂C₂-GATA, bZIP, C₂H₂, ARF, ARR-B, BSD, and C₂C₂-YABBY with the numbers of transcripts were 6, 6, 3, 3, 3, 1, 1, 1, and 1, respectively. These results indicated that the TFs above played critical roles in oats response to high altitude stresses.

Analyses of DEGs Related to Altitudinal Responses

To further identify the key genes related to responses at different altitudes, the DEGs in two pathways for the chlorophyll metabolism and the carotenoid biosynthesis were compared in detail between the low and the high altitudes. The results demonstrated that 23 genes in the porphyrin and the chlorophyll metabolism pathways were detected. Among them, 14 genes had been up-regulated and only 9 genes were down-regulated at high altitude. The DEGs played important roles during the chlorophyll biosynthesis and encoding chlorophyllase (Supplementary Table S6 and Table 3). Moreover, the comparison of the DEGs involved in the carotenoid biosynthesis showed that 32 genes were significantly down-regulated at high altitude (Supplementary Table S7). The results indicated that the photosynthetic pigment biosynthesis was an important process for oats under the high altitude stress.

TABLE 3

Table 3. DEGs associated with pigment metabolism.

High altitude stresses might have an impact on plant photosynthesis. In this study, we identified 15 DEGs were significantly regulated in photosynthesis which played a part in photosystem I and II. Meanwhile, 9 genes that mapped to the photosynthesis-antenna protein pathways exhibited differential expressions at high altitude, where the genes encoded light-harvesting complex II chlorophyll a/b binding protein (Supplementary Table S8 and Table 4). These results indicated that changes at the expression levels of the genes above might have blocked photosynthesis, thereby influencing the chlorophyll biosynthesis of oats at high altitude.

TABLE 4

Table 4. DEGs associated with photosynthesis.

Hormones are pivotal to plants in stress adaptive signaling cascades and act as central integrators to connect and reprogram different responses. In this experiment, we identified 273 hormone-related DEGs between the low and the high altitude conditions. Among them, 60 genes were identified to encode DELLA protein and acted as a part of gibberellin receptor and phytochrome-interacting factor. Forty-eight genes were identified in the auxin signal transduction pathway in which some encoded auxin-responsive protein IAA and others participated in the process of auxin synthesis and metabolism. In jasmonic acid signal transduction pathway, 41 genes encoded jasmonate ZIM domain-containing protein and transcription factor MYC2. The three pathways above were the most abundant phytohormone signal transduction pathways. In addition, 33 genes in cytokinin pathway, 31 genes in abscisic acid signal pathway, 24 genes in brassinosteroid pathway, 22 genes in ethylene pathway, and 14 genes in salicylic acid pathway were differentially expressed (Supplementary Table S9). The results above suggested that these hormone-related genes can be involved in cell division, stem elongation, stomata aperture, fruit maturation, and disease prevention processes in response to high altitude conditions.

Furthermore, DEGs pathway enrichment study determined that the mechanism of cutin, suberine, and wax biosynthesis was considerably enriched in several genes. Cuticle wax, outside the cell wall of the plant epidermis, is an important hydrophobic structure. Its central aim is to play an important role in plant resistance to abiotic stresses such as salt and heat, but in the cutin, suberine and wax biosynthesis pathways involved in high altitude stress, a total of 32 DEGs were annotated (Table 5). Among them, 2 down-regulated and 1 up-regulated omega-hydroxylase and peroxygenase fatty acid encoding genes were found. We noticed that 10 DEGs of aldehyde decarbonylase and alcohol-forming fatty acyl-CoA reductase were all highly regulated. In addition, 5 (3 up regulated and 2 down regulated) and 14 (4 up regulated and 10 down regulated) omega-hydroxypalmitate O-feruloyl transferase and wax-ester synthase/diacylglycerol O-acyltransferase genes differentially expressed were identified. In addition to high altitude stress in oats, the above findings showed that cutin, suberine and wax biosynthesis played important roles.

TABLE 5

Table 5. DEGs associated with cutin, suberine and wax biosynthesis.

Relative DEGs Measured by qRT-PCR

We randomly performed qRT-PCR of 20 genes involved in altitude stress to further check the accuracy of RNA-seq data and validate the extent of differential expression genes, as shown in Figure 6A. The transcript abundance of these genes (log 2 fold change) is shown in Supplementary Table S10. Real-time PCR revealed that the expression patterns of 20 low and high altitude genes were compatible with the sequencing findings of expression profiles. A strong correlation (R² = 0.88) was observed between the two methods, suggesting that the relevance of RNA-seq data and qRT-PCR showed good accuracy between the two methods (Figure 6B).

FIGURE 6

Figure 6. The relative expression levels of 20 DEGs identified in the comparison between RNA-Seq and qRT-PCR. (A) The gene relative expression levels were normalized to the expression level of UBC. Error bars represented the standard deviations of three PCR replicates. Significant differences from the low altitude were denoted by one or two asterisks corresponding to P < 0.05 and P < 0.01, respectively, by the Student’s t-test. (B) R-values are the correlation coefficients between qRT-PCR and RNA-seq.

Discussion

The Qinghai-Tibet Plateau itself has unique climate and geographical features which pose huge challenges to the survival of plants. Although plants cannot escape from unfavorable environments, they can change their morphological structures, physiological and molecular processes in which changes happen in endogenous secondary metabolite levels, ion homeostasis, gene expressions and protein activity regulations to adapt to different environmental conditions. Oats are an important source of food and feed in the Qinghai-Tibet Plateau and are beneficial to human health because they are rich in protein and minerals (Noga et al., 2016). Now the data available on the molecular basis of oats response to Qinghai-Tibet Plateau conditions stress is not limited any longer with the development of high-throughput sequencing.

The Adaptation of the Oat Physiological Characteristics to High Altitude Environments

The increased levels of protective antioxidant enzymes could maintain the homeostasis of ROS. In this study, oats could regulate many physiological processes to adapt to high altitude stresses, e.g., increasing the activities of SOD (Figure 1D) which suggested that oats can improve their ability to resist damages by establishing a better antioxidant enzyme system under high altitude stresses. The membrane lipid peroxidation product MDA was also one of the important indicators of cell damages which could bind to proteins and cause membrane protein denaturation (Yang et al., 2017). Our result showed that oats suffered more severe damages to the cell membrane and antioxidant system in high altitude environments than those at low altitude (Figure 1C). In an alpine environment, low temperatures could affect the osmotic pressure of plant cells. As shown in Figure 1I, high altitude stresses induced the accumulation of soluble sugars in oats. Therefore, soluble sugars might play important roles in oat growth and development and responses to abiotic stresses, and eventually affect yield and quality. In addition, sugars acted as not only an essential source of carbon but also signaling molecules in response to the integration of information from environments as well as developmental and metabolic cues in plants (Lilley et al., 2012; Lastdrager et al., 2014; Ma et al., 2017). Therefore, it could further improve the tolerance of oats under low temperatures at high altitude by increasing the content of soluble sugars to avoid protein coagulation. Meanwhile, oat plants showed sharp decreases in chlorophyll content and NPQ at high altitude (Figures 1J,K). The decreases of chlorophyll content might be because the enhancement of radiation could exacerbate the photooxidation process of photosynthetic pigments or restrain some synthetic processes to decrease the light energy absorbed by leaves at high altitude to avoid or reduce the radiation damages on plants (Lizana et al., 2009). The decrease in NPQ indicated that the plant did not dissipate excess excitation energy in the form of heat, thereby increased the damage on PSII by the accumulation of the excitation energy under high altitude stresses (Costa, 2003). Moreover, the chlorophyll content and NPQ as photosynthesis markers were closely linked to yield which was probably ascribed to the reductions of both plant height and the ratio of stem/leaf, subsequently led to marked decreases in hay yield index (Figure 2). In addition, high altitude conditions could regulate the development of oat plants, e.g., decreased the stomatal density to reduce water loss to adapt to high altitude stress (Mahajan and Tuteja, 2005). This mightbe due to the water shortage caused by low temperatures which was also the result of adaptation to low temperature and strong radiation in high mountains.

TFs Responding to High Altitude Stresses

When subjected to abiotic, TFs could play an important role by binding to cis-regulated DNA elements in gene promoters and activating or inhibiting their expressions (Singh and Laxmi, 2015). bHLH had been reported to directly bond the phytochrome action factor in the photoreceptor signaling network which specifically binded to phytochrome responsive gene promoter G-box regulatory elements in vitro (Huq, 2002). In addition, overexpression of Arabidopsis bHLH TFs enhanced chilling or osmotic tolerance by reducing malondialdehyde content and accumulating ROS (Liu et al., 2014; Zhao Q. et al., 2018). AP2-EREBP and ABI3 TF could participate in the process of plant growth, development and respond to stresses by regulating the expression of downstream target genes (Liu and Zhang, 2017; Zhao M. et al., 2018). In our study, bHLH, AP2-EREBP, and ABI3VP1 were the three most abundant transcription factors identified as DEGs in response to high altitude stresses. Moreover, many transcription factors containing C₂C₂ zinc finger domain were found including C₂C₂-CO-like, C₂C₂-Dof, C₂C₂-GATA, and C₂C₂-YABBY. Researches had shown that they played vital biological functions in vegetative tissue development and multiple stress responses (Ge et al., 2011). It could be inferred that bHLH, AP2-EREBP, ABI3VP1, and C₂C₂ TFs might play important roles in response to the high altitude stress. In addition, many other TF families were also identified to be differentially expressed such as bZIP, C₂H₂, ARF, ARR, and BSD. This also explained that numerous TFs participated actively in the regulation of oats in response to high altitude conditions through different pathways.

High Altitude Affected the Expression of Genes Related to Pigment Metabolism and Photosynthesis

Plant pigments determine color variations of leaves which are complex biological processes. Among them, chlorophyll and carotenoid are two important antenna pigments in the photosynthetic system. Chlorophyll plays a central role in the process of absorbing and transmitting energy (Tanaka and Tanaka, 2006). Carotenoids can influence chlorophyll content as a major component of the photosynthetic apparatus and photoprotection system. Plateau plants can adapt to UV-B radiation by changing the content of photosynthetic pigments in leaves such as reducing the chlorophyll and carotenoid content (Sicora, 2006). Chlorophyll and carotenoid metabolisms are complex processes involving multiple enzymes. Comparative transcriptome profiling data for the low altitude and the high altitude showed that the expression of numbers of DEGs were significantly regulated including uroporphyrinogen decarboxylase, copropor phyrinogen III oxidase, protochlorophyllide reductase, chlorophyll b reductase, and cytochrome c oxidase assembly protein subunit 15, which were key enzymes associated with the process of chlorophyll synthesis (Mazor et al., 2015). Meanwhile, our comparison of the DEGs involved in carotenoid biosynthesis showed that 32 genes were significantly down-regulated in high altitude (Supplementary Table S3). It might imply that oats can adapt to UV-B radiation stress by changing the content of chlorophyll and carotenoids in leaves.

Chloroplasts are the main sites for photosynthesis which are comprised of membranes, thylakoid system, and stroma with functional units in it. In higher plants, the multi-subunit pigment–protein complexes (i.e., photosystem (PS) I, PSII, light harvesting complexes, cytochrome b6/f, and ATP synthase) are embedded in the highly folded thylakoid membrane where they are responsible for light absorption and energy transfer (Bashir, 2011). Strong radiation in high-altitude environments might have been related to thylakoid development and photosynthesis at the two different altitudes. Based on transcriptome data it was found that genes encoding photosystem II CP47 chlorophyll apoprotein, photosystem II oxygen-evolving enhancer protein, photosystem I P700 chlorophyll a apoprotein and photosystem I subunit PsaO were significantly regulated at high altitude. In higher plants, light-harvesting complex II chlorophyll a/b binding proteins that accumulate in the thylakoid membranes, where their proposed function is in photoprotection (Hutin, 2003; Beck et al., 2017; Ozawa et al., 2018). The chlorophyll a/b-binding proteins are correlated with photodamage in the PSII reaction centers (Heddad et al., 2006). We found that the proteins were significantly enriched suggesting that strong photosynthetic apparatus damage occur in high altitude (Table 3). These results were consistent with our previous photosynthetic physiological index analysis (Figure 2B) which further demonstrated that the genes related to pigment metabolism and photosynthesis in oats could respond to high altitude stresses.

Identification of Signal Transduction Genes Responsible for High Altitude Stress

The different expression of genes involved in the signal cascade mechanism can affect the expression of genes participating in the formation of plant hormones such as auxin, ABA, ET, SA, and JA. Theses hormones may amplify the cascades or initiate some new signaling pathways (Sarwat et al., 2013; Wang et al., 2018). It is now known that these phytohormones extensively regulate all aspects of plant stress responses ranging from signal cascade transduction to modifications in plant developmental processes (Shu et al., 2018). In this experiment, GA, IAA, and JA were the top three responsive hormone signal transduction pathways in oats and there were a number of DEGs involved in them. Studies had shown that the gibberellin receptor sensing the GA signal could activate the signaling pathway, thereby regulating the expression of downstream genes to affect plant growth and morphogenesis (Ito et al., 2018). Totally 60 DEGs involved in the response to GA signal transduction pathway were identified as being differentially expressed, which encoded gibberellin receptor GID and DELLA proteins. This indicated that high altitude-induced expression of oats appeared to be GA dependent.

Auxin is a plant hormone extensively involved in regulating cell growth, cell division and cell-specific differentiations to participate in a series of growth and development processes (Tuan et al., 2019). ARF and SAUR family protein could bind to the cis-acting element of auxin-inducible gene and regulate its transcriptional expression (Spartz et al., 2014; Zhang et al., 2019). Forty-eight DEGs encoding three key proteins (transport inhibitor response, auxin-responsive and SAUR family protein) of IAA signaling pathway were identified in our study. It could infer that both ABF and SAUR can activate the down-stream gene expressions to respond to high altitude stresses.

JAs are a group of important hormones that play pivotal roles in a variety of plant growth, developmental and stress response processes (Sharma and Laxmi, 2015). In this study, 41 DEGs encoding jasmonate ZIM domain-containing protein and transcription factor MYC2 were enriched in the JA signaling pathway. It had been shown that JAZ gene and MYC2 were the major regulators of jasmonic acid signaling pathway in which JAZ protein degradation could release MYC2 transcription factor and activate the expression of downstream genes of jasmonic acid signaling pathway (Boter et al., 2004).

High Altitude Resistant Genes Associated With Cutin, Suberine, and Wax Biosynthesis

The cuticle protects the surface of plant organs and plays a crucial role in the contact between plants and the environment. High altitude tolerant genes associated with cutin, suberine, and wax biosynthesis In many plants under salt stress, cuticle wax content is strongly negatively associated with cutin transpiration, suggesting that a large volume of cutin wax accumulation forms a protective barrier for controlling cutin transpiration, maximizing water usage quality, and improving the resistance of plant salt (Kosma et al., 2009; Hasanuzzaman et al., 2017). Plants can respond by increasing or decreasing the content of wax to high temperature stress. If the wax content of the cuticle is increased or decreased by high temperature stress is still inconclusive, which is related to the form of plant (Awika et al., 2017; Huggins et al., 2018). A significant number of cutin, suberine and wax biosynthesis genes were expressed differentially at elevated altitudes in Table 5. The results of the study on cuticle wax involved in the abiotic stress mechanism are clear, but there are few studies on high stress at altitude. Therefore, in reaction to high altitude stress, we investigated the molecular function of cuticle wax. The deposition of wax may be caused by stomatal growth to increase plant tolerance, suggesting that the improvements in wax synthesis are closely connected to irregular epidermal stomata. We found in our research that the density of stomata has changed significantly (Figure 1L). We therefore speculate that oats will change the density of stomata to regulate wax biosynthesis, further change the loss of stomata and non-stomata water, and ensure that the water content adjust to the condition of the plateau. This presents new ideas for the study under high altitude stress of the physiological and molecular processes of oats.

Materials and Methods

Plant Materials and Growth Conditions

Seeds of Avena sativa L. Qingyin No. 1 were collected from the Qinghai with the permission of the Qinghai Academy of Animal Science and Veterinary Medicine in the Qinghai Province of China. Animal Science and Veterinary Medicine in the Qinghai Province of China undertook the formal identification of the samples and provided details of specimen deposited. In order to investigate the effects of altitudes on the physiology and the transcriptions, oat leaf samples of the cultivar (cv.) “Qingyin No.1” were sown and collected at heading stage, 2015 in two regions with significant differences in altitude in Qinghai Province, China, County of Minhe, (36°14′, 102°41′) at the altitude ca. 2,080 m (defined as the low altitude) and County of Huangzhong (36°21′,101°44′) at ca. 2,918 m (defined as the high altitude). During the heading period, the climate data in the two altitude regions were recorded (Supplementary Table S11). Compared with low-altitude areas, high-altitude areas have longer sunshine duration and lower average monthly temperatures. The field plots were arranged as a completely randomized block design and adopted line sowing method. The distance of row spacing and seeding depth of were 20 and 5 cm, respectively. The experiment was carried out without fertilization and irrigation and weeding manually. The whole flag leaf was collected randomly and immediately frozen in liquid nitrogen. Samples were kept at −80°C for further analyses.

The Measurements of Agronomic Traits

Plant heights were measured at heading. Each of the 5 oat plants in each plot was selected randomly to measure its height to calculate their mean. At the same stage, the ratio of stem to leaf was measured. About 0.5 kg of oat samples were weighed in each plot, and stems, leaves and inflorescences were separated into two parts. After air-dry, samples were weighed and the ratio of stem to leaf was calculated (repeated with above). The leaf sheath of the oat grasses was collected into stems and the ears into leaves. The fresh grass samples were naturally air-dried for the dry weights which were converted into hay yield per hectare. After weighed, the hay samples were crushed and sieved with 40-mesh then the crude fat (EE) was determined by the soxhlet extraction.

Measurements of the Stomatal Density and Physiological Analyses

The stomatal density of oats grown at the two altitudes were observed by SEM (Bello et al., 2017).

SOD activity was detected by recording the decrease in optical density of NBT dye by the enzyme. MAD content was measured using a modified TBA method. H₂O₂ levels, O₂^– activities, CAT, APX, GR, and Proline contents were tested using the related biological assay kits (Nanjing Jiancheng BioengineeringInstitute, China) following the manufacturer’s instructions with some modifications according to Tang et al. (2013). Approximate 0.1 g of each leaf sample was incubated on ice for 30 min and immersed in 10 ml of 95% ethanol in 25 ml brown flask at room temperature. Until the green color disappeared, the solution was set to volume, shaken, and stored in the dark. A series of 100 μl aliquots were taken every 10 min after initial immersion, and subjected to spectrophotometry (absorption measured at 647 and 664 nm) to quantify the amount of chlorophyll leached.

Before Fv/Fm measurement, whole oat plants were dark-adapted for 20 min, and Fv/Fm was measured using a steady-state gas-exchange system with an integrated fluorescence chamber head (Heinz Walz GmbH,Effeltrich,Germany). Fv/Fm values were calculated as Fv/Fm = (Fm – Fo)/Fm, NPQ = F m/F m ‘^–1.

All assays described above were repeated four times, with four biological replicates. The data, which are shown as the means ± SDs, were subjected to ANOVA to determine significant differences. The least significant differences (LSDs) of the means were determined via Duncan’s test at the level of significance defined as α = 0.05.

Library Preparation for Transcriptome Sequencing

RNA extraction was performed according to the manufacturer’s protocol of the RNA prep Pure Plant kit (KangWei Biotech Co., Ltd., Beijing, China). RNA purity and integrality were checked before the cDNA library construction. The mRNA was enriched from the total RNA through the magnetic beads with oligo (dT). Afterward, mRNA was randomly broken into fragments. Catalyzed by reverse transcriptase, double-stranded cDNA was synthesized through mRNA by combining with random primers. After ds-cDNA purified, the adapters and poly (A) tails were ligated at both ends of ds-cDNA. With the adapter sequences, ds-cDNA fragments were amplified through PCR. After ds-cDNA quantification, the cDNA libraries were ready for sequencing. using the Illumina HiSeq 2000 (single strands for sequencing).

Transcriptome Assembly and Gene Function Annotations

The clean data were obtained by removing reads containing adapters and low quality ones from raw data. Transcriptome assemblies were performed using Trinity (Grabherr et al., 2011) based on clean data. The sequences of unigenes were compared with Nr (NCBI non-redundant), Nt (Nucleotide Sequence Database), SwissProt (a manually annotated and reviewed protein sequence database), and GO (Gene Ontology), KOG (euKaryotic Orthologous Groups) databases by BLAST and the KEGG (Kyoto Encyclopedia of Genes and Genomes). GO enrichment analysis was implemented by the Goseq R package on the basis of Wallenius non-central hypergeometric distribution. KOBAS software was subsequently used to test the statistical enrichment of DEGs in the KEGG pathways (Altschul, 1997).

Differential Expression Analyses

The EdgeR program package was used to adjust read counts of each sequenced library and the RSEM package for differential expression level analyses of two samples. The p-value was adjusted using the q-value. In this study, |log2 (fold change)| > 2 were set as the thresholds for significantly differential expressions.

Validation of RNA-Seq Data by qRT-PCR

The RNA samples isolated above were used as templates and were reverse transcribed with a HiScript II Q Select RT SuperMix for qPCR (gDNA eraser) kit (Vazyme, China). The expression of the actin gene was used as an internal control. qRT-PCR was performed with ChamQ^TM Universal SYBR qPCR Master Mix(Vazyme, China) on a LightCycler 480 II58 device (Roche, Switzerland) according to the manufacturers’ protocol. Relative gene expression levels were evaluated according to the 2^–ΔΔCT method.

Conclusion

In order to understand the mechanisms of oats adapting to high altitudes, we accomplished physiological and agronomical characterizations and transcriptome sequence analysis. The physiological and agronomical results showed that oats displayed significant increases in MDA, SOD and soluble sugar contents, accompanied by decreases in chlorophyll, NPQ, stomatal density, and hay yield. The transcriptomic profiles showed apparent differences between low and high altitudes where many pathways associated with high altitude stresses were identified including photosynthetic pigment synthesis, photosynthesis, and plant hormone signaling transduction pathways. Moreover, large amount of differentially expressed TFs between low and high altitudes was identified to participate in the response to altitude stress of oats, which can serve as candidate genes for altitude stresses molecular breeding project of oat. Through molecular technology of gene cloning or gene silencing to further clarify the functions of candidate genes and identify potential target genes that increase oat yield and photosynthesis under high altitude stress in oats. Association analysis of candidate gene can find marker sites that are significantly associated with agronomic traits of oats to improve breeding efficiency. In general, this results provided the basis for understanding better the biological regulation on oats under two scenarios of altitudes and thus contributed to breeding efforts aimed at increasing yield under such situations.

Data Availability Statement

The raw sequencing data have been submitted to the NCBI Sequence Read Archive database with accession number PRJNA639320.

Author Contributions

YJ, CG, and CY conceived and designed the experiments. YJ, LB, ST, HJ, KZ, and LL performed the experiments. HW, HT, HX, and LZ analyzed and interpreted the sequence data. YJ wrote the manuscript. All authors read and approved the final manuscript.

Funding

This research was supported by National Natural Science Foundation of China (31872998). The funding body had no influence on the experimental design, data analysis and interpretation, or writing of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank the Biomarker Technologies (Beijing, China) for the RNA-seq and for the raw data analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.638683/full#supplementary-material

Supplementary Table 1 | Function annotation results of oat.

Supplementary Table 2 | Annotation of 11,639 significantly differentially expressed genes transcripts.

Supplementary Table 3 | Gene Ontology (GO) functional annotation of transcripts;

Supplementary Table 4 | KEGG pathway mapping.

Supplementary Table 5 | DEGs associated with TFs between low and high altitude.

Supplementary Table 6 | Differently expressed genes related to porphyrin and chlorophyll metabolism.

Supplementary Table 7 | Differently expressed genes related to carotenoid biosynthesis.

Supplementary Table 8 | Differently expressed genes related to photosynthesis.

Supplementary Table 9 | Differently expressed genes related to hormone signal transduction pathway.

Supplementary Table 10 | 20 differently expressed genes transcript abundance.

Supplementary Table 11 | Meteorological data of the growing season of Minhe and Huangzhong County (April–October 2015).

Abbreviations

NPQ, non-photochemical quenching coefficient; DEGs, differentially expressed genes; GO, Gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; Nr, NCBI non-redundant; Nt, Nucleotide Sequence Database; SwissProt, a manually annotated and reviewed protein sequence database; KOG, euKaryotic Orthologous Groups; TFs, transcription factors; Chl a, chlorophyll a; APX, ascorbate oxidase; CAT, catalase; GR, glutathione reductase; MDA, Malondialdehyde; SOD, superoxide dismutase; ROS, Reactive oxygen species; qRT-PCR, Quantificational real-time polymerase chain reaction; RNA-seq, RNA sequencing; SEM, scanning electron microscopy; IAA, Auxin; ABA, Abscisic acid; ET, Ethylene; SA, salicylicacid; JA, jasmonic acid; ARF, Auxin response factor; NBT, nitro-blue tetrazolium; TBA, thiobarbituric acid; FPKM, per Million mapped reads; RSEM, RNA sequencing by expectation maximization.

References

Altschul, S. F. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Anwar, A., She, M., Wang, K., Riaz, B., and Ye, X. (2018). Biological roles of Ornithine Aminotransferase (OAT) in plant stress tolerance: present progress and future perspectives. Int. J. Mol. Sci. 19:3681. doi: 10.3390/ijms19113681

PubMed Abstract | CrossRef Full Text | Google Scholar

Awika, H. O., Hays, D. B., Mullet, J. E., Rooney, W., and Weers, B. (2017). QTL mapping and loci dissection for leaf epicuticular wax load and canopy temperature depression and their association with QTL for staygreen in Sorghum bicolor under stress. Euphytica 213:207. doi: 10.1007/s10681-017-1990-5

CrossRef Full Text | Google Scholar

Bashir, H. (2011). Chloroplast and photosystems: impact of cadmium and iron deficiency. Trends Plant Sci. 53, 321–335. doi: 10.1007/s11099-015-0152-z

CrossRef Full Text | Google Scholar

Beck, J., Lohscheider, J. N., Albert, S., Andersson, U., Mendgen, K. W., Rojas-Stütz, M. C., et al. (2017). Small one-helix proteins are essential for photosynthesis in Arabidopsis. Front. Plant Sci. 8:7. doi: 10.3389/fpls.2017.00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Bello, M. A., Ruiz-León, Y., Sandoval-Sierra, J. V., Rezinciuc, S., and Diéguez-Uribeondo, J. (2017). Scanning Electron Microscopy (SEM) protocols for problematic plant, Oomycete, and fungal samples. J. Vis. Exp. 120:55031. doi: 10.3791/55031

PubMed Abstract | CrossRef Full Text | Google Scholar

Boter, M., Ruíz-Rivero, O., Abdeen, A., and Prat, S. (2004). Conserved MYC transcription factors play a key role in jasmonate signaling both in tomato and Arabidopsis. Genes Dev. 18, 1577–1591. doi: 10.1101/gad.297704

PubMed Abstract | CrossRef Full Text | Google Scholar

Ch, F. (2005). Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell 17, 1866–1875. doi: 10.1105/tpc.105.033589

PubMed Abstract | CrossRef Full Text | Google Scholar

Chi, X., Zhang, F., Gao, Q., Xing, R., and Chen, S. (2019). Genetic structure and eco-geographical differentiation of Lancea tibetica in the Qinghai-Tibetan Plateau. Genes 10:97. doi: 10.3390/genes10020097

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, E. S. (2003). Chlorophyll a fluorescence analysis in response to excitation irradiance in bean plants (Phaseolus vulgaris L. and Vigna unguiculata L. Walp) submitted to high temperature stress. Photosynthetica 41, 77–82. doi: 10.1023/A:1025860429593

CrossRef Full Text | Google Scholar

Ge, Y., Li, Y., Lv, D. K., Bai, X., Ji, W., Cai, H., et al. (2011). Alkaline-stress response in Glycine soja leaf identifies specific transcription factors and ABA-mediated signaling factors. Funct. Integr. Genom. 11, 369–379. doi: 10.1007/s10142-010-0191-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasanuzzaman, M. D., Davies, N. W., Shabala, L., Zhou, M., Brodribb, T. J., and Shabala, S. (2017). Residual transpiration as a component of salinity stress tolerance mechanism: a case study for barley. BMC Plant Biol. 17:107. doi: 10.1186/s12870-017-1054-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Heddad, M., Norén, H., Reiser, V., Dunaeva, M., Andersson, B., and Adamska, I. (2006). Differential expression and localization of early light-induced proteins in Arabidopsis. Plant Physiol. 142, 75–87. doi: 10.1104/pp.106.081489

PubMed Abstract | CrossRef Full Text | Google Scholar

Huggins, T. D., Mohammed, S., Sengodon, P., Ibrahim, A. M. H., Tilley, M., Hays, D. B., et al. (2018). Changes in leaf epicuticular wax load and its effect on leaf temperature and physiological traits in wheat cultivars (Triticum aestivumL.) exposed to high temperatures during anthesis. J. Agron. Crop Sci. 204, 49–61. doi: 10.1111/jac.12227

CrossRef Full Text | Google Scholar

Huq, E. (2002). PIF4, a phytochrome-interacting bHLH factor, functions as a negative regulator of phytochrome B signaling in Arabidopsis. EMBO J. 21, 2441–2450. doi: 10.1093/emboj/21.10.2441

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutin, C. (2003). Early light-induced proteins protect Arabidopsis from photooxidative stress. Proc. Natl. Acad. Sci. U.S.A. 100, 4921–4926. doi: 10.1073/pnas.0736939100

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, T., Okada, K., Fukazawa, J., and Takahashi, Y. (2018). DELLA-dependent and -independent gibberellin signaling. Plant Signal. Behav. 13:e1445933. doi: 10.1080/15592324.2018.1445933

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, Y., Liu, M. L., Yue, M., Zhao, Z., Zhao, G. F., Li, Z. H., et al. (2017). Comparative transcriptome analysis reveals adaptive evolution of Notopterygium incisum and Notopterygium franchetii, two high-alpine herbal species endemic to China. Molecules 22:1158. doi: 10.3390/molecules22071158

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosma, D. K., Bourdenx, B., Bernard, A., Parsons, E. P., Lü, S., Joubès, J., et al. (2009). The impact of water deficiency on leaf cuticle lipids of Arabidopsis. Plant Physiol. 151, 1918–1929. doi: 10.1104/pp.109.141911

PubMed Abstract | CrossRef Full Text | Google Scholar

Lastdrager, J., Hanson, J., and Smeekens, S. (2014). Sugar signals and the control of plant growth and development. J. Exp. Bot. 65, 799–807. doi: 10.1093/jxb/ert474

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Yang, Y., Ma, L., Sun, X., Yang, S., Kong, X., et al. (2014). Comparative proteomics analyses of Kobresia pygmaea adaptation to environment along an elevational gradient on the central Tibetan Plateau. PLoS One 9:e98410. doi: 10.1371/journal.pone.0098410

PubMed Abstract | CrossRef Full Text | Google Scholar

Lilley, J. L., Gee, C. W., Sairanen, I., Ljung, K., and Nemhauser, J. L. (2012). An endogenous carbon-sensing pathway triggers increased auxin flux and hypocotyl elongation. Plant Physiol. 160, 2261–2270. doi: 10.1104/pp.112.205575

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., and Zhang, T. (2017). Expansion and stress responses of the AP2/EREBP superfamily in cotton. BMC Genom. 18:118. doi: 10.1186/s12864-017-3517-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Tai, H., Li, S., Gao, W., Zhao, M., Xie, C., et al. (2014). bHLH122 is important for drought and osmotic stress resistance in Arabidopsis and in the repression of ABA catabolism. New Phytol. 201, 1192–1204. doi: 10.1111/nph.12607

PubMed Abstract | CrossRef Full Text | Google Scholar

Lizana, X. C., Hess, S., and Calderini, D. F. (2009). Crop phenology modifies wheat responses to increased UV-B radiation. Agric. For. Meteorol. 149, 1964–1974. doi: 10.1016/j.agrformet.2009.07.003

CrossRef Full Text | Google Scholar

Ma, L., Sun, X., Kong, X., Galvan, J. V., Li, X., Yang, S., et al. (2015a). Physiological, biochemical and proteomics analysis reveals the adaptation strategies of the alpine plant Potentilla saundersiana at altitude gradient of the Northwestern Tibetan Plateau. J. Proteom. 112, 62–83. doi: 10.1016/j.jprot.2014.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, L., Yang, L., Zhao, J., Wei, J., Kong, X., Wang, C., et al. (2015b). Comparative proteomic analysis reveals the role of hydrogen sulfide in the adaptation of the alpine plant Lamiophlomis rotata to altitude gradient in the Northern Tibetan Plateau. Planta 241, 887–906. doi: 10.1007/s00425-014-2209-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Q. J., Sun, M. H., Lu, J., Liu, Y. J., Hu, D. G., Hao, Y. J., et al. (2017). Transcription factor AREB2 is involved in soluble sugar accumulation by activating sugar transporter and amylase genes. Plant Physiol. 174, 2348–2362. doi: 10.1104/pp.17.00502

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahajan, S., and Tuteja, N. (2005). Cold, salinity and drought stresses: an overview. Arch. Biochem. Biophys. 444, 139–158. doi: 10.1016/j.abb.2005.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazor, Y., Borovikova, A., and Nelson, N. (2015). The structure of plant photosystem I super-complex at 2.8 A resolution. eLife 4:e07433. doi: 10.7554/eLife.07433

PubMed Abstract | CrossRef Full Text | Google Scholar

Noga, A., Skrzypek, E., Warchoł, M., Czyczyło-Mysza, I., Dziurka, K., Marciñska, I., et al. (2016). Conversion of oat (Avena sativa L.) haploid embryos into plants in relation to embryo developmental stage and regeneration media. Vitro Cell Dev. Biol. Plant 52, 590–597. doi: 10.1007/s11627-016-9788-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozawa, S. I., Bald, T., Onishi, T., Xue, H., Matsumura, T., Kubo, R., et al. (2018). Configuration of Ten light-harvesting chlorophyll a/b complex I subunits in Chlamydomonas reinhardtii Photosystem I. Plant Physiol. 178, 583–595. doi: 10.1104/pp.18.00749

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, Q., Wang, Q., Han, X., Guan, Y., Sun, H., Zhong, Y., et al. (2015). Transcriptome sequencing of Crucihimalaya himalaica (Brassicaceae) reveals how Arabidopsis close relative adapt to the Qinghai-Tibet Plateau. Sci. Rep. 6:21729. doi: 10.1038/srep21729

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarwat, M., Naqvi, A. R., Ahmad, P., Ashraf, M., and Akram, N. A. (2013). Phytohormones and microRNAs as sensors and regulators of leaf senescence: assigning macro roles to small molecules. Biotechnol. Adv. 31, 1153–1171. doi: 10.1016/j.biotechadv.2013.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, M., and Laxmi, A. (2015). Jasmonates: emerging players in controlling temperature stress tolerance. Front. Plant Sci. 6:1129. doi: 10.3389/fpls.2015.01129

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, P., Tiwari, S., Tripathi, N., and Mehta, A. K. (2016). Polymorphism analysis in advanced mutant population of oat (Avena sativa L.) using ISSR markers. Physiol. Mol. Biol. Plants 22, 115–120. doi: 10.1007/s12298-015-0333-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, Y., Li, W., Zhao, J., Liu, Y., and Guo, C. (2018). Transcriptome sequencing and expression profiling of genes involved in the response to abiotic stress in Medicago ruthenica. Genet. Mol. Biol. 41, 638–648. doi: 10.1590/1678-4685-gmb-2017-0284

PubMed Abstract | CrossRef Full Text | Google Scholar

Sicora, C. (2006). UV-B and UV-A radiation effects on photosynthesis at the molecular level. Earth Environ. Sci. 57, 121–135. doi: 10.1007/1-4020-3697-3_10

CrossRef Full Text | Google Scholar

Singh, D., and Laxmi, A. (2015). Transcriptional regulation of drought response: a tortuous network of transcriptional factors. Front. Plant Sci. 6:895. doi: 10.3389/fpls.2015.00895

PubMed Abstract | CrossRef Full Text | Google Scholar

Spartz, A. K., Ren, H., Park, M. Y., Grandt, K. N., Lee, S. H., Murphy, A. S., et al. (2014). SAUR inhibition of PP2C-D phosphatases activates plasma membrane H+-ATPases to promote cell expansion in Arabidopsis. Plant Cell 26, 2129–2142. doi: 10.1105/tpc.114.126037

PubMed Abstract | CrossRef Full Text | Google Scholar

Storz, J. F., Scott, G. R., and Cheviron, Z. A. (2010). Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates. J. Exp. Biol. 213(Pt 24), 4125–4136. doi: 10.1242/jeb.048181

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, A., and Tanaka, R. (2006). Chlorophyll metabolism. Curr. Opin. Plant Biol. 9, 248–255. doi: 10.1016/j.pbi.2006.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, L., Cai, H., Ji, W., Luo, X., Wang, Z., Wu, J., et al. (2013). Overexpression of GsZFP1 enhances salt and drought tolerance in transgenic alfalfa (Medicago sativa L.). Plant Physiol. Biochem. 71, 22–30. doi: 10.1016/j.plaphy.2013.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiebaut, F., Hemerly, A. S., and Ferreira, P. C. G. (2019). A role for epigenetic regulation in the adaptation and stress responses of non-model plants. Front. Plant Sci. 10:246. doi: 10.3389/fpls.2019.00246

PubMed Abstract | CrossRef Full Text | Google Scholar

Tuan, P. A., Yamasaki, Y., Kanno, Y., Seo, M., and Ayele, B. T. (2019). Transcriptomics of cytokinin and auxin metabolism and signaling genes during seed maturation in dormant and non-dormant wheat genotypes. Sci. Rep. 9:3983. doi: 10.1038/s41598-019-40657-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Vinocur, B., and Altman, A. (2005). Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr. Opin. Biotechnol. 16, 123–132. doi: 10.1016/j.copbio.2005.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Bai, Z. Y., Liang, Q. Y., Liu, Q. L., Zhang, L., Pan, Y. Z., et al. (2018). Transcriptome analysis of chrysanthemum (Dendranthema grandiflorum) in response to low temperature stress. BMC Genom. 19:319. doi: 10.1186/s12864-018-4706-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B., Huo, P., Zhang, Q., Chen, X., and Zhan, Z. (2017). Transcriptome analysis of hexaploid hulless oat in response to salinity stress. PLoS One 12:e0171451. doi: 10.1371/journal.pone.0171451

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, G. S., Long, R. J., Guo, X. S., Irvine, J., Ding, L. M., Ding, L. L., et al. (2011). Blood mineral status of grazing Tibetan sheep in the Northeast of the Qinghai–Tibetan Plateau. Livestock Sci. 136, 102–107. doi: 10.1016/j.livsci.2010.08.007

CrossRef Full Text | Google Scholar

Yang, H. L., Korivi, M., Lin, M. K., Chang, H. C., Wu, C. R., Lee, M. S., et al. (2017). Antihemolytic and antioxidant properties of pearl powder against 2,2’-azobis(2-amidinopropane) dihydrochloride-induced hemolysis and oxidative damage to erythrocyte membrane lipids and proteins. J. Food Drug Anal. 25, 898–907. doi: 10.1016/j.jfda.2016.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Wang, L., Han, J., Tang, X., Ma, M., Wang, K., et al. (2015). Comparative transcriptomic analysis revealed adaptation mechanism of Phrynocephalus erythrurus, the highest altitude Lizard living in the Qinghai-Tibet Plateau. BMC Evol. Biol. 15:101. doi: 10.1186/s12862-015-0371-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zeng, Z., Chen, C., Li, C., Xia, R., and Li, J. (2019). Genome-wide characterization of the auxin response factor (ARF) gene family of litchi (Litchi chinensis Sonn.): phylogenetic analysis, miRNA regulation and expression changes during fruit abscission. PeerJ 7:e6677. doi: 10.7717/peerj.6677

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Q., Xiang, X., Liu, D., Yang, A., and Wang, Y. (2018). Tobacco transcription factor NtbHLH123 confers tolerance to cold stress by regulating the NtCBF pathway and reactive oxygen species homeostasis. Front. Plant Sci. 9:381. doi: 10.3389/fpls.2018.00381

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, M., Qilong, Li, Chen, Z., Lv, Q., Bao, F., Wang, X., et al. (2018). Regulatory Mechanism of ABA and ABI3 on vegetative development in the moss Physcomitrella patens. Int. J. Mol. Sci. 19:2728. doi: 10.3390/ijms19092728

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, S. Y., Chen, L.-Y., Muchuku, J. K., Hu, G.-W., and Wang, Q.-F. (2016). Genetic adaptation of giant lobelias (Lobelia aberdarica and Lobelia telekii) to different altitudes in East African Mountains. Front. Plant Sci. 7:488. doi: 10.3389/fpls.2016.00488

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, P.-J., and Yang, L. (2015). Ambient UV-B radiation inhibits the growth and physiology of Brassica napus L. on the Qinghai-Tibetan plateau. Field Crops Res. 171, 79–85. doi: 10.1016/j.fcr.2014.11.006

CrossRef Full Text | Google Scholar