All relevant data are within the paper and its Supporting Information files.

Introduction {#sec001}
============

Artemisinin (AN), a sesquiterpene lactone produced mainly in *Artemisia annua* as a secondary metabolite, is a highly effective natural product against malaria and other parasitic diseases, as well as an anti-cancer agent. The World Health Organization (WHO) has recommended the Artemisinin-based Combination Therapy (ACT) for the treatment of malaria \[[@pone.0140266.ref001],[@pone.0140266.ref002]\]. Our group has previously reported artemisinin in 12 *Artemisia* species found in Pakistan \[[@pone.0140266.ref003]\] but at very low levels, i.e. 0.01--0.8% of the dry weight of the plant \[[@pone.0140266.ref004]\], reaching a maximum of 1.5% in some cases \[[@pone.0140266.ref005],[@pone.0140266.ref006]\]. There is therefore a need either to enhance the metabolite concentration in the natural plant or seek alternative sources of artemisinin.

Among different strategies used to improve secondary metabolite production, the recombinant DNA technology has made it possible to directly modify the expression of biosynthetic genes, and manipulate the pathways that lead to secondary plant compounds \[[@pone.0140266.ref007]\]. Various studies have shown *rol* genes to be powerful inducers of secondary metabolism in different plant families \[[@pone.0140266.ref008]\]. *Rol* A is a DNA-binding protein and stimulator of growth, while the tyrosine phosphatase activity of *rol B* regulates the signal transduction pathway of auxin \[[@pone.0140266.ref009]\]. *Rol B* has been used to increase the production of resveratrol in *Vitis amurensis* \[[@pone.0140266.ref010]\] and anthraquinones in *Rubia cardifolia* \[[@pone.0140266.ref011]\]. The *rol* C gene, which has cytokinin glucosidase activity, is capable of stimulating the production of many secondary compounds in transformed plants and cell cultures, such as tropane alkaloids, pyridine alkaloids, indole alkaloids, ginsenosides and anthraquinones \[[@pone.0140266.ref012],[@pone.0140266.ref013],[@pone.0140266.ref014]\]. Our previous work also showed that *Artemisia dubia*, when transformed with *rol ABC* genes, produces 10 times more artemisinin and its derivatives \[[@pone.0140266.ref015]\].

The aim of the present investigation, after the detection of artemisinin and its derivatives \[[@pone.0140266.ref016]\] ([Fig 1](#pone.0140266.g001){ref-type="fig"}) in *Artemisia carvifolia* Buch, was to enhance the levels of these metabolites by the expression of *rol* B and *rol* C genes. As the *A*. *carvifolia* Buch transgenics showed a significant increase in target compound levels, real time qPCR analysis of artemisinin biosynthetic genes including those encoding amorpha-4, 11 diene synthase (ADS), cytochrome P450 (CYP71AV1), and aldehyde dehydrogenase 1 (ALDH1) \[[@pone.0140266.ref017],[@pone.0140266.ref018],[@pone.0140266.ref019],[@pone.0140266.ref020]\] ([Fig 2](#pone.0140266.g002){ref-type="fig"}), was carried out to find a correlation between the level of gene expression and metabolite concentration. As artemisinin and other secondary metabolites are produced in 10-celled glandular trichomes located on leaves and other aerial parts \[[@pone.0140266.ref021],[@pone.0140266.ref022],[@pone.0140266.ref023]\] and sequestered in the epicuticular sac at the apex of the trichome \[[@pone.0140266.ref022]\], we also carried out RT-qPCR analysis of the TAFR1 (trichome-specific fatty acyl-CoA reductase 1) gene. TAFR1 is an enzyme involved in trichome development and sesquiterpenoid biosynthesis. We also calculated the trichome density to find a relationship between trichome development and artemisinin production.

![Artemisinin and derivatives.\
Structure of artemisinin and derivatives studied. Fig taken from literature \[[@pone.0140266.ref016]\].](pone.0140266.g001){#pone.0140266.g001}

![Artemisinin biosynthesis.\
A schematic diagram of artemisinin biosynthetic pathway.](pone.0140266.g002){#pone.0140266.g002}

Material and Methods {#sec002}
====================

Seed Germination and DNA Barcoding {#sec003}
----------------------------------

Seeds of *Artemisia carvifolia* were collected from Astore, in the Northern regions of Pakistan (35.3667° N, 74.8500° E; 8,500 ft elevation). No specific permissions were required for the mentioned location for collection of seeds. However, the project was approved by the Institutional Biosafety Committee (IBC) Quaid-i-Azam University Islamabad, Pakistan. After collection, seeds were surface sterilized with 70% ethanol and germinated on half-strength MS medium. Genomic DNA was extracted from the germinated plantlets according to our established lab protocol \[[@pone.0140266.ref024]\]. For the identification of *A*.*carvifolia* Buch, a non-coding spacer between the *psbA* and *trnH* genes of chloroplast DNA was amplified by PCR using primers of *psbA*: 5′ -GTTATGCATGAACGTAATGCTC-3′; *trnH*: 5\'-CGCGCATGGTGGATTCACAATC-3\'. The PCR reaction was carried out according to the reaction conditions reported previously \[[@pone.0140266.ref025]\]. Rapid PCR Purification System 9700 (Marligen Biosciences, Ijamsville, MD, USA) was used to purify the PCR product, which was then sequenced following the dideoxy-chain termination method using an ABI Prism 310 Automated DNA Sequencer (PE, Applied Biosystems, Foster City, CA, USA). Sequences were identified and aligned via the BioEdit sequence alignment tool (editor version 7.2.5.0).

Bacterial Strains and Plasmids {#sec004}
------------------------------

*Agrobacterium tumefaciens* strain GV3101 containing plasmids pPCV002-CaMVBT and pPCV002-CaMVC, kindly provided by Dr. A. Spena, Max-Planck-Institut fur Zuchtungsforschung, 5000 Koin 30, FRG \[[@pone.0140266.ref026]\], was used for transformation purposes. The coding sequence of the *rol* B and rol C genes was expressed in the T-DNA region of the plasmids pPCV002-CaMVBT and pPCV002-CaMVC, respectively, under the control of the CaM35S promoter. T-DNA of pPCV002-CaMVBT and pPCV002-CaMVC also contained the *neomycin phosphotransferase* (*NPTII*) gene with the nopaline synthase (NOS) promoter and NOS terminator sequences ([Fig 3](#pone.0140266.g003){ref-type="fig"}). *Agrobacterium tumefaciens* containing the plasmids pPCV002-CaMVBT and pPCV002-CaMVC were grown overnight in Luria broth (Sigma Cat \# L-1900). After inoculation, bacterial cultures were maintained at 28°C and 120 rpm in a shaking incubator. Growth was obtained in 24 hours and OD was checked by spectrophotometer, when in the range of 0.2--0.8, the bacterial suspension was used for the transformation.

![Detail of vectors used for transformation of *A*. *carvifolia*.\
Vectors (pPCV002-CaMVBT and pPCV002-CaMVC) used in the transformation of *Artemisia carvifolia*. Fig taken from literature \[[@pone.0140266.ref026]\].](pone.0140266.g003){#pone.0140266.g003}

Transformation and Regeneration {#sec005}
-------------------------------

One-month-old *in vitro* grown plants were used for transformation. Nodal, explants were prepared, and after 2 days of preculturing on shooting media (0.5 mg/l BAP and 0.1 mg/l NAA) supplemented with 200μM acetosyringone, they were infected with bacterial strains containing the desired constructs. After 10--15 minutes of bacterial infection, the explants were dried on autoclaved filter paper for 2--3 minutes and placed on MS shooting media supplemented with 200μM acetosyringone. After 2 days of incubation in darkness at 28°C, explants were washed with antibiotics and placed on selection medium (0.5mg/l BAP, 0.1mg/l NAA, 50mg/l kanamycin, 300mg/l cefotaxime). Regeneration occurred within one month, and explants were shifted to fresh selection media every 15 days. When shoots attained considerable length of 3-4cm, they were shifted to rooting media (NAA 0.1mg/l, kanamycin 50mg/l) to get complete plant with roots. After three to four selection cycles, the complete plants were regenerated on selection media.

Molecular Analysis {#sec006}
------------------

Molecular analysis was performed after extraction of genomic DNA from aerial parts of 2 month old transformed and wild type plants by the CTAB method \[[@pone.0140266.ref024]\]. The plasmid from GV3101 was isolated by the alkaline lysis method. PCR analysis was performed using a programmed DNA thermal cycler (Biometra, USA). The *rol* B gene forward 5'-GCTCTTGCAGTGCTAGATTT-3' and reverse primer 5'-GAAGGTGCAAGCTACCTCTC-3', the *rol* C gene forward 5'-GAAGACGACCTGTGTTCTC-3' and reverse primer 5'- CGTTCAAACGTTAGCCGATT-3' and the *nptII* gene forward 5'-AAGATGGATTGCACGCAGGTC-3' and reverse primer 5´GAAGAACTCGTCAAGAAGGCG-3´ were used for PCR analysis. Conditions applied for PCR were as described previously \[[@pone.0140266.ref015]\].

Southern Blotting {#sec007}
-----------------

Southern blot analysis was performed using the DIG High Prime DNA Labeling and Detection Starter Kit II (Roche Cat. No. 11585614910), following the manufacturer's instructions. For southern blotting, 1--5 μg of genomic DNA of *rol B* transformants was digested with X-baI and that of *rol C* gene transformants was digested with EcoRI and electrophoresed on 0.7% agarose gel overnight. The transfer of digested DNA to the positively charged nylon membranes was carried out according to the standard procedure \[[@pone.0140266.ref027]\]. The probe used for hybridization was the PCR products of *rol B* and *rol C* genes from plasmids, which was labeled using digoxigenin (DIG)-11-dUTP with DIG High Prime DNA Labeling reagents (Roche, Mannheim, Germany). Hybridization was done at 42--44°C and then immunological detection was carried out on X-ray film using CSPD substrate in accordance with the manufacturer's instructions.

Analysis of Artemisinin and Its Derivatives by HPLC-ESI/MS {#sec008}
----------------------------------------------------------

Extraction of artemisinin and its derivatives from shoots of 4 months old transformed and wild type plants was carried out using the reported methodology \[[@pone.0140266.ref028]\]. HPLC-MS analysis was performed with a Varian 500 Ion Trap Mass spectrometer (Varian, Spain) coupled to a 212-HPLC system equipped with an auto sampler 410 (Varian, Spain). The analytes were separated by a C18 (5μm) column (150mmx 4.6mm) (AKDAY Chromatografica, Spain) using a mobile phase consisting of water with 0.1% formic acid (A) and acetonitrile (B) in a gradient as follows: (t (min), % B) (0, 50) (10, 50) (15, 99) (23, 50) (26, 50). The flow rate was 1ml/min. Electro spray ionization (ESI) in positive mode was applied with a capillary voltage of 4000V and plate voltage of 600V. Other parameters include: spray chamber temperature 50 C, nebulizing gas (N~2~) pressure: 35psi, dry gas (N~2~) 10 l/min at 35psi and 400C. Mass analyzer scanned from 66--500 u. Fragmentation amplitude was 1.0 V. Three injections per sample were applied.

Expression Analysis of *rol* Genes and Artemisinin Biosynthetic Genes by Real Time Quantitative PCR {#sec009}
---------------------------------------------------------------------------------------------------

The young leaves of transformed and wild-type plants (4 months old) were selected for the RT-qPCR analysis. Total RNA was extracted from three replicates of each plant according to the procedure reported previously \[[@pone.0140266.ref029]\] and also by using the trizol reagent (Life Technologies). Turbo DNAse (Ambion) was used to ensure complete removal of DNA from RNA after extraction. Purity and quantity were checked by taking absorbance at 260 and 280nm on a Nanodrop ND-2000 spectrophotometer (Thermo scientific). The quality of RNA was also assessed by running RNA samples on 1.2% agarose gel. Reverse transcription of 1 μg of RNA was carried out using a 1^st^ strand cDNA synthesis kit (Invitrogen), following the manufacturer's instructions. In order to check the expression of *rol* B and *rol* C genes, a semi-quantitative reverse transcriptase-polymerase chain reaction was performed with *rol* B and *rol* C gene primers as previously, using 1 μl of cDNA reaction mixture as a template. The PCR reactions were run in triplicate. Gel images of PCR products were scanned by Kodak Molecular Imaging software (version 4.2) and integrated density values were found to be different for each band.

To evaluate the possible effects of *rol* genes on the expression of artemisinin biosynthetic genes, quantitative real time PCR of four selected genes was performed, namely those encoding amorpha-4, 11 diene synthase (ADS), cytochrome P450 (CYP71AV1) and aldehyde dehydrogenase 1 (ALDH1). The gene involved in trichome development and sesquiterpenoid biosynthesis (TFAR1) was also studied. The β-actin gene was used as reference gene \[[@pone.0140266.ref030]\]. For real time qPCR, a 1:4 dilution of cDNA was used. The amplification reaction was performed by gene specific primers i.e. ADS forward: 5'-ATTACTGGCGGTGCTAAC-3' and reverse: 5'-GTGCAGAGACAGCCCATT-3', CYP71AV1 forward: 5'-ATTTTGGATATGTTTGGAGCAGGC-3' and reverse: 5'-TCCGCTTGTACTTTCTCCATTGCT-3', ALDH1 forward: 5'-CAGGAGCTAATGGAAGTTCTAAGTCAG-3' and reverse: 5'-TTTCTTCCTTCGGCCACTGTTG-3', TAFR1 forward: 5'-CCTTGGAGATCCTGAAGCTG-3' and reverse: 5'-CGTTGGATTGTGCTGAACTG-3', β-actin forward: 5'-ATCAGCAATACCAGGGAACATAGT-3' and reverse: 5'-AGGTGCCCTGAGGTCTTGTTCC-3'. Primer amplification efficiency was determined as previously \[[@pone.0140266.ref031]\]. The qPCR was performed using iTAqTM Universal SYBR Green Supermix (BioRad, Hercules, CA, USA) in a 384-well platform system (ABI Prism® 7900HT sequence detection system, Applied Biosystems, Foster, CA, USA). The reaction conditions for real time qPCR were as follows: denaturation for 5 min at 95 ◦C, followed by 45 cycles each of denaturation for 10s at 95 ◦C, annealing for 10 s at 60 ◦C, followed by elongation for 10 s at 72 ◦C. The RT-qPCR reactions were run in triplicate with mean values and standard deviation calculated for all cycle thresholds (Ct). For each gene, the relative expression levels were normalized with respect to the wild-type plant (reference value = 1).

Calculations of Trichome Density {#sec010}
--------------------------------

Young leaves of wild-type and transformed *A*. *carvifolia* plants were selected to calculate trichome density. The ventral surfaces of the leaves were attached to the microscopic slides with the help of glue and maintained overnight at 4°C for dehydration. The next day, slides were examined under the fluorescent microscope (Leica) with a FITC green filter. Slides were studied under 10 x resolution and photographed.

Statistical Analysis {#sec011}
--------------------

All the experiments including, PCR, LC-MS and trichomes analysis were run in triplicate. Results of metabolites content, real-time quantitative PCR analysis and calculation of trichome density are presented as mean values ± S.E. The data obtained for the qualitative and quantitative analysis of artemisinin and derivatives, was analyzed statistically by ANOVA and Duncan's multiple range test using Mstat C software. Statistical significance of trichomes was determined by t-test (\*\*, P\<0.01; \*, P\<0.05).

Results {#sec012}
=======

Plant Identification through DNA Barcoding {#sec013}
------------------------------------------

We successfully amplified the 500bp *psbA-trnH* region of the chloroplast genome ([Fig 4](#pone.0140266.g004){ref-type="fig"}). DNA samples were sequenced in triplicate to confirm the authenticity of species-specific nucleotides and the same results were obtained. GenBank Accession number \[NCBI: FJ418751\] was used as the reference sequence to confirm our plant species. After performing the CLUSTAL W in BioEdit software (version 7.2.5.0) and BLAST in NCBI, our sequence was confirmed as *psbA-trnH* of *A*. *carvifolia*.

![Sequenced DNA for identification of *Artemisia carvifolia*.\
*psbA-trnH* sequence of *A*. *carvifolia* \[NCBI: FJ418751\].](pone.0140266.g004){#pone.0140266.g004}

Transformation and Regeneration {#sec014}
-------------------------------

Successful transformation of *A*. *carvifolia* with *A*. *tumefaciens* GV3101 harboring the *rol B* and *rol C* genes was carried out. Two independent transformation experiments were performed and 300 explants were used per transformation event. Transformation efficiency was found to be 20--30%, but only 4 *rol B* transformants and 3 *rol C* transformants survived to maturity on the selection media. Morphological variability was observed between wild type plants and transgenics bearing *rol B* and *rol C* genes. *Rol B* transgenics, which grew faster on the selection media, had wider leaves and more inflorescence, while *rol C* transgenics, which were recalcitrant to regeneration, showed a narrow leaf blade and decreased internode length ([Fig 5](#pone.0140266.g005){ref-type="fig"}).

![Vegetative propagation of *A*. *carvifolia*.\
*A*.*carvifolia*: wild type (a), tranformant of *rol* B (b), and transformant of *rol* C gene (c).](pone.0140266.g005){#pone.0140266.g005}

Molecular Analysis {#sec015}
------------------

PCR performed for the *rol* gene transformants of *A*.*carvifolia* showed the amplified products of 779bp for *rol B*, 540bp for *rol C* and 781bp for the *nptII* gene, as shown in [Fig 6A, 6B & 6C](#pone.0140266.g006){ref-type="fig"}. Similar amplified products were obtained from plasmid DNA of GV3101-CaMVBT and GV3101-CaMVC, respectively. Wild-type plants did not show the presence of these genes in their genome. Southern blot analysis confirmed the integration of *rol* genes in the plant genome and also gave an idea of the gene copy number in independent transgenic lines ([Fig 7A & 7B](#pone.0140266.g007){ref-type="fig"}). RT-PCR confirmed gene expression in all the regenerants, although it varied according to line. The *rol B* transgenic line TB4 and the *rol C* transgenic line TC1 showed the highest levels of expression ([Fig 6D & 6E](#pone.0140266.g006){ref-type="fig"}), and were found to harbor two copies of the respective genes.

![Gel pictures showing results of normal and semi quantitative reverse transcriptase PCR.\
PCR amplified products of *rol* B (a), *rol* C (b) and *nptII* (c) gene are shown in the figure. Semi quantitative RT-PCR showing the relative expression of *rol B* (d) and *rol C* gene (e). β-actin used as housekeeping gene (f). TB1-TB4 represents *rol* b transgenics whereas TC1-TC3 represent *rol* C transgenics. WT stands for wild type plant of *A*. *carvifolia*, lane "P" refers to the plasmid DNA and lane M corresponds to the marker DNA (Fermentas). NC = negative control.](pone.0140266.g006){#pone.0140266.g006}

![Southern blot analysis.\
Southern blot analysis of *rol* B (a) and *rol* C transgenic lines (b). TB1-TB4 represents *rol* B transgenics whereas TC1-TC3 represent *rol* C transgenics. Whereas \"P" stands for the plasmid DNA run as positive control. WT indicates the wild type plant.](pone.0140266.g007){#pone.0140266.g007}

Analysis of Artemisinin and Derivatives by LC-MS {#sec016}
------------------------------------------------

Analysis of artemisinin and its derivatives by an HPLC-UV system was not successful in *A*. *carvifolia* plant extracts, as these compounds lack UV properties or fluorescent chromophores. Thus, a method coupling HPLC with an Ion trap mass spectrometer was developed for the detection and quantification of the target compounds ([S1 Fig](#pone.0140266.s001){ref-type="supplementary-material"}), with the following results found in *A*.*carvifolia* wild-type extracts: artemisinin (8μg/g), artesunate (2.24μg/g), dihydroartemisinin (D.H.A. 13.6μg/g) and artemether (12.8μg/g). Differences in content between transformed and untransformed plants were observed. All the transgenic lines bearing *rol B* and *rol C* genes showed enhanced levels of artemisinin and its derivatives ([Fig 8](#pone.0140266.g008){ref-type="fig"}). In *rol B* gene transgenics, this increase was 3-7-fold for artemisinin, 3--10 fold for artesunate, and 2.6--4 fold for D.H.A and artemether, whereas in *rol C* transgenics it was 3.8--6 fold for artemisinin, 4.4--8.9 fold for artesunate, 2.3--3.2 fold for D.H.A and 2.3--5 fold for artemether.

![Quantitative analysis of artemisinin and derivatives.\
Comparative analysis of artemisinin content and its derivatives in wild type *A*. *carvifolia* and transgenics of *rol* B and *rol* C gene. Values with same alphabet did not differ significantly at 5% probability level using LSD. Error bars indicate the standard deviation.](pone.0140266.g008){#pone.0140266.g008}

Statistical analysis was conducted in factorial design (8 X 4 X 3) to see the effect of different transgenic lines on the production of artemisinin and derivatives, where wild type *A*. *carvifolia* and transgenics of *rol B* and *rol C* gene showed significant difference (P\<0.05) in the production of studied compounds ([Table 1](#pone.0140266.t001){ref-type="table"}).

10.1371/journal.pone.0140266.t001

###### Analysis of variance for quantitative analysis of artemisinin and derivatives.

![](pone.0140266.t001){#pone.0140266.t001g}

  Source                                             Degree of freedom   Sum of Squares   Mean Square   F-Value     Probability
  -------------------------------------------------- ------------------- ---------------- ------------- ----------- -------------
  Factor A (Transgenic lines)                        7                   12415.863        1773.695      1457.1584   0.0000
  Factor B (artemisinin standards)                   3                   11897.913        3965.971      3258.1977   0.0000
  A X B (Transgenic lines X artemisinin standards)   21                  3329.692         158.557       130.2605    0.0000
  Error                                              64                  77.903           1.217                     
  Total                                              95                  27721.371                                  

Coefficient of Variation: 3.38%

Expression Analysis of Artemisinin and Trichome Biosynthetic Genes by Real Time qPCR {#sec017}
------------------------------------------------------------------------------------

Higher transcript levels of all the studied genes were found in the transgenics expressing the *rol B* and *rol C* genes, although in a variable fashion ([Fig 9](#pone.0140266.g009){ref-type="fig"}). Among the artemisinin biosynthetic genes, CYP71AV1 was the most expressed in the transgenics, increasing 15-30-fold in *rol B* transformants and 10-26-fold in *rol C* transformants, compared to wild-type plants. The ALDH1 and ADS genes showed a respective 3-9-fold and 2-6-fold enhanced expression in *rol B* transformants, and a 5-9-fold and 3-6-fold increase in *rol C* gene transformants. The transgenic lines TB4 and TC1 exhibited the highest expression of these genes. Real time qPCR analysis of the TFAR1 gene revealed variable expression among the transgenics, with a 3-6-fold and 2-4-fold increase, respectively, in *rol B* and *rol C* transformants.

![Expression analysis of artemisinin biosynthetic pathway genes and TAFR1 gene.\
Expression level of Artemisinin biosynthetic genes, ADS (a), CYP71AV1 (b) and ALDH1 (c) in wild type plant of *A*. *carvifolia* and transgenics of *rol B* and *rol C* gene. TFAR1 is the gene involved in trichome development and sesquiterpenoid biosynthesis (d).](pone.0140266.g009){#pone.0140266.g009}

The Trichome Density {#sec018}
--------------------

The number of glandular trichomes was determined on the adaxial side of 3 random pieces of fresh leaf material from each sample, with an accurately determined leaf area of approximately 1cm^2^ per leaf, and significant differences were observed, as shown in [Fig 10A & 10B](#pone.0140266.g010){ref-type="fig"}. The leaves of *A*. *carvifolia* Buch transformed with *rol B* and *rol C* genes produced more trichomes (20--30 trichomes/cm^2^) compared to untransformed leaves (5--7 trichomes/ cm^2^).

![Trichome density.\
Comparison of trichome density of transformed and untransformed plants of *A*. *carvifolia* (a). Graphical representation of trichomes densitiy of wild type *A*. *carvifolia* and transformants of *rol B* and *rol C* gene (b). TB1-TB4 represent *rol B* transgenics whereas TC1-TC3 represent transgenics of *rol C* gene. WT represents the control plant of *A*. *carvifolia*.](pone.0140266.g010){#pone.0140266.g010}

Discussion {#sec019}
==========

In the *Artemisia* genus, many species have been evaluated for their artemisinin content, as mentioned previously, but not *A*. *carvifolia*. We not only detected this important phytochemical and its derivatives in this species, but also demonstrated the effectiveness of *rol* genes in the enhancement of their production.

Transformants harboring the *rol B* and *rol C* genes showed a significant increase in their artemisinin content compared to the wild-type plants. Data was obtained from four independent regenerants bearing *rol B* and three independent regenerants bearing *rol C* genes from two independent transformation events. Significant differences were observed in terms of growth, morphology and artemisinin accumulation. Similar findings were reported by another group who showed for the first time how *rol* genes induce ginsenoside overproduction in transformed cell cultures of *P*. *ginseng*, using plasmid DNA containing the individual *rol* genes from the TL-DNA of *A*. *rhizogenes* strain A4 \[[@pone.0140266.ref032]\]. *Rol C* cultures were found to accumulate 1.8-3-fold more ginsenoside than the control plant, although *rol B* lines were less productive \[[@pone.0140266.ref032]\]. In contrast, in another study, *rol B* gene transgenics of *Rubia cardifolia* showed enhanced production of anthraquinones. In fact, β-glucosidase, the product of the *rol B* gene, releases indole acetic acid (IAA) from its inactive glucose conjugates, thus increasing auxin sensitivity in *rol B* transgenics. Increased tyrosine phosphatase activity in *rol B* transformed cells distracts the signal transduction pathway of the hormone. This sensitivity and transduction alters the physiological configuration of transformed cells and eventually the whole plant \[[@pone.0140266.ref011],[@pone.0140266.ref033],[@pone.0140266.ref034],[@pone.0140266.ref035],[@pone.0140266.ref036],[@pone.0140266.ref037],[@pone.0140266.ref038]\]. Results of current study are also supported by our previous report, demonstrating the effect of *rol ABC* gene construct in over production of artemisinin in *A*. *dubia* \[[@pone.0140266.ref015]\]. It may be concluded that *rol* genes when combined show more increase in artemisinin i.e. up to 10 fold \[[@pone.0140266.ref015]\] than individual *rol* genes as we observed 6--7 fold increase.

Morphological alterations and differences in metabolite concentrations among different transgenic lines have been related to the variability in the expression level of *rol* genes and can be attributed to position effect \[[@pone.0140266.ref039]\]. Transgenic *A*. *carvifolia* plants showing increased secondary metabolism and morphological alterations were analyzed by Southern blot analysis and semi quantitative RT-PCR to check for correlation with *rol* gene expression. The results confirmed that the observed changes were due to the presence of *rol* gene transcripts. The lines TB4 and TC1, which produced the most artemisinin, also showed more *rol* B and *rol* C gene transcripts, respectively, as well as bearing two copies of these genes. On the other hand, the unimproved accumulation of artemisinin and derivatives observed in TB2 can be attributed to the less *rol B* transcripts in this line. Identical conclusions were reached by other group who correlated the altered morphology of tobacco plants with the detectable level of *rol* gene transcripts \[[@pone.0140266.ref039]\]. Our findings are also supported by Arshad et al. (2014), who found that the expression of the *rol* B gene of *A*. *rhizogenes* in tomato enhanced nutritional contents and foliar tolerance against fungal pathogens \[[@pone.0140266.ref040]\].

Variation in the target compound accumulation also reflected differential expression of artemisinin biosynthetic genes, with transcript levels being higher in the *rol B* and *rol C* gene transformants than in wild-type *A*. *carvifolia*. Key enzymes in artemisinin biosynthesis include ADS, CYP71AV1 and ALDH1, which are involved in the conversion of farensyl diphosphate (FDP) to artemisininc acid, a late precursor of artemisinin. Biochemical and molecular studies have shown that amorpha-4,11-diene synthase (ADS) catalyzes the first committed step of the artemisinin biosynthetic pathway by cyclization of FDP to amorpha-4,11-diene \[[@pone.0140266.ref041],[@pone.0140266.ref042],[@pone.0140266.ref043],[@pone.0140266.ref044],[@pone.0140266.ref045]\], which is then oxidized to artemisinic alcohol by the cytochrome P450 enzyme, CYP71AV \[[@pone.0140266.ref018],[@pone.0140266.ref046]\]. Conversion of artemisinic alcohol to artemisinin through formation of either artemisinic aldehyde or dihydroartemisinic aldehyde is catalyzed by the enzyme ALDH1. *A*. *carvifolia* transgenics bearing *rol B* and *rol C* genes showed an overexpression of the studied genes, particularly CYP71AV, in accordance with the accumulation of artemisinin and its derivatives. These results are supported by a previous study showing that overexpression of farensyl pyrophosphate synthase (FPS) in *A*. *annua* increases the accumulation of artemisinin through conversion of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMADP) into FDP \[[@pone.0140266.ref047]\]. Similarly, an increased amount of artemisinin was reported in *A*. *annua* transformed with *Agrobacterium tumefaciens* wild type nopaline strains \[[@pone.0140266.ref048]\], and in *A*. *annua* bearing the *Ipt* gene \[[@pone.0140266.ref049]\].

The TFAR1 gene, which stimulates trichome development and catalyzes sesquiterpenoid biosynthesis \[[@pone.0140266.ref047]\], also showed increased expression in all transgenic lines. Accordingly, the transformed lines of *A*. *carvifolia* showed higher trichome density compared to the untransformed plant ([Fig 10](#pone.0140266.g010){ref-type="fig"}), indicating a relationship between trichome development and artemisinin production. Glandular trichomes are the sites of production of important phytochemicals, including artemisinin, as well as many others that have found numerous applications in the pesticide, pharmaceutical and flavor industries \[[@pone.0140266.ref050],[@pone.0140266.ref051]\]. Artemisinin content has been directly related to the trichome index in previous studies \[[@pone.0140266.ref052],[@pone.0140266.ref053]\].

Conclusion {#sec020}
==========

Transformation of *A*. *carvifolia* with *rol B* and *rol C* genes resulted in the enhancement of its secondary metabolites, particularly artemisinin and its derivatives. Transcript levels of the *rol B* and *rol C* genes found in the transgenics also correlated with their artemisinin accumulation pattern. An altered expression of genes involved in artemisinin and trichome biosynthesis was observed, indicating that *rol* gene integration modified the plant metabolism. The relationship between trichome density and accumulation of artemisinin and its derivatives was confirmed.

Supporting Information {#sec021}
======================

###### LC-MS chromatograms.

LC-MS chromatograms for standard artemisinin, artesunate, D.H.A and artemether for C18, 5μm (150mm x 4.6mm) column along with their mass spectra.

(TIF)

###### 

Click here for additional data file.

We are grateful to Dr. A. Spena, Max-Planck-Institut fur Zuchtungsforschung, 5000 Koin 30, FRG for providing the vectors. Furthermore, we are appreciative to Dr. Muhammad Tahir Waheed for his intellectual support.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: MB. Performed the experiments: ED. Analyzed the data: KRE MB. Contributed reagents/materials/analysis tools: MB RMC. Wrote the paper: ED BM RMC.
