Enhanced photosynthesis and carbon metabolism favor arsenic tolerance in Artemisia annua, a medicinal plant as revealed by homology-based proteomics.
In recent years, metal/metalloid contamination has resulted in degradation of large areas of cultivable lands, posing a serious threat to sustainable agriculture and food security in developing countries like India . The frequent use of Asladen ground water for irrigation elevates As levels in the soils and makes land unproductive and toxic, thereby resulting in severe reduction of growth and yield of food crops [2,3]. This calls for judicious management of contaminated lands to put them back in cultivation. Although a number of As-tolerant noncommercial plants, including some species of the genera Pteris, Agrostis, Lepidium, Lemna, and Jasione, are known to accumulate high levels of As and have been successfully exploited for As phytoremediation [4, 5], none of these have any commercial use. Hence, it is highly desirable to identify nonedible plant species which may not only be cultivated for phytoremediation of As contaminated soils but also put back into commercial use.
Artemisia annua (Asteraceae) is an important dicotyledonous medicinal plant native from China and widely distributed in temperate and subtropical zones of the world, especially in Asia. It is a nonmodel plant with limited genomic information, although sequencing of limited number of randomly selected cDNA clones has been achieved . Despite the above, this plant has provoked wide interest for biotechnological studies because it produces an antimalarial drug, called artemisinin. Currently, artemisinin and its semisynthetic derivatives are extensively used in the treatment of malaria, mostly in combination therapies , as well as in the treatment of cancers and viral diseases [8, 9]. In addition to the therapeutic properties, A. annua possess important attributes such as a fast growth rate and a high biomass production, easy cultivation, and an unusual weedy habit, all of them important traits for plants employed for phytoextraction of metals and metalloids from contaminated soils . Our previous study demonstrated a significant induction of artemisinin production under As stress. The high tolerance of A. annua against As stress was coupled with an appreciable increase in biomass and As accumulation capacity, thus enabling this plant to be worth cultivating in As contaminated soils .
The enhanced artemisinin production and As tolerance of A. annua raises the question as which proteins could be potentially responsible for the above traits and how far the changes in different metabolic pathways, that is, primary and secondary, are interrelated with each other. While some mechanisms including, surface exclusion, intracellular detoxification, efflux, and production of phytochelatins have been proposed for the development of As tolerance in plants , the molecular mechanism underlying physiological changes still remains to be elucidated. Since proteins are directly involved in stress responses, any change in biochemical and physiological parameters is likely to be reflected at the proteome level. Proteomics has become a useful tool to unravel the relationships between protein abundance and stress acclimation . High resolution 2D-PAGE (2-DE) is the most direct approach for defining gene functions, and studies based on transcriptomics and proteomics can provide significant information on the processes involved in these stress responses .
Proteomics has been widely used to investigate the proteome of plants under abiotic stresses [15-17]. However, proteomic studies on As-induced changes are still few and limited to maize roots , maize shoots , rice roots , rice leaves , a pseudometallophyte Agrostis tenuis , and the hyperaccumulator species Pteris vittata fronds  and roots . With the exceptions of Agrostis and Pteris, all the above mentioned proteomic investigations primarily focus on plants which are agronomically important food crops. It is worth mentioning that country like India, where food-security is of foremost concern, to use such plants for phytoremediation of As-contaminated agricultural lands cannot be acceptable. In view of the fact that none of above plant species is well adapted to high As concentrations, A. annua emerged as the plant of choice. Despite potential application and tolerance the proteome of A. annua has remained unexplored due to unsequenced genome and inadequate representation in annotated, nonredundant (nr) protein databases. This study is the first attempt to decode the As tolerance strategy of A. annua combining physiological and proteomic analyses and to establish interactive network between its primary and secondary metabolism.
It requires mentioning that As (V) exposure (i) triggers ROS production causing alterations in the antioxidant pathways [11, 20], (ii) disturbs phosphate metabolism (e.g., ATP synthesis) and electron transport systems , and (iii) induces photoinhibitory conditions [25, 26]. In view of the above it was hypothesized that for survival under As stress, A. annua may (i) increase the accumulation of APX and DHAR to detoxify ROS, (ii) increase the electron transport to overcome the shortage of ATP and NADPH, and (iii) enhance phytochelatin synthase (pcs) and ABC transporter like protein genes for detoxification.
Taking recourse to the above the present study was undertaken to have an insight into (i) the metabolic proteins involved in As tolerance and to ascertain how far the physiological attributes go hand in hand with proteomic data and (ii) establish protein network between primary and secondary metabolism pathways. The results of this study enabled us to associate physiologically significant candidate proteins with As stress tolerance mechanism in A. annua.
2. Materials and Methods
2.1. Hydroponic Dose Screening. Artemisia annua L (family: Asteraceae) plants were cultured hydroponically as described by Rai et al. . Briefly, seeds of A. annua were collected from Central Institute of Medicinal and Aromatic Plants (CIMAP), Lucknow, India, and grown in botanical garden, Banaras Hindu University, Varanasi, India. Twenty one days after germination, seedlings of equal size (approximately 10 cm) were washed with deionized water and acclimated in a hydroponic system in 10% Hoagland's nutrient solution . The nutrient solution was aerated continuously using an aquarium pump and renewed twice a week. Seedlings were kept in a controlled room with a 14 h light period, 100 [micro]mol [m.sup.-2] [s.sup.-1] photosynthetic photon flux density (PPF), 28/26[degrees]C day/night temperature, and 60-70% relative humidity. After one week, the plants were transferred into 10% strength Hoagland's nutrient solution, spiked with different doses such as 50, 100, and 150 [micro]M of As. Arsenate was applied as [Na.sub.2]HAs[O.sub.4] x 7[H.sub.2]O. Experiments were conducted in triplicate in 150-mL Erlenmeyer flasks containing 100 mL of nutrient solution. Each replicate contained 3 plants of equal height (approximately 10 cm). Control seedlings were never treated with As. Three replicates, each including the two concentrations and the control, were placed separately in growth chamber and were harvested after 7 days of treatment. After 7 d of As exposure no visible toxicity symptoms were seen up to 100 [micro]M. Moreover, lower concentrations of As stimulated root growth and the biomass the maximum being 21.3 and 44.4%, respectively, at 100 [micro]M As over the control. However, toxicity symptoms appeared and the root length decreased by (6%) at 150 [micro]M. This clearly indicated that plants fail to tolerate As beyond 100 [micro]M (data not shown). Based on these data, the diagnostic concentrations for the proteome analysis were set at 100 and 150 [micro]M. Moreover, these doses represent the likely range of arsenic in contaminated lands of eastern Uttar Pradesh (U.P.), India .
2.2. Pot Culture Experiment. Based on the results of the hydroponic study another set of pot culture experiments were conducted in natural conditions (field) using soil treated with similar doses of As, 100 and 150 [micro]M As in the form of sodium arsenate. Experiments were done in triplicate. The soil used for the pot experiments was alluvial, collected from the premises of botanical garden, Banaras Hindu University, Varanasi, India. The soil was mixed with sand at the ratio of 3: 1 (v/v). As solution was uniformly mixed with air-dried soil, kept for one month to stabilize and filled in 16 kg pots. The soil without As treatment was used as control. Four seedlings of twenty-one day-old healthy A. annua were transplanted into each experimental pot. Soil moisture was maintained at 50-60% of field water holding capacity by adding water during the experimental period. Plants were grown under natural daylight and ambient temperature for eight months, after which they were harvested for protein isolation as well as to investigate the effects of As on flowering time, inflorescence, and other morphological as well as physiological characters.
2.3. Extraction and Estimation of Arsenic and Phosphate. Total As was estimated in an atomic absorption spectrometry (Perkin-Elmer; A Analyst 600) fitted with a graphite furnace . The As reference standard of 1000 mg [mL.sup.-1] (AA03N5) was supplied by Accustandard, USA. Total phosphate was extracted from shoots of control and As treated seedlings according to the method of Fiske and Subbarow . Absorbance was measured at 660 nm using K[H.sub.2]P[O.sub.4] as a phosphate standard.
2.4. Measurements of Photosynthesis and Chlorophyll Fluorescence. Photosynthetic rates ([P.sub.s]), stomatal conductance ([g.sub.s]), intercellular C[O.sub.2] (C;), water use efficiency (WUE), and transpiration rates (E) of test plants were determined using a LI-6400, (LICOR Inc., Neb., USA) portable gas exchange system on the uppermost fully expanded leaves from each individual during the period from 8:30 a.m. to 10:30 a.m. in six randomly selected plants of each treatment. During the experiments, system was calibrated with known C[O.sub.2] concentration of 509 ppm and PAR ranged between 1000 and 1200 [micro]mol [m.sup.-2] [s.sup.-1]; Chlorophyll fluorescence was determined with a portable plant efficiency analyzer (Model PEAMK2 9414, Hansatech Instrument Ltd., Norfolk, UK) on the same leaf where ([P.sub.s]) was measured. WUE was calculated as the ratio of [P.sub.s] and E. Leaf clips for dark adaptation were placed on the adaxial side for 15 min before measurement at excitation irradiance of 2000 [micro]mol[m.sup.-2] [s.sup.-1]. Minimum fluorescence ([F.sub.o]) and maximum fluorescence ([F.sub.m]) were measured. Variable fluorescence was calculated as ([F.sub.v] = [F.sub.m] - [F.sub.o]) and the ratio of variable to maximum fluorescence as ([F.sub.v]/[F.sub.m]).
2.5. Measurement of ATP, NADP, PSI, PSII, and Whole Chain Activity. ATP and NADPH contents as well as PSI and PSII activities were measured on the uppermost fully expanded leaves of each treatment and control (sampled in six). The ATP pool was measured by the method of Larsson and Olssonl . ATP was quantified using Luminometer (LKB, 1250). The NADPH/NADH level of the cell extract in TrisCl (pH 8.0) was measured by recording the absorbance at 340 nm . PSI- and PSII-mediated electron transport were measured using isolated thylakoids with an oxygen electrode (digital oxygen system model-10, Rank Brothers, UK). Thylakoid was isolated from leaves using a modified protocol of Mishra and Sabat . Leaves were crushed in extraction buffer containing 300 mM sucrose, 10 mM NaCl, 10 mM Ca[Cl.sub.2], 5mM Mg[Cl.sub.2], 10 mM ascorbate, and 20 mM Tricine-NaOH (pH 7.5). The homogenate was filtered and centrifuged at 6000 g for 10 min. The pellet represents the thylakoid fraction and was suspended in 1 mL of buffer containing 100 mM sucrose, 10 mM NaCl, 10 mM Ca[Cl.sub.2], 5mM Mg[Cl.sub.2], and 20 mM Tricine-NaOH (pH 7.5). The reaction mixture for assaying PSII activity consisted of 50 mM HEPES buffer (pH 7.5) containing 400 mM sucrose, 5mM Mg[Cl.sub.2], and 10 mM NaCl in a final volume of 1 mL. The activity was measured as oxygen evolution with water as electron donor and 500 [micro]M phenyl p-benzoquinone (pBQ) as the electron acceptor . The whole-chain ([H.sub.2]O-MV) electron transport was measured polarographically as [O.sub.2] uptake. The 1 mL reaction mixture contained 20 mM HePES (pH 7.5), 100 mM sucrose, 10 mM NaCl, 2 mM Mg[Cl.sub.2], 0.5 mM MV, and 1 mM sodium azide. In the third type of assay, the electron transport through PSI was also measured polarographically as above except that reduced DCPIP was used as the electron donor. Assay conditions were identical to those of whole chain assay except that 5 [micro]M DCMU (3-3,4-dichlorophenyl-1,1dimethylurea), 1mM sodium ascorbate, and 100 [micro]M DCPIP were added to the above reaction mixture. A thylakoid aliquot containing 10 [micro]g chlorophyll was used in the measurements. The chlorophyll content was estimated according to Arnon . All the photochemical assays were carried out at a saturating intensity of white light (1000 [micro]E[m.sup.-2] [s.sup.-1]) at 26 [+ or -] 2[degrees]C.
2.6. Determination of Artemisinin and Anthocyanin Contents. Artemisinin was extracted and quantified by the method of Zhao and Zeng . Concentration of artemisinin was measured using calibration curve of artemisinin standard. Anthocyanins were extracted from 1 g of fresh leaves with 20 mL of mixture: n-propanol: HCl: [H.sub.2]O (18:1:81, v/v/v) according to Lange et al.  and Bette and Kutschera  in modification. Samples were heated in a boiling water bath for 10 min and then incubated for 24 h in the dark at 4[degrees]C. The extract was filtered and its absorbance was measured at 529 and 650 nm. The red peak absorbance for degraded chlorophyll in the sample extracts occurs at 650 nm and has a tail that overlaps with the anthocyanin peak at 529 nm. Therefore, the absorbance of anthocyanin was corrected for the effect of chlorophyll using this equation: AA = [A.sub.529]-(0.288 [A.sub.650]) where AA is corrected anthocyanin absorbance. Total anthocyanin content was then calculated using this corrected absorbance and a molar absorbance coefficient for anthocyanin at 529 nm of 30,000 [mol.sup.-1] [cm.sup.-1] The content of total anthocyanins was expressed as [micro]mol per g FW .
2.7. Protein Extractionand 2-DE Separation. The experiments were conducted in triplicate. Three plants per treatment (both 100 & 150 [micro]M) and control were collected after 8 months of pot culture experiment. Two leaves per replicate (2nd pair and weighing 0.5 to 1.0 g fresh weight) were ground in liquid nitrogen to a fine powder. The powder was suspended completely in 10 mL of TCA extraction buffer (10% w/v TCA in acetone with 0.5% w/v [beta]-mercaptoethanol) at -20[degrees]C overnight. After centrifugation at 6000 xg, for 30 min at 4[degrees]C the pellets were washed 3 to 4 times with 10 mL ice cold acetone (with 0.5% w/v [beta]-mercaptoethanol). The collected protein pellets were dried in air. The dried powder was resuspended completely in 2 mL lysis buffer containing 7 M urea, 2 M thiourea, 4% CHAPS, 40 mM DTT, and 1.0% IPG buffer (4-7). Traces of bromophenol blue were added and the solution was centrifuged at 19000 xg for 10 min. Protein contents were quantified according to Bradford using BSA (bovine serum albumin) as a standard (2mg[mL.sup.-1]) . Sample entry was made through in-gel rehydration. A total of 300 [micro]L of solubilization buffer containing 250 [micro]g protein sample was incubated with the dry IPG strips (pH 4-7; 13 cm; GE healthcare, USA) at 20[degrees]C for 16 h. The first dimension separation was conducted at 20[degrees]C with an Ettan IPGphor system (GE Healthcare BIO-Science, USA). Focusing was performed in 7 steps: linear 30 V for 00:30 h, 150 V for 2:00 h, 300 V for 00:40 h, 500 V for 4 h, gradient 1000 V for 1h, gradient 8000 V for 2h, and finally 8000 V for 13:00 h. Focused IPG strips were then equilibrated by first incubating them in an equilibration solution (6 M urea, 30% v/v glycerol, 2% w/v SDS, 50 mM Tris-HCl, pH 8.8, and 1% w/v DTT and a trace amount of bromophenol blue) for 15 min, followed by incubation in 2.5% w/v iodoacetamide in the same equilibration solution instead of 1% DTT for 15 min. The strip was placed on the top of 12.5% SDS-PAGE (13 x 15 cm) and sealed with 0.5% agarose. Electrophoresis was carried out at 10mA/gel for 30 min then 25mA/gel for 5 h using a Hoefer SE 600 apparatus (Amersham Biosciences, USA). The gels were stained with CBB R-250 as described by Pandey et al. . Three replicates of CBB stained gels from three independent biological samples were used for analysis.
2.8. 2-DE Gel Data Analysis. Gels were scanned and protein spots were analyzed by using PDQuest software version 7.1 (Bio-Rad). After background subtraction and spot detection, spots were matched. Missing values imputation was done by K-Nearest Neighbor (KNN) algorithm [MATLAB]. After missing values imputation, total spot intensity per gel was used to normalize spot intensities (% of individual spot intensity/[SIGMA]% spot intensity of each gel) to compensate for variations between gel replicates. These data sets were also log-transformed to reduce the spot volume-spot deviation dependency. Only spots with statistically significant (Duncan Multiple Range Test, DMRT, P < 0.05) and reproducible changes were considered valid, and the protein spots with an abundance ratio of at least 1.5-fold were selected as differentially expressed proteins.
2.9. In-Gel Digestion and Peptide Extraction. 46 differentially expressed protein pots from CBB stained gels were manually excised and in-gel digestion was performed as per the protocol of Bruker Daltonics adapted from Schevchenko et al. . 100 mL of washing solution (1: 1 ratio of acetonitrile and 100 mM ammonium bicarbonate) was added to excised spots and incubated for 30 min for destaining. It was further incubated in absolute acetonitrile for 3 min and pellets were air-dried. Reduction of the protein spots was done using 10 mM DTT dissolved in 100 mM N[H.sub.4]C[O.sub.3] for 1h at 50[degrees]C. The samples were further alkylated using 55 mM iodoacetamide in 100 mM N[H.sub.4]C[O.sub.3] for 45 min in dark at room temperature. The gel pieces were dehydrated with 100% acetonitrile, air-dried, and then digestion was done by using sequence-grade modified trypsin (Promega gold) which was added (25 [micro]g/mL) in a digestion buffer containing 25 mM ammonium bicarbonate to the dehydrated gel pieces. Ingel digestion of protein was performed overnight at 37[degrees]C, followed by three extractions of peptides with 50% (v/v), acetonitrile 0.1% (v/v), and trifluoroacetic acid (TFA). Combined supernatants were dried in a vacuum centrifuge (Thermo electron corporation, UVS400A) for 2-4 h. Dry pellets were resuspended in a small volume (3-5 mL) of 0.1% (v/v) TFA in 50% (v/v) acetonitrile and stored at -20[degrees]C until analysis.
2.10. Mass Spectrometry Analysis and Database Searching. For MS analysis, one microliter of alpha cyanocinnamic acid (CCA) matrix (10mg/mL in 50% acetonitrile and 0.1% TFA) (Sigma Aldrich, St. Louis, Missouri, USA) was mixed with 1 [micro]L of digested protein samples. The mixture was spotted on MALDI target plate (MTP 384 ground steel, Bruker Daltonics, Germany) and allowed to form the crystals. The peptide spectra was acquired in an AUTOFLEX speed MALDI TOF/TOF instrument (Bruker Daltonics, Germany) having Nd: YAG smart Laser beam of 335 nm wavelength. External calibration was done with peptide calibration standard supplied by Bruker, with masses ranging from 1046 to 3147 Da. The obtained spectra were acquired using FlexControl version 3.3 software (Bruker Daltonics, Germany) in reflectron ion mode with an average of 2000 laser shots at the m/z detection range between 700 and 4000. The peak list was generated automatically and at least their most intense and discriminating peaks were subjected to further MS/MS fragmentation in LIFT mode.
Ion Source-1 19 kV and Reflectron (grid)-1 21 kV
Ion Source-2 16 kV and Reflectron (grid)-2 9.65 kV.
The data were analyzed using Flex Analysis software version 3.3 (Bruker Daltonics, Germany) and searched in MASCOT web server (Matrix Science; http://www.matrixscience.com) using Biotools version 3.2 software (Bruker Daltonics, Germany). The protein databases employed were NCBInr (The National Center for Biotechnology Information nonredundant) and Swiss-Prot (date of access, 5.11. 2012). A combined peptide mass fingerprinting (PMF) and tandem MS/MS were performed with the following MASCOT settings: taxonomy as Viridiplantae; peptide mass tolerance of [+ or -] 100 ppm for peptide mass fingerprinting and [+ or -] 1.2 Da for MS/MS, monoisotopic mass, alkylation of cysteine by carbamidomethylation as a fixed modification and oxidation of methionine as a variable modification. Proteins with threshold MASCOT score > 43 (P < 0.05) were considered as significant match. The peptide sequences with highest scores of each spot were validated by BLAST P analysis. Also, the borderline MASCOT hits that showed poor score and those annotated as unknown, predicted, or hypothetical proteins were further searched in BLAST P analysis against the NCBI nr database to query their homology.
2.11. RNA Isolation and DNase Treatment. Total RNA was isolated from the control and As treated leaf samples in TRIZOL reagent (GIBCO-BRL) using the instructions given in the manufacturer's protocol. Briefly, 100 mg of leaf samples wasgroundinliquidnitrogenand mixedwithTRIzol reagent. After centrifugation at 12,000 g for 15 min at 4[degrees]C, the upper aqueous layer containing RNA was removed and placed in a new centrifuge tube. About 0.5 mL of isopropyl alcohol per 1 mL of TRIzol was used to precipitate RNA. This suspension was centrifuged at 12,000 g for 8 min at 4[degrees]C and the pellet washed with 75% ethanol and resuspended in 50 [micro]L DEPC treated water. Forty units of RNase free DNase was added to remove any contaminating DNA followed by heating at 70[degrees]C for 5 min to denature the DNase. The quality of RNA was checked on a 1.0% agarose gel and the concentration determined by measuring the absorbance at 260 nm.
2.12. Primer Design and RT-PCR Analysis. The cDNAs of DXS, DXR, HMGR, GPPS, FPS, and CHS were synthesized using RT-PCR based on the A. annua nucleotide sequences published in database. The database employed was GenBank (GenBank accession numbers AF182286, AF182287, AF142473, and AF112881, AF112881, resp.). The sequences of the primers used for PCR amplifications are listed in Supplementary Table S2 in Supplementary Material available online at http://dx.doi.org/10.1155/2014/163962. Amplification was performed with initial denaturation at 94[degrees]C for 5 min, followed by 30 cycles of 30 s at 94[degrees]C, 30 s at 54[degrees]C, 90 s at 72[degrees]C, and a final extension at 72[degrees]C for 10 min. Ribosomal protein S9 (RPS9) was used as an internal reference gene. The amplification of each gene was repeated at least twice. Since the genome of A. annua is largely unsequenced, to design the gene specific primers, such as of GAPDH and PAL, the highly conserved region was picked up through multiple sequence alignment of different dicot sequences available in database by using Primer 3 software and these primers were used to amplify genomic DNA from A. annua. Specificity of all the primers was confirmed by sequence analysis of RTPCR amplicons derived from A. annua.
2.13. Western Blot Analysis. SDS-PAGE was carried out as per the method of Mishra et al. . 20 [micro]g protein was separated by SDS-PAGE and electroblotted onto PVDF membrane (Millipore Immobilon-P), using a dual Mini-electroblot system (Precision Instruments, Varanasi, India) as described . The gel cassette was kept in transfer buffer 3.03 g/L Tris base, 14.4 g/L glycine, and 200 mL/L methanol (99% v/v pure) for 12 h or overnight at 15 V at 4[degrees]C. Membrane was blocked for 4 h in TTBS (Tris buffer saline containing 0.1% Tween20) and 5% (w/v) nonfat dried milk. The primary antibody (anti-PCS antibody) was diluted as per the instructions of the donors. The membrane incubated overnight at 4[degrees]C in the diluted solution of the primary antibody was washed five times for 5 min each in TTBS. This was then incubated in a Goat anti-Rabbit IgG HRP (horseradish peroxidases) conjugated secondary antibody (Genei, India) for 4 h. Following four consecutive 5 min wash in TTBS the membrane was developed with DAB/Ni[Cl.sub.2] visualization solution. The reaction was terminated by washing the PVDF membrane with deionized distilled water. The blots were dried between filter paper, which considerably reduced the background staining. Polyclonal antibodies used for the detection of PCS by immunoblotting were obtained as generous gift from Dr. Stephan Cuine (CEA Cadarache, France). The immunoblotswere photographed using a gel documentation system (Bio-Rad, USA) and their quantitative densitometric analysis was performed using the Quantity One Software (Bio-Rad, USA).
2.14. Statistical Analysis. The statistical significance of the data was analyzed using one-way ANOVA, followed by Duncan's new multiple range test (DMRT) using the statistical software package SPSS 12.0. The results were accepted as significant at P < 0.05. All values shown in the figures and tables are means [+ or -] SD. For statistical treatment and cluster analysis of protein abundance values, the web-based software NIA array analysis tool used was Sharov et al.  available at http://lgsun.grc.nia.nih.gov/anova/index.html. This software tool selects statistically valid protein spots based on analysis of variance (ANOVA). The data (MS excel sheet with indication of biological replications) were statistically analyzed using the following settings: error model "max (average, actual)," 0.01 proportion of highest variance values to be removed before variance averaging, 10 degrees of freedom for the Bayesian error model, 0.05 FDR threshold, and zero permutations. Hierarchical clustering was performed to check the reproducibility of data set (based on 39 differentially expressed protein spots), and the results were represented in dendrograms using the Euclidean distance. PCA was done using the following settings: covariance matrix type, three principal components, 1-fold change threshold for clusters, and 0.6 correlation threshold for clusters. PCA results were represented as a biplot, with proteins with similarity in expression pattern located in the same area of the graph.
3.1. Arsenic and Phosphate Concentrations. A. annua accumulated As in a concentration dependent manner in both root and shoots. Maximum As content was found to be (191.01 [micro]g [g.sup.-1] DW) in 150 [micro]M As treatment as compared to the control (6.08 [micro]g [g.sup.-1] DW). In contrast, phosphate concentration showed a dose dependent decrease with lowest value in the 150 [micro]M As treated plants as compared to the control (Supplementary Table S1).
3.2. Morphological Responses of A. annua to As Stress. Upon visual inspection, A. annua plants did not show any phytotoxicity symptoms at 100 [micro]M As concentration. The plants looked healthy and green with an increased biomass, affirming their tolerant nature (Figure 1, Table 1). In contrast, 150 [micro]M As treated plants showed visible toxicity symptoms such as wilting, whitish chlorosis in young leaves and yellowing and necrosis in old leaves (Figure 1(c)). The most emblematic symptoms observed under 150 [micro]M As treatments were intense red coloration of leaflets particularly at the apex (Figure 1(c)). The pigment analyses from leaves of comparable developmental stages revealed a 3.0-fold higher anthocyanin content in 150 [micro]M As treated plants compared to the control (Table 1). The data compiled in Table 1 depicted an inhibitory effect of As on plant height, number of inflorescence, and inflorescence branches being more pronounced at 150 [micro]M As concentrations. Interestingly, however, the number of florets, the size of capitula, the ratio of capitula/inflorescence, and the number of oil glands were found to be increased significantly under As treatments conceivably contributing to the higher artemisinin and essential oils production in As treated plants (Table 1, Figures 1(h) and 1(i)). Another remarkable observation made was the early onset of flowering in As treated plants. As treated plants flowered about 4-5 weeks before than control plants after sowing (Table 1).
3.3. Physiological Responses of A. annua to As Stress. Table 2 compiles data on the various physiological traits such as photosynthetic rate ([P.sub.s]), stomatal conductance ([g.sub.s]), transpiration rate ([pounds sterling]), intercellular C[O.sub.2] concentration ([C.sub.i]), water use efficiency (WUE), and chlorophyll fluorescence emission patterns ([F.sub.v]/[F.sub.m]) of control and As treated A. annua plants. While in the 100 [micro]M As treatments [P.sub.s] was found to be slightly depressed by (4.7%) and was near the control, the 150 [micro]M As treatments resulted in a reduction of the leaf photosynthetic activity by 26% with respect to the control. With regard to [g.sub.s], the value of this parameter decreased by 18% and 23% in 100 and 150 [micro]M As treatments, respectively. However, in contrast, the intercellular C[O.sub.2] concentration showed a little increase under 100 [micro]MAs treatments, although was equivalent to the control Transpiration rate demonstrated a negligible decrease under As treatments. An increase in the water use efficiency (WUE) due to As stress was observed in both treatments, being 37% and 34% of the control in 100 [micro]M and 100 [micro]M As treatments, respectively. Considering photosynthetic electron transport chain activities, PSII demonstrated a minor decline of 11% and 34% in 100 [micro]M and 150 [micro]M As treatments, respectively, whereas PSI showed a significant increase of 70% in 100 [micro]M As treatments. In addition, the chlorophyll fluorescence parameters were monitored to determine the performance of photosystem 2 (PSII). The [F.sub.v]/[F.sub.m] ratio is correlated with the efficiency of leaf photosynthesis and a decline in this ratio is a good indicator of photoinhibitory damage under stress . The [F.sub.v]/[F.sub.m] ratios in the 100 [micro]M As treatments were slightly below 0.83, which is the theoretical optimum . In contrast it dropped down to 0.72 in the plants exposed to 150 [micro]M of As, indicating some damage in the PSII .
3.4. Biochemical Responses of A. annua to As Stress. Biochemical parameters including anthocyanin, artemisinin, ATP, and NADPH contents were measured from 8-month-old plants (Tables 1 and 2). Anthocyanin content was increased by 175% and 275% in 100 and 150 [micro]M As treatment, respectively. In tune with the above, the concentration of artemisinin was increased by 53.1 and 75% in 100 and 150 [micro]M As treatments, respectively, over the control (Table 1). While in 150 [micro]M As treatments, both NADPH and ATP contents decreased by 24% and 29%, NADPH registered an increase of 29% in 100 [micro]M As treatments, but insignificant increase in ATP content was observed (Table 2).
3.5. 2-DE Analysis of As-Responsive Leaf Proteins in A. annua. Figure 2 presents a comparative account of the CBB stained gels of untreated control and 100 and 150 [micro]M As treated A. annua seedlings after 8 months of treatment. Quantitative image analysis using PD Quest Version 7.1 software (Bio-Rad) from three biologically independent replicate experiments revealed a total of 391 [+ or -] 35, 412 [+ or -] 41, and 447 [+ or -] 27 spots in the control, 100 [micro]M and 150 [micro]M As treatments in the range of pH 4-7, and relative molecular masses 14-90 kDa. Supplementary Table S3 contains the description of the 2-DE proteomic experiment.
Fifty-one protein spots consistently changed their abundance (vol. %) by [greater than or equal to] 1.5-fold in at least one of the As treatments. However, a spot displaying control versus treatment ratio of at least 1.5-fold and consistently present or absent in all three replicates was considered valid. Fold change was calculated as a ratio of the averaged means of normalized spot volumes of control and the treatments. Essentially arbitrarily threshold values ranging from 1.3 to 2.0 fold have been used in earlier proteomic studies [48,49]. Nevertheless, a 1.5 threshold value was selected in order to focus on the most responsive proteins and for consistency with previous proteomic and microarray experiments [18, 41, 50]. While most spots showed quantitative changes, few protein spots demonstrated qualitative changes. For example, spots 29, 30, and 32 were absent in As treated gels at one or more time points, whereas spots 9,14,39, and 40 were absent in the control gels (Supplementary Table S3). Raw data was processed in three steps, missing value estimation, normalization, and transformation. Analysis of variance before and after data transformation is given as Supplementary Table S4.
3.6. MALDI-TOF/TOF Analysis, Functional Classification, and Identification of As-Responsive Proteins. Forty-six differentially expressed protein spots (based on P [less than or equal to] 0.05) with reproducible alterations were subjected to MALDI-TOF/TOF analysis (Supplementary Table S5). Figure 3 shows the Venn diagram of the differentially expressed proteins at different concentrations of As. From 46 spots excised, 37 were successfully identified; two were not hit (spots 7 and 15) whereas 7 proteins were annotated either as unknown, putative, or hypothetical proteins (Table 3, Supplementary Figure S1). To gain functional information about these proteins, BLASTP (http://www.ncbi.nlm.nih.gov/BLAST) was used to investigate their homologies with other proteins in the database. Of the seven proteins, three were successfully BLAST identified (spots 5, 20, and 23); however, spots 32, 33, 43, 45, and 46 showed a poor BLAST score. Since the genome of A. annua is still unsequenced therefore the MS/MS results (the peptide sequences with highest scores) of each spot were also validated by BLAST P analysis which revealed similarity thereby validating the MS/MS results. These proteins were sorted into 7 functional groups including energy metabolism, ROS scavenging and defense, secondary metabolism, transcriptional regulator, protein metabolism, transport, and unknown/unidentified proteins according to the ontological classification of Bevan et al.  (Supplementary Figure S2). BLAST search results are given as Supplementary Table S6.
3.7. Expression Profiles of Proteins Associated in Different Metabolic Pathways. Expression pattern of selected proteins of C, control, and 100, 150 [micro]M As treatments is depicted in Supplementary Figure S3. Few proteins were represented by more than one spot with slightly different Mr and pI values, although they were excised from the same gel. For example, Chl a/b type I was identified in four spots (2, 10, 11, and 20), NADP-ferredoxin oxidoreductase was identified twice (8 and 38), carbonic anhydrase was identified in three spots (17, 27, and 40), triose-phosphate-isomerase (spots 42 and 44), and glyceraldehyde-3-phosphate (spots 6, 23, and 37). Most of the protein matches and theoretical and experimental pI and Mr were in good agreement, supporting confidence in the identifications; however, some of the identified proteins showed discrepancy with their Mr and pI. These kinds of phenomena are commonly observed in 2-D gels for several reasons including protein migration, different isoforms derived from different genes, proteolytic cleavage, posttranslational modification variants of the same gene product, or artificial modification of proteins, such as carbamylation during protein extraction . According to the expression pattern, the identified proteins are categorized as proteins showing increased or decreased abundance under 100 or 150 [micro]M As treatments. More than any functional group, energy metabolism proteins (60%) were differentially expressed in 100 [micro]M As treatment and were divided into four functional subgroups: (1) chlorophyll a, b-binding proteins, (2) oxygen-evolving complex proteins of PSII, (3) proteins participating in the Calvin cycle, and (4) proteins involved in photosynthetic electron transport chain. The first subgroup harbored four protein species of chlorophyll a/b binding proteins (CP26 in PS II) displaying a mixed pattern of expression. While two isoforms of chlorophyll a/b binding protein (spot 10 and spot 20), an unknown protein that was BLAST identified as chlorophyll a/b binding protein, were found in higher amounts in either 100 or 150 [micro]M As treatments, spot (2) was decreased in relative abundance in response to As stress. Another protein (spot 11) identified as chloroplast PS I type III chlorophyll a/b binding protein was increased in relative abundance in 100 As treatment. The second subgroup consists of oxygen-evolving complex protein (OECP2) where one intact protein (spot 4) was increased in relative abundance (2.3-fold in 100 [micro]M and 2.2-fold in 150 [micro]M), whereas one degraded unnamed protein product (spot 5) BLAST identified as oxygen-evolving complex protein (OEC1) was decreased in amounts (1.3-fold and 1.2-fold in 100 and 150 [micro]M As, resp.). The third subgroup consisted of 15 proteins. Most of the enzymes involved in different metabolic pathways such as plastidic aldolases (ALD, spot 12, 2.69-fold, and spot 36, 1.5-fold), fructose bisphosphatase aldolase class I (FBA, spot 36, 1.5-fold), glyceraldehyde 3-phosphate dehydrogenase (GAPDH, spot 6, 1.8 fold, spot 23, 3.1 fold, and spot 37,1.6 fold), Ribose-5-phosphate,3-epimerase (spot 41, 1.7 fold), triose-phosphateisomerase (spot 42, 1.6 fold and spot 44, 1.9-fold), and carbonic anhydrase (CA, spot 17, 3.1 fold, spot 27, 3.22 fold, and spot 40, 1.64-fold) were increased in relative abundance in 100 [micro]M As treatment except for transketolase (TK; spot 1) and RuBisCO activase (spot 22). Consistent with changes in Calvin cycle, three other important proteins involved in electron transport chain (fourth subgroup) such as, ATP synthase a subunit (spot 21, 2.7-fold), rieske protein (FeS, spot 24, 2.4-fold), and ferredoxin-NADP reductase (FNR, spots 8 and 38, 2.4- and 1.7-fold) were significantly enhanced only under 100 [micro]M As treatments. A marked increase in P protein subunit of glycine decarboxylase (spot 14, 2.0-fold in 100 [micro]M and 2.5-fold in 150 [micro]M) was found in As stressed plants, this being more significant in case of 150 [micro]M As treatment. Cytoplasmic malate dehydrogenase (MDH, spot 39) depicted a significant upregulation (2.12-fold) under 150 [micro]M As but no change under 100 [micro]M As treatment. Two upregulated antioxidant proteins identified in this study, dehydroascorbate peroxidase (DHAR, spot 18) and ascorbate peroxidase (APX, spot 19), revealed significant accumulation of 3.4- and 3.1-fold in 100 [micro]M As treatments. Three proteins of the transcriptional regulator category, namely, glycine rich RNA binding protein (SGRP, spot 25), maturase K (spot26), and LEAFY like protein (spot 28), displayed similar abundance patterns in 100 and 150 [micro]M As treatments, with a slightly higher amount in 150 [micro]M As treatments. Four protein spots belonged to protein metabolism category. Of the four proteins, the rossmann NAD(P) fold (spot16) was found to be increased in both 100 and 150 [micro]M As treatments, whereas ribosomal L12 1a (spot 13) revealed a 2.5-fold increase in 100 [micro]M As treatments and decrease in abundance in 150 [micro]M As treatments. In contrast to the above, Hsp 70 (spot 29) and chloroplast protease (spot 30) revealed decrease in relative abundance under 100 and 150 [micro]M As treatments. Besides this, two proteins were identified belonging to secondary metabolism category. The spot number 31 was identified as S-adenosyl methionine synthase (SAMS), a signaling protein involved in biosynthesis of homocysteine and ethylene, and revealed a 1.5-fold induction under 100 [micro]M As treatment. The other protein, chalcone synthase (spot 9), revealed a 2.2- and 2.3-fold increase in the amount of this protein under 100 and 150 [micro]M As treatments, respectively. Spot 35 was identified as a multidrug resistance-associated protein (MRP), a subclass of ATP-binding cassette (ABC) transporters depicted a 1.5-fold upregulation under 100 [micro]M As treatment. There were five protein spots that did not find match to any of the characterized proteins in the database and hence were named as UK (unknown).
3.8. Multivariate Statistical Analysis of Protein Abundance. Principal component analysis (PCA) was done to obtain a more accurate grouping of the samples and to determine the most discriminant spots. In order to check out the reproducibility in the biological replicates, a hierarchical clustering was performed. Two main clusters were obtained in the dendrogram, namely, cluster 5 (control) and cluster 4 (100 and 150 [micro]M As treatments) (Figure 4(a)). This clustering indicated that 100 and 150 [micro]M As treatments had a somewhat similar protein abundance profiles but different from that of control. Moreover, the hierarchical clustering of biological repetitions confirmed that the data were reproducible for the experiments. PCA analysis indicated a distinct separation of 100 and 150 [micro]M As treatments (Figure 4(b)). Photosynthetic proteins were found in the categories of proteins correlated with PC1 (positive and negative direction) and PC3 (positive direction). Carbon metabolism proteins were found in the categories of proteins correlated with PC1 (negative direction), PC2 (positive and negative direction), and PC3 (positive direction). Proteins belonging to transcription as well as protein metabolism category were found in PC1 (positive and negative direction) and PC3 (positive direction). Defense and stress-related proteins were found in PC1 (negative direction). Transport proteins were found in PC3 (positive direction). A high proportion of photosynthetic, carbon metabolism and defense proteins were more abundant in 100 [micro]M As treatments with respect to the 150 [micro]M As treatments. The quantity of protein abundance change within a specific PC was calculated by the slope of regression of log-transformed protein abundance versus the corresponding eigenvector multiplied by the range of values within the eigenvector (Figure 4(c)).
3.9. RT-PCR Analysis of Isoprenoid/Artemisinin Biosynthetic Pathway Genes. The expression profile of eight genes encoding enzymes of isoprenoid (MEP) and artemisinin biosynthetic pathways were monitored in the control and As treated plants (Figure 5). While all the six genes were upregulated in response to 100 [micro]M As treatment, a significant overexpression of 4.4- and 1.2-fold was observed for 1-deoxyxylulose-5-phosphate synthase (DXS) and 1-deoxy-D-xylulose 5-phosphate reductoisomerase (DXR), respectively. A similar pattern of gene expression was observed under 150 [micro]M As treatment. GAPDH also depicted an induction of 1.0- and 1.1-fold under 100 and 150 [micro]M As treatments, respectively. In contrast, geranyl pyrophosphate synthase (GPPS) and 3hydroxy-3-methyl-glutaryl-CoA reductase (HMGR) showed minor down regulation in response to both As treatments (Figure 5).
An integrated morphological, physiological, and proteomic approach has been used for the first time to investigate the mechanism of As tolerance of A. annua, a medicinal plant whose genome has not yet been sequenced. Although 46 proteins identified by MALDI-TOF/MS/MS and MS-BLAST analyses represent only a very small fraction of A. annua leaf proteome, some novel As stress responsive proteins appear involved in cellular homeostasis and survival of A. annua plants under As stress.
4.1. As Regulates Plant Architecture and Flowering Time in A. annua. The alterations observed in inflorescence branching patterns, vegetative axillary shoot development, and flowering time in As treated A. annua clearly indicate that these changes in plant form may arise, at least in part, due to change in expression pattern of conserved regulators or proteins. It is quite likely that the early onset of flowering and the increased number of florets and capitulum size in As treated A. annua were due to the significant accumulation of LEAFY- (LFY) like protein, a plant-specific transcription factor involved in growth, development , and flowering in Arabidopsis , first time observed under As stress. Furthermore, the red coloration of the leaf apex and shoots of 150 [micro]M As treated A. annua plants could be due to transient accumulation of anthocyanin pigments below the leaf epidermis. This finding is supported by an increased accumulation of chalcone synthase, involved in the production of secondary metabolites and pigments such as flavonoids and anthocyanin as also observed in As-stressed plants .
4.2. Photosynthetic Apparatus under As Stress. The physiological traits monitored revealed contrastingly differences under 100 and 150 [micro]M As treatments. While PSI activity, transpiration rate, ATP, and NADPH contents of 100 [micro]M As treated A. annua plants showed an increase over the control values, a decrease for these parameters was recorded at 150 [micro]M As treatments. On the other hand, PSII activity registered decline under both the treatments, this being more pronounced at 150 [micro]M As treatments. The [F.sub.v]/[F.sub.m] ratio, considered as indicator of the photochemical processes in photosystem 2 (PSII), depicted only a minor inhibition at 100 [micro]M but severe reduction under 150 [micro]M As treatments. Above findings were supported by the protein data having increased abundance of ATP synthase a subunit, rieske protein, and ferredoxin-NADP reductase (FNR) under 100 [micro]M As treatments, suggesting their likely role in As tolerance. Moreover, this enhancement unswervingly indicates operation of Fd-dependent cyclic electron transport (CET) in the test organism. In general, plants exposed to abiotic stresses show downregulation of linear electron flow (LEF) and activation of cyclic electron transport (CET) when linear electron flow becomes saturated [56, 57]. It could be anticipated that under 100 [micro]M As treatment the LEF is partially replaced by CET around PSI and the upregulated CET in turn provides energy to the Calvin cycle  thereby meeting the energy demand. Appreciable increase in NADPH content seems justified in view of its requirement as reducing equivalents during C[O.sub.2] fixation on one hand and maintaining a ratio of reduced glutathione for scavenging ROS. By contrast, under 150 [micro]M As treatments, the down regulated FNR disrupts CET. The electrons so accumulated suppress the charge separation of PSII reaction centre leading to damage of PSII as attested by the results of decreased [F.sub.v]/[F.sub.m] ratios and decreased ATP content (Table 2, Figure 7). Four proteins of chlorophyll a/b binding category (CP26 in PS II) displayed a mixed expression under 100 and 150 [micro]M As treatments. Altered levels of multiple forms of this protein have been reported in rice seedling under high temperature  and [H.sub.2][O.sub.2] stress . Another protein (spot 11) identified as chloroplast PS I type III chlorophyll a/b binding protein was increased in 100 [micro]M As treatment. This may confer enhanced resistance of PSI as known for Hordeum vulgare against boron stress . While (OECP2) protein spot 4 (intact) depicted increased relative abundance, the protein (spot 5) BLAST identified as oxygen-evolving complex protein (OEC1) was decreased under both treatments. OECP2 are manganese stabilizing proteins playing crucial roles in photosystem II stability , their decreased accumulation supports the damaged PSII under As treatments. A marked increase in P protein subunit of glycine associated with photorespiration was found in 150 mM As-stressed plants. Its overexpression stimulated photorespiration in Anabaena sp. under As stress . These results clearly demonstrated that inhibition of photosynthetic rate in 150 [micro]M As treated A. annua was not the consequence of stomatal closure (reflected by high intercellular C[O.sub.2] value) hence C[O.sub.2] limitation may not be the real cause for reduced carbon assimilation rates. The decrease in assimilation rates could be due to low ATP and NADPH contents (Table 2) and consequent impairment of the photosynthesis [47, 62]. Furthermore, the photosynthetic rate is known to weaken in the presence of heavy metals due to alterations in the active site of rubisco subunits [63, 64]. By contrast, PSI complex machinery might be envisioned to support the damaged PSII and maintain NADPH pool in 100 [micro]M As treated plants as observed in As-tolerant cordgrass Spartina densiflora  and Ricinus communis .
4.3. Metabolic Networks under As Stress. Of the 15 identified proteins of the metabolic pathways like PLA, FBA class I, GAPDH, RuBP-5-P, TPI, and CA were increased in 100 [micro]M As treatment except for TK and RuBisCO activase. This finds support from Pandey et al. , where TK was decreased in abundance under As stress. Transketolase is a universal amphibolic enzyme catalyzing reactions in the Calvin cycle and the oxidative pentose pathway producing erythrose-4-phosphate, leading to phenylpropanoid metabolism. Thus, any change in its abundance is expected to affect both photosynthetic carbon assimilation and secondary metabolism in plants. Despite a decrease in the abundance of transketolase in As treated plants neither Calvin cycle nor secondary metabolite synthesis (artemisinin, carotenoid, etc.) reflected any downregulation. This indicates the likely complementation of its role by some other enzymes. Henkes et al.  reported that Arabidopsis TK shares 23% identity with DXS which depicted increased expression under decreased TK activity. Further, TK has been reported to be a DXS like protein  and hence speculated to play a role similar to transketolase. An in silico analysis of TK showed 38% homology with DXS in A. annua and RT-PCR analysis revealed a significant transcript in As treated A. annua, which together support our above speculation (Figure 5). Of the three protein species corresponding to chloroplastic CA, two appeared as products of the same gene (spots 17 and 27) whereas the third one (spot 40) was different. An increased accumulation of carbonic anhydrase may be required for interconversion of C[O.sub.2] and HC[O.sub.3.sup.2-] and is thought to play a role in photosynthetic carbon assimilation as known in Arabidopsis  and other C3 plants [ 70,71]. The increased abundance of various isoforms of FBA, GAPDH, and TPI clearly indicates the activation of the glycolytic pathway under As stress. FBA reversibly catalyzes the conversion of fructose 1, 6-bis phosphate (FBP) to GAP and dihydroxy acetone 3-phosphate (DHAP). During stress, GAP and FBP may be converted to glucose 6-phosphate for reentry into the PPP for NADPH synthesis. In the light of the reports of Pandey et al. , one is tempted to presume that enhanced expression of FBA maintained the glycolysis, PPP, and turnover rate of Calvin cycle. In contrast to the report of Ahsan et al. , significant accumulation of three isoforms of glyceraldehyde-3P-dehydrogenase (GAPDH) was noticed. GAPDH detoxifies As using arsenate instead of phosphate converting glyceraldehyde 3-phosphate into 1-arseno-3-phospho-glycerate. Enhanced expression of GAPDH has also been detected in Pteris under As stress  rice seedlings under drought, submergence, and ABA treatment . These findings imply that increased accumulation of GAPDH may contribute to abiotic stress tolerance in plant. Of particular interest was the cytoplasmic MDH depicting significant accumulation under 150 [micro]M but no change under 100 [micro]M As treatment. This is conceivable because As stress leads to phosphate (Pi) starvation in plants (Supplementary Table S1). Though arsenic substitutes for phosphorus in the plants, it cannot function like phosphate, hence plants realize phosphate deficiency . Under these situations, plants ought to cope up with Pi deficiency. Taken together, these results suggest an increased turnover of C[O.sub.2] assimilating pathway and redox energy production in the 100 [micro]M As treated plants which are required for survival under downregulated photosynthetic rate. In contrast, decrease in carbon metabolism proteins associated with breakdown of energy generation results in failure of the plants to cope up with 150 [micro]M As toxicity.
4.4. Regulation of Antioxidant Machinery, Protein Metabolism, and Transport under As Stress. Two antioxidant proteins identified in this study include dehydroascorbate peroxidase and ascorbate peroxidase in 100 [micro]M As treatments. While DHAR (spot 18) has been reported to offer tolerance against oxidative stress , APX catalyses the reduction of [H.sub.2][O.sub.2] to water in ascorbate-glutathione cycle . Increased abundance of APX under stress is nicely documented in earlier proteomic works [76, 77]. It is quite likely that the upregulated ascorbate-glutathione cycle regulates excess ROS generated at PSII thereby offering tolerance in 100 [micro]M As treated plants. This is in accordance with the findings of Rai et al. . Three proteins belonging to transcriptional regulator category, namely, glycine rich RNA binding protein, maturase K, and LEAFY like protein displayed analogous abundance in 100 and 150 [micro]M As treatments, with a slightly higher level in 150 [micro]M As. Glycine rich RNA binding proteins are implicated in osmotic, cold stress, and RNA splicing activity , and maturase K functions like a ribosomal protein with RNA chaperone activity . Increased accumulation of these proteins suggest their conceivable role in overcoming the As stress. Another interesting observation was upregulated protein synthesis in A. annua under As stress. Of the four proteins belonging to the protein metabolism category, the rossmann NAD(P) fold, a structural motif found in nucleotide binding proteins , was increased under both As treatments, whereas ribosomal L 12 was increased in 100 [micro]M but decreased in 150 [micro]M As treatments. In contrast to the above, two molecular chaperones, Hsp70 and chloroplast protease were decreased under both As treatments. Significant down regulation of Hsp70 suggests it to be the key target under As stress as observed in Poplar trees under Cd stress  and needs further investigation.
Of particular interest are the multidrug resistance-associated proteins (MRPs), a subclass of ATP-binding cassette (ABC) transporters with a hypothetical role in As extrusion in plant. Our previous finding  impelled to check the role of pcs gene in As detoxification. The upregulated (1.5-fold) mdr-like ABC transporter and pcs gene (as confirmed by Western blot analysis, Figure 6) under 100 [micro]M As suggest that detoxification was brought about by reduction of As(V) to As(III) and sequestration with phytochelatin. The increased abundance of (ABC) transporters has also been reported in Comamonas sp. under As stress . Another interesting observation was the increased accumulation of S-adenosyl methionine synthase (SAMS), known to act as a signaling molecule and protects cells from arsenic stress by forming Sadenosyl-methionine . Furthermore, in the light of the report of Noriega et al.  that SAM protects cells from deleterious effects of cadmium, its increased accumulation under 100 [micro]M As treatment seems justified.
4.5. Hypothetical Model Depicting As Tolerance and Interaction between Primary and Secondary Metabolism in A. annua. Taking recourse to 39 As-responsive proteins as well as the data available in the literature, a protein network model has been proposed. This network consists of several functional components, including ROS production and scavenging, electron transport chain, protein metabolism, photosynthesis, energy supply, and biosynthesis of signaling molecules (Figure 7). The model clearly shows that As(V) may be taken up by the test plants through phosphate transporter Pht1 and reduced to arsenite As(III) by arsenate reductase (ACR2/CDC 25) in the cell. As(III) is rapidly complexed with soluble thiols and phytochelatins. Upon binding to phytochelatins (PCs), PC2-As(III) complexes are deposited into the vacuole by ABC transporters. This is clearly reflected by increased phytochelatin content and transcript of pcs gene under As stress. To cope up with the As stress, A. annua plants in turn activate their antioxidative and photosynthetic machineries. Two groups of proteins identified in multiple isoforms are chlorophyll a, b binding protein as well as OEC of PSII. OEC appears to be involved in the oxidation of water molecules to form molecular oxygen and ATP production via proton pumps. The excess electrons are transferred to oxygen at PSI or via Mehler reaction. In order to protect PSII from ROS, As-stressed plants augment excess electron removal by activating ascorbate-glutathione cycle such as APX and DHAR. An increased abundance of FBA, Ru5P, ALD, and CA (Table 3) points not only enhanced turnover of C[O.sub.2], assimilating pathway but also augmentation in ATP and NADPH synthesis. This was clearly reflected by an increased cyclic electron flow around PSI through increased accumulation of Fd-NADP reductase, rieske-Fes, and ATP synthase. The increased ATP and redox energy support enhanced C[O.sub.2] assimilation. The enhanced glycolysis provides not only reducing power required for As detoxification and repair but also precursors for biosynthesis of aromatic amino acid, flavonoid, and anthocyanin (via shikimate pathway) and artemisinin (via MEP pathway). Interestingly, two intermediates of glycolysis, GAPDH and pyruvate, are direct substrates of the DXP synthase. The upregulation of GAPDH on 2DE gel indicates persistent triose supply for the synthesis of artemisinin, anthocyanin, and carotenoids. Furthermore, the upregulated Calvin cycle and PPP produce erythrose-4-phosphate, a precursor for the shikimate pathway leading to phenylpropanoid metabolism. The rapid consumption of glyceraldehyde-3-P and other primary metabolites requiring continuous replenishment of carbon intermediates may in turn be provided by ribulose-5-P (with 2-fold accumulation on 2DE gel). Presence of only one protein (CHS) belonging to a flavonoid biosynthetic pathway may be due probably to limited number of protein spots analyzed. The network connection between primary and secondary metabolism of A. annua measured by the transcript of enzymes of flavonoid and artemisinin biosynthetic pathway such as DXS, DXR, FPS, HMGR, GPPS, CHS, and PAL indicated their fine regulatory network.
Of the different proteins identified in this study, proteins of photosynthesis and carbon metabolism appear as possible candidates for augmenting As tolerance in A. annua, although the exact mechanism of tolerance remains to be investigated. Our results also suggest that ABC-transporter like protein maybe involved in As sequestration into vacuoles as the final detoxification step. Data obtained in this study provide clues to the molecular aspects of As tolerance, as well as in development of As resistant crop varieties through plant breeding and biotechnological interventions.
Abbreviations ANOVA: Analysis of variance CBB: Coomassie Brilliant Blue CHAPS: 3-(3-Cholamidopropyl) dimethylammonio-1-propane sulfonate FBA: Fructose bisphosphatase aldolase class I GAPDH: Glyceraldehyde-3-phosphate dehydrogenase Ru-5-P: Ribulose-5-phosphate, 3-epimerase CA: Carbonic anhydrase TK: Transketolase MDH: Malate dehydrogenase TPI: Triose phosphate isomerase GRP: Glycine-rich RNA binding protein DMAPP: Dimethylallyl diphosphate DXS: 1-Deoxy-d-xylulose 5-phosphate synthase DXR: 1-Deoxy-d-xylulose 5-phosphate reductoisomerase FPS: Farnesyl diphosphate synthase HMGR: 3-Hydroxy-3-methylglutaryl coenzyme A reductase GPPS: Geranyl pyrophosphate synthase IPP: isopentenyl diphosphate MEP: 2-C-methyl-d-erythritol-4-phosphate MVA: Mevalonic acid CHS: Chalcone synthase PAL: Phenyl ammonia lyase PCA: Principal component analysis FDR: False discovery rate.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Rashmi Rai, Sarita Pandey, and Alok Shrivastava thank the Council of Scientific and Industrial Research for Senior Research Fellowship. Professor L. C. Rai FNA, Coordinator Centre of Advanced Study in Botany, Banaras Hindu University, is gratefully acknowledged for liberally providing the facilities for proteomic work as well as for his critical reading of the paper and valuable suggestions. The authors thank Mr. Sanjeev Kumar of BioCOS Life Sciences Private Limited, Bangalore, India, for multivariate analysis.
 S. A. Adejump, A. O. Togun, J. A. Adediran, and M. B. Ogundiran, "In-situ remediation of heavy metal contaminated soil using Mexican Sunflower (Tithonia densiflora) and Cassava waste compost," World Journal of Agricultural Sciences, vol. 7, pp. 224-233, 2011.
 D. Chakraborti, S. C. Mukherjee, S. Pati et al., "Arsenic groundwater contamination in Middle Ganga Plain, Bihar, India: a future danger?" Environmental Health Perspectives, vol. 111, no. 9, pp. 1194-1201, 2003.
 J. F. Ma, N. Yamaji, N. Mitani et al., "Transporters of arsenite in rice and their role in arsenic accumulation in rice grain," Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 29, pp. 9931-9935, 2008.
 X. Cao, L. Q. Ma, and C. Tu, "Antioxidative responses to arsenic in the arsenic-hyperaccumulator Chinese brake fern (Pteris vittata L.)," Environmental Pollution, vol. 128, no. 3, pp. 317-325, 2004.
 J. R. Rofkar, D. F. Dwyer, and J. M. Frantz, "Analysis of arsenic uptake by plant species selected for growth in northwest Ohio by inductively coupled plasma-optical emission spectroscopy," Communications in Soil Science and Plant Analysis, vol. 38, no. 17-18, pp. 2505-2517, 2007.
 W. Wang, Y. Wang, Q. Zhang, Y. Qi, and D. Guo, "Global characterization of Artemisia annua glandular trichome transcriptome using 454 pyrosequencing," BMC Genomics, vol. 10, article 465, 2009.
 R. K. Haynes, "From artemisinin to new artemisinin antimalarials: biosynthesis, extraction, old and new derivatives, stereochemistry and medicinal chemistry requirements," Current Topics in Medicinal Chemistry, vol. 6, no. 5, pp. 509-537, 2006.
 T. Efferth, M. R. Romero, D. G. Wolf, T. Stamminger, J. J. G. Marin, and M. Marschall, "The antiviral activities of artemisinin and artesunate," Clinical Infectious Diseases, vol. 47, no. 6, pp. 804-811, 2008.
 T. Efferth, "Willmar Schwabe award 2006: antiplasmodial and antitumor activity of artemisinin--from bench to bedside," Planta Medica, vol. 73, no. 4, pp. 299-309, 2007
 WHO, WHO Monograph on Good Agricultural and Collection Practices (GACP) for Artemisia annua L., WHO, Geneva, Switzerland, 2006.
 R. Rai, S. Pandey, and S. P. Rai, "Arsenic-induced changes in morphological, physiological, and biochemical attributes and artemisinin biosynthesis in Artemisia annua, an antimalarial plant," Ecotoxicology, vol. 20, no. 8, pp. 1900-1913, 2011.
 R. D. Tripathi, P. Tripathi, S. Dwivedi et al., "Arsenomics: omics of arsenic metabolism in plants," Frontiers in Physiology, vol. 3, article 275, 2012.
 K. Kosova, P. Vitamvas, I. T. Prasil, and J. Renaut, "Plant proteome changes under abiotic stress--contribution of proteomics studies to understanding plant stress response," Journal ofProteomics, vol. 74, no. 8, pp. 1301-1322, 2011.
 S. Bohler, M. Bagard, M. Oufir et al., "A DIGE analysis of developing poplar leaves subjected to ozone reveals major changes in carbon metabolism," Proteomics, vol. 7, no. 10, pp. 1584-1599, 2007.
 X. Zang and S. Komatsu, "A proteomics approach for identifying osmotic-stress-related proteins in rice," Phytochemistry, vol. 68, no. 4, pp. 426-437, 2007.
 F. Han, H. Chen, X. J. Li, M. F. Yang, G. S. Liu, and S. H. Shen, "A comparative proteomic analysis of rice seedlings under various high-temperature stresses," Biochimica et Biophysica Acta--Proteins and Proteomics, vol. 1794, no. 11, pp. 1625-1634, 2009.
 P. Kieffer, J. Dommes, L. Hoffmann, J. F. Hausman, and J. Renaut, "Quantitative changes in protein expression of cadmium-exposed poplar plants," Proteomics, vol. 8, no. 12, pp. 2514-2530, 2008.
 R. Requejo and M. Tena, "Proteome analysis of maize roots reveals that oxidative stress is a main contributing factor to plant arsenic toxicity," Phytochemistry, vol. 66, no. 13, pp. 1519-1528, 2005.
 R. Requejo and M. Tena, "Maize response to acute arsenic toxicity as revealed by proteome analysis of plant shoots," Proteomics, vol. 6, pp. 156-162, 2006.
 N. Ahsan, D. G. Lee, I. Alam et al., "Comparative proteomic study of arsenic-induced differentially expressed proteins in rice roots reveals glutathione plays a central role during As stress," Proteomics, vol. 8, no. 17, pp. 3561-3576, 2008.
 N. Ahsan, D. G. Lee, K. H. Kim et al., "Analysis of arsenic stress-induced differentially expressed proteins in rice leaves by two-dimensional gel electrophoresis coupled with mass spectrometry," Chemosphere, vol. 78, no. 3, pp. 224-231, 2010.
 I. Duquesnoy, P. Goupil, I. Nadaud, G. Branlard, A. PiquetPissaloux, and G. Ledoigt, "Identification of Agrostis tenuis leaf proteins in response to As(V) and As(III) induced stress using a proteomics approach," Plant Science, vol. 176, no. 2, pp. 206-213, 2009.
 E. Bona, C. Cattaneo, P. Cesaro et al., "Proteomic analysis of Pteris vittata fronds: two arbuscular mycorrhizal fungi differentially modulate protein expression under arsenic contamination," Proteomics, vol. 10, no. 21, pp. 3811-3834, 2010.
 E. Bona, F. Marsano, N. Massa et al., "Proteomic analysis as a tool for investigating arsenic stress in Pteris vittata roots colonized or not by arbuscular mycorrhizal symbiosis," Journal of Proteomics, vol. 74, no. 8, pp. 1338-1350, 2011.
 N. Stoeva and T. Bineva, "Oxidative changes and photosynthesis in oat plants grown in As contaminated soil," Bulgarian Journal of Plant Physiology, vol. 29, pp. 87-95, 2003.
 E. Miteva and M. Merakchiyska, "Response of chloroplasts and photosynthetic mechanism of bean plants to excess arsenic in soil," Bulgarian Journal of Agricultural Science, vol. 8, pp. 151156, 2002.
 D. R. Hoagland and D. I. Arnon, "The water-culture method for growing plants without soil," California Agricultural Experiment Station, vol. 347, pp. 1-3,1950.
 S. Ahamed, M. Kumar Sengupta, A. Mukherjee et al., "Arsenic groundwater contamination and its health effects in the state of Uttar Pradesh (UP) in upper and middle Ganga plain, India: a severe danger," Science of the Total Environment, vol. 370, no. 2-3, pp. 310-322, 2006.
 T. Roychowdhury, T. Uchino, H. Tokunaga, and M. Ando, "Survey of arsenic in food composites from an arsenic-affected area of West Bengal, India," Food and Chemical Toxicology, vol. 40, no. 11, pp. 1611-1621, 2002.
 C. H. Fiske and Y. Subbarow, "The colorimetric determination of phosphorus," The Journal of Biological Chemistry, vol. 66, pp. 375-400, 1925.
 C. Larsson and T. Olssonl, "Firefly assay of adenine nucleotides from algae: comparison of extraction methods," Plant and Cell Physiology, vol. 20, no. 1, pp. 145-155,1979.
 D. A. Smyth and W. M. Dugger, "Cellular changes during boron-deficient culture of the diatom Cylindrotheca fusiformis" Physiologia Plantarum, vol. 51, no. 1, pp. 111-117,1981.
 S. R. Mishra and S. C. Sabat, "Effect of calcium ion on Hydrilla verticillata thylakoid membrane [O.sub.2] evolution," Indian Journal of Biochemistry and Biophysics, vol. 32, no. 2, pp. 94-99,1995.
 N. Atal, P. P. Saradhi, and P. Mohanty, "Inhibition of the chloroplast photochemical reactions by treatment of wheat seedlings with low concentrations of cadmium: analysis of electron transport activities and changes in fluorescence yield," Plant and Cell Physiology, vol. 32, no. 7, pp. 943-951, 1991.
 D. I. Arnon, "Copper enzymes in isolated chloroplast: polyphenol oxidase in beet (Beta vulgaris)" Plant Physiology, vol. 24, no. 1, pp. 1-15, 1949.
 S. S. Zhao and M. Y. Zeng, "Spectrometric high-pressure liquid chromatographic studies on analytical chemistry of Qinghaosu," Planta Medica, vol. 3, pp. 233-237, 1985.
 H. Lange, W Shropshire, and H. Mohr, "An analysis of phytochrome mediated anthocyanin synthesis," Plant Physiology, vol. 47, no. 5, pp. 649-655, 1971.
 A. Bette and U. Kutschera, "Pigment accumulation and photosynthesis in developing rye coleoptiles," Botanica Acta, vol. 109, no. 3, pp. 194-198, 1996.
 D. A. Sims and J. A. Gamon, "Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages," Remote Sensing of Environment, vol. 81, no. 2-3, pp. 337-354, 2002.
 M. M. Bradford, "A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein dye binding," Analytical Biochemistry, vol. 72, no. 1-2, pp. 248-254, 1976.
 S. Pandey, R. Rai, and L. C. Rai, "Proteomics combines morphological, physiological and biochemical attributes to unravel the survival strategy of Anabaena sp. PCC7120 under arsenic stress," Journal of Proteomics, vol. 75, no. 3, pp. 921-937, 2012.
 A. Shevchenko, M. Wilm, O. Vorm, and M. Mann, "Mass spectrometric sequencing of proteins from silver-stained polyacrylamide gels," Analytical Chemistry, vol. 68, no. 5, pp. 850858, 1996.
 Y. Mishra, P. Bhargava, N. Chaurasia, and L. C. Rai, "Proteomic evaluation of the non-survival of Anabaena doliolum (Cyanophyta) at elevated temperatures," European Journal of Phycology, vol. 44, no. 4, pp. 551-565, 2009.
 A. A. Sharov, D. B. Dudekula, and M. S. H. Ko, "A web-based tool for principal component and significance analysis of microarray data," Bioinformatics, vol. 21, no. 10, pp. 2548-2549, 2005.
 O. Bjorkman and B. Demmig, "Photon yield of [O.sub.2] evolution and chlorophyll fluorescence characteristics at 77 K among vascular plants of diverse origins," Planta, vol. 170, no. 4, pp. 489-504, 1987.
 K. Maxwell and G. N. Johnson, "Chlorophyll fluorescence--a practical guide," Journal of Experimental Botany, vol. 51, no. 345, pp. 659-668, 2000.
 Z. Krupa and M. Moniak, "The stage of leaf maturity implicates the response of the photosynthetic apparatus to cadmium toxicity," Plant Science, vol. 138, no. 2, pp. 149-156, 1998.
 S. P. Yan, Q. Y. Zhang, Z. C. Tang, W. A. Su, and W. N. Sun, "Comparative proteomic analysis provides new insights into chilling stress responses in rice," Molecular and Cellular Proteomics, vol. 5, no. 3, pp. 484-496, 2006.
 Q. Jiang, H. Chen, X. Pan et al., "Proteomic analysis of wheat (Triticum aestivum L.) hybrid necrosis," Plant Science, vol. 175, no. 3, pp. 394-401, 2008.
 J. M. Abercrombie, M. D. Halfhill, P. Ranjan et al., "Transcriptional responses of Arabidopsis thaliana plants to As (V) stress," BMC Plant Biology, vol. 8, article 87, 2008.
 M. Bevan, I. Bancroft, E. Bent et al., "Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana" Nature, vol. 391, no. 6666, pp. 485-488, 1998.
 C. Xu, W. M. Garrett, J. Sullivan, and T. J. Caperna, "Separation and identification of soybean leaf proteins by two-dimensional gel electrophoresis and mass spectrometry," Phytochemistry, vol. 67, no. 22, pp. 2431-2340, 2006.
 Q. Meng, C. Zhang, F. Huang, J. Gai, and D. Yu, "Molecular cloning and characterization of a LEAFY-like gene highly expressed in developing soybean seeds," Seed Science Research, vol. 17, no. 4, pp. 297-302, 2007
 P. Han, B. Garcia-Ponce, G. Fonseca-Salazar, E. R. AlvarezBuylla, and H. Yu, "AGAMOUS-LIKE 17, a novel flowering promoter, acts in a FT-independent photoperiod pathway," Plant Journal, vol. 55, no. 2, pp. 253-265, 2008.
 P. Catarecha, M. D. Segura, J. M. Franco-Zorrilla et al., "A mutant of the Arabidopsis phosphate transporter PHT1;1 displays enhanced arsenic accumulation," Plant Cell, vol. 19, no. 3, pp. 1123-1133, 2007.
 B. Dhir, P. Sharmila, and P. P. Saradhi, "Photosynthetic performance of Salvinia natans exposed to chromium and zinc rich wastewater," Brazilian Journal of Plant Physiology, vol. 20, no. 1, pp. 61-70, 2008.
 J. E. Backhausen, C. Kitzmann, P. Horton, and R. Scheibe, "Electron acceptors in isolated intact spinach chloroplasts act hierarchically to prevent over-reduction and competition for electrons," Photosynthesis Research, vol. 64, no. 1, pp. 1-13,2000.
 J. K. Hurley, R. Morales, M. Martinez-Julvez et al., "Structure-function relationships in Anabaena ferredoxin/ferredoxin:NADP+ reductase electron transfer: insights from site-directed mutagenesis, transient absorption spectroscopy and X-ray crystallography," Biochimica et Biophysica Acta--Bioenergetics, vol. 1554, no. 1-2, pp. 5-21, 2002.
 F. Han, H. Chen, X. J. Li, M. F. Yang, G. S. Liu, and S. H. Shen, "A comparative proteomic analysis of rice seedlings under various high-temperature stresses," Biochimica et Biophysica Acta--Proteins and Proteomics, vol. 1794, no. 11, pp. 1625-1634, 2009.
 X. Y. Wan and J. Y. Liu, "Comparative proteomics analysis reveals an intimate protein network provoked by hydrogen peroxide stress in rice seedling leaves," Molecular and Cellular Proteomics, vol. 7, no. 8, pp. 1469-1488, 2008.
 A. E. Atik, G. O. Bozdag, E. Akinci et al., "Proteomic changes during boron tolerance in barley (Hordeum vulgare) and the role of vacuolar proton-translocating ATPase subunit E," Turkish Journal of Botany, vol. 35, no. 4, pp. 379-388, 2011.
 M. Burzynski and G. Klobus, "Changes of photosynthetic parameters in cucumber leaves under Cu, Cd, and Pb stress," Photosynthetica, vol. 42, no. 4, pp. 505-510, 2004.
 A. Siedlecka and Z. Krupa, "Rubisco activity maintenance in environmental stress conditions-how many strategies," Cellular and Molecular Biology Letters, vol. 9, pp. 56-57, 2004.
 M. T. Dominguez, T Maradon, J. M. Murillo, and S. Redondo-Gomez, "Response of Holm oak (Quercus ilex subsp. ballota) and mastic shrub (Pistacia lentiscus L.) seedlings to high concentrations of Cd and Tl in the rhizosphere," Chemosphere, vol. 83, no. 8, pp. 1166-1174, 2011.
 E. Mateos-Naranjo, L. Andrades-Moreno, and S. Redondo-Gomez, "Tolerance to and accumulation of arsenic in the cordgrass Spartina densiflora Brongn," Bioresource Technology, vol. 104, pp. 187-194, 2012.
 R. Mahmud, N. Inoue, S. Y. Kasajima et al., "Response of common buckwheat and castor oil plant against different levels of soil arsenic concentration: a comparative study," Fagopyrum, vol. 23, pp. 45-51, 2006.
 S. Henkes, U. Sonnewald, R. Badur, R. Flachmann, and M. Stitt, "A small decrease of plastid transketolase activity in antisense tobacco transformants has dramatic effects on photosynthesis and phenylpropanoid metabolism," Plant Cell, vol. 13, no. 3, pp. 535-551, 2001.
 M. Rohmer, M. Seemann, S. Horbach, S. Bringer-Meyer, and H. Sahm, "Glyceraldehyde 3-phosphate and pyruvate as precursors of isoprenic units in an alternative non-mevalonate pathway for terpenoid biosynthesis," Journal of the American Chemical Society, vol. 118, no. 11, pp. 2564-2566, 1996.
 F. J. Ferreira, C. Guo, and J. R. Coleman, "Reduction of plastid-localized carbonic anhydrase activity results in reduced Arabidopsis seedling survivorship," Plant Physiology, vol. 147, no. 2, pp. 585-594, 2008.
 L. E. Fridlyand and V. L. Kaler, "Possible C[O.sub.2] concentrating mechanism in chloroplasts of C3 plants. Role of carbonic anhydrase," General Physiology and Biophysics, vol. 6, no. 6, pp. 617-636, 1987.
 P. Kieffer, S. Planchon, M. Oufir et al., "Combining proteomics and metabolite analyses to unravel cadmium stress-response in poplar leaves," Journal of Proteome Research, vol. 8, no. 1, pp. 400-417, 2009.
 M. A. Pillai, Z. Lihuang, and T. Akiyama, "Molecular cloning, characterization, expression and chromosomal location of OsGAPDH, a submergence responsive gene in rice (Oryza sativa L.)," Theoretical and Applied Genetics, vol. 105, no. 1, pp. 34-42, 2002.
 C. Tu and L. Q. Ma, "Effects of arsenic on concentration and distribution of nutrients in the fronds of the arsenic hyperaccumulator Pteris vittata L," Environmental Pollution, vol. 135, no. 2, pp. 333-340, 2005.
 Z. Wang, Y. Xiao, W. Chen, K. Tang, and L. Zhang, "Increased vitamin C content accompanied by an enhanced recycling pathway confers oxidative stress tolerance in Arabidopsis," Journal of Integrative Plant Biology, vol. 52, no. 4, pp. 400-409, 2010.
 Y. Nakano and K. Asada, "Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts," Plant and Cell Physiology, vol. 22, no. 5, pp. 867-880, 1981.
 K. Lee, D. W. Bae, S. H. Kim et al., "Comparative proteomic analysis of the short-term responses of rice roots and leaves to cadmium," Journal of Plant Physiology, vol. 167, no. 3, pp. 161168, 2010.
 R. Rakwal, G. K. Agrawal, and M. Yonekura, "Separation of proteins from stressed rice (Oryza sativa L.) leaf tissues by two-dimensional polyacrylamide gel electrophoresis: induction of pathogenesis-related and cellular protectant proteins by jasmonic acid, UV irradiation and copper chloride," Electrophoresis, vol. 20, no. 17, pp. 3472-3478, 1999.
 A. Mangeon, R. M. Junqueira, and G. Sachetto-Martins, "Functional diversity of the plant glycine-rich proteins superfamily," Plant Signaling and Behavior, vol. 5, no. 2, pp. 99-104, 2010.
 I. S. Minas, G. Tanou, M. Belghazi et al., "Physiological and proteomic approaches to address the active role of ozone in kiwifruit post-harvest ripening," Journal of Experimental Botany, vol. 637, pp. 2449-2464, 2012.
 N. S. Scrutton, A. Berry, and R. N. Perham, "Redesign of the coenzyme specificity of a dehydrogenase by protein engineering," Nature, vol. 343, no. 6253, pp. 38-43, 1990.
 Y. Zhang, Y. F. Ma, S. W. Qi et al., "Responses to arsenate stress by Comamonas sp. strain CNB-1 at genetic and proteomic levels," Microbiology, vol. 153, no. 11, pp. 3713-3721, 2007
 G. O. Noriega, K. B. Balestrasse, A. Batlle, and M. L. Tomaro, "Cadmium induced oxidative stress in soybean plants also by the accumulation of d-aminolevulinic acid," BioMetals, vol. 20, no. 6, pp. 841-851, 2007.
Rashmi Rai, (1) Sarita Pandey, (2) Alok Kumar Shrivastava, (2) and Shashi Pandey Rai (1)
(1) Laboratory of Morphogenesis, Center of Advanced Study in Botany, Banaras Hindu University, Varanasi 221005, India
(2) Laboratory of Algal Biology, Molecular Biology Section, Center of Advanced Study in Botany, Banaras Hindu University, Varanasi 221005, India
Correspondence should be addressed to Shashi Pandey Rai; email@example.com
Received 19 October 2013; Accepted 3 February 2014; Published 29 April 2014
Academic Editor: David Sheehan
TABLE 1: Effect of As on various morphological and physiological parameters including, plant height, flowering time, number of inflorescences-plant, number of branches-inflorescence, number of capitula-inflorescence, capitula size, florets distance in whorl, number of florets-capitulum, number of oil glands, and artemisinin and anthocyanin contents in A. annua plants. Measurements were conducted under natural conditions when plants were 8-month old. Characters Control 100 [micro]M Mature plant 1.4 [+ or -] 0.04 (a) 1.2 [+ or -] 0.02 (a) height (m) Time of first 30-31 weeks 28-29 weeks flower Number of 28 [+ or -] 1.02 (a) 27 [+ or -] 1.34 (ab) inflorescence/ plant Number of 23 [+ or -] 1.1 (a) 21 [+ or -] 1.1 (a) branches/ inflorescence Number of 20 [+ or -] 1.02 (b) 32 [+ or -] 0.98 (a) capitula/ inflorescence Capitula size (mm) 4.3 [+ or -] 0.95 (c) 5.1 [+ or -] 0.76 (b) Distance of florets 1.1 [+ or -] 0.004 (a) 1.0 [+ or -] 0.01 (a) in involucre whorl (mm) Florets/capitulum Ray = 13 Disc = 35 Ray = 16 Disc = 45 Artemisinin content 0.39 [+ or -] 0.02 (a) 0.46 [+ or -] 1.02 (a) (mg/100 g DW) Number of oil 5 [+ or -] 0.05 (c) 7 [+ or -] 0.09 (b) glands/florets Anthocyanin content 0.06 [+ or -] 0.01 (a) 0.07 [+ or -] 0.06 (a) ([micro]mol [g.sup.-1] FW) Biomass (DW g 241 [+ or -] 4.5 (b) 267.0 [+ or -] 5.2 (a) [plant.sup.-1]) Characters 150 [micro]M Mature plant 1.1 [+ or -] 0.09 (a) height (m) Time of first 22-24 weeks flower Number of 24 [+ or -] 0.98 (c) inflorescence/ plant Number of 18 [+ or -] 1.4 (b) branches/ inflorescence Number of 33 [+ or -] 1.04 (a) capitula/ inflorescence Capitula size (mm) 5.4 [+ or -] 0.01 (a) Distance of florets 0.8 [+ or -] 0.003 (b) in involucre whorl (mm) Florets/capitulum Ray = 17 Disc = 55 Artemisinin content 0.52 [+ or -] 0.02 (a) (mg/100 g DW) Number of oil 9.9 [+ or -] 0.01 (a) glands/florets Anthocyanin content 0.12 [+ or -] 0.09 (a) ([micro]mol [g.sup.-1] FW) Biomass (DW g 221 [+ or -] 5.0 (c) [plant.sup.-1]) All values are mean [+ or -] SD (n = 3). Different letters show significantly different values (P < 0.05, DMRT). TABLE 2: Photosynthetic rate ([P.sub.s]), stomatal conductance ([g.sub.s]), intercellular C[O.sub.2] ([C.sub.i]), transpiration rate (E), water use efficiency (WUE), [F.sub.i]-[F.sub.m] (the maximal photochemical efficiency), photosystem I activity (PS- I), photosystem II activity (PS-II), whole chain electron transport chain, ATP and NADPH contents of control, and As treated A. annua plants. Measurements were conducted when plants were 8-month old. Physiological parameters Control [P.sub.s] ([micro]mol 31.7 [+ or -] 0.43 (a) C[O.sub.2] [m.sup.-2] [s.sup.-1]) [g.sub.s] (mol 0.44 [+ or -] 0.04 (a) [H.sub.2]O[m.sup.-2] [s.sup.-1]) [C.sub.i] (ppm) 201.0 [+ or -] 8.9 (a) E (mmol [H.sub.2]O 6.1 [+ or -] 0.45 (a) [m.sup.-2] [s.sup.-1]) WUE ([micro]mol 3.2 [+ or -] 0.02 (b) [mol.sup.-1]) [F.sub.v]/[F.sub.m] 0.82 [+ or -] 0.01 (a) PS-I activity 250 [+ or -] 9.01 (c) ([micro]mol [O.sub.2] consumed [mg.sup.-1] Chl [h.sup.-1]) PS-II activity 225 [+ or -] 11.0 (a) ([micro]mol [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) Whole chain ([micro]mol 148 [+ or -] 12 (a) [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) ATP ([micro]mol 0.37 [+ or -] 0.07 (a) [g.sup.-1] FW) NADPH (nmol 5.17 [+ or -] 0.61 (b) [g.sup.-1] FW) Physiological parameters 100 [micro]M [P.sub.s] ([micro]mol 30.0 [+ or -] 0.78 (b) C[O.sub.2] [m.sup.-2] [s.sup.-1]) [g.sub.s] (mol 0.36 [+ or -] 0.07 (a) [H.sub.2]O[m.sup.-2] [s.sup.-1]) [C.sub.i] (ppm) 189.0 [+ or -] 9.7 (a) E (mmol [H.sub.2]O 6.08 [+ or -] 0.67 (a) [m.sup.-2] [s.sup.-1]) WUE ([micro]mol 5.1 [+ or -] 0.82 (a) [mol.sup.-1]) [F.sub.v]/[F.sub.m] 0.79 [+ or -] 0.02 (a) PS-I activity 375 [+ or -] 11.0 (a) ([micro]mol [O.sub.2] consumed [mg.sup.-1] Chl [h.sup.-1]) PS-II activity 199 [+ or -] 12.0 (b) ([micro]mol [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) Whole chain ([micro]mol 149 [+ or -] 22.3 (a) [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) ATP ([micro]mol 0.40 [+ or -] 0.12 (a) [g.sup.-1] FW) NADPH (nmol 7.12 [+ or -] 0.82 (a) [g.sup.-1] FW) Physiological parameters 150 [micro]M [P.sub.s] ([micro]mol 23.5 [+ or -] 0.99 (c) C[O.sub.2] [m.sup.-2] [s.sup.-1]) [g.sub.s] (mol 0.34 [+ or -] 0.08 (a) [H.sub.2]O[m.sup.-2] [s.sup.-1]) [C.sub.i] (ppm) 212.0 [+ or -] 5.7 (a) E (mmol [H.sub.2]O 6.02 [+ or -] 0.54 (a) [m.sup.-2] [s.sup.-1]) WUE ([micro]mol 4.9 [+ or -] 1.1 (a) [mol.sup.-1]) [F.sub.v]/[F.sub.m] 0.72 [+ or -] 0.04 (b) PS-I activity 279 [+ or -] 10.0 (b) ([micro]mol [O.sub.2] consumed [mg.sup.-1] Chl [h.sup.-1]) PS-II activity 154 [+ or -] 12.0 (c) ([micro]mol [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) Whole chain ([micro]mol 109 [+ or -] 11 (b) [O.sub.2] evolved [mg.sup.-1] Chl [h.sup.-1]) ATP ([micro]mol 0.29 [+ or -] 0.07 (a) [g.sup.-1] FW) NADPH (nmol 3.91 [+ or -] 0.71 (b) [g.sup.-1] FW) All values are mean [+ or -] SD (n = 6). Different letters show significantly different values (P < 0.05, DMRT). TABLE 3: Arsenic-induced differentially expressed proteins in Artemisia annua identified by MALDI-MS/MS analysis. Spot No. Protein name Plant species/ Mr/pZ accession Theoretical number (a) Energy metabolism (primary carbon metabolism) 1 Transketolase Spinacia 80.7/6.2 oleracea/ gi|68052991 6 Glyceraldehyde- Mikaniamicrantha/ 37.0/8.5 3-phosphate gi|269856436 dehydrogenase 37 Glyceraldehyde- Ricinus 49.1/7.5 3-phosphate communis/ dehydrogenase, gi|255539282 putative 23 Unknown Glycine mnx/ 48.8/6.7 gi|255641007 22 Ribulose-1,5- Glycine mnx/ 48.8/6.2 bisphosphate gi|358249078 carboxylase activase 36 Plastidic Solanum 38.6/5.8 aldolases tuberosum/ g|1781348 12 Plastidic Solanum 38.6/5.8 aldolases tuberosum/ g|1781348 17 Carbonic Flaveria 36.0./5.8 anhydrase bidentis/ gi|1089983 27 Carbonic Flaveria 36.0/5.8 anhydrase 3 bidentisi/ gi|1089983 40 Carbonic Dimocarpus 35.2/6.5 anhydrase ATP longan/ dependent gi|339958979 30 Chloroplast Ricinus 30.9/8.7 protease communis/ putative gi|255537123 34 Fuctose-bis- Spinacia 42.7/7.5 phospahte oleracea/ aldolase, gi|22633 39 Malate Solanum 36.4/8.8 dehydrogenase tuberosum/ gi|21388544 41 Ribose-5- Spinacia 30.9/6.5 phosphate- oleracea/ isomerase gi|18654317 44 Triosephosphate Medicago 27.7/6.5 isomerase tranculata/ gi|1351282 42 Triosephosphate Medicago 27.7/6.5 isomerase tranculata/ gi|1351282 4 PS II PsbP Pisum sativum/ 28.2/8.2 protein (Oygen gi|131390 Evolving enhancer 2) 3 Photosynthetic Panax ginseng/ 19.9/5.7 electron gi|89475526 transfer-like protein 5 Unnamed protein Spinacia 31.9/7.6 product oleracea/ gi|21283 8 Ferredoxin Solanum 5.4/77 NADP+ peruvianum/ oxidoreductase gi|61969078 38 Ferredoxin Solanum 35.4/7.7 NADP+ peruvianum/ oxidoreductase gi|61969078 2 Chlorophyll a, Solanum 29.3/8.9 b, binding lycopersicum/ protein gi|115813 11 Chloroplast PSI Helianthus 13.8/4.5 Type III annuus/ gi|159138839 10 Chlorophyll a, Pisum sativum/ 29.3/5.2 b, binding gi|20671 protein 20 Unknown protein Medicago 31.2/6.5 tranculata/ gi|217071344 21 ATP synthase Viburnum 22.7/4.3 [beta] chain 0opulus/ gi|7688419 24 Rieske-FeS Arabidopsis 24.6/8.8 protein thaliana/ gi|9843639 Energy metabolism (photorespiration) 14 P protein Flaveria 11.4/6.5 pringlei/ gi|438003 ROS scavenging and defence 18 Dehydro Lotus 29.1/7.7 ascorbate japonicas/ reductase gi|66732627 19 Ascorbate Tragopogan 13.1/4.5 peroxidase-2- dubius/ like protein gi|290796648 Transcriptional Regulator 25 SGRP glycine Daucus cnrota/ 15.7/5.3 rich binding gi|544426 protein 26 Maturase K Prinsepia 32.2/9.0 uniflora/ gi|290583766 28 LEAFY like Mimulus 9.1/9.6 protein guttatus/ gi|14573451 Protein metabolism 13 Ribosomal L12 Nicotiana 20.3/6.04 la tabacum/ gi|20020 16 NAD(P)-rossaman Arabidopsis 34.9/8.37 binding protein thaliana/ gi|18404496 29 Chloroplast Ipomoea nil/ 74.4/5.1 heat shock gi|166919370 protein 70-1 Secondary metabolism 9 Chalcone Leibnitzia 25.4/5.4 synthase anandria/ gi|1403057 31 S-adenosyl Arabidopsis 42.4/5.7 methionine thaliana/ synthase gi|15228048 Transport protein 35 ABC transporter Arabidopsis 66.6/9.6 like protein thaliana/ gi|224090097 Unknown proteins 32 Unknown protein Populus 83.0/5.9 trichorpa/ gi|224130876 33 Hypothetical Vitis vinifera/ 60.5/8.68 protein gi|147845220 43 Unknown protein Zea mnys/ 17.0/7.8 gi|302844664 45 Hypothetical Sorghum 46.2/6.5 protein bicolor/ gi|242061954 46 Unknown protein Vitis vinifera/ 47.9/9.4 gi|296081316 Spot No. Mr/pI Score (c) Matched Experimental peptides (b) Protein Peptide (d) Energy metabolism (primary carbon metabolism) 1 66.0/5.9 87 50 2 6 32.0/6.8 96 63 2 37 36.0/6.0 257 132 4 23 29.0/4.1 257 136 4 22 40.0/4.9 418 122 4 36 32.0/5.8 303 113 5 12 20.0/5.5 279 103 3 17 26.0/6.9 128 87 2 27 29.0/6.5 209 111 4 40 24.0/6.4 128 87 2 30 58.0/5.4 90 90 1 34 27.0/4.5 93 39 3 39 34.0/6.4 248 157 3 41 25.0/6.1 69 69 1 44 25.0/4.9 227 130 3 42 26.0/5.6 129 89 2 4 22.0/6.4 138 54 6 3 20.0/6.1 132 94 2 5 28.0/5.2 252 76 3 8 22.0/6.8 247 112 4 38 28.0/6.0 358 146 5 2 23.0/5.9 63 63 1 11 24.0/5.8 113 71 2 10 23.0/5.9 86 44 2 20 38.0/4.3 149 54 3 21 45.0/5.2 55 16 2 24 15.0/5.1 49 49 1 Energy metabolism (photorespiration) 14 8.0/4.5 79 41 2 ROS scavenging and defence 18 26.0/6.0 83 83 1 19 27.0/5.3 91 91 1 Transcriptional Regulator 25 15.0/5.3 157 107 3 26 26.0/5.8 56 47 1 28 25.0/4.8 20 20 1 Protein metabolism 13 20.0/4.6 102 102 1 16 15.0/6.3 114 61 2 29 63.0/4.8 308 124 4 Secondary metabolism 9 25.0/6.5 50 50 1 31 38.0/5.8 281 93 5 Transport protein 35 29.0/5.5 27 27 1 Unknown proteins 32 97.0/6.2 26 26 1 33 55.0/5.8 38 38 1 43 24.0/5.4 30 30 1 45 22.0/4.7 43 29 3 46 21.0/4.3 33 33 1 Spot No. Coverage Cellular Blast P (%) (e) location results (f) Energy metabolism (primary carbon metabolism) 1 12 Chloroplast -- 6 9 Cytosol 37 9 Cytosol -- 23 11 Cytosol Glyceraldehyde- 3-phosphate dehydrogenase 22 13 Chloroplast -- 36 28 Chloroplast -- 12 12 Chloroplast 17 13 Chloroplast -- 27 9 Chloroplast -- 40 10 Chloroplast -- 30 5 Chloroplast -- 34 10 Chloroplast -- 39 16 Mitochondria -- 41 10 Chloroplast -- 44 8 Cytosolic -- 42 11 Cytosolic -- 4 9 Chloroplast -- 3 14 Chloroplast PSII-O, 5 20 Chloroplast photosystem II, oxygen evolving enhancer 8 7 Chloroplast -- 38 16 Chloroplast -- 2 10 Chloroplast -- 11 19 Chloroplast -- 10 14 Chloroplast -- 20 10 Chloroplast Chlorophyll a, b, binding protein 21 76 Chloroplast -- 24 5 Chloroplast -- Energy metabolism (photorespiration) 14 10 Chloroplast -- ROS scavenging and defence 18 4 Cytosol -- 19 20 Cytosol -- Transcriptional Regulator 25 25 Cytoplasm -- 26 13 Cytoplasm -- 28 16 Cytoplasm -- Protein metabolism 13 7 Cytoplasm -- 16 10 Cytoplasm -- 29 10 Cytoplasm -- Secondary metabolism 9 8 Cytoplasm -- 31 16 Cytoplasm -- Transport protein 35 10 Cytoplasm Unknown proteins 32 2 -- Unknown protein 33 5 -- Hypothetical protein 43 4 -- Unknown protein 45 1 -- Hypothetical protein 46 2 -- Unknown protein Spot No. Fold change (g) Energy metabolism (primary carbon metabolism) 1 1.89 (b) -2.10 (a) 6 1.85 (a) 1.87 (a) 37 1.60 (b) 2.04 (a) 23 4.51 (a) 1.70 (b) 22 1.55 (a) -1.72 (a) 36 1.50 (a) --1.08 (b) 12 2.69 (a) 1.96 (a) 17 3.10 (a) 2.00 (b) 27 3.22 (a) 2.86 (b) 40 4.4 (a) 2.3 (b) 30 -3.02 (a) -3.63 (a) 34 2.12 (a) 1.04 (b) 39 -0.75 (b) 2.12 (a) 41 1.72 (a) 1.54 (a) 44 1.29 (a) 1.50 (a) 42 1.49 (a) 1.50 (a) 4 2.24 (a) 2.46 (a) 3 -1.04 (b) 1.39 (a) 5 1.39 (a) 1.49 (a) 8 2.41 (a) 1.50 (b) 38 1.70 (a) 1.18 (b) 2 1.93 (b) -1.50 (a) 11 2.83 (b) 2.36 (a) 10 1.99 (b) 2.17 (a) 20 3.12 (a) 1.35 (b) 21 2.68 (a) -0.56 (b) 24 2.35 (a) -1.06 (b) Energy metabolism (photorespiration) 14 2.24 (a) 2.46 (a) ROS scavenging and defence 18 3.07 (a) 1.98 (b) 19 3.44 (a) 2.67 (b) Transcriptional Regulator 25 2.83 (b) 3.00 (a) 26 2.90 (b) 3.16 (a) 28 2.32 (a) 2.38 (a) Protein metabolism 13 2.45 (b) -1.25 (a) 16 2.19 (a) 2.29 (a) 29 -2.43 (a) --1.7 9 (b) Secondary metabolism 9 2.24 (b) 2.30 (a) 31 1.61 (a) 1.41 (a) Transport protein 35 1.50 (a) 1.28 (a) Unknown proteins 32 -1.9 (a) -1.7 (a) 33 1.5 (b) 1.9 (a) 43 1.99 (a) 1.67 (b) 45 1.67 (b) 1.80 (a) 46 (a) Theoretical mass (kDa) and pI of identified proteins were calculated consulting databases NCBI nr and Swissprot. (b) Experimental mass (kDa) and pI of identified proteins. (c) Score: the protein score is derived from the ions scores from an MS/MS search. (d) Matched peptides derived from MS/MS data. (e) Coverage (%): it is the percent of the residues in each protein sequence that have been identified. (f) Blast P result similar to the result as obtained in NCBInr database is indicated by (-). (g) Fold change is calculated as a ratio of the averaged means of normalized spot volumes of control and the treatments.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Rai, Rashmi; Pandey, Sarita; Shrivastava, Alok Kumar; Rai, Shashi Pandey|
|Publication:||International Journal of Proteomics|
|Date:||Jan 1, 2014|
|Previous Article:||Dimerization of peptides by calcium ions: investigation of a calcium-binding motif.|
|Next Article:||The influence of flanking secondary structures on amino acid content and typical lengths of 3/10 helices.|