- Open Access
In silico characterization, molecular phylogeny, and expression profiling of genes encoding legume lectin-like proteins under various abiotic stresses in Arabidopsis thaliana
BMC Genomics volume 23, Article number: 480 (2022)
Lectin receptor-like kinases (Lec-RLKs), a subfamily of RLKs, have been demonstrated to play an important role in signal transduction from cell wall to the plasma membrane during biotic stresses. Lec-RLKs include legume lectin-like proteins (LLPs), an important group of apoplastic proteins that are expressed in regenerating cell walls and play a role in immune-related responses. However, it is unclear whether LLPs have a function in abiotic stress mitigation and related signaling pathways. Therefore, in this study, we examined the possible role of LLPs in Arabidopsis thaliana (AtLLPs) under various abiotic stresses.
The study was initiated by analyzing the chromosomal localization, gene structure, protein motif, peptide sequence, phylogeny, evolutionary divergence, and sub-cellular localization of AtLLPs. Furthermore, the expression profiling of these AtLLPs was performed using publicly accessible microarray datasets under various abiotic stresses, which indicated that all AtLLPs were differently expressed in both root and shoot tissues in response to abiotic stresses. The cis-regulatory elements (CREs) analysis in 500 bp promoter sequences of AtLLPs suggested the presence of multiple important CREs implicated for regulating abiotic stress responses, which was further supported by expressional correlation analysis between AtLLPs and their CREs cognate transcription factors (TFs). qRT-PCR analysis of these AtLLPs after 2, 6, and 12 h of cold, high light, oxidative (MV), UV-B, wound, and ozone stress revealed that all AtLLPs displayed differential expression patterns in most of the tested stresses, supporting their roles in abiotic stress response and signaling again. Out of these AtLLPs, AT1g53070 and AT5g03350 appeared to be important players. Furthermore, the mutant line of AT5g03350 exhibited higher levels of ROS than wild type plants till 12 h of exposure to high light, MV, UV-B, and wound, whereas its overexpression line exhibited comparatively lower levels of ROS, indicating a positive role of this gene in abiotic stress response in A. thaliana.
This study provides basic insights in the involvement of two important representative AtLLPs, AT1g53070 and AT5g03350, in abiotic stress response. However, further research is needed to determine the specific molecular mechanism of these AtLLPs in abiotic stress mitigation and related signaling pathways in A. thaliana.
Plants are sessile and frequently exposed to various external factors, including biotic and abiotic stresses that adversely affect their growth and productivity. The perception and transmission of these stresses are regulated by sensors and downstream signaling components, including receptor-like kinases (RLKs) [1,2,3,4]. More than 600 RLK genes have been identified in the Arabidopsis thaliana genome, with many of them being classified into subfamilies depending on their extracellular domains . RLKs extracellular domain diversity shows a wide range of activities and signal transduction modalities. To completely understand the functions of RLKs, it is important to know about their interacting ligands and their downstream targets. The adaptive responses of plants to extracellular ligands and stimuli are regulated by a functional continuity between the cell wall (CW) and plasma membrane (PM) in which RLKs with CW-bound extracellular domains presumably serve as cell wall integrity sensors . The lectin-receptor-like kinases (Lec-RLKs), a subfamily of RLKs with extracellular lectin motifs that bind various carbohydrates, are believed to be possible CW-PM linkers. Lec-RLKs have an important role in both plant development and stress responses. They have been linked to seed germination, lateral root growth, pollen development, hormone signaling, defenses against pathogens, and insect pests [2, 5,6,7,8,9,10,11]. Moreover, there is certain evidence on their probable functions in abiotic stresses [12, 13].
The lectin motif has been reported to be involved in stress signaling, development, and defense in plants . As per their lectin domain-type, Lec-RLKs are classified in three groups: GNA-related lectins (G-type Lec-RLKs), legume lectins (L-type Lec-RLKs), and Ca-dependent lectins (C-type Lec-RLKs) [15, 16].
L-type Lec-RLKs are one of the most important groups of Lec-RLKs in A. thaliana; moreover, they are known to play an important role in plant development and defense . The important features of L-type Lec-RLKs are the presence of signal peptide (SP), L-type lectin domain, transmembrane domain, and intracellular kinase domain. In A. thaliana, certain L-type Lec-RLKs are referred to as legume lectin-like proteins (LLPs) because they lack the transmembrane (TM) and kinase domains [1, 18]. As members of the L-type Lec-RLK group, LLPs may function as CW-PM linkers and cell wall integrity sensors during abiotic stress mitigation/signaling. However, their biological significance in A. thaliana has to yet be determined.
In A. thaliana and certain other Brassicaceae species, LLPs have been examined at the transcriptional and proteomic level in response to many biotic stresses. The proteomic analysis of cell wall regeneration in A. thaliana protoplast suspension culture demonstrated that certain LLPs, including AT1g53070, AT3g16530, and AT3g15356, may play an important role in cell wall organization or regeneration . Another LLP, AT5g03350, was reported to be among the top expressed genes in A. thaliana after treatment with salicylic acid . Armijo et al. (2013) reported its induction in A. thaliana against Pseudomonas syringae and named as salicylic acid induced legume lectin like protein 1 (SAI-LLP1) . Moreover, the up-regulation of this gene was reported in A. thaliana after infection with the cucumber mosaic virus (CMV(Y) . These results show that AT5g03350 provides protection against pathogen attack. Certain hormonal regulations such as cytokinins  and sritinol, an artificial inducer for auxin inducible genes , have been investigated on this LLP. In A. thaliana, AT3g15356 also demonstrated increased expression after infection with Alternaria brassicicola . In addition, expressional analysis of Brassica LLPs revealed that they are up-regulated in response to Sclerotinia sclerotiorum and Pieris rapae infection/attack [26, 27]. The increased expression of LLPs was detected in regenerated cell walls, indicating their significance in stress-induced cell wall rearrangement [19, 28, 29]. Another LLP, AT3g16530, demonstrated downregulation in A. thaliana after arsenic related stress , which is unlike the behavior of the aforementioned LLPs under biotic stresses.
To summarize, the abovementioned studies demonstrated that LLPs are related to various biotic stresses that are induced by pathogens and phytohormones [20, 21, 23, 25]. Since abiotic stresses may induce cell wall rearrangement and alter pH, redox, and osmotic balance, there is a possibility that LLPs are participating in abiotic stress mitigation and related signaling pathways. However, how involved LLPs are in abiotic stress mitigation/signaling is not very clear.
Therefore, in this study, we examined the possible role of LLPs in A. thaliana (hereafter AtLLPs) in response to several abiotic stresses. The gene sequences of AtLLPs were first retrieved, and then their chromosomal localization, gene structure, phylogeny, and sub-cellular localization was obtained. The evolutionary divergence of this gene family in the course of evolution was analyzed by comparing evolutionary rate ratio (dN/dS) of AtLLPs and its orthologous pairs in which dN corresponds to non-synonymous substitutions per non-synonymous site and dS corresponds to synonymous substitutions per synonymous site. Furthermore, the expression patterns of these genes were analyzed using publicly accessible microarray datasets from both root and shoot tissues under multiple abiotic stresses. The cis-regulatory elements (CREs) and their associated transcription factors (TFs) were identified in the promoter regions of AtLLPs to better understand their mechanism of transcriptional regulation and identify important regulatory TFs that may affect these gene expressions. Furthermore, the expression patterns of these AtLLP genes were examined using quantitative-reverse transcription (qRT-PCR) after 2, 6, and 12 h duration of cold, high light, oxidative (MV), UV-B, wound, and ozone stress. To characterize the function of AT5g03350, one of the key AtLLPs under tested abiotic stresses, we subsequently evaluated the abiotic stress response in terms of ROS generation in AT5g03350 mutant (T-DNA insertion SALK_036814 line), 35 S::AT5g03350 overexpression and wild type plants after 2, 6, and 12 h of cold, high light, oxidative, UV-B, wound, and ozone stress.
The results of this study will increase our understanding of the involvement of AtLLPs in abiotic stresses and related signaling transduction pathways; moreover, they will provide a solid base for the further functional characterization of AtLLP genes in A. thaliana.
Gene identification and genome-wide physical localization of AtLLPs
AtLLPs were identified on the basis of previous studies [1, 18] and the physical localization of this group of genes across the complete genome of A. thaliana was identified using TAIR chromosomal map tools (www.arabidopsis.org). In A. thaliana, seven AtLLPs, including AT1g07460, AT1g53060, AT1g53070, AT1g53080, AT3g16530, AT5g03350, and AT3g15356 have been identified in total .
Gene structure, protein motif, peptide sequence alignment, phylogenetic relationship, and physicochemical analyses of AtLLPs
The gene sequences of AtLLPs were retrieved from TAIR (https://www.arabidopsis.org/) and their gene structures were analyzed using TBTools 0.665 (https://github.com/CJ-Chen/TBtools) . To predict the functional motif of AtLLPs, protein sequences were downloaded from TAIR and submitted for online MEME suit analysis (https://meme-suite.org/meme/). The peptide sequence alignment of AtLLPs was performed in Jalview using the Clustal Omega program . To develop the phylogenetic tree, the orthologs of AtLLPs present in Marchantia polymorpha, Physcomitrella patens, Selaginella moellendorffii, Amborella trichopoda, Arabidopsis lyrata, Brassica oleracea, Brassica rapa, Daucus carota, Solanum lycopersicum, Solanum tuberosum, Oryza sativa ssp. japonica, Zea mays, Triticum aestivum, and Populus trichocarpa were identified using PLAZA 5.0 . Their peptide sequences were retrieved from NCBI and Phytozome . Peptide sequences were then aligned using the ClustalW program in MEGA 11. The neighbour-joining tree was constructed using 1000 bootstrap replicates, and the physicochemical characterization of protein sequences was performed using ProtParam (https://web.expasy.org/protparam/), a peptide sequence analysis tool .
Sub-cellular localization prediction of AtLLPs
Orthologs identification and evolutionary divergence analysis of AtLLPs
Orthologs of six AtLLPs (except At1G07460.1, which appears to be an outlier among the other AtLLPs) were retrieved from genome database of Arabidopsis lyrata, Brassica rapa, Solanum lycopersicum, Zea mays, Selaginella moellendorffii, Physcomitrium patens, and Chlamydomonas reinhardtii using NCBI genome database. Confirmation of orthologous relationship was done through phylogenetic tree construction using UPGMA method in MEGAX . dN (non-synonymous substitutions per non-synonymous site) and dS (synonymous substitutions per synonymous site) values of AtLLPs and its orthologous pairs were calculated using PAL2NAL (http://www.bork.embl.de/pal2nal/), which has codeml program of PAML package . The dS values of gene pairs so obtained were used to estimate their evolutionary divergence using formula T = dS/(2 × 6.5 × 10 − 9)×10 − 6 MYA. Divergence time of AtLLPs and its orthologs from plant family was compared with evolutionary time scale of speciation events (http://www.timetree.org/) .
In silico expression analysis of AtLLPs under various abiotic stress conditions
Publicly accessible microarray-based gene expression data were used for the in silico expression investigation of six AtLLPs under diverse abiotic stressors, with the exception of AT3g15356, which lacks expression data in the database. Nine abiotic stress conditions, including cold, osmotic, salt, drought, genotoxic, oxidative, UV-B, wounding, and heat stress, were selected to examine the expression of AtLLPs. Therefore, nine AtGeneExpress (Stress Series) samples were retrieved from Bio-Analytic Resource, e-Northerns w. Expression Browser (http://bar.utoronto.ca; , and the array was normalized with a TGT value of 100 using the Gene Chip Operating Software (GCOS). The datasets comprised an average of replicate treatments compared to the average of appropriate control and output values in the table and image were log2-transformed ratios. For differentially expressed genes (DEG) analysis, datasets comprised seven data points of nine samples that were extracted in notepad, and then the differential expression of six selected genes was depicted in a heatmap using default criteria of the TBtools software version 0.665 . The experimental conditions for the extracted expression data from Bio-Analytic Resource, e-Northerns w. Expression Browser (http://bar.utoronto.ca; ) were as follows: A. thaliana wild-type (col-0) seeds were spread in magenta boxes containing the Murashige and Skoog (MS) medium. After two days of incubation at 4 °C in dark, the samples were transferred to the growth chamber (photoperiod of 16/8 h, temperature of 24 °C, relative humidity of 50%, and light intensity of 150 µE/cm2 s). The above-mentioned stresses were applied to 16-day-old seedlings for 0.5, 1.0, 3.0, 6.0, 12.0, and 24.0 h, with control samples at 0 h. For tissue-specific expression, roots and shoots were separated, and all treatments were performed on the same batch of seedlings.
Promoter analysis, functional annotation of cis-regulatory elements (CREs), and identification of cognate TFs in AtLLPs
To understand the promoter structure and distribution of CREs, 500 bp non-coding sequences that were upstream of the transcription start site (TSS) of candidate genes were retrieved from TAIR and analyzed in PlantRegMap (http://plantregmap.gao-lab.org/) . The stress-related CREs was identified and displayed using TBtools 0.665 . Cognate TFs that bind to CREs at the promoter regions of AtLLPs were selected from PlantRegMap (http://plantregmap.gao-lab.org/) .
Expressional correlation analysis between AtLLPs and their CREs cognate TFs
The normalized expression data for AtLLPs and their CREs cognate TFs was retrieved from e-Northerns w. Expression Browser as described in the in silico expression analysis of AtLLPs under similar nine abiotic stress conditions. For six AtLLPs, the TFs were chosen based on the availability of their expression in the dataset. The seventh AtLLP, AT3g15356, is likewise not included in this analysis due to a lack of expression data in the database. The normalized expression data of each AtLLPs was plotted against respective TFs data, and the Pearson correlation coefficient (r) value was calculated. R values of > 0.8 and < − 0.8 were considered as significant. TFs that show significant expression correlation with corresponding AtLLPs at least in a single stress condition was plotted as heatmap using Morpheus (https://software.broadinstitute.org/morpheus, version 1.0).
Plant material, growth condition and abiotic stress treatment
The seeds of A. thaliana ecotype Columbia-0 (Col-0) were obtained from the Arabidopsis Biological Resource Centre (ABRC), The Ohio State University, USA. Seeds were sterilized for 10 min in a solution of freshly prepared 5% calcium hypochorite and 0.02% Triton-x 100. After being washed 4–5 times with sterilized water, seeds were plated on half-strength MS medium with 1% sucrose and 0.8% plant agar. The plates were then maintained at 4 °C for three days to achieve synchronized and higher seed germination. Then, plates were placed in a plant growth chamber (Conviron, Adaptis 1000) at 23 °C, bright fluorescent light of 150 µmol photons m− 2 s− 1, relative humidity of 65% and 16 h light: 8 h dark cycle.
For stress treatment, 10-day-old seedlings were subjected to six abiotic stresses, namely, cold (8 °C), high light (500 µmol photons m− 2 s− 1), oxidative (25 µM, MV), UV-B (0.99 W m− 2 s− 1), wound, and ozone (40 ppb) for a duration of 2, 6, and 12 h in triplicates, as previously described by Mondal et al. (2021) . In addition to each stress treatment, three similar sets of seedlings were maintained in their normal growing condition and used as a control.
RNA isolation and cDNA synthesis
Total RNA was isolated from stress-exposed seedlings at selected time points, including 0 (control), 2, 6, and 12 h using the RNeasy plant mini kit (Qiagen, USA) as per manufacturer’s instructions. The quantity and quality of isolated RNA were determined by spectrophotometry (Nanodrop 2000, Thermo Fisher Scientific, Waltham, MA, USA) and formaldehyde-based gel electrophoresis, respectively (Additional File 1: Fig S1). For cDNA synthesis, 1 µg of total RNA was reverse-transcribed in 20 µl using Revert Aid First Strand cDNA Synthesis Kit (Fermentas Life Sciences, USA) and oligo (dT) primers according to the manufacturer’s instructions.
Expression analysis of AtLLPs under six selected abiotic stresses
To examine the temporal expression patterns of AtLLPs under six different abiotic stresses, qRT-PCR was performed at 2, 6, and 12 h of exposure to cold, high light, oxidative, UV-B, wound, and ozone stress. qRT-PCR was performed in CFX-96 Real-time PCR Detection System (Bio-Rad, USA). Reactions were conducted in a total volume of 20 µl using 50 ng of cDNA, 0.5 to 1 pmol of forward and reverse primers, and 10 µl of 2x Sso Fast Eva Green qPCR Supermix (Bio-Rad, USA). The cycling conditions were as per the manufacturer’s protocol with a primer-specific annealing temperature. The threshold cycle (Ct) was automatically determined for each reaction using the system set with default parameters. The transcript levels were normalized to AtUBQ5 and AtAPT1 transcript and the fold differences of each amplified product in the samples was calculated using the 2−ΔΔCt method. All primer sequences used in this study were designed using the NCBI-Primer blast online tool (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) and are listed in Additional File 2: Table S1.
Vector construction, generation of At5g03350::overexpression plants and selection of homozygous lines
Gateway cloning strategy was employed to clone the CDS of AT5g03350. Primers were design with Vector NTI Software (ThermoFischer scientific) with flanking attB1 and attB2 sequences (Additional File 2: Table S1). attB-flanked PCR products were amplified in 25 µl PCR reactions that contained 5 µl of 5X Phusion HF Buffer (Thermo scientific), 0.5 µl of 10 mM dNTP mix (Thermo-R0191), 0.5 µl of 10 mM forward and reverse primers, 100 ng (1 µl) cDNA, 0.25 µl of Phusion™ DNA Polymerase (2 U µl− 1), and 17.25 µl nuclease free water (Sigma). Amplified PCR product of 885 bp was purified using QIAquick gel extraction kit (QIAGEN, cat. no. 28,704) as per the manufacturer instruction (Additional File 3: Fig. S2 A, B).
To generate entry clone, BP reaction was performed, in which 2 µl (75 ng) attB-flanked PCR product, 1 µl (150 ng) pDONR207 (Invitrogen, USA), 5 µl of TE buffer, and 1 µl Gateway 2 µl of BP Clonase™ II enzyme (Invitrogen, USA) were mixed and incubation at 25 °C for 1 h. The reaction was terminated with the addition of 1 µl Protease K and incubation at 37 °C for 15 min. Subsequently 1 µl of BP reaction was transformed in DH5-alpha competent cells. Transformed colonies were selected using gentamycin (10 µg ml− 1) on LB-plates. The transformed colonies were PCR verified for AT5g03350:CDS using same primers used for CDS amplification (Additional File 3: Fig. S2 C, D). The entry clone AT5g03350 CDS- pDONR207 plasmid was isolated using the Qiagen Plasmid mini kit (catalogue no. 12,123) according to the manufacturer’s instructions and further confirmed by DNA sequencing using the ATTL1 and ATLL2 primers (Additional File 2: Table S1; Additional File 4 and 5). Subsequently expression clone was generated by LR reaction which include 1 µl (100 ng µl− 1) entry clone (AT5g03350 CDS-pDONR207), 1 µl (150 ng µl− 1) pMDC32 destination vector  4 µl TE buffer and 2 µl LR Clonase II Enzyme mix (Invitrogen). The reaction mixture was incubated as described for BP reaction. The LR reaction was transformed in DH5-alpha competent cells and positive colonies were selected on 50 µg ml− 1 kanamycin. The selected colonies were cross-validated using PCR as described for BP reaction (Additional File 3: Fig S2 C, D). The 35 S::AT5g03350 CDS-pMDC32 clone was isolated using Qiagen Plasmid mini kit (cat. No. 12,123) and used for transformation of Agrobacterium tumefaciens GV3101 cells.
Electroporation strategy was used to transform the Agrobacterium cells. 100 µL electrocompetent Agrobacterium GV3101 (helper plasmid pMP90GR) cells and 50–100 ng of 35 S::AT5g03350 CDS-pMDC32 clone were mixed and placed in electroporation cuvette (Gene Pulser Cuvettes, Bio-Rad) on ice for 30 min. Gene pulser (Bio-Rad) was set to 2.5 kV, 200 Ω, and 25 µFD. Cuvette containing Agrobacterium cells and expression clone was placed in the Gene pulser and a short pulse was applied. Immediately 800 µL LB was added to cuvette. Electroporated Agrobacterium cells was transferred to a microcentrifuge tube incubated in shaker at 28 °C for 3 h. 200 µL aliquot was spread an on a LB plate containing 10 µg ml− 1 rifampicin, 10 µg ml− 1 gentamycin, and 50 µg ml− 1 kanamycin and positive transformed colonies were selected. Colony-PCR was conducted to cross-validate the successful transformation.
A. thaliana Col-0 wild type plants were transformed with Agrobacterium cells carrying the 35 S::AT5g03350 CDS-pMDC32 expression clone (Additional File 3: Fig S2 E) using floral-dip method as described earlier . Antibiotic selection of transformed plants was done using modified method of Harrison et al. (2006) . The collected seeds of floral-dipped plants were surface sterilized with 5% calcium hypochlorite for 10 min and spread on ½ MS plates having 1% sucrose and 0.8% plant agar with 10 µg ml− 1 hygromycin (B) The plates were then incubated in 4 °C for stratification and synchronized germination. After 48 h of stratification, the plates were placed in plant growth chamber (Conviron, Adaptis 1000) at 23 °C, 150 µmol photons m− 2 s− 1 for 10–12 h and then incubated in dark at 23 °C for 3 days. The plates are then transferred to plant growth chamber at 23 °C temperature, 150 µmol photons m− 2s− 1, 60% relative humidity and 16 h light/ 8 h dark cycle. After another 2–3 days of growth, the hygromycin resistant seedlings were distinguished from the susceptible ones by longer hypocotyls and greener cotyledons (Additional File 6: Fig. S3A). Following that, the selected seedlings were transferred to ½ MS plates with 1% sucrose and 0.8% plant agar and grown normally in plant growth chamber for another 14 days. The seedlings with true leaves were then transferred to potted soil and grown for next generation seed collection. Transformed T0 plants were PCR verified using hygromycin phosphotransferase (hpt) gene specific primers (Additional File 6: Fig. S3A). Genomic DNA form each T0 plant was isolated using modified protocol of Murray et al. (1980) . For copy number identification and homozygosity test in T1, T2, and T3 generation, the modified method of Kihara et al. (2006) was used . Amplification of endogenous single-copy gene 4-hydroxyphenylpyruvate dioxygenase (4hppd) is used as a reference against the amplification of hpt gene. Briefly, PCR amplification of both the genes were done separately using DreamTaq PCR kit (DreamTaq, Thermo scientific) as per the manufacturer’s instrcutions from genomic DNA of each plant with amplification cycle of 28 in a 25 µl of reaction volume with same concentration of primers and templates (~ 150 ng DNA). 12 µl of each reaction was run in 1.5% (w/v) agarose gel (UltraPure™ Agarose, Invitrogen, 16500-500) and image was captured using chemidoc (ImageQuant LAS 500). Band intensity was quantified using ImageJ software (https://imagej.nih.gov/ij/download.html). Band intensity ratio of 1:1 between HPT and 4HPPD amplification for a line in T1, T2, and T3 generation confirms homozygosity with single insertion, 0.5:1 ratio confirms heterozygosity with single insertion and 2:1 suggests homozygosity with double insertion (Additional File 6: Fig S3 B, C, and D) For the following experiments, seeds of T3 homozygous single insertion line1 was used.
Screening of At5g03350 gene specific SALK_036814 T-DNA insertion line
AT5G03350 locus specific t-DNA insertion lines, N536814 (SALK_036814) were procured from the Nottingham Arabidopsis Stock Centre (NASC) (http://arabidopsis.info/). The seeds were surface sterilized with 5% calcium hypochlorite for 10 min and spread on ½ MS plates with 1% sucrose and 0.8% plant agar. The plates were then incubated at 4 °C for stratification and synchronized germination. After 72 h of stratification, the plates were placed in plant growth chamber (Conviron, Adaptis 1000) at 23 °C, 150 µmol photons m− 2 s− 1, 60% relative humidity, and 16 h light/8 h dark cycle. 14 days old seedlings were transferred to soil and grown for another 7 days. The 21 days old plant leaves were used for genomic DNA isolation using modified cTAB method . The primers for gene specific amplification (SALK_036814.56.00.x forward and reverse) were designed by T-DNA primer design tool (http://signal.salk.edu/tdnaprimers.2.html) and the LBP1.3 was used for T-DNA plus genomic DNA hybrid amplification (Additional File 2: Table S1). Predicted product length of PCR amplification of genomic DNA of this mutant line with SALK_036814.56.00.x forward and reverse primers was 1190 bp whereas, in wild type or heterozygous mutant plants LBP1.3 and SALK_036814.56.00.x RP primer specific product size was 598–898 bp. We have set two paired reactions for each genotypic study because of difference in Tm of LBP1.3 and gene specific primers (Additional File 7: Fig. S4).
Confirmation of At5g03350 overexpression and At5g03350 mutant (SALK_036814) lines using qRT-PCR
To confirm the 35 S::At5g03350 overexpression and At5g03350 mutant lines, the transcript levels of the At5g03350 gene were examined and compared with Col-wild type plants using qRT-PCR, as described in the section on expression analysis of AtLLPs under six selected abiotic stresses.
Detection of ROS production in At5g03350 mutant, At5g03350 overexpression and wild type plants under six chosen abiotic stresses
The selected abiotic stresses, cold, high light, oxidative, UV-B, wound, and ozone are known to enhance the production of ROS, such as singlet oxygen (1O2), superoxide radical (O2•−), hydrogen peroxide (H2O2), and hydroxyl radical (OH•) in plants [47,48,49,50,51]. To assess differences in the abiotic stress response of At5g03350 mutant (SALK_036814), At5g03350 overexpression, and wild type plants under the aforesaid abiotic stresses, the presence of two classical stress markers i.e. superoxide radical and hydrogen peroxide were detected using nitroblue tetrazolium (NBT) and 3, 3’diaminobenzidine (DAB), respectively were monitored in the control and stress exposed seedlings at selected time points, 0 (control), 2, 6, and 12 h as previously described by Mondal et al. (2021) .
Data were expressed as mean ± standard deviation (SD) of at least three biological replicates. The results of expression data were statistically examined using one-way analysis of variance (ANOVA), followed by Student–Newman–Keuls test using SigmaPlot 12 (p < 0.05) to determine significant differences in control- and stress-treated samples .
Chromosomal localization, gene structure, protein motif, peptide sequence, and phylogeny analyses
Chromosomal localization investigation revealed that out of seven, four AtLLPs, AT1g53060, AT1g53070, AT1g07460, and AT1g53080, are located on chromosome 1 while two others, AT3g16530 and AT3g15356 are located on chromosomes 3 and one last, AT5g03350 is found on chromosome number 5 alone (Fig. 1 A). The localization study suggested that the AT1g53060, AT1g53070, and AT1g53080 genes on chromosome 1 were tandem-duplicated. Gene structure analysis demonstrated that all AtLLPs have a similar basic gene structure and are intronless (Fig. 1B). Although there is considerable variation in both the 5′ and 3′ UTR structure, this indicates that they are more flexible in terms of their functional diversity (Fig. 1B). Motif analysis revealed the presence of eight conserved motifs in most of the AtLLPs, all of which are arranged in the same order (Fig. 1B, C). The sequence alignments of peptide sequences in AtLLPs showed that lectin domains are substantially conserved while signal peptides (SP) are less conserved (Additional File 8: Fig. S5). AT1g53060 and AT1g07460 lack SP, whereas the other five AtLLPs include SP of varying lengths that are partially similar in sequence (Additional File 8: Fig. S5).
To understand the evolutionary relationship of AtLLPs with other plant groups, phylogenetic analysis was conducted and presented in the form of phylogenetic tree (Fig. 2). In the phylogenetic tree, lower non-seed plants made Group I, LLPs of Brassicaceae made a separate group as Group II, dicots made Group III, and monocots on Group IV (Fig. 2). Interestingly, AtLLPs are closely related to the non-seed group of plants rather than other dicots (Fig. 2). Although most of the LLPs are found to be closely related to non-seed plants, AT1g07460 and its orthlogs from members of other Brassicaceae species made a small group in Group III (highlighted in red bracket) and are more closely related to monocots (Fig. 2).
To have an improved understanding of AtLLPs, the in silico physicochemical analysis was conducted (Additional File 9: Table S2). With the exception of AT1g53060 and AT1g07460, which has a length of 242 and 258 amino acids (AA) respectively, AtLLPs are stable and have a length of 272–283 AA (Additional File 9: Table S2). AtLLPs have aliphatic indexes ranging from 69.26 to 80.50; however, their grand average hydropathicity is extremely diverse (Additional File 9: Table S2). AT3g16530 is nearly neutral in charge with a pI of 6.98 and AT1g07460 is acidic in nature with pI of 4.80, whereas the others are basic in nature.
Sub-cellular localization prediction suggested AtLLPs could be located in intracellular compartments
Based on the SUBAcon algorithm, subcellular localization analysis of AtLLPs revealed that all investigated AtLLPs, with the exception of AT1g53060 and AT1g07460, are found in the extracellular spaces (Additional File 10: Table S3). This may be due to a lack of signal peptide in AT1g53060 and AT1g07460, as shown in Additional File 8: Fig. S5. Moreover, the locations of different AtLLPs have been predicted in different organelles such as, nucleus, mitochondria, and plastid (Additional File 10: Table S3). The diversity in the locations of AtLLPs shows the versatility in their functions.
Evolutionary divergence analysis of AtLLPs gene family and its orthologs
The evolutionary divergence of six AtLLPs gene family with seven additional model plant species, including Arabidopsis lyrata, Brassica rapa, Solanum lycopersicum, Zea mays, Selaginella moellendorffii, Physcomitrium patens, and Chlamydomonas reinhardtii, was analyzed. The dN/dS ratios of LLP genes and their orthologous pairs were calculated to detect evolutionary pressure acting on these genes (Additional File 11: Table S4) according to Kryazhimskiy and Plotkint, (2008) . A total 42 orthologs of six AtLLP genes were found from Arabidopsis lyrata, Brassica rapa, Solanum lycopersicum, Zea mays, Selaginella moellendorffii, Physcomitrium patens, and Chlamydomonas reinhardtii. Chlamydomonas reinhardtii showed one orthologous gene in its genome. The other two, Physcomitrium patens and Selaginella moellendorffii showed three and four orthologs, respectively. In higher plant species such as Arabidopsis lyrata, Brassica rapa, Solanum lycopersicum, Zea mays 5, 4, 6 and 5 orthologous genes were found, respectively.
For evolutionary pressure quantification, we estimated values of dN/dS ratios of the gene pairs and this was followed by estimation of evolutionary divergence time (Table 1, Additional File 11: Table S4). The dN/dS ratio of the LLP gene orthologs from Angiospermic members was found to be in the range of 0.0109–0.4634, whereas the ratio values ranged 0.0114–0.0463 when compared with lower plant species, Selaginella moellendorffii, Physcomitrium patens and Chlamydomonas reinhardtii (Additional File 11: Table S4). This analysis shows that the values of dN/dS ratio falling less than 1 (Additional File 11: Table S4), suggesting that LLP genes are under purifying (negative) selection, i.e. natural selection here suppresses protein changes.
In silico expression analysis revealed that AtLLPs differently express under various abiotic stresses
Based on publicly available microarray datasets, AtLLPs were subjected to in silico expression analysis under various abiotic stresses, including cold, osmotic, salt, drought, genotoxic, oxidative, UV-B, wound, and heat stress (Fig. 3). The comprehensive analysis of expression profile from the microarray datasets demonstrated that a majority of tested AtLLPs expressed at high/intermediate levels in both root and shoot tissues under studied abiotic stresses (Fig. 3). As per this differential expression pattern, AtLLPs may be regulated by abiotic stresses and implicated in stress management.
The promoter region of AtLLPs contains CREs involved in growth and development, hormone, and multi-stress response
To understand transcriptional regulation of seven AtLLPs under various abiotic stresses, we analyzed the abundance of potential CREs using 500 bp upstream sequences of these genes (Fig. 4). The promoter structure analysis suggests that wide classes of CREs spanned across the 500 bp AtLLPs promoter. Dot plot analysis with functional annotation revealed that ~ 44% of CREs were involved in growth and development including WOX, AP2, bHLH, bZIP, C2H2, C3H, CPP, Dof, GRAS, LBD, RAV, SBP, SRS, TCP, and ZF-HD and ~ 27% modulate multiple stress (VOZ, CAMTA, HD-ZIP, MYB, MYB related, and NAC) (Fig. 4 A). Furthermore, different CREs regulated by plant hormones such as auxin (ARF and B3), cytokinin (GeBP), ethylene (EIL and ERF), salicylic acid (WRKY), and brassinosteroids (BES1) and certain important TFs regulates underlying response to light (FAR1, GATA), heat (HSF), cell cycle (E2F/DP), nitrate signaling (Nin-like), and identity of floral organ (MIKC_MADS) were also identified. Additionally, we have predicted the physical localization of motif of CREs across the promoter length (Fig. 4B). The proportional distribution pattern of different identified CREs and their functional annotations are represented in the pie chart (Fig. 4 C).
Expressional correlation coefficients analysis between AtLLPs and their CREs cognate TFs supported the role of AtLLPs in different abiotic stress conditions
To strengthen our CREs-based results and validate the role of AtLLPs in multiple abiotic stresses, we predicted the cognate TFs that can bind to the CREs present in 500 bp upstream region of TSS. Additional File 12: Table S5 lists the predicted cognate TFs from major families mentioned in the previous section that can bind to the identified CREs in each of the AtLLPs. Based on publicly available in silico expression data of different abiotic stress series, used for in silico expression profiling of AtLLPs (Fig. 3), we calculated the expressional correlation of coefficients (r-value) between AtLLPs and their CREs specific TFs (Additional File 13: Table S6), to identify probable TFs that regulate the expression of AtLLPs under different abiotic stresses in both root and shoot tissues. This analysis showed that under different abiotic stresses, the tested AtLLPs have a positive or negative correlation with their CREs cognate TFs in both root and shoot tissue (Additional File 13: Table S6). Among these TFs, only a few of them, including abiotic stress-responsive TFs, exhibited significant positive or negative expressional correlation [± 0.8 (p < 0.05)] (Additional File 14: Table S7). The TFs that have demonstrated a significant expressional correlation with AtLLPs at least in a single tested abiotic stress are presented in the form of heatmap in shoot and root tissues (Figs. 5 and 6). The analysis demonstrated that the primary regulatory TFs are different for roots and shoots and for different abiotic stresses. The presence of ~ 27% abiotic stress-responsive CREs in the promoter of tested AtLLPs and their significant expressional correlation with CREs cognate TFs suggested the tight regulation of AtLLPs under different abiotic stresses in both root and shoot tissue. These results suggest that AtLLP genes in A. thaliana may play various functions, including abiotic stress response.
Transcript analysis confirmed that AtLLPs are responsive to various abiotic stresses
Six abiotic stresses, including cold, high light, oxidative, UV-B, wound, and ozone stress, were selected to examine the temporal gene expression of seven AtLLPs. Most AtLLPs exhibited differential expression in all the tested stresses at 2, 6, and 12 h of exposure compared to control (Fig. 7; Additional File 15: Table S8).
The relative expression of AT1g53060 was significantly upregulated by > 1.5 fold at 2 h of cold, ozone, and UV-B. However, after 6 h of treatment, the expression was downregulated and reached to 0.27-, 0.85-, and 0.27-fold in cold, ozone, and UV-B treatment, respectively. For MV (oxidative) treatment, the expressional changes of this gene at all tested time points were insignificant. Moreover, AT1g53060 demonstrated significant downregulation under high light stress condition with the relative expression of 0.53, 0.63, and 0.42 fold at 2, 6, and 12 h of treatment, respectively. Furthermore, after wounding, the relative expression of this gene was insignificantly downregulated for the first 6 h; however, after 12 h of treatment, it showed a significant downregulation of 0.56-fold.
For AT1g53070, it shows a common pattern of expression under MV, ozone, and wounding. Changes in the relative expression of AT1g53070 was insignificant after 2 and 6 h of treatment; however, after 12 h of treatment, the gene expression was significantly upregulated to 1.81-, 2.73- and 1.53-fold in MV, ozone, and wounding, respectively. UV-B and high light stress initially downregulated the expression of this gene but reached a control level after 12 h of treatment. Under UV-B, the relative expression of AT1g53070 was 0.55-fold at 2 h; however, after 6 h, it reached the control level and remained the same till 12 h of treatment. Moreover, under high light, the relative expression of this gene was 0.55- and 0.64-fold at 2 and 6 h, respectively, and reached the control level after 12 h of treatment. Unlike the abovementioned expression patterns, this gene shows continuous downregulation under cold stress where the relative expression reached 0.51-fold after 12 h of treatment.
The relative expression of AT1g53080 and AT3g16530 did not show any distinct patterns under the studied abiotic stresses. Although after 12 h of MV exposure, AT1g53080 demonstrated significant upregulation (2.0-fold) in its expression. Furthermore, AT3g16530 showed a substantial increase in expression of 10.0- and 6.0-fold after 2 and 6 h of cold exposure, respectively. No significant and consistent changes of expression were observed for these two genes under other tested conditions.
Unlike these genes, AT5g03350 demonstrated the highest expressional upregulation after 6 and 12 h of treatments for most stresses. UV-B treatment did not significantly change the gene expression of this gene. Moreover, under ozone stress, 1.57 fold of expression was noticed at 6 h of treatment and then the declined expression reached to 0.26 fold after 12 h of treatment. The expression pattern of this gene was demonstrated to be similar under light, MV, and wounding. In these three stresses, the maximum increase of 2.22-, 2.01- and 14.19-fold were noticed after 12 h of stress treatment.
AT1g07460 did not display a clear pattern of expression during cold, high light, MV, and UV-B. However, under ozone and wounding stress, it demonstrated maximal expression of 18.71- and 3.83-fold after 12 and 6 h, respectively. The seventh gene in this family, AT3g15356, showed a substantial increase in expression of 10.57-, 6.96-, 10.91-, 4.19-, and 2.67-fold after 2 h of exposure to cold, high light, ozone, UV-B, and wounding, respectively. Its expression reduced as the duration of stress exposures increased, with the exception of wounding, when it rose 5.85-fold after 2 h.
The critical evaluation of the expression pattern of all the AtLLPs suggested that out of these seven, AT1g53070 and AT5g03350 are important candidates that may have a prominent role in abiotic stress response and related signaling pathways.
At5g03350 mutant (SALK_036814) plants produced more ROS than wild type plants, whereas 35 S::At5g03350 overexpression plants produced less under six abiotic stresses
At5g03350, one of the key genes under evaluated abiotic stresses among seven AtLLPs (Fig. 7; Additional File 15: Table S8), was selected for further functional characterization to validate AtLLPs involvement in abiotic stress tolerance. For that, ROS generation in terms of two stress indicators, superoxide radical (O2•−), hydrogen peroxide (H2O2), was examined in At5g03350: mutant (SALK_036814), 35 S::At5g03350 overexpression and Col-0 wild type plants subjected to the same six abiotic stresses, cold, high light, oxidative, UV-B, wound, and ozone, which were chosen for the transcript analysis of AtLLPs (Fig. 7; Additional File 15: Table S8). However, before examining abiotic stress-induced ROS generation in At5g03350: mutant (SALK_036814), 35 S::At5g03350 overexpression and Col-0 wild type plants, the expression levels of the At5g03350 gene were examined by qRT-PCR to evaluate relative expressional fold difference (Fig. 8). Data showed that the expression level of the At5g03350 gene was significantly reduced to ~ 4.0-fold in At5g03350 mutant plants, but significantly increased to ~ 5.0-fold in At5g03350 overexpression plants compared to wild type plants (Fig. 8). Moreover, expressional difference between mutant and overexpression plants is ~ 20 fold.
The production of superoxide radical (O2•−), hydrogen peroxide (H2O2) was detected using nitroblue tetrazolium (NBT) and 3,3- diaminobenzidine (DAB) staining, respectively in At5g03350: mutant (SALK_036814), 35 S::At5g03350 overexpression and Col-0 wild type plant seedlings after 0, 2, 6, and 12 h of aforesaid abiotic stresses (Fig. 9 A,B). Among the tested stresses, under high light, MV, UV-B, and wound, mutant seedlings showed higher stain colour intensity than wild type till the 12 h of treatment, indicating higher ROS accumulation, whereas overexpression seedlings showed lower stain colour intensity than wild type till the 12 h of treatment, indicating lower ROS accumulation (Fig. 9 A,B). In the case of cold and ozone stress, however, the observable differences in stain colour intensity are quite minimal (Fig. 9 A,B). According to these findings, At5g03350 gene has a protective role in A. thaliana against high light, MV, UV-B, and wound stress.
The transmission of environmental cues from the CW to the PM is an important aspect of stress perception, mitigation, and signaling. Apoplastic regions play an important function in information transmission. Environmental signals are received and recognized by the plant signaling pathways that may transmit these signals to the nucleus and finally result in the induction of stress-specific transcriptional and physiological responses. Plants have a large number of CW and PM receptors that detect environmental stimuli and initiate adaptive responses via downstream signal cascades. Among them, Lec-RLKs are known to be important for transmitting signals from the CW to the PM under stressful conditions. A group of Lec-RLKs, known as LLP, lack TM and kinase domains, indicating their full extracellular nature; moreover, they may act as a labile sensor/carrier for signal transmission from CW to PM. Because LLPs share structural similarity with the ectodomains of RLKs, which are well known to play a role for biotic and abiotic stresses [4, 54]. Therefore, it could be possible that the LLP subfamily of Lec-RLKs in A. thaliana may have certain role in abiotic stress mitigation/signaling. Therefore, the primary focus of this study is to examine the possible role of LLPs in A. thaliana (AtLLPs) under different abiotic stresses using in silico and wet lab experiments.
Since the molecular regulation and function of AtLLPs to multiple abiotic stresses is unclear, we investigated structural and functional aspects of AtLLPs in terms of chromosomal localization gene and protein structure, biochemical properties, phylogenetic relationship, sub-cellular localization, evolutionary divergence, in silico expression, promoter sequences, and transcript analyses under several abiotic stresses.
Most plant gene families are attributed to gene duplication either as an individual chromosomal duplication or segmental duplication in the same chromosome . The chromosomal distribution patterns of AtLLPs show that the expansion of these genes in A. thaliana may have occurred because of simultaneous tandem gene duplications on chromosome 1 (Fig. 1 A). Such duplicated genes are considered to facilitate the rapid neo-functionalization that helps plants acclimate to harsh environmental conditions. Furthermore, the gene structure analysis of AtLLPs shows that introns are completely absent in them (Fig. 1B). The lack of introns may cause them to be less stable; however, it makes them capable of providing rapid functioning capabilities by avoiding splicing. Similar explanations have been proposed for other genes, in which intronless genes have been reported to have lower RNA stability . Peptide sequence alignments and protein motif analyses demonstrated that AtLLPs have a highly conserved arrangement of motifs; however, their sequences are not similar (Fig. 1B,C, Additional File 8: Fig. S5). The signal peptides of these proteins were different in length and sequence (Additional File 8: Fig. S5). These results show a common mode of action for all AtLLPs, which might be spatially distant and distributed throughout the cell wall, intercellular space, apoplastic fluid, and other cellular organelles. This result is corroborated by sub-cellular localization analysis where AtLLPs are expected to be present in chloroplast, mitochondria, and nucleus (Additional File 10: Table S3).
Phylogenetic analysis revealed that AtLLPs except AT1g07460 are more closely related to their orthologs present in lower plants, i.e., bryophytes and pteridophytes rather than to other dicots and monocots (Fig. 2). The presence of orthologs of AtLLPs in early land plants such as bryophytes suggest that the related genes played an important role in combating against biotic and abiotic stresses. Subsequently, in the due course of evolution, functional diversification had taken place in these LLPs, which lead to the separation into dicot and monocots, which is nicely explained by the results of our phylogenetic analysis (Fig. 2).
Protein sequences and structures diverge from one another as they evolve from a common ancestor [57, 58]. The evolutionary rate ratio (dN/dS) is the oldest and most extensively used approach for determining selection pressure for protein-coding genes. Commonly, selective pressures on proteins are determined by comparing nucleotide sequences. Consequently, we compared the dN/dS ratios of AtLLP genes, as well as its orthologous pairs to examine the evolutionary divergences of the AtLLPs gene family over time (Additional File 11: Table S4, Table 1). As per evolutionary pressure analyses, AtLLPs are under negative selection, which indicates that natural selection suppresses protein modifications. Negative selection is effective for multiple housekeeping genes because mutations cause death or reduced fitness . From this analysis, nature retained the roles of LLPs as housekeeping, indicating that they have been important for plant fitness/survival and stress mitigation/signaling from very early days.
Since abiotic stresses alter transcriptional regulation and activate stress-responsive genes for stress mitigation, the analysis of gene expression in response to stresses is an important and frequently investigated aspect. Therefore, we initiated the expressional analysis of AtLLPs using publicly available expression data (Fig. 3). Differential expression pattern of six selected AtLLPs by in silico expression profiling across a wide range of abiotic stresses, namely, cold, osmotic, salt, drought, genotoxicicty, oxidative, UV-B, wounding, and heat strengthened the idea of the abiotic stress regulation of AtLLPs and their role in abiotic stress and related signaling pathways.
CREs are short DNA sequences that TFs use to regulate gene expression . Therefore, we performed structural analysis of AtLLP promoter sequences and CRE prediction for different TF binding sites (TFBS) (Fig. 4). Although computational identification and distribution of CREs do not confirm tissue-/time-specific regulatory mechanism, the presence of CREs provides an indication of the probable regulatory cues of gene expression. The promoter sequence analysis of AtLLPs identified multiple important consensus sequences under the functional categories of biotic and abiotic stress responsiveness, hormone signaling pathways, and development (Fig. 4). The presence of these abiotic stress-responsive CREs or TFBS in the promoter region of AtLLPs shows that the expression of AtLLPs may be influenced by multiple environmental signals such light, cold, and wounding.
Transcription is regulated by a number of proteins known as TFs, which form complexes involving protein–protein and protein–DNA interactions. Individual TFs bind to their cognate CREs or TFBS; therefore, we identified CREs cognate TFs in AtLLPs and analyzed the expressional correlation between tested AtLLPs and their CREs cognate TFs under cold, osmotic, salt, drought, genotoxic, oxidative, UV-B, wound, and heat stress (Figs. 5 and 6). We reported that that in each AtLLPs, the CREs cognate TFs differ but the expression of at least few of them was significantly correlated with their respective AtLLPs in each tested abiotic stresses in both root and shoot tissues (Figs. 5 and 6). This suggests that AtLLPs might be differentially regulated spatially and conditionally by the different sets of TFs under different abiotic stresses. Thus, as shown in Fig. 3, the results of expressional correlation analysis between tested AtLLPs and their CREs-specific TFs supported the differential expression patterns of AtLLPs under different abiotic stress conditions in both shoot and root tissues.
To experimentally verify the results of computational analysis, A. thaliana seedlings were subjected to abiotic stresses for 2, 6, and 12 h, and the expression profiling of seven AtLLPs was performed using qRT-PCR (Fig. 7).
The qRT-PCR result showed that the expression of AT1g53080 and AT3g16530 was inconsistent under the stress conditions that were tested (Fig. 7), suggesting they may not be an important player against abiotic stresses. The expression patterns of AT1g53060 were significantly downregulated in all tested stresses except MV (Fig. 7). Because AT1g53060 demonstrated initial upregulation in cold, ozone and UV-B stresses but significantly downregulated after 12 h except MV (Fig. 7), suggesting its role in abiotic early stress responses.
The AT1g07460 gene demonstrated maximal expression of 18.71- fold after 12 h in ozone and 3.83-fold after 6 h in wound stress (Fig. 7), indicating a possible role in these abiotic stresses. Further, another gene of this family, AT3g15356, showed a substantial increase in expression after 2 h of exposure to cold, high light, ozone, UV-B, and wound (Fig. 7), suggesting its role in early abiotic stress responses similar to AT1g53060.
However, the expression of other two genes, AT1g53070 and AT5g03350, was upregulated in most stresses tested and retained high transcript levels till 12 h of stress exposure (Fig. 7), suggesting their potential to be the important representative AtLLP genes under abiotic stresses (Fig. 7). The first key representative gene, AT1g53070, registered significant accumulations in its transcripts after exposures to MV and wound stress for 12 h (Fig. 7). Because these stresses are extensively known to increase reactive oxygen species (ROS) production in A. thaliana [49, 51, 61], the upregulation of its transcripts suggests that this gene may be regulated by ROS. The presence of an oxidative stress-induced CREs in the promoter region of this gene (Fig. 4) supports our results. The AT1g53070 expression was similarly up-regulated in response to ozone stress after 12 h (Fig. 7). The up-regulation of this gene can be supported by the fact that ozone induces a burst of ROS in apoplast and can induce extensive changes in gene expression . Similar to this gene, the other members of RLK family have been demonstrated with their increased gene expression under ozone stress [4, 63]. Temperature is one of the most important environmental factors that influence plant growth and development. Plants are affected by both high and low temperatures, which disturbs their normal cellular machinery. Unlike other abiotic stresses, AT1g53070 transcript level was maintained till 6 h but dropped after 12 h of cold stress (Fig. 7), indicating its thermo-sensitive nature.
The expression of the second key representative AtLLP, AT5g03350 demonstrated significant change in all of the assessed abiotic stresses (Fig. 7). The up-regulation of its transcripts was detected during cold stress from 6 to 12 h (Fig. 7), and this change in expression could be supported by recent study which showed that Lec-RLKs have a role in cold stress in antarctic mosses [64, 65]. AT5g03350 transcripts accumulated the most after 6 or 12 h of exposure to high Light, MV, and ozone (Fig. 7), all of which are well-known oxidative stressors in plants [50, 66, 67]. This expression pattern shows that it, like AT1g53070, is regulated by ROS. Other Lec-RLKs have been shown to have ROS-induced expressions when exposed to a ROS-inducing stress [12, 68]. Similarly, AT5g03350 expression showed maximal expression after 12 h in response to the wound (Fig. 7). Another Lec-RLK was demonstrated to be up-regulated in Populus nigra in response to wounding, which supports our results .
These expression results suggest that out of seven, AT1g53070 and AT5g03350 genes may be involved in Arabidopsis responses to multiple environmental stresses. To support this hypothesis, AT5g03350 was chosen as a key representative for further functional characterization using its overexpression and mutant lines (Figs. 8 and 9). The visualization of ROS generation in wild type, At5g03350 overexpression, and At5g03350 mutant (SALK_036814), plants subjected to the same six abiotic stresses as in the transcript study (Fig. 7) revealed that At5g03350 mutant had higher ROS than wild type while At5g03350 overexpression had lower ROS under high light, MV, UV-B, and wound (Fig. 9). This finding shows the essentiality of At5g03350 gene in combating these abiotic stresses. The findings were further corroborated by our qRT-PCR data, which showed that the transcript level of At5g03350 gene increased after 12 h of exposure to these stresses (Fig. 7).
Based on the above results, the possible roles of AtLLPs in response to abiotic stresses in A. thaliana were discussed and summarized in the form of a hypothetical model (Fig. 10).
Based on the results of this study, we propose that out of seven, two important representative AtLLPs, AT1g53070 and AT5g03350, have strong potential to play a role in abiotic stress mitigation and related signaling in A. thaliana. Additional investigation is required to establish the exact role of these genes in terms of ligands binding and signaling pathways in which they may be involved.
Availability of data and materials
Cucumber mosaic virus
General feature format
Legume lectin-like proteins
Plant growth regulator
- qRT PCR:
Quantitative-reverse transcription polymerase chain reaction
Salicylic acid induced-Legume Lectin like protein 1
Systemically acquired resistance
The Arabidopsis information resource
Translation start site
Shiu SH, Bleecker AB. Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 2003;132(2):530–43.
Wang Y, Nsibo DL, Juhar HM, Govers F, Bouwmeester K. Ectopic expression of Arabidopsis L-type lectin receptor kinase genes LecRK-I.9 and LecRK-IX.1 in Nicotiana benthamiana confers Phytophthora resistance. Plant Cell Rep. 2016;35:845–55.
Jose J, Ghantasala S, Roy Choudhury S. Arabidopsis transmembrane receptor-like kinases (RLKs): A bridge between extracellular signal and intracellular regulatory machinery. Int J Mol Sci. 2020;21(11):4000.
Mondal R, Biswas S, Srivastava A, Basu S, Trivedi M, Singh SK, Mishra Y. In silico analysis and expression profiling of S-domain receptor-like kinases (SD-RLKs) under different abiotic stresses in Arabidopsis thaliana. BMC Genom. 2021;22(1):1–5.
Wan J, Patel A, Mathieu M, Kim SY, Xu D, Stacey G. A lectin receptor-like kinase is required for pollen development in Arabidopsis. Plant Mol Biol. 2008;67:469–82.
Deng K, Wang Q, Zeng J, Guo X, Zhao X, Tang D, Liu X. A Lectin Receptor Kinase Positively Regulates ABA Response During Seed Germination and Is Involved in Salt and Osmotic Stress Response. J. Plant Biol. 2009;52:493–500.
Xin Z, Wang A, Yang G, Gao P, Zheng ZL. The Arabidopsis A4 subfamily of lectin receptor kinases negatively regulates abscisic acid response in seed germination. Plant Physiol. 2009;149:434–44.
Gilardoni PA, Hettenhausen C, Baldwin IT, Bonaventure G. Nicotiana attenuata LECTIN RECEPTOR KINASE1 suppresses the insect-mediated inhibition of induced defense responses during Manduca sexta Herbivory. Plant Cell. 2011;23:3512–32.
Cheng X, Wu Y, Guo J, Du B, Chen R, Zhu L, He G. A rice lectin receptor-like kinase that is involved in innate immune responses also contributes to seed germination. Plant J. 2013;76:687–98.
Singh P, Zimmerli, L. Lectin receptor kinases in plant innate immunity. Front Plant Sci. 2013;4:124.
Deb S, Sankaranarayanan S, Wewala G, Widdup E, Samuel MA. Te S-Domain Receptor Kinase Arabidopsis Receptor Kinase2 and the U Box/Armadillo Repeat-Containing E3 Ubiquitin Ligase9 Module Mediates Lateral Root Development under Phosphate Starvation in Arabidopsis. Plant Physiol. 2014;165:1647–56.
He XJ, Zhang ZG, Yan DQ, Zhang JS, Chen SY. A salt-responsive receptor-like kinase gene regulated by the ethylene signaling pathway encodes a plasma membrane serine/threonine kinase. Theor Appl Genet. 2004;109:377–83.
Joshi A, Dang HQ, Vaid N, Tuteja N. Pea lectin receptor-like kinase promotes high salinity stress tolerance in bacteria and expresses in response to stress in planta. Glycoconj J. 2010;27:133–50.
Eggermont L, Verstraeten B, Van Damme EJ. Genome wide screening for lectin motifs in Arabidopsis thaliana. Plant Genome. 2017;10(2). https://doi.org/10.3835/plantgenome2017.02.0010.
Barre A, Hervé C, Lescure B, Rougé P. Lectin receptor kinases in plants. Crit Rev Plant Sci. 2002;21(4):379–99.
Bouwmeester K, Govers F. Arabidopsis L-type lectin receptor kinases: phylogeny, classification, and expression profiles. J Exp Bot. 2009;60(15):4383–96.
Dievart A, Gottin C, Périn C, Ranwez V, Chantret N. Origin and diversity of plant receptor-like kinases. Ann Rev Plant Biol. 2020;71:31–156.
Bellande K, Bono JJ, Savelli B, Jamet E, Canut H. Plant lectins and lectin receptor-like kinases: how do they sense the outside? Int J Mol Sci. 2017;18(6):1164.
Kwon HK, Yokoyama R, Nishitani K. A proteomic approach to apoplastic proteins involved in cell wall regeneration in protoplasts of Arabidopsis suspension-cultured cells. Plant Cell Physiol. 2005;46(6):843–57.
Krinke O, Ruelland E, Valentová O, Vergnolle C, Renou JP, Taconnat L, Zachowski A. Phosphatidylinositol 4-kinase activation is an early response to salicylic acid in Arabidopsis suspension cells. Plant Physiol. 2007;144(3):1347–59.
Armijo G, Salinas P, Monteoliva MI, Seguel A, García C, Villarroel-Candia E, Holuigue L. A salicylic acid–induced lectin-like protein plays a positive role in the effector-triggered immunity response of Arabidopsis thaliana to Pseudomonas syringae Avr-Rpm1. Mol Plant Microbe Interact. 2013;26(12):1395–406.
Ishihara T, Sakurai N, Sekine KT, Hase S, Ikegami M, Shibata D, Takahashi H. Comparative analysis of expressed sequence tags in resistant and susceptible ecotypes of Arabidopsis thaliana infected with Cucumber mosaic virus. Plant Cell Physiol. 2004;45(4):470–80.
Hoth S, Ikeda Y, Morgante M, Wang X, Zuo J, Hanafey MK, Chua NH. Monitoring genome-wide changes in gene expression in response to endogenous cytokinin reveals targets in Arabidopsis thaliana. Febs Lett. 2003;554(3):373–80.
Zhao Y, Dai X, Blackwell HE, Schreiber SL, Chory J. SIR1, an upstream component in auxin signaling identified by chemical genetics. Science. 2003;301(5636):1107–10.
Mukherjee AK, Carp MJ, Zuchman R, Ziv T, Horwitz BA, Gepstein S. Proteomics of the response of Arabidopsis thaliana to infection with Alternaria brassicicola. J proteomics. 2010;73(4):709–20.
Zhao J, Wang J, An L, Doerge RW, Chen ZJ, Grau CR, Osborn TC. Analysis of gene expression profiles in response to Sclerotinia sclerotiorum in Brassica napus. Planta. 2007;227(1):13–24.
Broekgaarden C, Poelman EH, Steenhuis G, Voorrips RE, Dicke M, Vosman B. Genotypic variation in genome-wide transcription profiles induced by insect feeding: Brassica oleracea-Pieris rapae interactions. BMC Genom. 2007;8(1):1–13.
Hori C, Yu X, Mortimer JC, Sano R, Matsumoto T, Kikuchi J, Ohtani M. Impact of abiotic stress on the regulation of cell wall biosynthesis in Populus trichocarpa. Plant Biotechnol. 2020;37(3):273–83.
Rui Y, Dinneny JR. A wall with integrity: surveillance and maintenance of the plant cell wall under stress. New Phytol. 2020;225(4):1428–39.
Abercrombie JM, Halfhill MD, Ranjan P, Rao MR, Saxton AM, Yuan JS, Stewart CN. Transcriptional responses of Arabidopsis thaliana plants to As (V) stress. BMC Plant Biol. 2008;8(1):1–15.
Chen C, ChenH, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.
Van Bel M, Silvestri F, Weitz E M, Kreft L, Botzki A, Coppens F, Vandepoele K. PLAZA 5.0: extending the scope and power of comparative and functional genomics in plants. Nucleic Acids Res. 2022;50(D1):D1468–D1474.
Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N, Rokhsar DS. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012;40(D1):1178–86.
Gasteiger E, Hoogland C, Gattiker A, Wilkins MR, Appel RD, Bairoch A. Protein identification and analysis tools on the ExPASy server. In: Walker JM, editor. The proteomics protocols handbook. Totowa: Humana Press; 2005. p. 571–607.
Hooper CM, Castleden IR, Tanz SK, Aryamanesh N, Millar AH. SUBA4: the interactive data analysis centre for Arabidopsis subcellular protein locations. Nucleic Acids Res. 2017;45(D1):1064–74.
Kumar S, Stecher G, Li M, Knyaz C,Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35(6):1547–9.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
Kumar S, Stecher G, Suleski M, Hedges SB. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol. 2017;34(7):1812–9.
Toufighi K, Brady SM, Austin R, Ly E,Provart NJ.The Botany Array Resource: e-Northerns, expression angling, and promoter analyses. Plant J. 2005;43(1):153–63.
Jin J, Tian F, Yang DC, Meng YQ, Kong L, Luo J, Gao G. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2016;45(D1):D1040–5.
Curtis MD, Grossniklaus U. A gateway cloning vector set for high-throughput functional analysis of genes in planta. Plant Physiol. 2003;133(2):462–69.
Zhang X, Henriques R, Lin SS, Niu QW, Chua NH. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nat Protoc. 2006;1(2):641–6.
Harrison SJ, Mott EK, Parsley K, Aspinall S, Gray JC, Cottage A. A rapid and robust method of identifying transformed Arabidopsis thaliana seedlings following floral dip transformation. Plant methods. 2006;2(1):1–7.
Murray MG, Thompson W. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–6.
Kihara T, Zhao CR, Kobayashi Y, Takita E, Kawazu T, Koyama H. Simple identification of transgenic Arabidopsis plants carrying a single copy of the integrated gene. Bioscience, biotechnology, and biochemistry. 2006;70(7):1780–3.
Sharma Y K, León J, Raskin I, Davis KR. Ozone-induced responses in Arabidopsis thaliana: the role of salicylic acid in the accumulation of defense-related transcripts and induced resistance. Proc Natl Acad Sci. 1996;93(10):5099–104.
Chassot C, Buchala A, Schoonbeek HJ, Métraux JP, Lamotte O. Wounding of Arabidopsis leaves causes a powerful but transient protection against Botrytis infection. Plant J. 2008;55(4):555–67.
Zeller G, Henz SR, Widmer CK, Sachsenberg T, Rätsch G, Weigel D, Laubinger S. Stress-induced changes in the Arabidopsis thaliana transcriptome analyzed using whole‐genome tiling arrays. Plant J. 2009;58(6):1068–82.
Galvez-Valdivieso G, Fryer MJ, Lawson T, Slattery K, Truman W, Smirnoff N, Asami T, Davies WJ, Jones AM, Baker NR, Mullineaux PM. The high light response in Arabidopsis involves ABA signaling between vascular and bundle sheath cells. Plant Cell. 2009;21(7):2143–62.
Benina M, Ribeiro DM, Gechev TS, Mueller-Roeber B, Schippers JH. A cell type-specific view on the translation of mRNAs from ROS-responsive genes upon paraquat treatment of Arabidopsis thaliana leaves. Plant Cell Environ. 2015;38(2):349–63.
Kornbrot D. Statistical software for microcomputers: SigmaPlot 2000 and SigmaStat2. Br J Math Stat Psychol. 2000;53:335–7.
Kryazhimskiy S, Plotkin JB. The population genetics of dN/dS. PLoS Genet. 2008;4.
Jamieson PA, Shan L, He P. Plant cell surface molecular cypher: Receptor-like proteins and their roles in immunity and development. Plant Sci. 2018;274:242–51.
Panchy N, Lehti-Shiu M, Shiu SH. Evolution of gene duplication in plants. Plant Physiol. 2016;171:2294–316.
Wang HF, Feng L, Niu DK. Relationship between mRNA stability and intron presence. Biochem Biophysic Res Commun. 2007;354(1):203–8.
Chothia C, Lesk AM. The relation between the divergence of sequence and structure in proteins. EMBO J. 1986;5(4):823–6.
Povolotskaya IS, Kondrashov FA. Sequence space and the ongoing expansion of the protein universe. Nature. 2010;465(7300):922–6.
Loewe L. Negative selection. Nat. Educ. 2008;1(1):59.
Heidari P, Ahmadizadeh M, Najafi-Zarrini H. In Silico Analysis of Cis-Regulatory Elements on Co-Expressed Genes. J Biol Environ Sci. 2015;9(25):1–9.
Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24.
Overmyer K, Brosché M, Pellinen R, Kuittinen T, Tuominen H, Ahlfors R, Keinänen M, Saarma M, Scheel D, Kangasjärvi J. Ozone-induced programmed cell death in the Arabidopsis radical-induced cell death1 mutant. Plant Physiol. 2005;137(3):1092–104.
Wrzaczek M, Brosché M, Salojärvi J, Kangasjärvi S, Idänheimo N, Mersmann S, Robatzek S, Karpiński S, Karpińska B, Kangasjärvi J. Transcriptional regulation of the CRK/DUF26 group of receptor-like protein kinases by ozone and plant hormones in Arabidopsis. BMC Plant Biol. 2010;10(1):1–19.
Liu S, Wang J, Chen K, Zhang Z, Zhang P. The L-type lectin receptor-like kinase (PnLecRLK1) from the Antarctic moss Pohlia nutans enhances chilling-stress tolerance and abscisic acid sensitivity in Arabidopsis. Plant Growth Regul. 2017;81(3):409–18.
Minozzo MM, Metz GF, de Vargas MVM, Pereira AB, de Carvalho Victoria F, Phylogenetic and selection pressure analyses of cold stress-associated PAL-Like and Lec-RLK genes in antarctic mosses. Curr Plant Biol. 2020;24:100178.
Wei Z, Gao T, Liang B, Zhao Q, Ma F, Li C. Effects of Exogenous Melatonin on Methyl Viologen-Mediated Oxidative Stress in Apple Leaf. Int J Mol Sci. 2018;19(1):316.
Sharma YK, León J, Raskin I, Davis KR. Ozone-induced responses in Arabidopsis thaliana: the role of salicylic acid in the accumulation of defense-related transcripts and induced resistance. Proc Natl Acad Sci. 1996;93(10):5099–104.
Vaid N, Pandey P, Srivastava VK, Tuteja N. Pea lectin receptor-like kinase functions in salinity adaptation without yield penalty, by alleviating osmotic and ionic stresses and upregulating stress-responsive genes. Plant Mol Biol. 2015;88:193–206.
Nishiguchi M, Yoshida K, Sumizono T, Tazaki K. A receptor-like protein kinase with a lectin-like domain from lombardy poplar: gene expression in response to wounding and characterization of phosphorylation activity. Mol Genet Genom. 2002;267:506–14.
Mr. Subhankar Biswas is thankful to Council of Scientific & Industrial Research (CSIR), New Delhi, India, for the Junior and Senior Research Fellowship (JRF & SRF); Ms. Akanksha Srivastava is thankful to the Department of Science and Technology-Innovation of Science Pursuit for Inspire Research (DST–INSPIRE), New Delhi, India, for the fellowship.
Dr. Yogesh Mishra is thankful to DST-SERB for project grants (ECR/2016/000186). Moreover, Dr. Mishra is thankful to Institute of eminence (IOE) incentive grant, Banaras Hindu University (R/Dev/D/IOE/Incentive/2021-22/32401). We would like to thank the Head and the Programme Coordinator (CAS) in Botany, Banaras Hindu University, Varanasi, India, for providing instrumentation-related facilities.
Ethics approval and consent to participate
The seeds of A. thaliana ecotype Columbia-0 (Col-0) were obtained from the Arabidopsis Biological Resource Centre (ABRC). In this study, the authors report that experimental research works on the plant described comply with institutional, national and international guidelines.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fig. S1. A denaturing formaldehyde agarose gel electrophoresis image for RNA.
List of primers used in the current study.
Fig. S2. Represents the steps of overexpression line (At5g03350::overexpression) creation using gateway cloning technology followed by Agrobacterium mediated transformation.
Sequencing information of AT5g03350CDS-pDONR207 entry clone (forward).
Sequencing information of AT5g03350CDS-pDONR207 entry clone (reverse).
Fig. S3. Represents the steps for PCR-based genotyping for At5g03350:: overexpression line.
Fig. S4. Steps of PCR-based genotypic analysis of AT5g03350 gene-specific T-DNA insertion line.
Fig. S5. Peptide sequence alignment of AtLLPs showing the presence of highly conserved lectin-binding domain and less conserved signal peptide.
Physicochemical properties of AtLLPs.
Predicted and experimentally confirmed subcellular location of AtLLPs.
Estimation of divergence time of Arabidopsis thaliana LLP gene family and its orthologs in Arabidopsis lyrata, Brassica rapa, Solanum lycopersicum, Zea mays, Selaginella moellendorffii, Physcomitrium patens, and Chlamydomonas reinhardtii.
List of transcription factors (TFs), which are predicted to bind cis-regulatory elements (CREs) present in 500 bp upstream promoter region of AtLLPs.
Table S6. Expressional correlation of coefficients (r-value) analysis between AtLLPs and their cis-regulatory elements (CREs) cognate transcription factors (TFs) in root and shoot tissues under different abiotic stress conditions.
Number of predicted transcription factors (TFs), which can bind cognate cis-regulatory elements (CREs) in the promoter region of AtLLPs, as well as the list of TFs not included in this study and the number of TFs that are significantly correlated under different stress conditions in both root and shoot tissues.
Raw data of qRT-PCR.
About this article
Cite this article
Biswas, S., Mondal, R., Srivastava, A. et al. In silico characterization, molecular phylogeny, and expression profiling of genes encoding legume lectin-like proteins under various abiotic stresses in Arabidopsis thaliana. BMC Genomics 23, 480 (2022). https://doi.org/10.1186/s12864-022-08708-0
- Abiotic stress
- Arabidopsis thaliana
- Legume lectin-like proteins
- Receptor-like kinases