Impact Factor 4.402 | CiteScore 7.8
More on impact ›


Front. Plant Sci., 07 July 2020 |

Transgenic Cotton (Gossypium hirsutum L.) to Combat Weed Vagaries: Utility of an Apical Meristem-Targeted in planta Transformation Strategy to Introgress a Modified CP4-EPSPS Gene for Glyphosate Tolerance

  • 1ICAR-National Institute for Plant Biotechnology, Pusa Campus, New Delhi, India
  • 2Amity Institute of Biotechnology, Amity University Haryana, Gurugram, India

Weeds burden plant growth as they compete for space, sunlight, and soil nutrients leading to 25–80% yield losses. Glyphosate [N-(phosphonomethyl)glycine] is a widely used broad spectrum non-selective herbicide that controls weeds by inhibiting 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) enzyme and interfering with the shikimate biosynthesis pathway. Cotton (Gossypium hirsutum L.) is one of the most important commercial crops grown worldwide for its fiber. We have developed herbicide tolerant transgenic cotton (cv. P8-6) by introgression of a codon-optimized and modified EPSPS gene (CP4-EPSPS) possessing an N-terminal chloroplast targeting peptide from Petunia hybrida. Because of the recalcitrant nature of cotton, a genotype-independent non-tissue culture-based apical meristem-targeted in planta transformation approach was used to develop transformants. Although in planta transformation methodologies are advantageous in developing a large number of transgenic plants, effective screening strategies are essential for initial identification of transformants. In the present study, the use of a two-level rigorous screening strategy identified 2.27% of T1 generation plants as tolerant to 800 and 1,500 mg/L of commercially available glyphosate (Roundup). Precise molecular characterization revealed stable integration, expression, and inheritance of CP4-EPSPS in advanced generations of the promising transgenic events. Further, superiority of selected transgenic plants in tolerating increasing levels of glyphosate (500–4,000 mg/L) was ascertained through reduced accumulation of shikimate. This report is the first of its kind where cotton transformants tolerating high levels of glyphosate (up to 4,000 mg/L) and accumulating low levels of shikimate have been identified. This study not only reiterated the genotype-independent nature of the transformation strategy but also reiterated the translational utility of the CP4-EPSPS gene in management of weeds.


Cotton (Gossypium hirsutum L.) is one of the most important commercial crop plants grown worldwide for fiber and seeds. India is the second largest cotton producer contributing to approximately 22% of the world production (Directorate of Cotton Development Government of India, 2017; Mohapatra and Saha, 2019). Weeds are one of the very many factors that negatively challenge productivity of cotton accounting to approximately 50–85% of yield losses (Jain et al., 1981; Frisvold et al., 2009; Ayyadurai and Poonguzhalan, 2011; Awan et al., 2015; Ramachandra et al., 2016). During the early stages of growth, high rate of competition exists between weeds and crop plants for nutrition, space, and amenability, leading to a high percentage of plant mortality in the field (Janaki et al., 2015). Chemical herbicides are an effective way to control weeds, but have raised concerns of environmental pollution. Hence, to overcome these weeds, use of integrated weed management system in the field is an option (Sharma, 2008).

Crop improvement through genetic engineering has given the feasibility to engineer plants to combat herbicides being sprayed and thus help in the safety of the crops without major yield losses (Benbrook, 2016; Bohn and Lovei, 2017). Various genes that confer herbicide tolerance have been identified, validated, and introgressed into a multitude of crops including cotton (Awan et al., 2015; Liang et al., 2017; Zhang et al., 2017).

Glyphosate [N-(phosphonomethyl) glycine] is a widely used broad-spectrum non-selective herbicide (Funke et al., 2006; Beckie, 2011; Heap, 2014) that controls weeds by interfering with the shikimate biosynthesis pathway and inhibiting the rate-limiting enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS; Steinrucken and Amrhein, 1980; Gianessi, 2005; Cao et al., 2012; Gong et al., 2016). The EPSPS gene derived from the soil bacterium Agrobacterium tumefaciens strain CP4 is widely known for its tolerance toward glyphosate in an array of transgenic crops (Stalker et al., 1985; Kishore et al., 1986; Padgette et al., 1991; Berry et al., 1992; Padgette et al., 1996; Zhou et al., 2006; Dun et al., 2014; Ma et al., 2016; Fartyal et al., 2018; Yi et al., 2018).

In the present study, we provide a meticulous and in-depth utility of a codon-optimized and modified CP4-EPSPS gene in cotton. The gene was assembled with varying GC content for plant codon usage, fused with N-terminal Petunia hybrida chloroplast targeting signal peptide (CPT) as the shikimate pathway takes place in chloroplast (Della-Cioppa et al., 1986; Chhapekar et al., 2015). Effectiveness of the gene was further validated in rice for its tolerance to glyphosate (Chhapekar et al., 2015).

Non-tissue culture-based transformation strategies using apical meristem facilitated in the development of stable transgenics in several recalcitrant crop species (Kumar et al., 2011; Keshavareddy et al., 2013; Singh et al., 2018). One such apical meristem-targeted in planta transformation approach was exploited in the present study for the development of transformants in cotton (Rohini and Rao, 2002; Keshamma et al., 2008; Karthik and Rohini, 2017). Despite various efforts to introgress herbicide tolerance in cotton by different transformation strategies (Bayley et al., 1992; Lyon et al., 1993; Nida et al., 1996; Rajasekaran et al., 1996b; Keller et al., 1997; Awan et al., 2015; Liang et al., 2017; Zhang et al., 2017), this study is the first of its kind to demonstrate development of stable transgenics with the ability to tolerate high glyphosate using a non-conventional transformation strategy. In this direction, the present study demonstrates the utility of a modified CP4-EPSPS gene in transgenic cotton and also reiterates the applicability of the in planta transformation approach for cotton improvement.

Materials and Methods

In planta Transformation of Cotton

Agrobacterium tumefaciens strain EHA 105 harboring the binary vector pBinAR was used for transformation (Figure 1A). The binary vector contained codon optimized and modified CP4-EPSPS gene (GenBank accession number is KJ787649) driven by CaMV35s promoter and neomycin phosphotransferase II (nptII) gene driven by the nopaline synthase promoter (Figure 1A). For in planta transformation, seeds of cotton cultivar P8-6 were surface sterilized and maintained at 32°C for germination. Two-day-old seedlings were used for transformation by following the standardized protocol (Rohini and Rao, 2002; Keshamma et al., 2008; Karthik and Rohini, 2017).


Figure 1. Agrobacterium tumefaciencs–mediated in planta transformation in cotton. (A) Vector map of pBinAR/EPSPS used in the study. Codon-optimized and modified CP4-EPSPS gene (1,561 bp) was assembled with varying GC content for plant codon usage and fused at the N-terminal (1–216 bp) with Petunia hybrida chloroplast targeting signal peptide (CPT). (B–D) Various stages of recovery and establishment of the primary transformants in the greenhouse.

Glyphosate Tolerance Assays

Standardization of Glyphosate Concentration for Screening at Seedling Stage

Seeds of cotton cultivar P8-6 were grown in plastic trays (equal in number) for 8–10 days and at the fourth leaf stage; the plants in each tray were sprayed with glyphosate (200, 400, 600, 800, and 1,000 mg/L) equivalent to 2.5% (vol/vol) Roundup (Monsanto, St. Louis, Missouri, United States). Along with these trays, untreated control was also maintained with plants that did not receive the glyphosate spray. The plants were later scored for their recovery, and the concentration of glyphosate that resulted in the death of wild-type plants was used to select putative transformants.

Standardization of Glyphosate Concentration for Screening at Plant Stage

Screening at plant stage was conducted with leaves of 45-day-old wild-type plants using different glyphosate concentrations (50, 100, 150, 200, 250, 300, 500, 1,000, 1,500, 2,000, 2,500, and 3,000 mg/L). Cut leaf discs were placed in Petri plates in three replicates, suspended with 10 mL of glyphosate and maintained for 96 h at room temperature. The leaf discs were later scored for chlorosis.

Estimation of Chlorophyll Content in Glyphosate-Treated Leaves

Well-expanded leaves of 45-days-old transgenic and wild-type cotton plants grown under similar conditions were selected for the experiment. Six leaf discs (nearly 200 mg of fresh weight) each of transgenic and wild-type plants were treated with 1,500 mg/L of glyphosate and incubated in Petri plates for 48 h at room temperature in dark. They were later transferred to Eppendorf tubes (Eppendorf, Hamburg, Germany) containing 10 mL of 80% acetone: dimethyl sulfoxide (DMSO) (1:1 vol/vol) solution and incubated in dark for 72 h. Chlorophyll a and b that leached out of the discs were measured spectrophotometrically at 665 and 650 nm, respectively (Hiscox and Israelstam, 1979). Because the assay was conducted with individual T1 generation plants, chlorophyll was estimated in three technical replicates, and the average was presented. The data were analyzed using R Statistical Language version 3.4.3 (2019-2107-30) and Tukey honestly significant difference test to assess significant difference within selected samples.

Shikimate Assay

Shikimate assay was performed with leaves (six leaf discs nearly 200 mg of fresh weight for each concentration of glyphosate) of 45-days-old transgenic plants and wild type. Leaf discs were treated with 10 mL of glyphosate of different concentrations (500–4,000 mg/L) for 96 h. Leaf discs of wild-type treated with water were taken as untreated control, and wild-type leaf discs treated with different concentrations of glyphosate were taken as treated control. Further, the leaf discs were collected into an Eppendorf tube; 0.4 mL of 0.25 N HCl was added and incubated for 60 min. It was further ground using a micro pestle, mixed with 0.4 mL 0.25% periodic acid and 0.25% metaperiodate solution for 60 min. After the periodic acid–metaperiodate incubation, a 0.4 mL aliquot of 0.6 M sodium hydroxide with 0.22 M sodium sulfite solution was added (Gaitonde and Gordon, 1958; Koger and Reddy, 2005; Shaner et al., 2005). A shikimate standard curve was developed by adding known amounts of shikimate to vials containing water. Optical density of the solution at 382 nm was determined spectrophotometrically for all the samples taken in the study. The assay was performed with three biological and two technical replicates of each transgenic event. The amount of shikimic acid accumulation in each of the samples was measured, and average values were plotted to evaluate the performance of transgenic events.

Molecular Analysis of Transgenics for T-DNA Integration

Polymerase Chain Reaction Analysis

Genomic DNA from leaves of transgenic and wild-type plants was isolated following a modified cetyl trimethyl ammonium bromide method (Porebski et al., 1997). Polymerase chain reaction (PCR) analysis was performed to amplify a 500-bp T-DNA-specific right border region and 750-bp nptII gene using gene specific primers (Table 1). In case of T2 generation plants, PCR analysis was also performed with primers to amplify 420-bp CP4-EPSPS gene fragment (Table 1). Polymerase chain reaction amplification was carried out in a thermal cycler (Eppendorf) using the following program: one cycle of initial denaturation at 95°C for 5 min followed by 35 cycles of denaturation at 95°C for 1 min, annealing at 58°C for 1 min, extension at 72°C for 1 min, and one cycle of final extension of 10 min at 72°C. Each PCR reaction (25 μL) consisted of 100 ng of genomic DNA, 2.5 μL of 10× Taq buffer, 10 pM each of forward and reverse primer, 200 μM dNTPs, and 1 U of Taq DNA polymerase (Bangalore Genei, Bengaluru, India) and made up to a final volume of 25 μL with nuclease-free water. While “Blank” had nuclease-free water instead of genomic DNA, negative control contained 100 ng of genomic DNA of wild-type and positive control contained 25 ng of modified CP4-EPSPS plasmid.


Table 1. List of primers used in the study.

Genomic Southern Analysis

For genomic Southern analysis to detect T-DNA copy number, genomic DNA from transgenic and wild-type cotton plants were isolated and purified, and 15 μg of DNA was digested separately with HindIII (NEB high fidelity, NEB) and BamHI (NEB high fidelity, NEB). The digested DNA samples were then separated by electrophoresis on a 0.8% agarose gel in 1X TAE buffer at constant voltage of 40 V, blotted onto a positively charged nylon membrane (Bio-Rad, Hercules, California, United States) and UV cross-linked. While the membrane with HindIII-digested genomic DNA samples was probed with DIG PCR-labeled 750-bp nptII fragment, the membrane with BamHI-digested genomic DNA samples was probed with DIG PCR-labeled 420-bp CP4-EPSPS fragment (Table 1). Hybridization, washing, blocking, and development were carried out according to manufacturer’s instructions (Roche Holding AG, Basel, Switzerland). The membranes were exposed to X-ray film for 1 h in dark and later observed for bands.

Western Blot Analysis

Total protein was isolated from 100 mg leaf samples of both transgenic and wild-type plants by finely grinding with liquid nitrogen and later reconstituting in 500 μL of extraction buffer [0.1 M Tris-HCl, 0.5 M EDTA, 30% sucrose, 1% PVPP, 0.1 M KCl, 2% sodium dodecyl sulfate (SDS), 1 mM phenylmethanesulfonyl fluoride, 5% 2-mercaptoethanol; pH 8.8] (Wu et al., 2014). Isolated protein was quantified with Bradford’s assay (Bio-Rad, Hercules, California, United States), and 20 μg of protein was run on SDS–polyacrylamide gel electrophoresis and blotted onto nitrocellulose membrane (Millipore, Burlington, Massachusetts, United States) at a constant voltage (40 V) for 3 h in transfer buffer (1.5 g Tris, 7.2 g glycine and 100 mL methanol made up to 500 mL with distilled water). Initiation of development was performed by blocking the membrane with commercially available NAP blocker (G-Biosciences, St. Louis, Missouri, United States) followed by hybridization of the membrane with protein-specific primary antibodies (Amar Immunodiagnostics, Jubilee Hills, Hyderabad, India) against EPSPS with a dilution of 1:800. Subsequently, the membrane was hybridized with horseradish peroxidase–conjugated secondary antibody (Bangalore Genei) with 1:500 dilution. After blocking with secondary antibodies, the membrane was washed four times with 1X phosphate-buffered saline for 10 min each, followed by the addition of TMB substrate (Promega, Madison, Wisconsin, United States) for color development.

Real-Time PCR

Total RNA was isolated from both transgenic and wild-type cotton plants (SpectrumTM; Sigma-Aldrich, St. Louis, Missouri, United States); quantity and integrity of total isolated RNA were assessed using a Nanodrop Micro Photometer (Thermo Fisher Scientific, Waltham, Massachusetts, United States). cDNA was synthesized from 2.5 μg total RNA according to manufacturer’s instructions (SuperScript® VILOTM; Invitrogen, Carlsbad, California, United States). Gene-specific CP4-EPSPS primers and ubiquitin primers (internal control) were used to set up the reaction for quantitative PCR (qPCR; Table 1). Polymerase chain reaction program of 5-min initial denaturation at 95°C followed by 40 cycles of denaturation at 95°C for 20 s, annealing at 60°C for 20 s, and extension at 72°C for 15 s was employed. A final melting curve consisting of 95°C for 30 s, 65°C for 30 s, and 95°C for 30 s per cycle was used. Expression analysis was performed with three biological and four technical replicates of each of the events in T2 generation. The Ct values of samples were normalized with the Ct value of ubiquitin (internal control) to calculate Δ Ct value.


Development of Transformants in Cotton Harboring the Modified CP4-EPSPS Gene

Apical meristem-targeted in planta transformation strategy was employed for the development of transgenic cotton. Initially, 180 seedlings were infected with Agrobacterium, of which 65 plants recovered. Primary transformants were initially subjected to recovery for 2 weeks under controlled conditions and later transferred to pots in the greenhouse (Figures 1B–D), where they grew normally, flowered, and set seeds. As in planta transformation results in the development of chimeric T0 plants, analysis for the selection of putative transformants was carried out in T1 generation.

Identification of the Suitable Glyphosate Concentration for the Selection of Putative Transformants

In an effort to identify putative transformants tolerant to glyphosate, the concentration of glyphosate to be used for screening of transgenic seedlings was initially standardized (Figures 2Ai–vi). Approximately 7–8 days after spraying different concentrations of glyphosate, the plants were scored for wilting and chlorosis. Although lower concentrations of glyphosate (200 and 400 mg/L; Figures 2Aii,iii) did not seem to be effective, symptoms were seen in plants that were sprayed with 600 mg/L glyphosate and above (Figures 2Aiv–vi). Higher concentrations of glyphosate (800 and 1,000 mg/L; Figures 2Bi–iii) showed chlorosis and necrosis, eventually leading to death of plants. Based on the response of the wild-type plants, 800 mg/L glyphosate was used for the identification of putative transformants.


Figure 2. Standardization of glyphosate tolerance assays in cotton cv. P8-6 at seedling and plant level. (A) Standardization of glyphosate concentration to be sprayed at seedling level: (i) untreated control; seedlings sprayed with different concentrations of glyphosate, (ii) 200 mg/L, (iii) 400 mg/L, (iv) 600 mg/L, (v) 800 mg/L, (vi) 1,000 mg/L. (B) (i–iii) Response of the cotton seedlings to different glyphosate concentrations 7–8 days after the spray. (C) Standardization of glyphosate concentration in leaf disc assay (i) 50 mg/L, (ii) 100 mg/L, (iii) 150 mg/L, (iv) 200 mg/L, (v) 250 mg/L, (vi) 300 mg/L, (vii) 500 mg/L, (viii) 1,000 mg/L, (ix) 1,500 mg/L, (x) 2,000 mg/L, (xi) 2,500 mg/L, (xii) 3,000 mg/L.

Similarly, leaf disc assay with wild-type plants provided insights about the glyphosate tolerance ability at plant level. The assay demonstrated that the leaf discs could not resist glyphosate concentration greater than 150 mg/L (Figures 2Ci–iii). Although symptoms of glyphosate effect was seen with 200 mg/L; 500 mg/L glyphosate and above turned out to be lethal (Figures 2Civ–xii). However, for a stringent plant level selection, 1,500 mg/L was used to identify tolerant transformants.

Identification and Molecular Characterization of Transgenic Plants in T1 Generation

For the identification of putative transformants, three to four leaf stage T1 generation plants were sprayed with 800 mg/L glyphosate. After 7–8 days, as in the treated control (Figure 3Aii; treated control), some T1 generation plants showed necrotic symptoms, turned yellow initially, eventually to brown, and died because of susceptibility toward glyphosate (Figures 3Aiii–vi). However, a set of plants was healthy and on par with the untreated control (Figure 3Ai). Accordingly, a total of 35 plants out of 660 (6.06%) that endured glyphosate spraying were transferred to soil in pots and grown under greenhouse conditions.


Figure 3. Glyphosate tolerance assays for the identification of putative transformants. (A) Seedling level screening: T1 generation and wild-type seedlings were sprayed with 800 mg/L glyphosate (i) untreated control, (ii) treated control, (iii–vi) representative trays containing putative transformants. (B) Plant level screening: Response of leaf discs of 45-day-old representative transgenic and wild-type plants to 1,500 mg/L glyphosate. (C) Estimation of chlorophyll in the glyphosate treated transgenic plants vis-à-vis wild type.

Further, to reduce the number of false positives, a second plant level screening was performed with the leaf discs of 35 T1 plants that recovered after glyphosate spraying. Treatment of leaf discs of transgenics and wild-type plants with 1,500 mg/L glyphosate identified 16 tolerant transgenic plants (1-1, 2-1, 3-4, 4-1, 4-4, 4-5, 4-7, 6-5, 17-2, 21-1, 22-5, 39-8, 39-9, 42-1, 42-3, and 46-1) that remained green when compared to other transgenic plants that showed symptoms of bleaching and necrosis (Figure 3B). There was corroboration between the phenotype and the amount of chlorophyll, indicating the superiority of the 16 transgenic plants (Figure 3C). Nearly 2.27% of the putative transformants (16 of 660 plants) were finally selected for molecular characterization on the basis of two stringent glyphosate tolerance assays, as well as normal growth and phenotype.

To prove the presence and integration of the transgenes in the selected cotton transformants, PCR analysis performed using primers for T-DNA-specific right border and nptII gene-specific primers demonstrated the desired amplicons of size of 500 and 750 bp, respectively (Figures 4A,B). Evidence for the expression of modified CP4-EPSPS gene in nine selected transgenic events (1-1, 2-1, 3-4, 4-4, 22-5, 39-9, 42-1, 42-3, and 4-5) was further supported by Western blot analysis (Figures 4Ci,ii). Transgenic plants developed a single intense band of 45 kDa at the expected position as in the positive control (purified CP4-EPSPS), whereas there was no band observed in the wild type.


Figure 4. Molecular analysis of transformants in T1 generation. Polymerase chain reaction analysis of T1 generation plants for the amplification of panel (A) 500-bp right border region of T-DNA. (B) 750-bp nptII fragment. Lane M: 1 kb Ladder (Thermo Fisher Scientific); Lane B: Water Blank (all PCR components but without template DNA); Lane WT: wild type; Lanes 1–24 of panel (A) and 1–26 of panel (B): T1 generation transgenic plants; Lane P: positive control (modified CP4-EPSPS binary vector). (C) Western blot analysis of selected transgenic plants along with wild type. (i) Lanes 1–8: transgenic plants (1-1, 2-1, 3-4, 4-4, 22-5, 39-9, 42-1, 42-3). Lane L: Prestained protein ladder; Lane WT: wild type; (ii) Lane L: prestained protein ladder; Lane P: purified protein of CP4-EPSPS (30 ng); Lane WT: wild type; Lane 4: transgenic plant 4-5. (D) Genomic Southern analysis of T1 generation transgenic plants. (i) HindIII digested genomic DNA probed with DIG-labeled 750-bp nptII fragment. Lane L: Lambda DNA HindIII digest; Lane WT: wild-type DNA; Lanes 3–14: DNA of transgenic plants (46-1, 42-1, 42-3, 39-9, 22-5, 21-1, 17-2, 4-4, 4-5, 3-4, 2-1, 1-1); Lane P: positive control (10 pg DNA of modified CP4-EPSPS binary vector). (ii) BamHI digested genomic DNA probed with 420-bp CP4-EPSPS gene specific fragment. Lane L: Lambda DNA HindIII digest; Lane 2: 4-4; Lane 3: 4-5; Lane 4: positive control (10-pg modified CP4-EPSPS vector DNA); Lane 5: wild-type DNA.

To assess the integration of T-DNA, 12 T1 generation plants (46-1, 42-1, 42-3, 39-9, 22-5, 21-1, 17-2, 4-4, 4-5, 3-4, 2-1, and 1-1) with robust phenotype that had performed well in both the glyphosate tolerance assays were selected. Genomic Southern analysis confirmed the integration of transgene in three transgenic plants and demonstrated the independent nature of the transgene (Figure 4Di). Two events, 4-4 and 4-5, showed single-copy insertions, whereas plant 42-3 had two copies of the transgene integrated into its genome. Corroboratory results were observed when genomic Southern analysis was performed with the two positive samples (events, 4-4 and 4-5) and probed with 420-bp CP4-EPSPS probe (Figure 4Dii), indicating their transformed nature.

Characterization of Transformants for Inheritance of Transgene in T2 Generation

To assess the inheritance of the transgenes, analysis of the transgenics in T2 generation was carried out with Southern positive T1 generation events viz., 4-4, 4-5, 42-3, and another event 22-5 that performed well in glyphosate tolerance assay. Equal numbers of seeds (approximately 30) of both transgenics and wild-type plants were germinated in individual plastic trays, and fourth leaf stage (approximately 8–10 days after germination) healthy plants (Figures 5Ai–iv) were sprayed with 1,000 mg/L glyphosate. It was observed that after 7 days of spraying, the transgenic plants could apparently withstand the treatment (Figures 5Biii,iv). The treated control plants turned necrotic and eventually died (Figure 5Bii). After 15 days of glyphosate spraying, transgenic events 4-4, 4-5, and 42-3 showed approximately 80% recovery (Figures 5Ciii,iv), whereas line 22-5 had fewer number of plants that has recovered the stress (only three plants), and none recovered among treated wild type (Figure 5Cii). Untreated control plants remained healthy as they were not sprayed with glyphosate (Figure 5Ci). Further proof for the stable inheritance of the T-DNA (Table 2) was provided by segregation analysis carried out with the selected events based on their response to glyphosate spray.


Table 2. Transgene inheritance analysis in the selected transgenic events.


Figure 5. Glyphosate tolerance assay to assess inheritance in T2 generation. (A) Wild-type [(i) untreated control, (ii) treated control] and selected T2 generation transgenic cotton plants [(iii) Event 4-4, (iv) Event 4-5] maintained under greenhouse conditions. (B) Response of transgenic plants vis-à-vis wild-type to 1,000 mg/L glyphosate 7–8 days after spray [(i) untreated control, (ii) treated control, (iii) Event 4-4, (iv) Event 4-5]. (C) Recovery of the transgenic plants vis-à-vis wild-type 15 days after glyphosate spray [(i) untreated control, (ii) treated wild type, (iii) Event 4-4, (iv) Event 4-5].

The stable inheritance of the transgene was further confirmed by PCR in the plants that survived glyphosate spray using primers specific for CP4-EPSPS, T-DNA right border and nptII gene (Figures 6Ai–iii). To corroborate integration, expression, and efficacy, the selected transgenic plants, when further assessed for transcript accumulation by qPCR (Figure 6B), showed varied levels of expression of the modified CP4-EPSPS gene between 8- and 12-fold.


Figure 6. Characterization of cotton transformants in T2 generation. (A) Polymerase chain reaction analysis of T2 generation plants for the amplification of panel (i) 420-bp CP4-EPSPS fragment, (ii) 500-bp right border region of T-DNA, and (iii) 750-bp nptII fragment. Lane L: 1 kb ladder (Thermo Fisher Scientific); Lane B: template blank; Lane WT: wild type; Lane 1-12: transgenic plants (three lanes each of events 4-4, 4-5, 42-3, and 22-5); Lane P: positive control (modified CP4-EPSPS vector). (B) Quantitative RT-PCR analysis of transgenic plants for transcript accumulation of modified CP4-EPSPS. (C,D) Estimation of shikimate accumulation in the transgenic plants vis-à-vis wild type. (C) Representative picture showing the response of leaf discs treated with 3,000 mg/L glyphosate taken for shikimate assay. (D) Graphical representation of shikimate accumulation in transgenic events vis-à-vis wild type incubated with different concentration of glyphosate for 96 h.

Quantification of shikimic acid in T2 generation plants vis-à-vis wild type showed that the leaf discs of selected transgenic lines seemed to be performing consistently under the imposed stress (Figure 6C). Transgenic event 4-4 had the lowest amount of shikimate accumulation in all glyphosate concentrations (500–4,000 mg/L) in accordance with the shikimate standard. Transgenic events 4-5 and 43-2 showed similarity in the accumulation of shikimate (ranging between 1.0 and 1.2 μg/mL) in leaf samples under all the concentrations of glyphosate (Figure 6D). Although a higher amount of shikimate accumulation was observed in event 22-5, it was considerably lower than that of the treated control, which accumulated maximum amount of shikimate (>1.5 μg/mL). Untreated control had the lowest amount of shikimic acid accumulation in the plant (1.0 μg/mL).


Use of herbicides in modern agriculture has emerged as an integral approach to tackle the problem of weeds as they provide cost-effective increase in agricultural productivity. Because of herbicide usage, there has been at least 12% increase in yield owing to reduced weed competition (Jabran, 2016). Despite the availability of cheap agricultural labor force, large quantities of herbicides are being used in developing countries that could destabilize the functions of soil ecosystem due to high sensitivity to chemicals (Fletcher, 1960).

Glyphosate, one of the widely used herbicides effective in controlling weed populations, is not synthesized in animals as they do not produce aromatic amino acids (Steinrucken and Amrhein, 1980; Gianessi, 2005; Cao et al., 2012; Gong et al., 2016). Hence, glyphosate is less toxic to humans and highly toxic to plants. Resistance to glyphosate has therefore been a pertinent endeavor in crop improvement programs with conscious inputs from breeders and biotechnologists (Peters and Strek, 2018). Resistance to herbicides was among the first traits for which biotechnological approaches were applied (Wilkins et al., 2000). Biotechnological intervention in cotton improvement programs has been long standing with continued inputs toward many traits, including herbicide tolerance using varied genes (Rajasekaran et al., 1996b; Awan et al., 2015; Liang et al., 2017; Zhang et al., 2017). Nevertheless, this study is the first ever comprehensive report demonstrating development, analysis and efficacy of transgenic cotton with a modified CP4 EPSPS gene.

Considering the recalcitrant nature of cotton toward tissue culture regeneration, the demonstration of genotype-independent nature of transformation protocols assumes significance (Keshamma et al., 2008). As the initial transformants are chimeric, stringent screening using glyphosate was opted to identify putative transformants in T1 generation. The crux of the study was in the use of a stringent two-level screening of T1 generation plants, at seedling as well as at plant level. This strategy was deliberately designed to select the transformants that could withstand high levels of glyphosate, thus avoiding the advancement of false positives. Although approximately 6% of the seedlings survived the initial glyphosate (1,000 mg/L) challenge, 2.27% of the plants were finally selected for molecular characterization. This success in the identification of transformants in a difficult to transform crop such as cotton is highly commendable (Keshamma et al., 2008).

The foolproof efficacy of transgenics could be concluded by their consistent performance as shown by four events in T2 generation when treated with 1,000 mg/L of glyphosate and their subsequent recovery from the imposed stress. This further reiterates the superiority of the gene, CP4 EPSPS in conferring tolerance to high levels of glyphosate (Chhapekar et al., 2015) due to the presence of N-terminal P. hybrida CPT. In addition to the molecular characterization to confirm the transgene integration and inheritance, the study also demonstrated the superiority of the selected transformants through the ability to accumulate shikimate. It is known that EPSPS is a rate-limiting biochemical in the shikimate pathway that occurs in the chloroplasts, which is competitively inhibited by glyphosate. Treatment of plants with glyphosate inhibits the EPSPS activity leading to increased accumulation of shikimic acid (Berry et al., 1992). Studies have demonstrated this phenomenon in glyphosate sensitive wild-type cotton plants but reduced shikimate accumulation in transgenic plants (Tian et al., 2010). In the present study, despite high levels of glyphosate (up to 4,000 mg/L), shikimate assay precisely confirmed that the accumulation of shikimate was considerably low in the transgenics as evidenced by their recovery. Three transgenic events were very efficient in regulating shikimate accumulation even at higher concentrations of 4,000 mg/L as depicted by chlorophyll retention. This study is the first of its kind that has shown regulated accumulation of shikimate even at 4,000 mg/L in transgenic events (Rajasekaran et al., 1996a, b; Awan et al., 2015; Liang et al., 2017; Zhang et al., 2017).

The study explored the utility of an Indian cotton cultivar, P8-6, in the development of transgenics by employing an apical meristem-targeted in planta transformation strategy. Further, the findings of the study also ascertained the role of the modified CP4-EPSPS gene possessing the chloroplast transit peptide (CTP) to target the transgene to chloroplast for tolerating increased glyphosate concentrations in transgenic plants. The transgenics thus developed can form an important pool of resource for introgression of herbicide tolerance into other cotton genotypes and aid in weed management.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Author Contributions

KK has developed the transgenic plants. KK, MN, and AT performed the molecular analyses. KK has written the manuscript. SS, PM, and MR helped KK in preparing the manuscript. NS, PD, and MS have critically edited the manuscript. RS is responsible for acquisition of funds, designing experiments, editing, and revising the manuscript. All the authors have read the manuscript and declare no conflict of interest.


This research was funded by Network project for transgenic crops (NPTC; ICAR-NPTC/EFC/XII plan/3075), Indian Council of Agricultural Research (ICAR).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


Awan, M. F., Abass, A., Muzaffar, A., Ali, A., Tabassum, B., and Rao, A. Q. (2015). Transformation of insect and herbicide resistance genes in cotton (Gossypium hirsutum L.). J. Agric. Sci. Technol. 17, 287–298.

Google Scholar

Ayyadurai, P., and Poonguzhalan, R. (2011). Critical period of crop weed competition in zero-till cotton. Indian. J. Weed Sci. 43, 228–230.

Google Scholar

Bayley, C., Trolinder, N., Ray, C., Morgan, M., Quisenberry, J. E., and Ow, D. W. (1992). Engineering 2, 4-D resistance into cotton. Theor. Appl. Genet. 83, 645–649.

Google Scholar

Beckie, H. J. (2011). Herbicide resistant weed management: focus on glyphosate. Pest Manag. Sci. 67, 1037–1048.

Google Scholar

Benbrook, C. M. (2016). Trends in glyphosate herbicide use in the United States and globally. Environ. Sci. Eur. 28:3. doi: 10.1186/s12302-016-0070-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Berry, G., Kishore, G., Padgette, S., Taylor, M., Kolacz, K., et al. (1992). “Inhibitors of amino acid biosynthesis: strategies for imparting glyphosate tolerance to crop plants,” in Biosynthesis and Molecular Regulation of Amino Acids in Plants, eds D. Sinch, B. K. Flores, and J. C. Shannon (Rockville: American Society of Plant Physiologists), 139–145.

Google Scholar

Bohn, T., and Lovei, G. L. (2017). Complex outcomes from insect and weed control with transgenic plants: ecological surprises. Front. Environ. Sci. 5, 1–8. doi: 10.3389/fenvs.2017.00060

CrossRef Full Text | Google Scholar

Cao, G., Liu, Y., Zhang, S., Yang, X., Chen, R., Zhang, Y., et al. (2012). A novel 5-enolpyruvylshikimate-3-phosphate synthase shows high glyphosate tolerance in Escherichia coli and tobacco plants. PLoS One 7:e38718. doi: 10.1371/journal.pone.0038718

PubMed Abstract | CrossRef Full Text | Google Scholar

Chhapekar, S., Raghavendrarao, S., Pavan, G., Ramakrishna, C., Singh, V. K., Phanindra, M. L. V., et al. (2015). Transgenic rice expressing a codon-modified synthetic CP4-EPSPS confers tolerance to broad-spectrum herbicide, glyphosate. Plant Cell Rep. 34, 721–731. doi: 10.1007/s00299-014-1732-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Della-Cioppa, G., Bauer, S. C., Klein, B. K., Shah, D. M., Fraley, R. T., and Kishore, G. M. (1986). Translocation of the precursor of 5-enolpyruvylshikimate-3-phosphate synthase into chloroplasts of higher plants in vitro. Proc. Natl. Acad. Sci. U.S.A. 83, 6873–6877. doi: 10.1073/pnas.83.18.6873

PubMed Abstract | CrossRef Full Text | Google Scholar

Directorate of Cotton Development Government of India (2017). Status Paper of Indian Cotton. Nagpur: Directorate of Cotton Development Government of India.

Google Scholar

Dun, B., Wang, X., Lu, W., Chen, M., Zhang, W., Ping, S., et al. (2014). Development of highly glyphosate-tolerant tobacco by co-expression of glyphosate acetyltransferase gat and EPSPS G2-aroA genes. Crop. J. 2, 164–169. doi: 10.1016/j.cj.2014.03.003

CrossRef Full Text | Google Scholar

Fartyal, D., Agarwal, A., James, D., Borphukan, B., Ram, B., Sheri, V., et al. (2018). Co-expression of P173S mutant rice EPSPS and igrA genes results in higher glyphosate tolerance in transgenic rice. Front. Plant Sci. 9:144. doi: 10.3389/fpls.2018.00144

PubMed Abstract | CrossRef Full Text | Google Scholar

Fletcher, W. W. (1960). The Effect of Herbicides on Soil Micro-Organisms: Symposium Herbicides and the Soil. Oxford: Blackwell’s Science Publication, 20–62.

Google Scholar

Frisvold, G. B., Hurley, T. M., and Mitchell, P. D. (2009). Adoption of best management practices to control weed resistance by corn, cotton, and soybean growers. Aging Bioforum. 12, 370–381.

Google Scholar

Funke, T., Han, H., Healy-Fried, M. L., Fischer, M., and Schönbrunn, E. (2006). Molecular basis for the herbicide resistance of roundup ready crops. Proc. Natl. Acad. Sci. U.S.A. 103, 13010–13015. doi: 10.1073/pnas.0603638103

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaitonde, M. K., and Gordon, M. W. (1958). A microchemical method for the detection and determination of shikimic acid. J. Biol. Chem. 230, 1043–1050.

Google Scholar

Gianessi, L. P. (2005). Economic and herbicide use impacts of glyphosate resistant crops. Pest Manag. Sci. 61, 241–245. doi: 10.1002/ps.1013

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, Y. Y., Xu, Z. Z., Guo, S. Q., Shu, H. M., Jiang, L., and Ni, W. C. (2016). Genome-wide analysis of the EPSPS genes in upland cotton. Acta. Agric. Bor. Sin. 31, 15–21.

Google Scholar

Heap, I. (2014). Global perspective of herbicide resistant weeds. Pest Manag. Sci. 70, 1306–1315. doi: 10.1002/ps.3696

PubMed Abstract | CrossRef Full Text | Google Scholar

Hiscox, J. D., and Israelstam, G. F. (1979). A method for the extraction of chlorophyll from leaf tissue without maceration. Can. J. Bot. 57, 1332–1334. doi: 10.1139/b79-163

CrossRef Full Text | Google Scholar

Jabran, K. (2016). Weed flora, yield losses and weed control in cotton crop. Julius Kühn Archiv. 452, 177–182.

Google Scholar

Jain, S. C., Iyer, B. G., Jain, H. C., and Jain, N. K. (1981). “Weed management and nutrient losses in upland cotton under different ecosystems of Madhya Pradesh,” in Proceedings of 8th Asian – Pacific Weed Science Society, Bangalore, 131–135.

Google Scholar

Janaki, P., Sharma, N., Chinnusamy, C., Sakthivel, N., and Nithya, C. (2015). Herbicide residues and their management strategies. Indian. J. Weed Sci. 47, 329–344.

Google Scholar

Karthik, K., and Rohini, S. (2017). Apical meristem-targeted In planta transformation strategy: an overview on its utility in crop improvement. Agri. Res. Tech Open Access J. 8:555734. doi: 10.19080/ARTOAJ.2017.08.555734

CrossRef Full Text | Google Scholar

Keller, G., Spatola, L., Cabe, M. C., Martinell, D. B., Swain, W. B., and John, M. (1997). Transgenic cotton resistant to herbicide bialaphos. Transgenic. Res. 6, 385–392.

Google Scholar

Keshamma, E., Rohini, S., Rao, K. S., Madhusudhan, B., and Kumar, M. U. (2008). Tissue culture-independent in planta transformation strategy: an Agrobacterium tumefaciens-mediated gene transfer method to overcome recalcitrance in cotton (Gossypium hirsutum L.). J. Cotton Sci. 12, 264–272.

Google Scholar

Keshavareddy, G., Rohini, S., Ramu, S. V., Sundaresha, S., Kumar, A. R. V., Ananda Kumar, P., et al. (2013). Transgenics in groundnut (Arachis hypogaea L.) expressing cry1AcF gene for resistance to Spodoptera litura (F.). Physiol. Mol. Biol. Plants 19, 343–352. doi: 10.1007/s12298-013-0182-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kishore, G. M., Brundage, L., Kolk, K., Padgette, S. R., Rochester, D., Huynh, Q. K., et al. (1986). Isolation, purification and characterization of a glyphosate tolerant mutant E. Coli EPSP synthase. Proc. Fed. Am. Soc. Exp. Biol. 45:1506.

Google Scholar

Koger, C. H., and Reddy, K. N. (2005). Role of absorption and translocation in the mechanism of glyphosate resistance in horseweed (Conyza canadensis). Weed Sci. 53, 84–89. doi: 10.1614/ws-04-102r

CrossRef Full Text | Google Scholar

Kumar, A. M., Sreevathsa, R., Reddy, K. N., Ganesh, P. T., and Udayakumar, M. (2011). Amenability of castor to an Agrobacterium-mediated in planta transformation strategy using a cry1AcF gene for insect tolerance. J. Crop Sci. Biotechnol. 14, 125–132. doi: 10.1007/s12892-010-0013-z

CrossRef Full Text | Google Scholar

Liang, C. Z., Sun, B., Meng, Z. G., Meng, Z. H., Wang, Y., Sun, G. Q., et al. (2017). Co expression GR79 EPSPS and GAT yields herbicide-resistant cotton with low glyphosate residues. Plant Biotechnol. J. 15, 1622–1629. doi: 10.1111/pbi.12744

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyon, B. R., Cousins, Y. L., Llewellyn, D. J., and Dennis, E. S. (1993). Cotton plants transformed with a bacterial degradation gene are protected from accidental spray drift damage by the herbicide 2, 4-dichlorophenoxyacetic acid. Transgenic Res. 2, 162–169. doi: 10.1007/bf01972610

CrossRef Full Text | Google Scholar

Ma, X. Y., Wu, H. W., Jiang, W. L., Ma, Y. J., and Ma, Y. (2016). Weed and insect control affected by mixing insecticides with glyphosate in cotton. J. Integr. Agr. 15, 373–380. doi: 10.1016/s2095-3119(15)61188-1

CrossRef Full Text | Google Scholar

Mohapatra, L., and Saha, G. (2019). “Cotton Farming in India: alternative perspectives and paradigms,” in Transition Strategies for Sustainable Community Systems, ed. D. A. Nayak (Berlin: Springer), 195–213. doi: 10.1007/978-3-030-00356-2_17

CrossRef Full Text | Google Scholar

Nida, D. L., Kolacz, K. H., Buehler, R. E., Deaton, W. R., Schuler, W. R., Armstrong, T. A., et al. (1996). Glyphosate-tolerant cotton: genetic characterization and protein expression. J. Agr. Food Chem. 44, 1960–1966. doi: 10.1021/jf9505640

CrossRef Full Text | Google Scholar

Padgette, S. R., Re, D. B., Berry, G. F., Eichholtz, D. E., Delannay, X., Fuchs, R. L., et al. (1996). “New weed control opportunities: development of soybeans with a Roundup Ready gene,” in Herbicide-Resistant Crops: Agricultural, Economic, Environmental, Regulatory and Technological Aspects, ed. D. Duke (Boca Raton, FL: CRC Press), 53–84. doi: 10.1201/9781351073196-4

CrossRef Full Text | Google Scholar

Padgette, S. R., Re, D. B., Gasser, C. S., Eichholtz, D. A., Frazier, R. B., Hironaka, C. M., et al. (1991). Site directed mutagenesis of a conserved region of the 5-enolpyruvylshikimic acid-3-phosphate synthase active site. J. Biol. Chem. 266, 22361–22369.

Google Scholar

Peters, B., and Strek, H. J. (2018). Herbicide discovery in light of rapidly spreading resistance and ever−increasing regulatory hurdles. Pest Manag. Sci. 74, 2211–2215. doi: 10.1002/ps.4768

PubMed Abstract | CrossRef Full Text | Google Scholar

Porebski, S., Bailey, L. G., and Baum, B. R. (1997). Modification of a CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components. Plant Mol. Biol. Rep. 15, 8–15. doi: 10.1007/bf02772108

CrossRef Full Text | Google Scholar

Rajasekaran, K., Grula, J. W., and Anderson, D. M. (1996a). Selection and characterization of mutant cotton (Gossypium hirsutum L.) cell lines resistant to sulfonylurea and imidazolinone herbicides. Plant Sci. 119, 115–124. doi: 10.1016/0168-9452(96)04444-5

CrossRef Full Text | Google Scholar

Rajasekaran, K., Grula, J. W., Hudspeth, R. L., Pofelis, S., and Anderson, D. M. (1996b). Herbicide-resistant Acala and Coker cottons transformed with a native gene encoding mutant forms of acetohydroxyacid synthase. Mol. Breed. 2, 307–319. doi: 10.1007/bf00437909

CrossRef Full Text | Google Scholar

Ramachandra, D., Ramamohan, G., Bhan, A., and Suresh, P. J. (2016). Weed management in cotton: The potential of GM crops. Indian. J. Weed Sci. 48, 136–143.

Google Scholar

Rohini, V. K., and Rao, K. S. (2002). In planta strategy for gene transfer into plants: embryo transformation. Physiol. Mol. Biol. Plants. 8, 161–169.

Google Scholar

Shaner, D. L., Nadler-Hassar, T., Henry, W. B., and Koger, C. H. (2005). A rapid in vivo shikimate accumulation assay with excised leaf discs. Weed Sci. 53, 769–774. doi: 10.1614/ws-05-009r.1

CrossRef Full Text | Google Scholar

Sharma, R. (2008). Integrated weed management in field crops. Crop. Care. 35, 41–46.

Google Scholar

Singh, S., Kumar, N. R., Maniraj, R., Lakshmikanth, R., Rao, K. Y. S., Muralimohan, N., et al. (2018). Expression of Cry2Aa, a Bacillus thuringiensis insecticidal protein in transgenic pigeon pea confers resistance to gram pod borer, Helicoverpa armigera. Sci. Rep. 8:8820.

Google Scholar

Stalker, D. M., Hiatt, W. R., and Comai, L. D. (1985). A single amino acid substitution in the enzyme 5-enolpyruvylshikimic acid-3-phosphate synthase confers resistance to herbicide glyphosate. J. Biol. Chem. 260, 4724–4728.

Google Scholar

Steinrucken, H. C., and Amrhein, N. (1980). The herbicide glyphosate is a potent inhibitor of 5-enolpyruvyl-shikimic acid-3-phosphate synthase. Biochem. Biophys. Res. Commun. 94, 1207–1212. doi: 10.1016/0006-291x(80)90547-1

CrossRef Full Text | Google Scholar

Tian, Y. S., Xiong, A. S., Xu, J., Zhao, W., Gao, F., and Fu, X. Y. (2010). Isolation from Ochrobactrum anthropi of a novel class II 5-enopyruvylshikimate-3phosphate synthase with high tolerance to glyphosate. Appl. Environ. Microbiol. 76, 6001–6005. doi: 10.1128/aem.00770-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, T. A., Rajasekaran, K., and Anderson, D. M. (2000). Cotton biotechnology. Crit. Rev. Plant Sci. 19, 511–550.

Google Scholar

Wu, X., Gong, F., and Wang, W. (2014). Protein extraction from plant tissues for 2DE and its application in proteomic analysis. Proteomics 14, 645–658. doi: 10.1002/pmic.201300239

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, D., Ma, L., Lin, M., and Li, C. (2018). Development of glyphosate-resistant alfalfa (Medicagosativa L.) upon transformation with the GR79Ms gene encoding 5-enolpyruvylshikimate-3-phosphate synthase. Planta 248, 211–219. doi: 10.1007/s00425-018-2898-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. B., Tang, Q. L., Wang, X. J., and Wang, Z. X. (2017). Development of glyphosate-tolerant transgenic cotton plants harboring the G2-aroA gene. J. Integr. Agric. 16, 551–558. doi: 10.1016/s2095-3119(16)61458-2

CrossRef Full Text | Google Scholar

Zhou, M., Xu, H., Wei, X., Ye, Z., Wei, L., Gong, W., et al. (2006). Identification of a glyphosate-resistant mutant of rice 5-enolpyruvylshikimate 3-phosphate synthase using a directed evolution strategy. Plant Physiol. 140, 184–195. doi: 10.1104/pp.105.068577

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: transgenic cotton, weeds, CP4-EPSPS, glyphosate, N-phosphonomethyl glycine, herbicide tolerance, in planta transformation

Citation: Karthik K, Nandiganti M, Thangaraj A, Singh S, Mishra P, Rathinam M, Sharma M, Singh NK, Dash PK and Sreevathsa R (2020) Transgenic Cotton (Gossypium hirsutum L.) to Combat Weed Vagaries: Utility of an Apical Meristem-Targeted in planta Transformation Strategy to Introgress a Modified CP4-EPSPS Gene for Glyphosate Tolerance. Front. Plant Sci. 11:768. doi: 10.3389/fpls.2020.00768

Received: 30 January 2020; Accepted: 14 May 2020;
Published: 07 July 2020.

Edited by:

Yule Liu, Tsinghua University, China

Reviewed by:

Arun Jagannath, University of Delhi, India
Jia-He Wu, Institute of Microbiology, Chinese Academy of Sciences, China

Copyright © 2020 Karthik, Nandiganti, Thangaraj, Singh, Mishra, Rathinam, Sharma, Singh, Dash and Sreevathsa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rohini Sreevathsa,; Prasanta K. Dash,