1 Introduction
Programmed cell death (PCD) is an active process that results in selective elimination of damaged or unwanted cells [1]. It is also essential for growth and development of multicellular organisms as well as for proper response to environment [2,3]. Although PCD is extensively studied in animals, very little is known about genetic mechanisms underlying PCD in plant. Cell death machinery in plants differs genetically from that in animals [1]. This conclusion was supported by the detection of unique plant-specific proteases and nucleases that act during PCD. However, PCD in plant and animal share common morphological and biochemical characteristics, including cytoplasm shrinkage, chromatin aggregation, and nuclear DNA cleavage into shorter fragments [1,4].
Plant PCD occurs in a number of developmental processes, including embryo formation, differentiation of treachery elements in xylem, modelling of leaf shape, reinforcement of the cell wall, and leaf senescence [2,5,6]. Several plant adaptation processes, such as the hypersensitive response (HR) to pathogens, require PCD; however, many unfavorable abiotic stresses, such as heat or cold and salt [7–9] as well as necrotrophic pathogens [1], induce unwanted PCD. These abiotic stresses can result in oxidative stress and the consequent reactive oxygen species- or ROS-dependent PCD, in which ROS act as destructive toxic agents and as signals unlocking the PCD program [10]. In other words, PCD can serve as a positive or negative aspect of plant adaptation to the environment [1]. These processes reflect the unique aspects of plant architecture and physiology.
Signs of PCD accompany salt stress in several plant species [8,11]. In Arabidopsis, these signs were shown in seedling primary roots of both wild-type (WT) and plasma membrane Na/H antiporter salt overly sensitive (sos1) mutant. The WT could stand salt stress due to the formation of secondary roots, while the mutant died [11]. Salt stress also induces PCD in rice root tip cells and tobacco protoplasts by disturbing equilibrium of Caþ, Kþ, and Hþ ion and by overproduction of hydrogen peroxide [11,12]. A high level of hydrogen peroxide acts in conjunction with the open state of the mitochondrial permeability pores in inducing analogies between cell death program in plants and apoptosis in animals [8].
Virus-induced gene silencing (VIGS) is a method to induce gene knockdown that takes advantage of the plant RNAi-mediated antiviral defence mechanism [13]. The process is adapted for high-throughput functional genomics in plant via the use of the Agrobacterium tumefaciens Ti plasmid, a plant pathogen, to recover a recombinant virus (e.g., tobacco rattle virus or TRV) carrying part of the gene sequence targeted for silencing via the plant endogenous RNAi machinery [14–16]. The aim of the present study is to detect the influence of PCD-related genes in conferring tolerance to salt stress in Arabidopsis.
2 Materials and methods
2.1 Construction of virus-induced gene silencing (VIGS) lines
Wild-type tobacco plants were grown from seed in a growth room at 25 °C with 16-h light periods, as previously described [17]. A. tumefaciens and Escherichia coli strains were grown in a Luria–Bertani medium at 30 and 37 °C, respectively. Antibiotics were used at the following concentrations: ampicillin (100 μg/mL), kanamycin (50 μg/mL), and rifampicin (100 μg/mL). VIGS of selected PCD-related genes were done in 4-week-old tobacco (N. benthamiana) seedlings as described by Velásquez et al. [13]. Primers used in constructing pTRV2 vectors of different gene are shown in Table S1. A. tumefaciens harboring the gateway compatible pTRV2 [16] and pTRV2::PDS [18] were used as negative and positive controls, respectively. Silencing of phytoene desaturase (or PDS) causes the plants to photo-bleach, hence, it can be used as a reporter gene to detect silencing efficiency.
2.2 PCD assays
Leaf discs were obtained from 7-week-old plants of WT or VIGS lines using a 10-mm diameter cork borer, then, submerged in OA (Sigma–Aldrich) at 20 mM for 0, 12, 24 and 48 h for the initial experiment to detect the time of PCD occurrence in the WT, and for 0 and 24 h for the VIGS lines experiment. Cell death was determined in each experiment and quantified by using Evans blue as previously described [19]. Measurements were expressed as relative values, with “one” corresponding to a maximum of a sample tested. DNA laddering was detected across time for the WT experiment, where total plant DNAs were extracted from leaf discs across time using a modified method [20]. DNA samples were separated on a 1.5% agarose gel in 0.5 × Tris–borate–EDTA, stained with ethidium bromide and visualized using the Gel Doc XR from Bio-Rad.
2.3 RNA-Seq in tobacco WT plants
Leaf discs of 7-week-old wild-type tobacco (N. benthamiana) plants were submerged in OA at 20 mM for 2, 6, 12 and 24 h. Then, discs were frozen in liquid nitrogen to extract total RNAs from similar-sized leaf samples using Trizol (Invitrogen, Life Tech, Grand Island, NY, USA) and treated with RNase-free DNase (Promega Corporation, Madison, WI, USA) in the presence of 1 U/μL of RNasin® Plus RNase Inhibitor (Promega) for 2 h at 37 °C. Untreated leaf discs were used as a negative control. RNAs were quantified, and 30 μg (300 ng/μL) were used for RNA-Seq. Purified RNA samples were shipped to Genomics Unit, MSU, USA in three replicates for deep sequencing and generation of datasets (at least 100 million reads per sample). Filtered reads were aligned with up to two mismatches to the tobacco (N. benthamiana) draft genome sequences (v0.4.2 Scaffolds + NrContigs) of the Sol Genomics Network: (ftp.solgenomics.net/genomes/Nicotiana_benthamiana/annotation/Niben.genome.v0.4.4.transcripts.annotated.fasta). RSEM v1.1.6, an RNA-Seq quantification tool, was used to estimate the relative abundances and expected read counts for the transcripts. RSEM uses the Bowtie aligner (Bowtie v0.12.1) to map the reads against the transcripts. Transcript quantification of the reference-aligned reads was performed with RSEM, which allowed for the assessment of transcript abundances based on the mapping of RNA-Seq reads to the assembled transcriptome. Expected read counts were used as input to differential expression analysis by EdgeR (version 3.0.0, R version 2.1.5). The median was used as the common dispersion factor for differential expression (DE) analyses. The remaining unmapped sequences were re-aligned against the contigs collectively de novo assembled using the Trinity RNA-Seq Assembly package (r2013-02-25) from total unmapped sequences of all treatments and replicates [21]. Blastx was performed against the NCBI non-redundant protein database with an E-value cut off of 1e−5. Fold change of DE transcripts was calculated against the published tobacco actin gene that exhibited relatively constant expression levels. Accordingly, FC of ≥ 5 was selected to detect DE transcripts that might significantly have a role during PCD. To identify clusters with functional enrichment, we determined a significant Pearson correlation through permutation analysis [22]. The resulting clusters were refined by visual inspection and analyzed for GO term enrichment using Blast2GO (http://www.blast2go.org/). We also clustered the RPKM data to provide a representation of absolute abundance of the transcripts.
2.4 sqRT-PCR
This procedure was done to validate RNA-Seq data analysis and to validate the occurrence of VIGS in tobacco. In the first experiment, eight differentially expressed (DE) transcripts representing selected clusters were validated. In the second experiment, VIGS lines of the studied PCD-related genes were validated with the WT and VIGS of empty pTRV2 used as controls. Leaf discs of 7-week-old N. benthamiana plants of the different genotypes were submerged in OA at 20 mM and samples were taken after 0, 2, 6, 12 and 24 h for the first validation, while 0 and 2 h for the second validation. The latter indicates the time point where the selected genes were highly expressed. Then, discs were frozen in liquid nitrogen to extract RNAs. Total RNA was extracted from three similar-sized (10 mm2) leaf discs per plant (approximately 50 mg tissue) collected from upper leaves using Trizol (Invitrogen) and treated with RNase-free DNase (Promega Inc.). First-strand cDNA was synthesized using 2.5 μg of total RNA, 0.5 μg oligo (dT) primer and Superscript II reverse transcriptase (Invitrogen) to a final volume of 20 μL. sqRT-PCR was performed in 20-μL reactions using 1 μL cDNA, 1 × PCR buffer (with 1.5 mM MgCl2), 200 μM dNTPs, 200 nM of each gene-specific primers (Table S1) and 0.2 U of Taq DNA polymerase (Promega Inc.). Forty PCR cycles for each gene product include denaturation at 94 °C for 15 s, annealing at 55 °C for 30 s, and extension at 72 °C for 45 s. Amplicons were analyzed on a 1.2% agarose gel stained with ethidium bromide and visualized using the Gel Doc XR from Bio-Rad Laboratories (Hercules, CA, USA). The actin gene (372 bp) was used in either experiment as a housekeeping control. Primers were designed using Netprimer software (http://www.premierbiosoft.com/netprimer/index.html) with the following criteria: length ∼ 20 bases, GC content ∼ 50%, minimal secondary structures, comparable annealing temperatures (55 °C) of the primer pairs, and PCR products of 275–375 bp.
For detecting the TRV virus in the newly emerged leaves [23], a first-strand from TRV RNA in leaf RNA extracts of VIGS lines was synthesized using TRV1 primer (5′-CAGTCTATACACAGAAACAGA-3′). Then, PCR amplification was performed using primers of TRV1 and TRV2 (5′-GACGTGTGTACTCAAGGGTT-3′). VIGS of empty TRV2 was used as a positive control. Initial denaturation (2 min at 94 °C) was first done, followed by 40 cycles (denaturation for 15 s at 94 °C, annealing for 2 min at 58 °C, extension for 1 min at 72 °C; and a final extension step for 10 min at 72 °C).
2.5 Salt stress experiment for Arabidopsis T-DNA insertion SALK lines
An experiment was conducted to detect the performances of five T-DNA insertion KO SALK lines at 0, 100, 150 or 200 mM NaCl as compared to the WT plants (Col). Mutant seeds (Col) were provided by the SALK Institute, Genomic Analysis Laboratory (SIGnAL) (http://signal.salk.edu/tdnaprimers.2.html). The lines were screened for homozygous segregants by PCR using different PCD-related gene-specific primers and left border (LB) T-DNA primers. To score the germination rate at day 6, 30 seeds/replicate/treatment/line were surface sterilized, sown in Petri dishes containing MS medium [24] with different concentrations of NaCl (0, 100, 150 or 200 mM). Experiments were conducted in three replicates and analyzed using a randomized complete block design. The plates were kept in the dark at 4 °C for 2 days and then shifted to 21 ± 2 °C (day/night) under a 16-h light/8-h dark cycle. Emergence of the radicle from seed coat was considered as completion of germination. To score other characteristics under salt stress, two-week-old plants of different mutants, grown at 0 mM NaCl under a 16-h light/8-h dark cycle at 21 ± 2 °C (day/night), were transferred to the soil (4 plants/pot/replicate) with different concentrations of NaCl (0, 100, 150 or 200 mM) to grow for 2 more weeks. Light intensity was ∼ 175 μmoles/m2/s1 for the 16-h photoperiod. Then, the individual plant root length, the number of leaves per plant, and the rosette area were scored for 4-week-old plants of different SALK lines as well as the control WT plants. The experiment was conducted using a randomized complete block design with three replicates. Statistical analyses of different experiments were performed following the procedure outlined by Gomez and Gomez [25], and multiple comparisons were performed based on the Duncan's New Multiple Range test [26].
3 Results and discussion
3.1 Detection of PCD as triggered by OA treatment
Qualitative and quantitative time-course analysis of cell death upon treatment with OA in tobacco WT leaf discs were measured using Evans blue staining. Gradual cell death was clearly observed across time of OA treatment (Fig. S1). Our results indicated that treatment beyond 24 h is not necessary. Longer treatment time (48 h) resulted in recovering poor quality RNAs for further gene expression analysis. Kim et al. [19] indicated substantial increases in tobacco cell death 48 through 96 h after OA treatment with concentration (e.g. 20 mM) similar to this in the present study. However, they indicated that a 96-h treatment resulted in extensive tissue collapse. This reflects the difficulty of isolating good quality RNA for gene expression analysis. They also indicated that lower concentration of OA (10 mM) yielded similar but delayed results.
Based on the kinetics of the cell death in response to OA, tobacco gDNAs were isolated at several time points. DNA cleavage resulting in a characteristic DNA laddering pattern associated with apoptotic cell death started to appear 12 h after OA treatment (Fig. S2). OA was reported by Kim et al. [19] as a non-specific multifunctional phytotoxin in plant tissue that enhances fungal pathogenicity and induces PCD by inhibiting the oxidative burst of the plant.
3.2 Analysis of RNA-Seq dataset
This analysis was conducted in order to detect new PCD-related genes induced by OA treatment in tobacco (N. benthamiana). Involvement of these genes in conferring salt stress tolerance in Arabidopsis was done. These genes will be candidates for future improvement of plant salt tolerance via genetic transgenesis. Sequencing of cDNAs recovered from RNA samples of tobacco collected during the course of PCD yielded between 51 and 55 million reads corresponding to over 5 billion nucleotides per sample (Table 1). The raw sequencing reads will eventually be deposited in the Small Read Archive (SRA) at GenBank. Between ∼ 60 and 65% of tobacco reads were mapped to the draft genome sequences in exonic regions (v0.4.2 Scaffolds + NrContigs) of the Sol Genomics Network (http://solgenomics.net/), while the remaining ∼ 35–40% matched no sequences (Table 1). The unmapped reads of all treatments and replicates were collectively de novo assembled and the resulting contigs were used as a new reference genome. Then, the unmapped reads were aligned to this new reference genome and the resulted transcripts were added to those mapped on the draft genome for subsequent analysis.
Statistics of tobacco (N. benthamiana) RNA-Seq numerical data analysis.
Treatment time | Total no. readsa | Aligned readsb | % readsc |
Unaligned readsd | % readse |
Control | 55200292 | 35439431 | 64.20 | 19760861 | 35.80 |
2 h | 55133896 | 35399120 | 64.21 | 19734776 | 35.79 |
6 h | 54501770 | 34394296 | 63.11 | 20107474 | 36.89 |
12 h | 51761784 | 33333310 | 64.40 | 18428474 | 35.60 |
24 | 51688534 | 30835697 | 59.66 | 20852837 | 40.34 |
a Total number of reads recovered from RNA-Seq of OA-untreated and treated tobacco discs.
b Number of reads aligned with tobacco (N. benthamiana) draft genome sequences (v0.4.2 Scaffolds + NrContigs) of the Sol Genomics Network (http://solgenomics.net/).
c Percentage of tobacco reads aligned with tobacco (N. benthamiana) draft genome.
d Number of reads unaligned with tobacco (N. benthamiana) draft genome.
e Percentage of tobacco reads unaligned with tobacco (N. benthamiana) draft genome.
RPKM-derived read counts were compared using a likelihood ratio test [21]. Statistical analysis was reliable when applied to genes with an RPKM value ≥ 2. To determine DE transcripts, a two-fold (or greater) change in expression and false discovery rate (FDR) of 10−3 or less was required. The resulting number of DE transcripts from the alignment with tobacco (N. benthamiana) reference genome as well as the de novo assembled transcripts was ∼ 2650 (data provided upon request). DE transcripts were, then, utilized in generating the expression clusters across time of OA treatment. Expression profiles of the DE transcripts were determined by a cluster analysis based on the k-means method using Pearson's correlation distance so that the similarity in relative change for each transcript or among transcripts across time of OA treatment was determined. These data were, then, subjected to hierarchical clustering using the Pearson correlation as the distance metric (Fig. 1). Two-dimensional hierarchical clustering classified DE profiles into 43 expression clusters according to the similarity of their expression profiles in which 16 of them (Table S2) fitted within the required expression profiles of upregulation after 2 h. The selected clusters might include genes with major effects on triggering or retarding PCD. A number of 23 PCD-related genes were selected for further analysis (Fig. 2) in which eight of them successfully validated the RNA-Seq data via sqRT-PCR (Fig. 3) across different time points of OA treatment. Data across the different selected clusters indicated the regulation of a number of 23 PCD-related genes (Table 2 and Fig. 2). These genes encode Bax inhibitor-1 proteins, Mlo, metacaspase 1, programmed cell death protein 2-like, apoptosis-inducing factors, plasminogen activator inhibitor, Stm1, xylosyltransferase, and exostosin.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr1.jpg)
Hierarchical cluster analysis (heat map) based on log ratio RPKM data for transcripts resulted from tobacco leaf discs exposed to OA across time (0, 2, 6, 12 and 24 h) to trigger PCD-related genes.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr2.jpg)
Expression pattern of the selected upregulated PCD-related genes (see Table 2) 2 h after exposure to 20 mM oxalic acid (OA) in tobacco. The proteins encoded by these genes harbor PCD-related domains.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr3.jpg)
Semi-quantitative RT-PCR and profiles of fold change values resulted from RNA-Seq analysis of PCD-related genes in selected clusters (b–d) used for validating RNA-Seq data of tobacco leaves at different time points of OA (20 mM) treatment (0, 2, 6, 12, and 24 h). Gene codes (Gs) refer to PCD-related genes described in Table 2. The “actin” gene (a) was used as the unregulated housekeeping gene.
List of PCD-related genes or genes with PCD domains detected from RNA-Seq analysis and selected for further analysis.
Gene code | Transcript ID | Description | PCD-related domain |
G1 | GG79261|c1_g1_i1 | u3 small nucleolar RNA-associated protein 4-like | WD, G-beta repeat |
G2 | GG63750|c4_g2_i2 | Katanin p80 wd40 repeat-containing subunit b1 homolog | WD, G-beta repeat |
G3 | GG79836|c0_g1_i1 | Metacaspase-3-like isoform × 2 | Caspase |
G4 | GG3719|c1_g1_i2 | BI-1-like partial | Inhibitor of apoptosis-promoting Bax1 |
G5 | GG59017|c0_g1_i3 | Late blight resistance protein homolog r1a-10-like | NB-ARC |
G6 | GG64878|c6_g1_i2 | Ras-related protein raba2a-like | Miro-like protein |
G7 | GG5937|c1_g1_i1 | Metacaspase 1 isoform 1 | Caspase |
G8 | GG40942|c0_g1_i3 | Probable poly(ADP-ribose) polymerase | RCD1-SRO-TAF4 (RST) |
G9 | GG59097|c1_g1_i3 | Programmed cell death protein 2-like | Programmed cell death protein 2 |
G10 | GG38146|c5_g1_i5 | mlo-like protein 11-like (promoting cell death) | Mlo family |
G11 | GG70956|c4_g1_i3 | sec12-like protein 1-like | WD, G-beta repeat |
G12 | GG28140|c1_g1_i1 | Uncharacterized protein LOC102589284 | WD, G-beta repeat |
G13 | GG3719|c1_g1_i4 | Bax inhibitor-1 family protein | Inhibitor of apoptosis-promoting Bax1 |
G14 | GG84247|c3_g1_i2 | gtp-binding protein yptm2-like | Miro-like protein |
G15 | GG4666|c0_g1_i7 | Apoptosis-inducing factor (AIF) homolog a-like | Apoptosis-inducing factor |
G16 | GG85010|c2_g1_i1 | wd repeat-containing protein rup2-like | WD, G-beta repeat |
G17 | GG61596|c2_g1_i1 | Late blight resistance protein homolog r1b-16-like | NB-ARC |
G18 | GG68778|c0_g1_i6 | Mlo family | mlo-like protein 3-like |
G19 | GG19032|c1_g1_i2 | MACPF domain-containing protein nsl1-like | MAC/Perforin |
G20 | GG43397|c1_g1_i2 | Programmed cell death protein 2-like | Programmed cell death protein 2-like |
G21 | GG56912|c0_g2_i3 | Plasminogen activator inhibitor-1 RNA-binding | Stm1 |
G22 | GG6585|c1_g1_i2 | tmv resistance protein n-like | NB-ARC |
G23 | GG58302|c1_g1_i3 | Exostosin family | Beta-xylosyltransferase irx10l-like |
3.3 Performance of tobacco VIGS lines for selected PCD-related genes under OA treatment
Effectiveness of VIGS in tobacco was proven via the use of pTRV2 plasmid carrying a fragment of the PDS gene, e.g., TRV2::PDS. Silencing of the PDS gene has resulted in photo-bleaching of the new leaves 10 days after infiltration (Fig. 4). TRV virus was successfully traced in the 23 VIGS lines via sqRT-PCR with primers specific to the 16-kDa gene (463 bp) of TRV (Fig. 5). sqRT-PCR was also used to confirm occurrence of knockdown in different genes via VIGS (Fig. 6). The results of PCD assay, in terms of Evans blue staining after 24 h, in VIGS of different PCD-related genes are shown in Table 3. Two controls were used in the study, e.g., WT and VIGS with empty pTRV2. The results indicated that VIGS of five PCD-related genes (lines VIGS-G8, VIGS-G13, VIGS-G15, VIGS-G18 and VIGS-G21) responded differently from WT and VIGS with empty pTRV2. One of them, namely Bax Inhibitor-1 (BI-1), has an antiapoptotic effect where the mean relative cell death value due to the exposure OA in its knockdown VIGS lines was significantly higher than the two controls. The other four genes have an apoptotic effect as the mean relative cell death values under treatment in their knockdown VIGS lines were significantly lower than for the two controls. These four genes encode Mlo-like protein, apoptosis-inducing factor (AIF) homolog a, Stm1 and RCD1-SRO-TAF4 (RST). Analogs of these five genes in Arabidopsis were further analyzed for their role during salt stress.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr4.jpg)
Silencing of the PDS control gene that causes photo-bleaching in tobacco plants. The gene was used as a model gene to confirm the occurrence of gene silencing due to VIGS.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr5.jpg)
Semi-quantitative RT-PCR to detect TRV virus in the different 23 VIGS lines with PCD-related knocked down genes as well as the wild-type (WT) and TRV2-transformed (V2) plants. Gene codes refer to those indicated in Table 2.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr6.jpg)
Semi-quantitative RT-PCR for VIGS lines with PCD-related knocked down genes triggered by oxalic acid (OA) treatment for 2 h as compared to the tobacco wild-type (WT) and VIGS line with empty pTRV2 (V2) plants. Amplicon sizes of different genes and primers used are shown in Table S1. The actin gene of different VIGS was used as a housekeeping control. Gene codes refer to those indicated in Table 2.
Multiple comparisons of the mean relative cell death as responses of tobacco leaf discs of the 23 VIGS lines following OA (20 mM, pH 7.0) treatment for 24 h as determined by Evans blue staining. Dye released from dead cells was measured at absorbance at 600 nm. Measurements were expressed as relative values with 1 corresponding to the maximum of a sample. Data are presented as means from two independent experiments with three replicates each.
Gene code | Genotype | OA treatment | |
0 mM | 20 mM | ||
WT | 0.18DE | 0.60B | |
VIGS line with empty TRV2 | 0.22D | 0.61B | |
G1 | VIGS-G1 line | 0.20DE | 0.72B |
G2 | VIGS-G2 line | 0.27D | 0.62B |
G3 | VIGS-G 3 line | 0.19DE | 0.64B |
G4 | VIGS-G 4 line | 0.22D | 0.59B |
G5 | VIGS-G 5 line | 0.26D | 0.65B |
G6 | VIGS-G 6 line | 0.27D | 0.60B |
G7 | VIGS-G 7 line | 0.20DE | 0.68AB |
G8a | VIGS-G 8 line | 0.23D | 0.38C |
G9 | VIGS-G 9 line | 0.28D | 0.58B |
G10 | VIGS-G 10 line | 0.29D | 0.61B |
G11 | VIGS-G 11 line | 0.24D | 0.68AB |
G12 | VIGS-G 12 line | 0.17E | 0.62B |
G13b | VIGS-G 13 line | 0.23D | 1.00A |
G14 | VIGS-G 14 line | 0.33D | 0.67B |
G15a | VIGS-G 15 line | 0.23D | 0.42C |
G16 | VIGS-G 16 line | 0.19DE | 0.60B |
G17 | VIGS-G 17 line | 0.18DE | 0.61B |
G18a | VIGS-G 18 line | 0.26D | 0.41C |
G19 | VIGS-G 19 line | 0.17E | 0.64B |
G20 | VIGS-G 20 line | 0.20DE | 0.63B |
G21a | VIGS-G 21 line | 0.28D | 0.35C |
G22 | VIGS-G 22 line | 0.18DE | 0.67B |
G23 | VIGS-G 23 line | 0.22D | 0.60B |
a Lower level of cell death in leaf discs.
b Higher level of cell death upon knock down of transcript as compared to leaf discs of the wild-type non-transformed plant (WT) or those of VIGS lines transformed with empty pRTV2.
3.4 Performance of Arabidopsis knockout T-DNA insertion mutants under salt stress
We speculated to detect the gene analogues in Arabidopsis and use their knockouts (Table S3) available at the SALK Institute, Genomic Analysis Laboratory (SIGnAL, http://signal.salk.edu/tdnaprimers.2.html) in order to detect the influence of these genes under salt stress. The SALK lines were cultivated in order to recover the mutants in homozygous state on an MS medium with kanamycin (100 mg/L) used as a selectable marker for T-DNA insertions. It is worth noting that the overall performance of mutant plants under normal growth conditions was similar to the WT plants. Seeds of five lines of T-DNA insertion homozygous KO mutants as well as those of the WT were first examined at different concentrations of NaCl (0, 100, 150 and 200 mM) in order to detect the day where shoot emergence starts. This sign is considered as seed germination. We observed that seed germination was initiated at day 3 under different salt concentrations, while day 4 for seeds put to germinate under control condition. We observed no further changes in germination rate across genotypes and treatments after day 5 of seed germination. Therefore, we decided to score germination percentages of different lines under different salt concentrations at day 6 (Fig. 7a). The results indicated almost no significant differences among genotypes up to 150 mM NaCl. However, significant differences in terms of germination percentages across genotypes were found under 200 mM NaCl. The germination percentage of KO mutant of Mlo gene (G18) was significantly higher than the other mutants as well as the WT under 200 mM NaCl. On the other hand, the germination percentage of BI-1 mutant (G13) was significantly lower than the other genotypes under the highest concentration.
![](/biologies/article/CRBIOL_2016__339_3-4_105_0/jpg/src/tex/figures/gr7.jpg)
a: effect of salt stress on germination percentage (at day 6) of seeds germinated on MS medium supplemented with different concentrations of NaCl (0, 100, 150 or 200 mM) for five selected Arabidopsis T-DNA insertion knockout mutants (G8, G13, G15, G18 and G21) as well as WT (Col). Concurrently, seeds of different genotypes were germinated at 0 mM NaCl, left for 2 weeks on MS agar plates, then transferred to the soil for two more weeks with different salt concentrations. Growth parameters, e.g., (b) root length (mm), (c) number of leaves per plant and (d) rosette area (cm2), were scored two weeks after salt treatment at different concentrations. Gene codes refer to those indicated in Table S3.
Other sets of plants of different genotypes were grown on the MS medium for two weeks under normal condition and transferred to the soil to grow for two more weeks under different salt treatments. Performance of 4-week-old plants for two mutants with contrasting genotypes (G8 referring to the apoptotic gene with RST domain and G13 referring to the antiapoptotic gene BI-1) based on PCD assay in tobacco, grown under 0 and 200 mM NaCl as compared to the WT control (Col) is shown in Fig. S3. The characteristics statistically analyzed for the 4-week-old plants of different genotypes are shown in Fig. 7b–d. The results for the mean root length (mm) and rosette area (cm2) indicated that BI-1 mutant had the highest performance under 200 mM NaCl as compared to the WT control, while Mlo mutant has the lowest. There are no significant differences among mutant and WT plants under different salt concentrations for the number of leaves/plant. We can conclude that the results of the salt stress experiment in Arabidopsis are in harmony to those of PCD assay for two gene analogs in tobacco.
The gene encoding the inhibitor of apoptosis-promoting Bax1 (BI-1) protein interacts with Bcl-2 protein to block the action of Bax in inducing PCD. In plant, BI-1 gene confers resistance against stress-induced metabolic perturbations or PCD [27,28] by blocking Bax gene product in orchestration with Bcl-2 gene. The gene encoding the functional barley MLO protein is a negative control element of cell death and of defence responses to the fungus Blumeria graminis [29]. Interestingly, it was reported that the gene encoding Bax inhibitor-1 induces breakdown of Mlo-mediated penetration resistance to the fungus in barley. The recessive mutant mlo2 alleles resulted in spontaneous leaf cell death [30] in pepper plants and enhanced resistance against the virulent fungus Xanthomonas campestris pv. vesicatoria (Xcv) [31]. In other words, the gene encoding BI-1 might block PCD either indirectly or directly. The first approach involves the interaction with Bcl-2 gene to suppress Bax encoding gene and the second involves cross talking with Mlo gene to block PCD. Accordingly, the reduced levels of PCD of knockdown line of Mlo encoding gene in the present study (Table 3) might be due to the accumulative effects of BI-1 gene (natural) and VIGS (artificial).
In conclusion, we can either retard apoptotic genes or overexpress antiapoptotic genes in plant breeding programs as a useful approach to improve the salt or drought stress tolerance of agronomically important field crops. This conclusion does not apply to all the studied genes, as their function is much more complicated to be adopted without experimentation.
Acknowledgements
The research team expresses its gratitude and appreciation to King Abdulaziz City for Science and Technology (KACST) for providing the research grant no. AT-32-92 (A-K-11-159). The research team also expresses its thanks to Deanship of Scientific Research (DSR) and Department of Biological Sciences, Faculty of Science, King Abdulaziz University, Jeddah, for providing services and assistance during the course of the project.