Next Article in Journal
Effects of Invasive Alien Plant Species on Native Plant Diversity and Crop Yield
Previous Article in Journal
Species Identification and Genetic Diversity Analysis of Medicinal Plants Aconitum pendulum Busch and Aconitum flavum Hand.-Mazz.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 6
10.3390/plants13060887
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum)

by
Zeyu Li
1,†,
Jiyao Liu
1,†,
Yan Chen
1,
Aihua Liang
2,3,
Wei He
1,
Xiaoya Qin
4,
Ken Qin
4,* and
Zixin Mu
2,3,*
1
College of Life Sciences, Northwest A&F University, Yangling 712100, China
2
College of Life Sciences & Technology, Tarim University, Alaer 843300, China
3
State Key Laboratory Breeding Base for the Protection and Utilization of Biological Resources in Tarim Basin Co–Funded by Xinjiang Corps and the Ministry of Science and Technology, Alaer 843300, China
4
National Wolfberry Engineering Research Center, Ningxia Academy of Agriculture and Forestry Sciences, Yinchuan 750002, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2024, 13(6), 887; https://doi.org/10.3390/plants13060887
Submission received: 12 February 2024 / Revised: 13 March 2024 / Accepted: 15 March 2024 / Published: 20 March 2024
(This article belongs to the Section Plant Response to Abiotic Stress and Climate Change)
Browse Figures
Review Reports Versions Notes

Abstract

:
The characterization of the PYL/RCAR ABA receptors in a great deal of plant species has dramatically advanced the study of ABA functions involved in key physiological processes. However, the genes in this family are still unclear in Lycium (Goji) plants, one of the well–known economically, medicinally, and ecologically valuable fruit crops. In the present work, 12 homologs of Arabidopsis PYL/RCAR ABA receptors were first identified and characterized from Lycium (L.) barbarum (LbPYLs). The quantitative real–time PCR (qRT–PCR) analysis showed that these genes had clear tissue–specific expression patterns, and most of them were transcribed in the root with the largest amount. Among the three subfamilies, while the Group I and Group III members were down–regulated by extraneous ABA, the Group II members were up–regulated. At 42 °C, most transcripts showed a rapid and violent up–regulation response to higher temperature, especially members of Group II. One of the genes in the Group II members, LbPYL10, was further functionally validated by virus–induced gene silencing (VIGS) technology. LbPYL10 positively regulates heat stress tolerance in L. barbarum by alleviating chlorophyll degradation, thus maintaining chlorophyll stability. Integrating the endogenous ABA level increase following heat stress, it may be concluded that LbPYL–mediated ABA signaling plays a vital role in the thermotolerance of L. barbarum plants. Our results highlight the strong potential of LbPYL genes in breeding genetically modified L. barbarum crops that acclimate to climate change.

1. Introduction

The phytohormone abscisic acid (ABA) plays a vital role in nearly all plant physiological processes under both normal conditions and stress conditions [1]. In comparison to its biosynthesis and catabolism, ABA signaling sense and signaling transduction may be more fundamental because there are more complicated regulation mechanisms that have evolved for the latter [2]. The most characteristic ABA receptor to date is a family of soluble proteins called PYR (pyrabactin resistant), PYL (PYR–like), or RCAR (regulatory components of ABA receptors), which form the beginning of a “core ABA signaling pathway” [3,4]. The pathway includes a PYL/RCAR ABA receptor, a negative regulator of clade A protein phosphatases of type 2C (PP2C), and a positive regulator of sucrose nonfermenting–1 (SNF1)–related protein kinase 2s (SnRK2s), which have been identified and characterized in a variety of plant species [5,6,7]. In this model, ABA binding stabilizes the receptor complex; as a result, the active site of the phosphatase is blocked and the associated protein kinase is released from PP2C inhibition. The activated protein kinase further activates downstream transcription factors and plasma membrane proteins to regulate gene expression and stomatal closure, respectively [8,9]. Genetic modification the coding genes of PYL/RCAR ABA receptors can improve plant water use efficiency in Arabidopsis, wheat, and maize, presenting a prospective target for engineering drought–tolerant crops under changing environmental conditions [10,11,12]. Moreover, chemical manipulation of ABA signaling by ABA agonists [13,14] and direct modification of PYL/RCAR ABA receptors by protein engineering [15] all have significant agricultural and ecological applications.
The ambient air temperature caused by global warming continues to rise, which makes heat stress (HS) even more serious [16,17,18]. When the temperature is higher than 38 °C, it will inhibit the growth of many crops, including rice, corn, wheat, etc., especially the pollen fertility of many crops will be reduced, and then the seed setting rate will decline, resulting in a significant reduction in crop yield. It will also cause grain filling not to be full, so that the grain quality is significantly reduced. HS severely inhibits photosynthesis and reduces carbohydrate synthesis by reducing the enzyme activity of the chloroplast protein complex. At the same time, the respiration of the plant will become strong, decompose, and consume a lot of organic nutrients, resulting in the “chlorina” and premature aging of the leaves [16,17,18]. The heat stress response (HSR) of plants is initiated by increased plasma membrane mobility and the activation of phospholipase (PLC), which hydrolyzes phosphatidylinositol 4, 5–diphosphate (PIP2) to inositol triphosphate (IP3) and diacylglycerol (DAG). The controlled entry of Ca2+ from the endoplasmic reticulum (ER) activates calmodulin (CaM) and heat shock signaling pathways [18]. The downstream signaling components, including protein kinases, e.g., calcium–dependent protein kinases (CDPKs) and mitogen activated protein kinases (MAPKs), activate heat shock transcription factors (HSFs) by protein phosphorylation. HSFs then activate heat shock proteins (HSPs) and other target proteins that are related to heat tolerance. HSF/HSP molecular modules are at the core of plant transcriptional regulatory networks upon HS.
It is well known that in the plant hormone–mediated regulation of heat tolerance in response to global climate change [19], the stress phytohormone ABA plays a key role [20]. Up–regulation of HSFA2c and HSPs by ABA contributes to the improved heat tolerance of various plant species [21]. The ABA–HSF/HSP signaling pathway was acted on in grafted cucumber plants [22], tall fescue and Arabidopsis [23], and developing grains in wheat [24]. ABA also exerted functions by reprogramming metabolism upon high–temperature stress, e.g., by mediating sugar metabolism in rice spikelets to prevent pollen abortion [25], through regulation of osmolytes and antioxidants to protect photosynthesis and growth in wheat plants [26], as well as to improve grain filling capacity in rice [27]. However, whether PYL/RCAR–mediated ABA signaling was involved in the heat stress response remains unclear so far. The functional characterization of the orthologous genes of PYL/RCAR ABA receptors in more crop plants, as well as exploration of their transcriptional regulation network, holds great promise for both sustainable agriculture and ecological systems, especially under global climate change [13].
Most Lycium (Goji) plants not only have medicinal and edible value but also have good resistance to adverse environments [28,29,30]. Lycium plays an important role in the utilization of marginal land due to its characteristics of drought, salt, and alkali resistance and barren resistance. However, global climate change–derived heat stress has greatly impeded the yield and quality of L. barbarum, one of the main popular cultivars around the world, in recent years [31]. Continuous high–temperature weather resulted in fruit falling, accelerated ripening, shortened picking intervals, advanced aestivation (a physiological phenomenon under high–temperature conditions in summer, that is, growth stops and leaves fall off and enter a dormancy state), and a delayed growth period of autumn twigs. At the same time, it also causes the occurrence of wolfberry diseases and pests [31]. In this case, exploring the molecular mechanism underlying the heat stress response becomes a vital issue for Lycium biology. The role of PYL/RCAR–mediated ABA signaling in the heat stress response of Lycium plants is still unclear, and 12 homologs of Arabidopsis PYL/RCAR ABA receptors were first identified from L. barbarum at a genome–wide scale in the present work. After a systematic bioinformatics analysis, a transcript profiling assay was conducted in response to ABA or heat stress. It is found that all of the genes were up–regulated by heat stress in a time–dependent manner. One of the genes, LbPYL10, which was transcriptionally activated both by ABA and heat stress, has been shown positively mediate the heat stress tolerance of L. barbarum seedlings by virus–induced gene silencing (VIGS) analysis. All these findings imply that the PYL/RCAR–mediated ABA signaling could play important roles in the heat stress response in L. barbarum, which provides important gene resources for the genetic improvement of the heat resistance of L. barbarum in a changing environment.

2. Results

2.1. Heat Stress–Induced ABA Accumulation in L. barbarum

In view of the role of ABA exerted in the heat tolerance of L. barbarum, we conducted a time–course ABA assay under 42 °C (Figure 1). It is shown that heat stress significantly enhanced ABA accumulation, though the response time was relatively slow. For the time course, there was not a significant change in the one–hour treatment, while the level was greatly raised after the 3 h treatment and then remain unchanged within the 24 h treatment time. This result highlights the potential role of ABA function in the heat stress tolerance of L. barbarum and promotes that we further explore its core signaling components.

2.2. Characterization of the LbPYL Gene Family

Twelve different genes encoding complete LbPYLs were identified both by BLAST searches in the L. barbarum reference genome using 12 Arabidopsis PYL amino acid sequences as queries (Table 1) and by conserved domain MIP PF00230 validation. The CDS, TMDs, pI, and MW of 12 LbPYLs were analyzed. As can be seen from the Table, most LbPYLs contain CDS between 534 and 747 base pairs in length. The predicted proteins range in length from 177 to 248 amino acids, with MWs of 20.229 to 27.636 kDa. The pI values vary between 5.34 and 8.65. These genes were unevenly distributed in all 12 chromosomes (Chr), except Chr02, Chr10, and Chr11, and detailed information is presented in Figure S1.

2.3. Phylogeny Analysis

To systematically classify L. barbarum PYL genes and uncover the evolutionary relationship with PYL genes from other plants, a phylogenetic tree was constructed with MEGA11 (11.0.13) using neighbor–joining (NJ) analysis (Figure 2). To study the evolutionary relationships between LbPYLs and PYLs from Arabidopsis thaliana, Solanum lycopersicum, Populus trichocarpa., and Zea mays L., the 12 LbPYLs with 14 AtPYLs, 14 SlPYLs, 14 PtPYLs, and 13 ZmPYLs were clustered into three groups, designated Group I, Group II, and Group III, on an unrooted phylogenetic tree. In this tree, with the exception of LbPYL12, LbPYL1 and LbPYL2 belonged to Group I, LbPYL3, LbPYL4, LbPYL5, LbPYL10, and LbPYL11 belonged to Group II, and LbPYL6, LbPYL7, LbPYL8, and LbPYL9 belonged to Group III (Figure 2). LbPYL12 has high homology with SlPYL13, and previous reports have also indicated that SLPYL13 (Sl02g076770) is independent of other subfamilies [32]. Therefore, we believe that LbPYL12 is also independent of other subfamilies.

2.4. Gene Structure, Conserved Motifs, and Domain Analysis

In order to fully understand the protein structure of LbPYLs, MEME (5.5.4) software was used to analyze the motif distribution of all 12 LbPYLs, and 15 conserved motifs were found (Figure 3B). Motifs 1, 2, and 3 were highly conserved across all LbPYLs. Whereas motifs 6 and 15 were only present in the Group I subfamily, motifs 4, 8, 9, and 10 were only present in the Group II subfamily, and motifs 5, 7, and 12 were only present in the Group III subfamily, except LbPYL9. Motif 11 was present in both LbPYL2 and LbPYL7, motif 13 was present in both LbPYL11 and LbPYL12, and motif 14 was present in both LbPYL9 and LbPYL5 (Figure 3A,B).
The intron–exon structure of LbPYLs was analyzed by TBtools (Figure 3D). The LbPYLs can be distinctly grouped into a clade with two introns and two other clades. Groups I and II have one intron or no intron, and all members of Group III have two introns. This result indicates that members of the subfamily have a highly similar genetic structure, which is consistent with their phylogenetic relationships (Figure 3C,D).
Similar to the PYR/PYL/RCAR proteins in Arabidopsis, all LbPYL proteins are composed of a highly similar helix–grip structure, characterized by two α–helixes (red) on either side of a seven–strand β–sheet (green) (Figure 4).

2.5. Analysis of Cis–Elements in the LbPYL Promoters

To better understand the transcriptional regulation and potential function of the LbPYL genes, we isolated sequences within 2000 bp upstream of the LbPYLs’ start codon and identified the cis–elements in these promoter sequences using the PlantCARE database. The major cis–elements included stress–related elements, such as anaerobic induction regulatory elements (AREs), dehydration stress elements (DRE cores), light responsive elements (G–boxes), MYBs, MYCs, and AAGAA motifs, as well as hormone–related elements, such as abscisic acid–responsive elements (ABREs) and ethylene–responsive elements (EREs). The promoter region of 12 LbPYL genes contains at least six cis–elements, and all LbPYL genes contain MYB and MYC elements. As can be seen from the statistical table, LbPYL10 has the largest number of cis–acting elements, followed by LbPYL11 (Figure 5).

2.6. Tissue–Specific Expression

To characterize the expression of the LbPYL gene family, we analyzed five different tissues of L. barbarum based on qRT–PCR (Figure 6). The different expression levels of LbPYL members in different tissues indicate that different tissues require different functions. Notably, most LbPYL genes are expressed at higher levels in the roots than in other organs. The genes LbPYL3, LbPYL4, LbPYL5, LbPYL6, LbPYL7, LbPYL8, LbPYL9, LbPYL10, and LbPYL11 had maximal transcripts in roots rather than in other tissues. The LbPYL1 and LbPYL2 genes had the highest transcript levels in leaves, and LbPYL12 had the highest transcript levels in fruits. No genes were preferentially expressed in flowers. Among the three subfamilies, subfamily II and subfamily III were preferentially expressed in roots, and subfamily I was expressed not only in roots but also in leaves (Figure 6). Taken together, these overlapping and preferential expression patterns of LbPYLs may allow L. barbarum to perform different ABA responses in specific tissues.

2.7. Transcriptional Response to Exogenous ABA

Under 10 μM ABA treatment, the transcript content decreased with the extension of the treatment time in 7 of the 12 members. Compared with the control, the expression levels of LbPYL2, LbPYL4, LbPYL5, and LbPYL10 were up–regulated for a short time and then decreased (LbPYL5 and LbPYL10) or remained unchanged (LbPYL2 and LbPYL4) during the three–hour treatment (Figure 7). The expression of the LbPYL3 gene continued to increase within three hours. It is clear that the up–regulated genes mainly belonged to the Group II subfamily, while those genes belonging to the Group I and Group III subfamilies were consistently down–regulated by exogenous ABA (Figure 7).

2.8. Transcriptional Response to Heat Stress

Within 24 h, the expression of the LbPYL gene first increased and then decreased in response to heat stress (Figure 8). The time dynamics of the responses are different for different subfamilies. While one–half of them (LbPYL1, LbPYL2, LbPYL11, LbPYL4, and LbPYL9) were up–regulated by 2~5–fold, the remaining genes (LbPYL5, LbPYL3, LbPYL10, LbPYL7, LbPYL8, LbPYL6, and LbPYL12) were up–regulated over 5–fold. In the latter, two genes, LbPYL10 and LbPYL12, were up–regulated over 30–fold (variations in scale require special attention). In the time course, three genes (LbPYL5, LbPYL3 and LbPYL10) obtained maximal expression at 1 h, while three genes (LbPYL1, LbPYL2 and LbPYL11) obtained expression peaks at 3 h. The remaining genes (LbPYL4, LbPYL7, LbPYL8, LbPYL6, LbPYL9, and LbPYL12) had relatively slower response speeds (Figure 8). As you can see, members of subfamily II respond the fastest, followed by subfamily I, and subfamily III responds the slowest. These findings suggest that LbPYL members might play diverse roles in sensing the ABA signal for L. barbarum to adapt to heat stress.

2.9. Functional Validation of LbPYL10 in Heat Stress Tolerance

Considering that there is a significant transient up–regulation pattern for the LbPYL10 gene in comparison to its homologues in L. barbarum under heat stress, we further validated its function in heat tolerance by VIGS technology (Figure 9). The 329 bp fragment of the LbPYL10 gene was cloned into the TRV RNA2 vector (Figure 9A and Figure S2). After detecting LbPYL10 expression in TRV–infected plants by RT–qPCR technology, the significantly LbPYL10–silenced lines were used for subsequent experiments (Figure 9B). Compared with the control, LbPYL10 silencing significantly stimulated leaf disc bleaching and senescence following two hours of 42 °C treatment (Figure 9C).
In order to explore the mechanism underlying LbPYL10 silencing–mediated leaf disc bleaching, the chlorophyll degradation–related genes were identified from L. barbarum at a genome–scale (Figure S3), and their expression pattens were analyzed in LbPYL10–silenced plants upon heat stress. It is shown that LbPYL10 silencing significantly up–regulated pheophorbide a oxygenase (LbPAO) and stay–green proteins like LbSGRL expression (Figure 9D,E). LbPAO and LbSGRL were homologs of Arabidopsis NtPAO (EU294211.1) and AtSGRL (AT1G44000), respectively (Figures S4 and S5). Our results may highlight that LbPYL10 actively mediates the heat stress tolerance of L. barbarum by inhibiting chlorophyll degradation, thus delaying leaf senescence.

3. Discussion

The development of PYLs/RCARs, as an initial factor inducing ABA signaling, is considered to be a key event in the evolution of aquatic plants into terrestrial plants. No members of the PYL/RCAR family are present in the aquatic algae Chlamydomonas, while the number of PYL/RCAR family members increases from mosses and ferns to angiosperms. By far the largest numbers have been reported in Musa acuminata and Brassica rapa (24), and the lower plant Marchantia polymorpha had the least (four) [2]. A total of 14 and 23 members of PYL/RCAR ABA receptors have been identified in two Solanaceae families, tomato and tobacco, respectively [6,32]. In the present work, we found that 12 members of this family of proteins exist in one of the most well–known Lycium plants, L. barbarum, highlighting the high conservation of this family of genes. It is speculated that heat stress may indirectly regulate LbPYL gene expression by activating the upstream transcriptional factors, e.g., MYB, DRE, etc., due to the fact that no potential heat shock elements (HSEs) were identified in these genes. It should be noted that L. barbarum is only one of the eight species, and four varieties of Lycium plants exist in China [28,29,30]. Along with the pangenome published for more Lycium plants in the future, the number of PYL/RCAR ABA receptors may also vary within this plant species.
ABA takes part in the heat stress response and has also been reported in some plant species [19,20,21,23,33,34,35], in which the exogenous application of this hormone prior or parallel to heat stress renders plants more thermotolerant, although the signaling mechanism is not well documented. In this study, we found that most of the LbPYL transcripts showed rapid and sharp up–regulation under a high–temperature stress of 42 °C. Integrating the time course and amount of the transcriptional response, it is speculated that the genes LbPYL3, LbPYL5, and LbPYL10, which belong to subfamily II, may play a major role in the heat acclimation of L. barbarum seedlings. It is speculated that this subfamily’s members conservatively mediate plants’ acclimation to global climate changes, as the overexpression of its orthologs TaPYL4 and PePYL4 enhanced draught tolerance in wheat and populus, respectively [11,36]. Considering that the level of endogenous ABA increased following heat stress, and LbPYL10 silencing significantly stimulated leaf disc bleaching and senescence, it may be concluded that LbPYL–mediated ABA signaling plays a vital role in the heat stress tolerance of L. barbarum plants.
As is widely acknowledged, photosynthesis is significantly impacted by aging or abiotic stress, and plants enhance their survival rate through nutrient recycling from damaged organelles [37,38]. It has been reported that high temperature can disrupt the stability of the light–harvesting chlorophyll–protein complex and result in the deactivation of photosynthetic enzymes. Additionally, under high–temperature stress, there was a significant increase in the expression of genes associated with chlorophyll degradation (PPH, ACD2, NYE, and HCAR) in Arabidopsis thaliana [39,40]. The crucial role of AtPAO in chlorophyll degradation during leaf aging has previously been identified in Arabidopsis thaliana; it is believed that Pheide a degradation ultimately leads to leaf bleaching [41]. Furthermore, studies have demonstrated that Arabidopsis plants overexpressing AtSGRL exhibit earlier leaf bleaching under abiotic stress conditions while leaves of atsgrl–1 mutants remain green for an extended period [42]. In our current study on L. barbarum, we examined the relative expression levels of LbPAO and LbSGRL. We observed higher expression levels in LbPYL10–silenced plants under high–temperature stress. Silencing LbPYL10 strengthened the chlorophyll degradation pathway and accelerated leaf bleaching after an equivalent duration of high–temperature stress. This promoted the process of leaf senescence. In conclusion, LbPYL10 positively regulates heat stress tolerance in L. barbarum by suppressing the chlorophyll degradation pathway and maintaining stability within the chlorophyll–protein complex.
In the future, the further characterization of the downstream signaling molecules of LbPP2Cs and LbSnRK2s, followed by exploration of the LbPYL–LbPP2C–LbSnRK2 combinations, will greatly advance our understanding of LbPYL–mediated ABA signaling in Lycium plants, as complicated combinations have been shown in Arabidopsis, maize, and tomato [9,43,44,45,46]. Moreover, establishing an efficient CRISPR/Cas9 gene editing system for L. barbarum will greatly accelerate the genetic improvement of L. barbarum for climate adaptation, such as LbPYL10. Finally, it should be noted that the reproductive growth stage is usually the most sensitive stage to heat stress; therefore, exploring LbPYL–mediated ABA signaling in this stage would provide more valuable genetic information.

4. Materials and Methods

4.1. Plant Materials

The different tissues, including young roots, stems, leaves, and flowers, as well as the ripening fruits of Lycium (L.) barbarum, were collected from three of the 5–year–old trees at the Wolfberry (Lycium) Germplasm Repository of Ningxia, Academy of Agriculture and Forestry. Our field studies were conducted in accordance with local legislation and appropriate permissions.

4.2. Heat Stress Assay

Two–month–old uniform clonal seedlings of L. barbarum, which were grown in pots in a greenhouse, were transported to a climate chamber. After acclimation to the artificial environment (25 °C, 16 h day/20 °C, 8 h night) for one week, the seedlings were divided into two parts. Whereas one part was still left in the same climate chamber as the control, the other part was transported to a 42 °C growth chamber for heat stress. Other environmental factors, such as light intensity (150 μmol∙m−2∙s−1) and humidity (65%), remained consistent between the two climate chambers. The seedlings were subjected to heat stress for 0, 1, 3, 6, 12, and 24 h. The leaves were sampled at each time point and immediately frozen by liquid nitrogen for RNA extraction. Three biological replicates were set for each treatment.
The function of LbPYL10 in heat stress tolerance was validated by leaf disc assay, as described by Lee [47] with some modifications. Discs (5 mm in diameter) were punched from the second to fifth fully expanded leaves of 25–day–old soil–grown plants and then floated on 2 mL of 10 mM MES–KOH buffer, pH 6.8, on 12– or 24–well microplates. Plates were incubated at 42 °C for 2 h. After heat stress treatment, leaf discs were returned to 25 °C under 16 h of light/8 h of dark and photographed 5 d later.

4.3. ABA Treatment

One–month–old uniform seed seedlings of L. barbarum, which were cultured in 1/2 Hoagland nutrient solution in a climate chamber, were subjected to 10 μM ABA for 0 (Mock), 0.5, 1, and 3 h, respectively. The climate chamber conditions were as follows: light intensity (150 μmol·m−2·s−1), humidity (65%), 25 °C, 16 h day/20 °C, and 8 h night. The ABA storage solution (dissolved in methanol) was diluted in 1/2 Hoagland nutrient solution. The leaves were sampled at each time point and immediately frozen by liquid nitrogen for RNA extraction. Three biological replicates were set for each treatment.

4.4. ABA Content Determination

The ABA content was determined as described by Guo [48]. Fifty milligrams (fresh weight) of the sample powder was extracted with 1 mL methanol/water/formic acid (15:4:1, v/v/v). The combined extracts were evaporated to dryness under a nitrogen gas stream, followed by being resuspended in 100 μL 80% methanol (v/v), as well as filtered through a 0.22 μm filter for further LC–MS analysis. The sample extracts were analyzed using a UPLC–ESI–MS/MS system (UPLC ExionLC™ AD https://sciex.com.cn/, accessed on 27 August 2023; MS Applied Biosystems 6500 Triple Quadrupole, https://sciex.com.cn/, accessed on 27 August 2023).

4.5. Virus–Induced LbPYL10 Gene Silencing

The CDS fragment of the LbPYL10 gene carrying the adaptor sequence was obtained by a pair of amplification primers (Forward: 5′–AAGAGCTGCCACATCATCGG; Reverse: 5′–GCGATCTGGGAAAGGGATTG). The pTRV2 vector was digested with EcoR I and BamH I to not only linearize the vector but also expose the adaptor sequence. The PCR product was connected to the sticky end of the vector using a DNA OK Clon DNA Ligation Kit (Accurate Biotechnology, Changsha, China) and then transformed into E. coli to screen for pTRV2–derivative colonies with the silenced target gene (pTRV2–LbPYL10). After pTRV1 and pTRV2 or pTRV2–LbPYL10 were rapidly frozen in liquid nitrogen, they were transfected with the Agrobacterium tumefaciens strain GV3101 by heat shock.
The infection method was performed as previously described [49]. Agrobacterium strains of pTRV1, pTRV2, and pTRV2–LbPYL10 were collected and resuspended in infiltration buffer (10 mM MgCl2, 10 mM MES, and 200 μM acetosyringone) at an optical density of 0.3 at 600 nm, and then suspensions were maintained at room temperature for at least 3 h. Finally, pTRV2 or pTRV2–LbPYL10 was mixed at a 1:1 ratio with pTRV1 before injection. One–month–old seedlings were prepared and injected using the pressure permeation method. After 2 min of injection, the surface was cleaned with ddH2O. After 2–3 weeks, the silencing efficiency in pTRV2–LbPYL10 plants was analyzed using qRT–PCR.

4.6. Identification and Chromosomal Location of LbPYL Genes

To identify the putative PYL proteins in L. barbarum (LbPYLs), the 14 AtPYR/PYLs in the A. thaliana genome were downloaded from the TAIR database (https://www.arabidopsis.org/, accessed on 5 February 2023) and were then used as queries to BLAST search the L. barbarum genome (https://www.ncbi.nlm.nih.gov/genome/81199?genome_assembly_id=1656998, accessed on 23 December 2023) with an E–value of e−10. Then, preliminary amino acid sequences that may function as LbPYLs were obtained according to the homology of AtPYR/PYLs. On this basis, conservative PYL domains obtained from the PFAM database (http://pfam-legacy.xfam.org/, accessed on 1 May 2023) were blast searched against these candidate sequences. Finally, amino acid sequences without conserved PYL domains and redundant sequences were manually removed.
The values of the molecular weights (MWs) and isoelectric points (pIs) of the LbPYL proteins were analyzed with ProtParam (http://web.expasy.org/protparam/, accessed on 20 January 2023).
The chromosomal locations of LbPYL genes were displayed with TBtools (v2.057) software [50] based on the genome annotation files downloaded from the Genome Database for L. barbarum, in which the chromosome numbers and positions of each sequence in the genome are indicated.

4.7. Phylogenetic Tree Analysis of LbPYLs in L. barbarum

Sequences of PYLs from A. thaliana, S. lycopersicum, Populus L., and Zea mays L. were acquired from Ensemble Plants (http://plants.ensemble.org/index.html, accessed on 5 February 2023). The maximum likelihood tree (ML tree) was constructed by the neighbor–joining (NJ) method with MEGA11 (11.0.13) software derived from a ClustalW2 alignment of related amino acid sequences (bootstrap replicates = 1000). The ML tree was beautified by the ChiPlot website (https://www.chiplot.online/#, accessed on 7 February 2023) [51].

4.8. Analysis of Gene Exon–Intron Structures and Protein Conserved Motifs

TBtools (v2.057) software was used to draw the structure map. The conserved motif analysis was generated by the classic mode of MEME (Multiple Em for Motif Elicitation, https://meme-suite.org/meme/, accessed on 2 May 2023), where the number of motifs was set to 15 and other settings were consistent with the default parameters.

4.9. Analysis of Cis–Elements in the LbPYL Promoters

The promoter regions 2 kb upstream of the corresponding gene were analyzed, and the cis–elements were predicted by the PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 8 March 2023). Subsequently, according to the results of the PlantCARE calculation, we recorded the numbers of each cis–element in these sequences and summarized this information in a figure for subsequent analysis.

4.10. Quantitative Real–Time PCR (qRT–PCR) Assay

Total RNA was isolated from different tissues using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), and its quality was examined by NanoDrop (Thermo Scientific, Waltham, MA, USA). First–strand cDNAs were synthesized from 1 μg of RNA using the FastKing RT Kit (TIANGEN, Beijing, China). qRT–PCR was conducted using 0.1 μg of cDNA on a CFX96 Touch (BIO–RAD) with SuperReal PreMix Color (SYBR Green) and with an initial denaturing step at 95 °C for 15 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s, and 72 °C for 32 s. The PCR primers used for qRT–PCR are listed in Table S1.

4.11. Statistical Analysis

Statistical analyses were performed using SPSS version 19.0. Parameter differences were determined using one–way ANOVA (Duncan’s test) with appropriate post hoc analysis. In some cases, significant differences in mean values were determined by Student’s t–test and different letters indicate significant differences at p < 0.05. The column and line figures were drawn by Origin 2023b.

5. Conclusions

This study presented a comprehensive genome–wide analysis of the PYR1/PYL/RCAR ABA receptor gene family in L. barbarum. The LbPYL gene expression profiling analyzed by qRT–PCR revealed high expression in roots versus other tissues and differential sensitivity to exogenous ABA or heat stress. It is speculated that the genes LbPYL3, LbPYL5, and LbPYL10, which belong to subfamily II, may play a major role in the heat stress acclimation of L. barbarum seedlings. LbPYL10, which was instantaneously up–regulated by heat stress or ABA treatment, positively regulated heat stress tolerance by alleviating chlorophyll degradation, thus maintaining chlorophyll stability. Taken together, our results highlight a vital role of LbPYL ABA receptor–mediated ABA signaling in L. barbarum’s response to heat stress.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants13060887/s1, Figure S1: Chromosomal distribution of LbPYL genes; Figure S2: Analysis of silent regions in LbPYL10; Figure S3: Phylogenetic tree of genes related to chlorophyll degradation; Figure S4: Conserved domain; Figure S5: Conserved motifs; Table S1: The primer sequences in qRT–PCR; Table S2: PYL genes in five species examined; Table S3: The sequences of motifs and cis–acting elements; Table S4: The sequences of cis–acting elements.

Author Contributions

Conceptualization, Z.M.; methodology, Z.L.; validation, Z.L. and Y.C.; formal analysis, J.L. and W.H.; investigation, Z.L. and J.L.; resources, K.Q.; writing—original draft preparation, A.L. and Z.M.; writing—review and editing, Z.M. and Z.L.; project administration, X.Q. and K.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Research and Development Project foundation of Ningxia Province of China (No. 2020BFH03005) and the Foreign Science and Technology Cooperation Project of Ningxia Academy of Agriculture and Forestry Sciences (No. DW–X–2020009).

Data Availability Statement

Data are contained within the article.

Acknowledgments

We thank Yule Liu at Tsinghua University for generously providing the pTRV1 and pTRV2 vectors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yoshida, T.; Christmann, A.; Yamaguchi-Shinozaki, K.; Grill, E.; Fernie, A.R. Revisiting the basal role of ABA-Roles outside of stress. Trends Plant Sci. 2019, 24, 625–635. [Google Scholar] [CrossRef]
  2. Sun, Y.; Pri-Tal, O.; Michaeli, D.; Mosquna, A. Evolution of abscisic acid signaling module and its perception. Front. Plant Sci. 2020, 11, 934. [Google Scholar] [CrossRef]
  3. Ma, Y.; Szostkiewicz, I.; Korte, A.; Moes, D.; Yang, Y.; Christmann, A.; Grill, E. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science 2009, 324, 1064–1068. [Google Scholar] [CrossRef]
  4. Park, S.; Fung, P.; Nishimura, N.; Jensen, D.R.; Fujii, H.; Zhao, Y.; Lumba, S.; Santiago, J.; Rodrigues, A.; Chow, T.F.; et al. Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science 2009, 324, 1068–1071. [Google Scholar] [CrossRef]
  5. Fan, W.; Zhao, M.; Li, S.; Bai, X.; Li, J.; Meng, H.; Mu, Z. Contrasting transcriptional responses of PYR1/PYL/RCAR ABA receptors to ABA or dehydration stress between maize seedling leaves and roots. BMC Plant Biol. 2016, 16, 99–112. [Google Scholar] [CrossRef]
  6. Bai, G.; Xie, H.; Yao, H.; Li, F.; Chen, X.; Zhang, Y.; Xiao, B.; Yang, J.; Li, Y.; Yang, D.H. Genome-wide identification and characterization of ABA receptor PYL/RCAR gene family reveals evolution and roles in drought stress in Nicotiana Tabacum. BMC Genom. 2019, 20, 575. [Google Scholar] [CrossRef]
  7. Nian, L.; Zhang, X.; Yi, X.; Liu, X.; Ain, N.U.; Yang, Y.; Li, X.; Haider, F.U.; Zhu, X. Genome-wide identification of ABA receptor PYL/RCAR gene family and their response to cold stress in Medicago sativa L. Physiol. Mol. Biol. Plants 2021, 27, 1979–1995. [Google Scholar] [CrossRef]
  8. Fujii, H.; Chinnusamy, V.; Rodrigues, A.; Rubio, S.; Antoni, R.; Park, S.Y.; Cutler, S.R.; Sheen, J.; Rodriguez, P.L.; Zhu, J.K. In vitro reconstitution of an abscisic acid signalling pathway. Nature 2009, 462, 660–664. [Google Scholar] [CrossRef]
  9. Dittrich, M.; Mueller, H.M.; Bauer, H.; Peirats-Llobet, M.; Rodriguez, P.L.; Geilfus, C.M.; Carpentier, S.C.; Al Rasheid, K.A.S.; Kollist, H.; Merilo, E.; et al. The role of Arabidopsis ABA receptors from the PYR/PYL/RCAR family in stomatal acclimation and closure signal integration. Nat. Plants 2019, 5, 1002–1011. [Google Scholar] [CrossRef]
  10. Yang, Z.; Liu, J.; Poree, F.; Schaeufele, R.; Helmke, H.; Frackenpohl, J.; Lehr, S.; von Koskull-Döring, P.; Christmann, A.; Schnyder, H.; et al. Abscisic acid receptors and coreceptors modulate plant water use efficiency and water productivity. Plant Physiol. 2019, 180, 1066–1080. [Google Scholar] [CrossRef]
  11. Mega, R.; Abe, F.; Kim, J.S.; Tsuboi, Y.; Tanaka, K.; Kobayashi, H.; Sakata, Y.; Hanada, K.; Tsujimoto, H.; Kikuchi, J.; et al. Tuning water-use efficiency and drought tolerance in wheat using abscisic acid receptors. Nat. Plants 2019, 5, 153–159. [Google Scholar] [CrossRef]
  12. Mao, H.; Jian, C.; Cheng, X.; Chen, B.; Mei, F.; Li, F.; Zhang, Y.; Li, S.; Du, L.; Li, T.; et al. The wheat ABA receptor gene TaPYL1-1B contributes to drought tolerance and grain yield by increasing water-use efficiency. Plant Biotechnol. J. 2022, 20, 846–861. [Google Scholar] [CrossRef]
  13. Hewage, K.A.H.; Yang, J.F.; Wang, D.; Hao, G.F.; Yang, G.F.; Zhu, J.K. Chemical manipulation of abscisic acid signaling: A new approach to abiotic and biotic stress management in agriculture. Adv. Sci. 2020, 7, 2001265. [Google Scholar] [CrossRef]
  14. Lozano-Juste, J.; Infantes, L.; Garcia-Maquilon, I.; Ruiz-Partida, R.; Merilo, E.; Benavente, J.L.; Velazquez-Campoy, A.; Coego, A.; Bono, M.; Forment, J.; et al. Structure-guided engineering of a receptor-agonist pair for inducible activation of the ABA adaptive response to drought. Sci. Adv. 2023, 9, e9948. [Google Scholar] [CrossRef] [PubMed]
  15. Park, S.Y.; Peterson, F.C.; Mosquna, A.; Yao, J.; Volkman, B.F.; Cutler, S.R. Agrochemical control of plant water use using engineered abscisic acid receptors. Nature 2015, 520, 545–548. [Google Scholar] [CrossRef] [PubMed]
  16. Lesk, C.; Rowhani, P.; Ramankutty, N. Influence of extreme weather disasters on global crop production. Nature 2016, 529, 84–87. [Google Scholar] [CrossRef] [PubMed]
  17. Zhao, C.; Liu, B.; Piao, S.; Wang, X.; Lobell, D.B.; Huang, Y. Temperature increase reduces global yields of major crops in four independent estimates. Proc. Natl. Acad. Sci. USA 2017, 20, 17–62. [Google Scholar] [CrossRef] [PubMed]
  18. Guihur, A.; Rebeaud, M.E.; Goloubinoff, P. How do plants feel the heat and survive? Trends Biochem. Sci. 2022, 47, 824–838. [Google Scholar] [CrossRef] [PubMed]
  19. Li, N.; Euring, D.; Cha, J.Y.; Lin, Z.; Lu, M.; Huang, L.J.; Kim, W.Y. Plant hormone-mediated regulation of heat tolerance in response to global climate change. Front. Plant Sci. 2021, 11, 627969. [Google Scholar] [CrossRef] [PubMed]
  20. Islam, M.R.; Feng, B.H.; Chen, T.T.; Tao, L.X.; Fu, G.F. Role of abscisic acid in thermal acclimation of plants. J. Plant Biol. 2018, 61, 255–264. [Google Scholar] [CrossRef]
  21. Huang, Y.C.; Niu, C.Y.; Yang, C.R.; Jinn, T.L. The heat stress factor HSFA6b connects ABA signaling and ABA-mediated heat responses. Plant Physiol. 2016, 172, 1182–1199. [Google Scholar] [CrossRef]
  22. Li, H.; Liu, S.S.; Yi, C.Y.; Wang, F.; Zhou, J.; Xia, X.J.; Shi, K.; Zhou, Y.H.; Yu, J.Q. Hydrogen peroxide mediates abscisic acid-induced HSP70 accumulation and heat tolerance in grafted cucumber plants. Plant Cell Environ. 2014, 37, 2768–2780. [Google Scholar] [CrossRef]
  23. Wang, X.; Zhuang, L.; Shi, Y.; Huang, B. Up-regulation of HSFA2c and HSPs by ABA contributing to improved heat tolerance in tall fescue and Arabidopsis. Int. J. Mol. Sci. 2017, 18, 1981. [Google Scholar] [CrossRef]
  24. Hu, X.; Carver, B.F.; Powers, C.; Yan, L.; Zhu, L.; Chen, C. Effectiveness of Genomic Selection by Response to Selection for Winter Wheat Variety Improvement. Plant Genome 2019, 12, 1–15. [Google Scholar] [CrossRef] [PubMed]
  25. Rezaul, I.M.; Baohua, F.; Tingting, C.; Weimeng, F.; Caixia, Z.; Longxing, T.; Guanfu, F. Abscisic acid prevents pollen abortion under high-temperature stress by mediating sugar metabolism in rice spikelets. Physiol. Plant. 2019, 165, 644–663. [Google Scholar] [CrossRef] [PubMed]
  26. Iqbal, N.; Sehar, Z.; Fatma, M.; Umar, S.; Sofo, A.; Khan, N.A. Nitric oxide and abscisic acid mediate heats tress tolerance through regulation of osmolytes and antioxidants to protect photosynthesis and growth in wheat plants. Antioxidants 2022, 11, 372. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, X.; Zhong, X.; Liao, J.; Ji, P.; Yang, J.; Cao, Z.; Duan, X.; Xiong, J.; Wang, Y.; Xu, C.; et al. Exogenous abscisic acid improves grain filling capacity under heat stress by enhancing antioxidative defense capability in rice. BMC Plant Biol. 2023, 23, 619. [Google Scholar] [CrossRef]
  28. Yao, R.; Heinrich, M.; Weckerle, C.S. The genus Lycium as food and medicine: A botanical, ethnobotanical and historical review. J. Ethnopharmacol. 2018, 212, 50–66. [Google Scholar] [CrossRef] [PubMed]
  29. Cao, Y.L.; Li, Y.L.; Fan, Y.F.; Li, Z.; Yoshida, K.; Wang, J.Y.; Ma, X.K.; Wang, N.; Mitsuda, N.; Kotake, T.; et al. Wolfberry genomes and the evolution of Lycium (Solanaceae). Commun. Biol. 2021, 4, 671. [Google Scholar] [CrossRef]
  30. Gong, H.; Rehman, F.; Ma, Y.; Zeng, S.; Yang, T.; Huang, J.; Li, Z.; Wu, D.; Wang, Y. Germplasm resources and strategy for genetic breeding of Lycium species: A Review. Front. Plant Sci. 2022, 13, 802936. [Google Scholar] [CrossRef]
  31. He, W.; Liu, M.; Qin, X.; Liang, A.; Chen, Y.; Yin, Y.; Qin, K.; Mu, Z. Genome-Wide Identification and Expression Analysis of the Aquaporin Gene Family in Lycium barbarum during Fruit Ripening and Seedling Response to Heat Stress. Curr. Issues Mol. Biol. 2022, 44, 5933–5948. [Google Scholar] [CrossRef] [PubMed]
  32. González-Guzmán, M.; Rodríguez, L.; Lorenzo-Orts, L.; Pons, C.; Sarrión-Perdigones, A.; Fernández, M.A.; Peirats-Llobet, M.; Forment, J.; Moreno-Alvero, M.; Cutler, S.R.; et al. Tomato PYR/PYL/RCAR abscisic acid receptors show high expression in root, differential sensitivity to the abscisic acid agonist quinabactin, and the capability to enhance plant drought resistance. J. Exp. Bot. 2014, 65, 4451–4464. [Google Scholar] [CrossRef] [PubMed]
  33. Waadt, R.; Seller, C.A.; Hsu, P.K.; Takahashi, Y.; Munemasa, S.; Schroeder, J.I. Plant hormone regulation of abiotic stress responses. Nat. Rev. Mol. Cell Biol. 2022, 23, 680–694. [Google Scholar] [CrossRef] [PubMed]
  34. Neill, E.M.; Byrd, M.C.R.; Billman, T.; Brandizzi, F.; Stapleton, A.E. Plant growth regulators interact with elevated temperature to alter heat stress signaling via the unfolded protein response in maize. Sci. Rep. 2019, 9, 10392. [Google Scholar] [CrossRef]
  35. Sharma, L.; Priya, M.; Kaushal, N.; Bhandhari, K.; Chaudhary, S.; Dhankher, O.P.; Prasad, P.V.V.; Siddique, K.H.M.; Nayyar, H. Plant growth-regulating molecules as thermoprotectants: Functional relevance and prospects for improving heat tolerance in food crops. J. Exp. Bot. 2020, 71, 569–594. [Google Scholar] [CrossRef]
  36. Li, Q.; Shen, C.; Zhang, Y.; Zhou, Y.; Niu, M.; Wang, H.L.; Lian, C.; Tian, Q.; Mao, W.; Wang, X.; et al. PePYL4 enhances drought tolerance by modulating water-use efficiency and ROS scavenging in Populus. Tree Physiol. 2023, 43, 102–117. [Google Scholar] [CrossRef]
  37. Barros, J.A.; Cavalcanti, J.H.F.; Medeiros, D.; Nunes-Nesi, A.; Avin-Wittenberg, T.; Fernie, A.R.; Araújo, W.L. Autophagy deficiency compromises alternative pathways of respiration following energy deprivation in Arabidopsis thaliana. Plant Physiol. 2017, 175, 62–76. [Google Scholar] [CrossRef]
  38. Wada, S.; Ishida, H.; Izumi, M.; Yoshimoto, K.; Ohsumi, Y.; Mae, T.; Makino, A. Autophagy plays a role in chloroplast degradation during senescence in individually darkened leaves. Plant Physiol. 2008, 149, 885–893. [Google Scholar] [CrossRef]
  39. Kóta, Z.; Horváth, L.I.; Droppa, M.; Horváth, G.; Farkas, T.; Páli, T. Protein assembly and heat stability in developing thylakoid membranes during greening. Proc. Natl. Acad. Sci. USA 2002, 99, 12149–12154. [Google Scholar] [CrossRef]
  40. Jiang, Z.; Zhu, L.; Wang, Q.; Hou, X. Autophagy-Related 2 Regulates Chlorophyll Degradation under Abiotic Stress Conditions in Arabidopsis. Int. J. Mol. Sci. 2020, 21, 4515. [Google Scholar] [CrossRef] [PubMed]
  41. Pruzinská, A.; Tanner, G.; Anders, I.; Roca, M.; Hörtensteiner, S. Chlorophyll breakdown: Pheophorbide a oxygenase is a Rieske-type iron-sulfur protein, encoded by the accelerated cell death 1 gene. Proc. Natl. Acad. Sci. USA 2003, 100, 15259–15264. [Google Scholar] [CrossRef]
  42. Sakuraba, Y.; Kim, D.; Kim, Y.S.; Hörtensteiner, S.; Paek, N.C. Arabidopsis STAYGREEN-LIKE (SGRL) promotes abiotic stress-induced leaf yellowing during vegetative growth. FEBS Lett. 2014, 588, 3830–3837. [Google Scholar] [CrossRef]
  43. Chen, P.; Sun, Y.F.; Kai, W.B.; Liang, B.; Zhang, Y.S.; Zhai, X.W.; Jiang, L.; Du, Y.W.; Leng, P. Interactions of ABA signaling core components (SlPYLs, SlPP2Cs, and SlSnRK2s) in tomato (Solanum lycopersicon). J. Plant Physiol. 2016, 205, 67–74. [Google Scholar] [CrossRef]
  44. Sun, S.; Fan, W.; Mu, Z. The spatio-temporal specificity of PYR1/PYL/RCAR ABA receptors in response to developmental and environmental cues. Plant Signal. Behav. 2017, 12, e1214793. [Google Scholar] [CrossRef] [PubMed]
  45. Tischer, S.V.; Wunschel, C.; Papacek, M.; Kleigrewe, K.; Hofmann, T.; Christmann, A.; Grill, E. Combinatorial interaction network of abscisic acid receptors and coreceptors from Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2017, 114, 10280–10285. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, Y.G.; Fu, F.L.; Yu, H.Q.; Hu, T.; Zhang, Y.Y.; Tao, Y.; Zhu, J.K.; Zhao, Y.; Li, W.C. Interaction network of core ABA signaling components in maize. Plant Mol. Biol. 2018, 96, 245–263. [Google Scholar] [CrossRef] [PubMed]
  47. Lee, U.; Wie, C.; Fernandez, B.O.; Feelisch, M.; Vierling, E. Modulation of nitrosative stress by S-nitrosoglutathione reductase is critical for thermotolerance and plant growth in Arabidopsis. Plant Cell 2008, 20, 786–802. [Google Scholar] [CrossRef] [PubMed]
  48. Guo, Q.Q.; Li, X.; Niu, L.; Jameson, P.; Zhou, W.B. Transcription-associated metabolomic adjustments in maize occur during combined drought and cold stress. Plant Physiol. 2021, 186, 677–695. [Google Scholar] [CrossRef]
  49. Liu, Y.L.; Sun, W.; Zeng, S.H.; Huang, W.J.; Liu, D.; Hu, W.M.; Shen, X.F.; Wang, Y. Virus-induced gene silencing in two novel functional plants, Lycium barbarum L. and Lycium ruthenicum Murr. Sci. Hortic. 2014, 170, 267–274. [Google Scholar] [CrossRef]
  50. Chen, C.J.; Wu, Y.; Li, J.W.; Wang, X.; Zeng, Z.H.; Xu, J.; Liu, Y.L.; Feng, J.T.; Chen, H.; He, Y.H.; et al. TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Mol. Plant 2023, 16, 1733–1742. [Google Scholar] [CrossRef] [PubMed]
  51. Xie, J.M.; Chen, Y.R.; Cai, G.J.; Cai, R.L.; Hu, Z.; Wang, H. Tree Visualization by One Table (tvBOT): A web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res. 2023, 51, 587–592. [Google Scholar] [CrossRef] [PubMed]
Plants 13 00887 g001
Figure 1. ABA level response to heat stress. The ABA content was measured in the leaves of 2–month–old L. barbarum Ningqi No.7, which was subjected to an artificial environment and 42 °C for 1, 3, 6, 12, or 24 h. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Figure 1. ABA level response to heat stress. The ABA content was measured in the leaves of 2–month–old L. barbarum Ningqi No.7, which was subjected to an artificial environment and 42 °C for 1, 3, 6, 12, or 24 h. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Plants 13 00887 g001
Plants 13 00887 g002
Figure 2. Phylogenetic analysis of the 12 LbPYLs with Arabidopsis and tomato homologs. Deduced amino acid sequences were aligned using ClustalX, and the phylogenetic tree was constructed using the bootstrap maximum likelihood tree (1000 replicates) method and MEGA11 (11.0.13) software. The full–length PYL protein sequences, including 12 members from L. barbarum (Lb), 14 members from A. thaliana (At), 14 members from S. lycopersicum (Sl), 14 members from Populus trichocarpa (Pt), and 13 members from Zea mays L. (Zm), were classified into the I, II and III subfamilies. The branches of different classes have altered colors, and each represents a different subfamily. The sequence file information is presented in Table S2.
Figure 2. Phylogenetic analysis of the 12 LbPYLs with Arabidopsis and tomato homologs. Deduced amino acid sequences were aligned using ClustalX, and the phylogenetic tree was constructed using the bootstrap maximum likelihood tree (1000 replicates) method and MEGA11 (11.0.13) software. The full–length PYL protein sequences, including 12 members from L. barbarum (Lb), 14 members from A. thaliana (At), 14 members from S. lycopersicum (Sl), 14 members from Populus trichocarpa (Pt), and 13 members from Zea mays L. (Zm), were classified into the I, II and III subfamilies. The branches of different classes have altered colors, and each represents a different subfamily. The sequence file information is presented in Table S2.
Plants 13 00887 g002
Plants 13 00887 g003
Figure 3. Phylogenic relationship, protein structure conserved domains, and conserved motif analyses. (A) A phylogenetic tree of LbPYL proteins was constructed with MEGA11 (11.0.13) software. (B) Conserved motif distribution of LbPYL proteins. The conserved motifs identified with MEME are displayed with boxes in different colors. A total of 15 motifs were identified and the motifs’ sequences are described in Table S3. The scale at the bottom shows the length of the protein. (C) Predicted conserved structural domains of LbPYL proteins. Gray lines represent the length of each protein sequence, and conserved domains are indicated by colored boxes. (D) Exon–intron structures of the LbPYL genes. Gene structures were analyzed with Tbtools (v2.057) software. Exons and introns are indicated with green boxes and gray lines, respectively. The lengths of the exons and introns (bp) are indicated on the x–axis. The combined figure was illustrated with Tbtools (v2.057) software.
Figure 3. Phylogenic relationship, protein structure conserved domains, and conserved motif analyses. (A) A phylogenetic tree of LbPYL proteins was constructed with MEGA11 (11.0.13) software. (B) Conserved motif distribution of LbPYL proteins. The conserved motifs identified with MEME are displayed with boxes in different colors. A total of 15 motifs were identified and the motifs’ sequences are described in Table S3. The scale at the bottom shows the length of the protein. (C) Predicted conserved structural domains of LbPYL proteins. Gray lines represent the length of each protein sequence, and conserved domains are indicated by colored boxes. (D) Exon–intron structures of the LbPYL genes. Gene structures were analyzed with Tbtools (v2.057) software. Exons and introns are indicated with green boxes and gray lines, respectively. The lengths of the exons and introns (bp) are indicated on the x–axis. The combined figure was illustrated with Tbtools (v2.057) software.
Plants 13 00887 g003
Plants 13 00887 g004
Figure 4. Protein sequence alignment of LbPYLs. Conserved transmembrane domains (TM1–6) and amino acids at NPA domains, ar/R selectivity filters, and Froger’s residues identified in LbPYLs.
Figure 4. Protein sequence alignment of LbPYLs. Conserved transmembrane domains (TM1–6) and amino acids at NPA domains, ar/R selectivity filters, and Froger’s residues identified in LbPYLs.
Plants 13 00887 g004
Plants 13 00887 g005
Figure 5. The cis–acting element on the putative promoter of the LbPYL genes. The identified cis–acting elements were mainly divided into three major categories: phytohormone responsiveness, biotic/abiotic stress, and others. The cis–acting elements sequences are described in Table S4.
Figure 5. The cis–acting element on the putative promoter of the LbPYL genes. The identified cis–acting elements were mainly divided into three major categories: phytohormone responsiveness, biotic/abiotic stress, and others. The cis–acting elements sequences are described in Table S4.
Plants 13 00887 g005
Plants 13 00887 g006
Figure 6. Tissue–specific expression of LbPYL genes. Tissue–specific expression of LbPYLs was determined by quantitative real–time PCR (qRT–PCR) in leaves, young roots, stems, ripening fruits, and flowers with gene–specific primers. qRT–PCR was performed in triplicate, and the fold change was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments.
Figure 6. Tissue–specific expression of LbPYL genes. Tissue–specific expression of LbPYLs was determined by quantitative real–time PCR (qRT–PCR) in leaves, young roots, stems, ripening fruits, and flowers with gene–specific primers. qRT–PCR was performed in triplicate, and the fold change was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments.
Plants 13 00887 g006
Plants 13 00887 g007
Figure 7. Expression profiling response to ABA. The expression levels of LbPYLs were measured in the leaves of one–month–old L. barbarum subjected to 10 μM ABA for 0, 0.5, 1, or 3 h. The fold change in the relative expression level was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Figure 7. Expression profiling response to ABA. The expression levels of LbPYLs were measured in the leaves of one–month–old L. barbarum subjected to 10 μM ABA for 0, 0.5, 1, or 3 h. The fold change in the relative expression level was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Plants 13 00887 g007
Plants 13 00887 g008
Figure 8. Expression profiling response to heat stress. The expression levels of LbPYLs were measured in the leaves of one–month–old L. barbarum subjected to 42 °C for 0, 1, 3, 12, or 24 h. The fold change in the relative expression level was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Figure 8. Expression profiling response to heat stress. The expression levels of LbPYLs were measured in the leaves of one–month–old L. barbarum subjected to 42 °C for 0, 1, 3, 12, or 24 h. The fold change in the relative expression level was analyzed via the 2−ΔΔCT method using the LbACTIN1 gene as an internal control. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Plants 13 00887 g008
Plants 13 00887 g009
Figure 9. Validation in LbPYL10 function by VIGS technology. (A) LbPYL10 gene structure in L. barbarum and pTRV2– LbPYL10 vector map used in this study. White boxes represent 5′– and 3′–UTRs, black boxes represent CDS, and red boxes represent target silent fragment. (B) Expression of LbPYL10 in LbPYL10–silenced (TRV2::LbPYL10) and control (TRV2::00) plants, respectively. (C) Leaf disc phenotype difference between the LbPYL10–silenced seedlings and the empty vector–infected seedlings after heat stress. (D,E) Effect of LbPYL10 silencing on pheophorbide a oxygenase (LbPAO, XM_060313688.1) and STAY–GREEN–like (LbSGRL, XM_060342386.1) after heat stress. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Figure 9. Validation in LbPYL10 function by VIGS technology. (A) LbPYL10 gene structure in L. barbarum and pTRV2– LbPYL10 vector map used in this study. White boxes represent 5′– and 3′–UTRs, black boxes represent CDS, and red boxes represent target silent fragment. (B) Expression of LbPYL10 in LbPYL10–silenced (TRV2::LbPYL10) and control (TRV2::00) plants, respectively. (C) Leaf disc phenotype difference between the LbPYL10–silenced seedlings and the empty vector–infected seedlings after heat stress. (D,E) Effect of LbPYL10 silencing on pheophorbide a oxygenase (LbPAO, XM_060313688.1) and STAY–GREEN–like (LbSGRL, XM_060342386.1) after heat stress. Values are the means of three independent experiments. Different letters indicate significant differences at p < 0.05 (Duncan’s test).
Plants 13 00887 g009
Table 1. Detailed information on 12 PYL genes of L. barbarum and their encoded proteins.
Table 1. Detailed information on 12 PYL genes of L. barbarum and their encoded proteins.
Gene NameVersionGene Length (bp)CDS Length (bp)Deduced Protein
Size (aa)MW (kDa)PI
LbPYL1XM_060355288.1135656418721.0456.16
LbPYL2XM_060315403.1164057319021.3055.74
LbPYL3XM_060316745.197061520422.6466.23
LbPYL4XM_060357010.199263621124.4016.29
LbPYL5XM_060322965.1131572324026.6515.52
LbPYL6XM_060317057.1551353417720.2296.38
LbPYL7XM_060347829.1348955818521.0785.75
LbPYL8XM_060350890.1367056718821.4596.30
LbPYL9XM_060362463.1209256118620.7625.96
LbPYL10XM_060346502.1117862720822.8748.18
LbPYL11XM_060320401.1110374724827.6368.65
LbPYL12XM_060337877.1205259719822.2765.34
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Z.; Liu, J.; Chen, Y.; Liang, A.; He, W.; Qin, X.; Qin, K.; Mu, Z. Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum). Plants 2024, 13, 887. https://doi.org/10.3390/plants13060887

AMA Style

Li Z, Liu J, Chen Y, Liang A, He W, Qin X, Qin K, Mu Z. Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum). Plants. 2024; 13(6):887. https://doi.org/10.3390/plants13060887

Chicago/Turabian Style

Li, Zeyu, Jiyao Liu, Yan Chen, Aihua Liang, Wei He, Xiaoya Qin, Ken Qin, and Zixin Mu. 2024. "Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum)" Plants 13, no. 6: 887. https://doi.org/10.3390/plants13060887

APA Style

Li, Z., Liu, J., Chen, Y., Liang, A., He, W., Qin, X., Qin, K., & Mu, Z. (2024). Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum). Plants, 13(6), 887. https://doi.org/10.3390/plants13060887

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop